VirSorter2: a multi-classifier, expert-guided approach to detect diverse DNA and RNA viruses

Metagenomic perspectives on antibiotic resistance genes in tap water: The environmental characteristic, potential mobility and health threat

As an emerging environmental contaminant, antibiotic resistance genes (ARGs) in tap water have attracted great attention. Although studies have provided ARG profiles in tap water, research on their abundance levels, composition characteristics, and potential threat is still insufficient. Here, 9 household tap water samples were collected from the Guangdong-Hong Kong-Macao Greater Bay Area (GBA) in China. Additionally, 75 sets of environmental sample data (9 types) were downloaded from the public database. Metagenomics was then performed to explore the differences in the abundance and composition of ARGs. 221 ARG subtypes consisting of 17 types were detected in tap water. Although the ARG abundance in tap water was not significantly different from that found in drinking water plants and reservoirs, their composition varied. In tap water samples, the three most abundant classes of resistance genes were multidrug, fosfomycin and MLS (macrolide-lincosamide-streptogramin) ARGs, and their corresponding subtypes ompR, fosX and macB were also the most abundant ARG subtypes. Regarding the potential mobility, vanS had the highest abundance on plasmids and viruses, but the absence of key genes rendered resistance to vancomycin ineffective. Generally, the majority of ARGs present in tap water were those that have not been assessed and are currently not listed as high-threat level ARG families based on the World Health Organization Guideline. Although the current potential threat to human health posed by ARGs in tap water is limited, with persistent transfer and accumulation, especially in pathogens, the potential danger to human health posed by ARGs should not be ignored.

Unveiling the overlooked small-sized microbiome in river ecosystems

Enriching microorganisms using a 0.22-μm pore size is a general pretreatment procedure in river microbiome research. However, it remains unclear the extent to which this method loses microbiome information. Here, we conducted a comparative metagenomics-based study on microbiomes with sizes over 0.22 μm (large-sized) and between 0.22 μm and 0.1 μm (small-sized) in a subtropical river. Although the absolute concentration of small-sized microbiome was about two orders of magnitude lower than that of large-sized microbiome, sequencing only large-sized microbiome resulted in a significant loss of microbiome diversity. Specifically, the microbial community was different between two sizes, and 347 genera were only detected in small-sized microbiome. Small-sized microbiome had much more diverse viral community than large-sized fraction. The viruses had abundant ecological functions and were hosted by 825 species of 169 families, including pathogen-related families. Small-sized microbiome had distinct antimicrobial resistance risks from large-sized microbiome, showing an enrichment of eight antibiotic resistance gene (ARG) types as well as the detection of 140 unique ARG subtypes and five enriched risk rank I ARGs. Draft genomes of five major resistant pathogens having diverse ecological and pollutant-degrading functions were only assembled in small-sized microbiome. These findings provide novel insights into river ecosystems, and highlight the overlooked small-sized microbiome in the environment.

Decrypting the viral community in aerobic activated sludge reactors treating antibiotic production wastewater

Viruses are the most abundant yet understudied members that may influence microbial metabolism in activated sludge treating antibiotic production wastewater. This study comprehensively investigated virome community characteristics under the selection pressure of nine types and different concentrations of antibiotics using a metagenomics approach. Of the 15,514 total viral operational taxonomic units (tOTUs) recovered, only 37.5 % were annotated. Antibiotics altered the original viral community structure in activated sludge. The proportion of some pathogenic viral families, including Herpesviridae_like, increased significantly in reactors treating erythromycin production wastewater. In total, 16.5 % of the tOTUs were associated with two or more hosts. tOTUs rarely carried antibiotic resistance genes (ARGs), and the ARG types in the tOTUs did not match the ARGs carried by the bacterial hosts. This suggests that transduction contributes little to the horizontal ARG transfer. Auxiliary metabolic genes (AMGs) were prevalent in tOTUs, and those involved in folate biosynthesis were particularly abundant, indicating their potential to mitigate antibiotic-induced host damage. This study provides comprehensive insights into the virome community in activated sludge treating antibiotic production wastewater and sheds light on the potential role of viral AMGs in mitigating antibiotic-induced stress.

Investigating potential auxiliary anaerobic digestion activity of phage under polyvinyl chloride microplastic stress

Polyvinyl chloride (PVC) microplastics present in sewage were trapped in sludge, thereby hindering anaerobic digestion performance of waste active sludge (WAS). Phages regulate virocell metabolism by encoding auxiliary metabolic genes (AMGs) related to energy acquisition and material degradation, supporting hosts survive in harsh environments and play a crucial role in biogeochemical cycles. This study investigated the potential effects of phages on the recovery of WAS anaerobic digestion under PVC stress. We observed a significant alteration in the phage community induced by PVC microplastics. Phages encoded AMGs related to anaerobic digestion and cell growth probably alleviate PVC microplastics inhibition on WAS anaerobic digestion, and 54.2 % of hydrolysis-related GHs and 40.8 % of acidification-related AMGs were actively transcribed in the PVC-exposed group. Additionally, the degradation of chitin and peptidoglycan during hydrolysis and the conversion of glucose to pyruvate during acidification were more susceptible to phages. Prediction of phage-host relationship indicated that the phyla Pseudomonadota were predominantly targeted hosts by hydrolysis-related and acidification-related phages, and PVC toxicity had minimal impact on phage-host interaction. Our findings highlight the importance of phages in anaerobic digestion and provide a novel strategy for using phages in the functional recovery of microplastic-exposed sludge.

Host population structure and species resolution reveal prophage transmission dynamics

Much knowledge about bacteriophages has been obtained via genomics and metagenomics over the last decades. However, most studies dealing with prophage diversity have rarely conducted phage species delimitation (aspect 1) and have hardly integrated the population structure of the host (aspect 2). Yet, these two aspects are essential in assessing phage diversity. Here, we implemented an operational definition of phage species (clustering at 95% identity, 90% coverage) and integrated the host's population structure to understand prophage diversity better. Gathering the most extensive data set of Acinetobacter baumannii phages (4,152 prophages + 122 virulent phages, distributed in 46 countries in the world), we show that 91% (875 out of 963) of the prophage species have four or fewer prophages per species, and just five prophage species have more than 100 prophages. Most prophage species have a narrow host range and are geographically restricted; yet, very few have a broad host range being well spread in distant lineages of A. baumannii. These few broad host range prophage species are not only cosmopolitan but also the most abundant species. We also noted that polylysogens had very divergent prophages, belonging to different prophage species, and prophages can easily be gained and lost within the bacterial lineages. Finally, even with this extensive data set, the prophage diversity has not been fully grasped. Our study highlights how integrating the host population structure and a solid operational definition of phage species allows us to better appreciate phage diversity and its transmission dynamics.

IMPORTANCE

Much knowledge about bacteriophages has been obtained via genomics and metagenomics over the last decades. However, most studies dealing with prophage diversity have rarely conducted phage species delimitation (aspect 1) and have hardly integrated the population structure of the host (aspect 2). Yet, these two aspects are essential in assessing phage diversity. Here, we implemented an operational definition of phage species (clustering at 95% identity, 90% coverage) and integrated the host's population structure to understand prophage diversity better. Gathering the most extensive data set of Acinetobacter baumannii phages, we show that most prophage species have four or fewer prophages per species, and just five prophage species have more than 100 prophages. Most prophage species have a narrow host range and are geographically restricted; yet, very few have a broad host range being well spread in distant lineages of A. baumannii. These few broad host range prophage species are cosmopolitan and the most abundant species. Prophages in the same bacterial genome are very divergent, and prophages can easily be gained and lost within the bacterial lineages. Finally, even with this extensive data set, the prophage diversity has not been fully grasped. This study shows how integrating the host population structure and clustering at the species level allows us to better appreciate phage diversity and its transmission dynamics.

Decoupling of strain- and intrastrain-level interactions of microbiomes in a sponge holobiont

Holobionts are highly organized assemblages of eukaryotic hosts, cellular microbial symbionts, and viruses, whose interactions and evolution involve complex biological processes. It is largely unknown which specific determinants drive similarity or individuality in genetic diversity between holobionts. Here, we combine short- and long-read sequencing and DNA-proximity-linkage technologies to investigate intraspecific diversity of the microbiomes, including host-resolved viruses, in individuals of a model marine sponge. We find strong impacts of the sponge host and the cellular hosts of viruses on strain-level organization of the holobiont, whereas substantial overlap in nucleotide diversity between holobionts suggests frequent exchanges of microbial cells and viruses at intrastrain level in the local sponge population. Immune-evasive arms races likely restricted virus-host co-evolution at the intrastrain level, generated holobiont-specific genome variations, and linked virus-host genetics through recombination. Our work shows that a decoupling of strain- and intrastrain-level interactions is a key factor in the genetic diversification of holobionts.

4CAC: 4-class classifier of metagenome contigs using machine learning and assembly graphs

Microbial communities usually harbor a mix of bacteria, archaea, plasmids, viruses and microeukaryotes. Within these communities, viruses, plasmids, and microeukaryotes coexist in relatively low abundance, yet they engage in intricate interactions with bacteria. Moreover, viruses and plasmids, as mobile genetic elements, play important roles in horizontal gene transfer and the development of antibiotic resistance within microbial populations. However, due to the difficulty of identifying viruses, plasmids, and microeukaryotes in microbial communities, our understanding of these minor classes lags behind that of bacteria and archaea. Recently, several classifiers have been developed to separate one or more minor classes from bacteria and archaea in metagenome assemblies. However, these classifiers often overlook the issue of class imbalance, leading to low precision in identifying the minor classes. Here, we developed a classifier called 4CAC that is able to identify viruses, plasmids, microeukaryotes, and prokaryotes simultaneously from metagenome assemblies. 4CAC generates an initial four-way classification using several sequence length-adjusted XGBoost models and further improves the classification using the assembly graph. Evaluation on simulated and real metagenome datasets demonstrates that 4CAC substantially outperforms existing classifiers and combinations thereof on short reads. On long reads, it also shows an advantage unless the abundance of the minor classes is very low. 4CAC runs 1-2 orders of magnitude faster than the other classifiers. The 4CAC software is available at https://github.com/Shamir-Lab/4CAC.

An expanded database and analytical toolkit for identifying bacterial virulence factors and their associations with chronic diseases

Virulence factor genes (VFGs) play pivotal roles in bacterial infections and have been identified within the human gut microbiota. However, their involvement in chronic diseases remains poorly understood. Here, we establish an expanded VFG database (VFDB 2.0) consisting of 62,332 nonredundant orthologues and alleles of VFGs using species-specific average nucleotide identity ( https://github.com/Wanting-Dong/MetaVF_toolkit/tree/main/databases ). We further develop the MetaVF toolkit, facilitating the precise identification of pathobiont-carried VFGs at the species level. A thorough characterization of VFGs for 5452 commensal isolates from healthy individuals reveals that only 11 of 301 species harbour these factors. Further analyses of VFGs within the gut microbiomes of nine chronic diseases reveal both common and disease-specific VFG features. Notably, in type 2 diabetes patients, long HiFi sequencing confirms that shared VF features are carried by pathobiont strains of Escherichia coli and Klebsiella pneumoniae. These findings underscore the critical importance of identifying and understanding VFGs in microbiome-associated diseases.

Virome and metagenomic sequencing reveal the impact of microbial inoculants on suppressions of antibiotic resistome and viruses during co-composting

Co-composting with exogenous microbial inoculant, presents an effective approach for the harmless utilization of livestock manure and agroforestry wastes. However, the impact of inoculant application on the variations of viral and antibiotic resistance genes (ARGs) remains poorly understood, particularly under varying manure quantity (low 10 % vs. high 20 % w/w). Thus, employing virome and metagenomic sequencing, we examined the influence of Streptomyces-Bacillus Inoculants (SBI) on viral communities, phytopathogen, ARGs, mobile genetic elements, and their interrelations. Our results indicate that SBI shifted dominant bacterial species from Phenylobacterium to thermotropic Bordetella, and the quantity of manure mediates the effect of SBI on whole bacterial community. Major ARGs and genetic elements experienced substantial changes with SBI addition. There was a higher ARGs elimination rate in the composts with low (∼76 %) than those with high manure (∼70 %) application. Virus emerged as a critical factor influencing ARG dynamics. We observed a significant variation in virus community, transitioning from Gemycircularvirus- (∼95 %) to Chlamydiamicrovirus-dominance. RDA analysis revealed that Gemycircularvirus was the most influential taxon in shaping ARGs, with its abundance decreased approximately 80 % after composting. Collectively, these findings underscore the role of microbial inoculants in modulating virus communities and ARGs during biowaste co-composting.

Unveiling the genomic diversity and ecological impact of phage communities in hospital wastewater

Phages are pivotal in shaping microbial communities and biogeochemical cycles, while our understanding of the diversity, functions potential, and resistance gene carriage of phages in hospital wastewater (HWW) remains limited. We collected influent and effluent samples from the 3 hospital wastewater treatment plants (HWTPs) to assess the diversity and fate of phages, the interactions between phages and hosts, and the presence of resistance genes and auxiliary metabolic genes (AMGs) encoded by phages. Compared to influent, effluent showed reduced phage abundance and altered composition, with decreases in Microviridae and Inoviridae. The gene-sharing network highlights that many phages in HWW are not classified in known viral genera, suggesting HWW as a rich source of new viruses. There was a significant association between phages and microorganisms, with approximately 32.57 % of phages expected to be capable of infecting microbial hosts, characterized primarily by lytic activity. A total of 8 unique antibiotic resistance genes, 13 unique metal resistance genes, and 5 mobile genetic elements were detected in 3 HWTPs phageomes. Phage AMGs have the potential to influence carbon, nitrogen, phosphorus, and sulfur metabolism, impacting biogeochemical cycles. This study reveals the genomic diversity and ecological role of phages in HWTPs, highlighting their environmental and ecosystem impact.

Viral diversity within marine biofilms and interactions with corrosive microbes

In marine environments, a wide variety of microbes like bacteria, and archaea influence on the corrosion of materials. Viruses are widely distributed in biofilms among these microbes and may affect the corrosion process through interactions with key corrosive prokaryotes. However, understanding of the viral communities within biofilms and their interactions with corrosive microbes remains is limited. To improve this knowledge gap, 53 metagenomes were utilized to investigate the diversity of viruses within biofilms on 8 different materials and their interactions with corrosive microbes. Notably, the viruses within biofilms predominantly belonged to Caudoviricetes, and phylogenetic analysis of Caudoviricetes and protein-sharing networks with other environments revealed the presence of numerous novel viral clades in biofilms. The virus‒host linkages revealed a close association between viruses and corrosive microbes in biofilms. This means that viruses may modulate host corrosion-related metabolism through auxiliary metabolic genes. It was observed that the virus could enhance host resistance to metals and antibiotics via horizontal gene transfer. Interestingly, viruses could protect themselves from host antiviral systems through anti-defense systems. This study illustrates the diversity of viruses within biofilms formed on materials and the intricate interactions between viruses and corrosive microbes, showing the potential roles of viruses in corrosive biofilms.

Genomic insights into the antimicrobial resistance and virulence of Helicobacter pylori isolates from gastritis patients in Pereira, Colombia

BACKGROUND

Helicobacter pylori infects the stomach and/or small intestines in more than half of the human population. Infection with H. pylori is the most common cause of chronic gastritis, which can lead to more severe gastroduodenal pathologies such as peptic ulcer, mucosa-associated lymphoid tissue lymphoma, and gastric cancer. H. pylori infection is particularly concerning in Colombia in South America, where > 80% of the population is estimated to be infected with H. pylori and the rate of stomach cancer is one of the highest in the continent.

RESULTS

We compared the antimicrobial susceptibility profiles and short-read genome sequences of five H. pylori isolates obtained from patients diagnosed with gastritis of varying severity (chronic gastritis, antral erosive gastritis, superficial gastritis) in Pereira, Colombia sampled in 2015. Antimicrobial susceptibility tests revealed the isolates to be resistant to at least one of the five antimicrobials tested: four isolates were resistant to metronidazole, two to clarithromycin, two to levofloxacin, and one to rifampin. All isolates were susceptible to tetracycline and amoxicillin. Comparative genome analyses revealed the presence of genes associated with efflux pump, restriction modification systems, phages and insertion sequences, and virulence genes including the cytotoxin genes cagA and vacA. The five genomes represent three novel sequence types. In the context of the Colombian and global populations, the five H. pylori isolates from Pereira were phylogenetically distant to each other but were closely related to other lineages circulating in the country.

CONCLUSIONS

H. pylori from gastritis of different severity varied in their antimicrobial susceptibility profiles and genome content. This knowledge will be useful in implementing appropriate eradication treatment regimens for specific types of gastritis. Understanding the genetic and phenotypic heterogeneity in H. pylori across the geographical landscape is critical in informing health policies for effective disease prevention and management that is most effective at local and country-wide scales. This is especially important in Colombia and other South American countries that are poorly represented in global genomic surveillance studies of bacterial pathogens.

Metagenomic analysis reveals high diversity of gut viromes in yaks (Bos grunniens) from the Qinghai-Tibet Plateau

The Qinghai-Tibet Plateau (QTP), renowned for its exceptional biological diversity, is home to numerous endemic species. However, research on the virology of vulnerable vertebrates like yaks remains limited. In this study, our objective was to use metagenomics to provide a comprehensive understanding of the diversity and evolution of the gut virome in yak populations across different regions of the QTP. Our findings revealed a remarkably diverse array of viruses in the gut of yaks, including those associated with vertebrates and bacteriophages. Notably, some vertebrate-associated viruses, such as astrovirus and picornavirus, showed significant sequence identity across diverse yak populations. Additionally, we observed differences in the functional profiles of genes carried by the yak gut virome across different regions. Moreover, the virus-bacterium symbiotic network that we discovered holds potential significance in maintaining the health of yaks. Overall, this research expands our understanding of the viral communities in the gut of yaks and highlights the importance of further investigating the interactions between viruses and their hosts. These data will be beneficial for revealing the crucial role that viruses play in the yak gut ecology in future studies.

Integrated analysis of facial microbiome and skin physio-optical properties unveils cutotype-dependent aging effects

BACKGROUND

Our facial skin hosts millions of microorganisms, primarily bacteria, crucial for skin health by maintaining the physical barrier, modulating immune response, and metabolizing bioactive materials. Aging significantly influences the composition and function of the facial microbiome, impacting skin immunity, hydration, and inflammation, highlighting potential avenues for interventions targeting aging-related facial microbes amidst changes in skin physiological properties.

RESULTS

We conducted a multi-center and deep sequencing survey to investigate the intricate interplay of aging, skin physio-optical conditions, and facial microbiome. Leveraging a newly-generated dataset of 2737 species-level metagenome-assembled genomes (MAGs), our integrative analysis highlighted aging as the primary driver, influencing both facial microbiome composition and key skin characteristics, including moisture, sebum production, gloss, pH, elasticity, and sensitivity. Further mediation analysis revealed that skin characteristics significantly impacted the microbiome, mostly as a mediator of aging. Utilizing this dataset, we uncovered two consistent cutotypes across sampling cities and identified aging-related microbial MAGs. Additionally, a Facial Aging Index (FAI) was formulated based on the microbiome, uncovering the cutotype-dependent effects of unhealthy lifestyles on skin aging. Finally, we distinguished aging related microbial pathways influenced by lifestyles with cutotype-dependent effect.

CONCLUSIONS

Together, our findings emphasize aging's central role in facial microbiome dynamics, and support personalized skin microbiome interventions by targeting lifestyle, skin properties, and aging-related microbial factors. Video Abstract.

Diverse RNA viruses of parasitic nematodes can elicit antibody responses in vertebrate hosts

Parasitic nematodes have an intimate, chronic and lifelong exposure to vertebrate tissues. Here we mined 41 published parasitic nematode transcriptomes from vertebrate hosts and identified 91 RNA viruses across 13 virus orders from 24 families in ~70% (28 out of 41) of parasitic nematode species, which include only 5 previously reported viruses. We observe widespread distribution of virus-nematode associations across multiple continents, suggesting an ancestral acquisition event and host-virus co-evolution. Characterization of viruses of Brugia malayi (BMRV1) and Onchocerca volvulus (OVRV1) shows that these viruses are abundant in reproductive tissues of adult parasites. Importantly, the presence of BMRV1 RNA in B. malayi parasites mounts an RNA interference response against BMRV1 suggesting active viral replication. Finally, BMRV1 and OVRV1 were found to elicit antibody responses in serum samples from infected jirds and infected or exposed humans, indicating direct exposure to the immune system.

Long-term persistence of diverse clones shapes the transmission landscape of invasive Listeria monocytogenes

BACKGROUND

The foodborne bacterium Listeria monocytogenes (Lm) causes a range of diseases, from mild gastroenteritis to invasive infections that have high fatality rate in vulnerable individuals. Understanding the population genomic structure of invasive Lm is critical to informing public health interventions and infection control policies that will be most effective especially in local and regional communities.

METHODS

We sequenced the whole draft genomes of 936 Lm isolates from human clinical samples obtained in a two-decade active surveillance program across 58 counties in New York State, USA. Samples came mostly from blood and cerebrospinal fluid. We characterized the phylogenetic relationships, population structure, antimicrobial resistance genes, virulence genes, and mobile genetic elements.

RESULTS

The population is genetically heterogenous, consisting of lineages I-IV, 89 clonal complexes, 200 sequence types, and six known serogroups. In addition to intrinsic antimicrobial resistance genes (fosX, lin, norB, and sul), other resistance genes tetM, tetS, ermG, msrD, and mefA were sparsely distributed in the population. Within each lineage, we identified clusters of isolates with ≤ 20 single nucleotide polymorphisms in the core genome alignment. These clusters may represent isolates that share a most recent common ancestor, e.g., they are derived from the same contamination source or demonstrate evidence of transmission or outbreak. We identified 38 epidemiologically linked clusters of isolates, confirming eight previously reported disease outbreaks and the discovery of cryptic outbreaks and undetected chains of transmission, even in the rarely reported Lm lineage III (ST3171). The presence of animal-associated lineages III and IV may suggest a possible spillover of animal-restricted strains to humans. Many transmissible clones persisted over several years and traversed distant sites across the state.

CONCLUSIONS

Our findings revealed the bacterial determinants of invasive listeriosis, driven mainly by the diversity of locally circulating lineages, intrinsic and mobile antimicrobial resistance and virulence genes, and persistence across geographical and temporal scales. Our findings will inform public health efforts to reduce the burden of invasive listeriosis, including the design of food safety measures, source traceback, and outbreak detection.

Improving viral annotation with artificial intelligence

Viruses of bacteria, "phages," are fundamental, poorly understood components of microbial community structure and function. Additionally, their dependence on hosts for replication positions phages as unique sensors of ecosystem features and environmental pressures. High-throughput sequencing approaches have begun to give us access to the diversity and range of phage populations in complex microbial community samples, and metagenomics is currently the primary tool with which we study phage populations. The study of phages by metagenomic sequencing, however, is fundamentally limited by viral diversity, which results in the vast majority of viral genomes and metagenome-annotated genomes lacking annotation. To harness bacteriophages for applications in human and environmental health and disease, we need new methods to organize and annotate viral sequence diversity. We recently demonstrated that methods that leverage self-supervised representation learning can supplement statistical sequence representations for remote viral protein homology detection in the ocean virome and propose that consideration of the functional content of viral sequences allows for the identification of similarity in otherwise sequence-diverse viruses and viral-like elements for biological discovery. In this review, we describe the potential and pitfalls of large language models for viral annotation. We describe the need for new approaches to annotate viral sequences in metagenomes, the fundamentals of what protein language models are and how one can use them for sequence annotation, the strengths and weaknesses of these models, and future directions toward developing better models for viral annotation more broadly.

Epidemiological pattern and genomic insights into multidrug-resistant ST491 Acinetobacter baumannii BD20 isolated from an infected wound in Bangladesh: Concerning co-occurrence of three classes of beta lactamase genes

OBJECTIVE

Multidrug-resistant (MDR) Acinetobacter baumannii is a major issue within healthcare facilities in Bangladesh due to its frequent association with hospital-acquired infections. In this study we report on a carbapenem-resistant draft genome sequence of an A. baumannii BD20 sample isolated from an infected wound in Bangladesh.

METHODS

A. baumannii BD20 was isolated from an infected burn wound. Whole-genome sequencing was carried out and annotated using PGAP and Prokka. Sequence type, antimicrobial resistance genes, virulence factor genes, and metal resistance genes were investigated. Core genome multilocus sequence typing-based phylogenomic analysis between A. baumannii BD20 and 213 A. baumannii strains retrieved from the NCBI GenBank database was performed using the BacWGSTdb 2.0 server.

RESULTS

A. baumannii BD20 (MLST 491) was resistant to all the antibiotics tested, except for colistin and polymyxin B. Along with many other antibiotic resistance genes, the isolate harbored three classes of beta lactamase-producing genes: blaGES-11 (class A), blaOXA-69 (class D), blaADC-10 (class C), and blaADC-11 (class C). Additionally, the strain carried several virulence genes and metal resistance determinants, which may contribute to its increased virulence. Core genome MLST-based phylogenomic analysis revealed that A. baumannii BD20 was closely related to another ST491 strain isolated from Singapore.

CONCLUSIONS

The findings of this study underscore the growing challenge of MDR A. baumannii, emphasizing the need for vigilant surveillance and infection-control measures in healthcare settings in order to address these emerging threats effectively.

Soil keystone viruses are regulators of ecosystem multifunctionality

Ecosystem multifunctionality reflects the capacity of ecosystems to simultaneously maintain multiple functions which are essential bases for human sustainable development. Whereas viruses are a major component of the soil microbiome that drive ecosystem functions across biomes, the relationships between soil viral diversity and ecosystem multifunctionality remain under-studied. To address this critical knowledge gap, we employed a combination of amplicon and metagenomic sequencing to assess prokaryotic, fungal and viral diversity, and to link viruses to putative hosts. We described the features of viruses and their potential hosts in 154 soil samples from 29 farmlands and 25 forests distributed across China. Although 4,460 and 5,207 viral populations (vOTUs) were found in the farmlands and forests respectively, the diversity of specific vOTUs rather than overall soil viral diversity was positively correlated with ecosystem multifunctionality in both ecosystem types. Furthermore, the diversity of these keystone vOTUs, despite being 10-100 times lower than prokaryotic or fungal diversity, was a better predictor of ecosystem multifunctionality and more strongly associated with the relative abundances of prokaryotic genes related to soil nutrient cycling. Gemmatimonadota and Actinobacteria dominated the host community of soil keystone viruses in the farmlands and forests respectively, but were either absent or showed a significantly lower relative abundance in that of soil non-keystone viruses. These findings provide novel insights into the regulators of ecosystem multifunctionality and have important implications for the management of ecosystem functioning.

Massive expansion of the pig gut virome based on global metagenomic mining

The pig gut virome plays a vital role in the gut microbial ecosystem of pigs. However, a comprehensive understanding of their diversity and a reference database for the virome are currently lacking. To address this gap, we established a Pig Virome Database (PVD) that comprised of 5,566,804 viral contig sequences from 4650 publicly available gut metagenomic samples using a pipeline designated "metav". By clustering sequences, we identified 48,299 viral operational taxonomic units (vOTUs) genomes of at least medium quality, of which 92.83% of which were not found in existing major databases. The majority of vOTUs were identified as Caudoviricetes (72.21%). The PVD database contained a total of 2,362,631 protein-coding genes across the above medium-quality vOTUs genomes that can be used to explore the functional potential of the pig gut virome. These findings highlight the extensive diversity of viruses in the pig gut and provide a pivotal reference dataset for forthcoming research concerning the pig gut virome.

Prokaryotic-virus-encoded auxiliary metabolic genes throughout the global oceans

BACKGROUND

Prokaryotic microbes have impacted marine biogeochemical cycles for billions of years. Viruses also impact these cycles, through lysis, horizontal gene transfer, and encoding and expressing genes that contribute to metabolic reprogramming of prokaryotic cells. While this impact is difficult to quantify in nature, we hypothesized that it can be examined by surveying virus-encoded auxiliary metabolic genes (AMGs) and assessing their ecological context.

RESULTS

We systematically developed a global ocean AMG catalog by integrating previously described and newly identified AMGs and then placed this catalog into ecological and metabolic contexts relevant to ocean biogeochemistry. From 7.6 terabases of Tara Oceans paired prokaryote- and virus-enriched metagenomic sequence data, we increased known ocean virus populations to 579,904 (up 16%). From these virus populations, we then conservatively identified 86,913 AMGs that grouped into 22,779 sequence-based gene clusters, 7248 (~ 32%) of which were not previously reported. Using our catalog and modeled data from mock communities, we estimate that ~ 19% of ocean virus populations carry at least one AMG. To understand AMGs in their metabolic context, we identified 340 metabolic pathways encoded by ocean microbes and showed that AMGs map to 128 of them. Furthermore, we identified metabolic "hot spots" targeted by virus AMGs, including nine pathways where most steps (≥ 0.75) were AMG-targeted (involved in carbohydrate, amino acid, fatty acid, and nucleotide metabolism), as well as other pathways where virus-encoded AMGs outnumbered cellular homologs (involved in lipid A phosphates, phosphatidylethanolamine, creatine biosynthesis, phosphoribosylamine-glycine ligase, and carbamoyl-phosphate synthase pathways).

CONCLUSIONS

Together, this systematically curated, global ocean AMG catalog and analyses provide a valuable resource and foundational observations to understand the role of viruses in modulating global ocean metabolisms and their biogeochemical implications. Video Abstract.

Host specificity and cophylogeny in the "animal-gut bacteria-phage" tripartite system

Cophylogeny has been identified between gut bacteria and their animal host and is highly relevant to host health, but little research has extended to gut bacteriophages. Here we use bee model to investigate host specificity and cophylogeny in the "animal-gut bacteria-phage" tripartite system. Through metagenomic sequencing upon different bee species, the gut phageome revealed a more variable composition than the gut bacteriome. Nevertheless, the bacteriome and the phageome showed a significant association of their dissimilarity matrices, indicating a reciprocal interaction between the two kinds of communities. Most of the gut phages were host generalist at the viral cluster level but host specialist at the viral OTU level. While the dominant gut bacteria Gilliamella and Snodgrassella exhibited matched phylogeny with bee hosts, most of their phages showed a diminished level of cophylogeny. The evolutionary rates of the bee, the gut bacteria and the gut phages showed a remarkably increasing trend, including synonymous and non-synonymous substitution and gene content variation. For all of the three codiversified tripartite members, however, their genes under positive selection and genes involving gain/loss during evolution simultaneously enriched the functions into metabolism of nutrients, therefore highlighting the tripartite coevolution that results in an enhanced ecological fitness for the whole holobiont.

Viral Diversity and Ecological Impact of DNA Viruses in Dominant Tick Species in China

Ticks are blood-feeding ectoparasites that also transmit various pathogens, posing severe risks to human and animal health. DNA viruses play a crucial role in the microbial ecology of ticks, but their distribution and ecological significance remain largely undetermined. Here, we assembled an extensive catalog encompassing 4320 viral operational taxonomic units (vOTUs) from six main dominant tick species in China, of which 94.8% have not been found in any other environment. To bridge the knowledge gap in tick DNA virus research and provide a crucial resource platform, we developed the Tick DNA Virus Database. This database includes the vOTUs that are known to cause diseases. Most of the predicted vOTUs are associated with dominant bacterial and archaeal phyla. We identified 105 virus-encoded putative auxiliary metabolic genes (AMGs) that are involved in host metabolism and environmental adaptation, potentially influencing ticks through both top-down and bottom-up mechanisms. The identification of microbial communities and antibiotic resistance in wild tick species suggests that wild ticks are reservoirs of antibiotic resistance and potential spreaders of antibiotic resistance. These findings reveal the potential role of tick viruses in ecosystems, highlighting the importance of monitoring tick microbiomes to address global public health challenges.

Genetic diversity and cross-species transmissibility of bat-associated picornaviruses from Spain

BACKGROUND

Emerging zoonotic diseases arise from cross-species transmission events between wild or domesticated animals and humans, with bats being one of the major reservoirs of zoonotic viruses. Viral metagenomics has led to the discovery of many viruses, but efforts have mainly been focused on some areas of the world and on certain viral families.

METHODS

We set out to describe full-length genomes of new picorna-like viruses by collecting feces from hundreds of bats captured in different regions of Spain. Viral sequences were obtained by high-throughput Illumina sequencing and analyzed phylogenetically to classify them in the context of known viruses. Linear discriminant analysis (LDA) was performed to infer likely hosts based on genome composition.

RESULTS

We found five complete or nearly complete genomes belonging to the family Picornaviridae, including a new species of the subfamily Ensavirinae. LDA suggested that these were true vertebrate viruses, rather than viruses from the bat diet. Some of these viruses were related to picornaviruses previously found in other bat species from distant geographical regions. We also found a calhevirus genome that most likely belongs to a proposed new family within the order Picornavirales, and for which genome composition analysis suggested a plant host.

CONCLUSIONS

Our findings describe new picorna-like viral species and variants circulating in the Iberian Peninsula, illustrate the wide geographical distribution and interspecies transmissibility of picornaviruses, and suggest new hosts for calheviruses.

Exploring the diversity of anti-defense systems across prokaryotes, phages, and mobile genetic elements

The co-evolution of prokaryotes, phages, and mobile genetic elements (MGEs) over the past billions of years has driven the emergence and diversification of defense and anti-defense systems alike. Anti-defense proteins have diverse functional domains, sequences, and are typically small, creating a challenge to detect anti-defense homologs across the prokaryotic genomes. To date, no tools comprehensively annotate anti-defense proteins within a desired genome or MGE. Here, we developed "AntiDefenseFinder" - a free open-source tool and web service that detects 156 anti-defense systems (of one or more proteins) in any genomic sequence. Using this dataset, we identified 47,981 anti-defense systems distributed across prokaryotes, phage, and MGEs. We found that some genes co-localize in "anti-defense islands", including E. coli T4 and Lambda phages, although many are standalone. Out of the 112 systems detected in bacteria, 100 systems localize only or preferentially in prophages, plasmids, phage satellites, integrons, and integrative and conjugative elements. However, over 80% of anti-Pycsar protein 1 (Apyc1) resides in non-mobile regions of bacteria. Evolutionary and functional analyses revealed that Apyc1 likely originated in bacteria to regulate cNMP signaling, but was co-opted multiple times by phages to overcome cNMP-utilizing defenses. With the AntiDefenseFinder tool, we hope to facilitate the identification of the full repertoire of anti-defense systems in MGEs, the discovery of new protein functions, and a deeper understanding of host-pathogen arms race.

Viral communities in a pH>10 serpentinite-like environment: insight into diversity and potential roles in modulating the microbiomes by bioactive vitamin B9 synthesis

Viral communities exist in a variety of ecosystems and play significant roles in mediating biogeochemical processes, whereas viruses inhabiting strongly alkaline geochemical systems remain underexplored. In this study, the viral diversity, potential functionalities, and virus-host interactions in a strongly alkaline environment (pH = 10.4-12.4) exposed to the leachates derived from the serpentinization-like reactions of smelting slags were investigated. The viral populations (e.g., Herelleviridae, Queuovirinae, and Inoviridae) were closely associated with the dominating prokaryotic hosts (e.g., Meiothermus, Trueperaceae, and Serpentinomonas) in this ultrabasic environment. Auxiliary metabolic genes (AMGs) suggested that viruses may enhance hosts' fitness by facilitating cofactor biosynthesis, hydrogen metabolism, and carbon cycling. To evaluate the activity of synthesis of essential cofactor vitamin B9 by the viruses, a viral folA (vfolA) gene encoding dihydrofolate reductase (DHFR) was introduced into a thymidine-auxotrophic strain Escherichia coli MG1655 ΔfolA mutant, which restored the growth of the latter in the absence of thymidine. Notably, the homologs of the validated vDHFR were globally distributed in the viromes across various ecosystems. The present study sheds new light on the unique viral communities in hyperalkaline ecosystems and their potential beneficial impacts on the coexisting microbial consortia by supplying essential cofactors.

IMPORTANCE

This study presents a comprehensive investigation into the diversity, potential functionalities, and virus-microbe interactions in an artificially induced strongly alkaline environment. Functional validation of the detected viral folA genes encoding dihydrofolate reductase substantiated the synthesis of essential cofactors by viruses, which may be ubiquitous, considering the broad distribution of the viral genes associated with folate cycling.

Viral and thermal lysis facilitates transmission of antibiotic resistance genes during composting

While the distribution of extracellular ARGs (eARGs) in the environment has been widely reported, the factors governing their release remain poorly understood. Here, we combined multi-omics and direct experimentation to test whether the release and transmission of eARGs are associated with viral lysis and heat during cow manure composting. Our results reveal that the proportion of eARGs increased 2.7-fold during composting, despite a significant and concomitant reduction in intracellular ARG abundances. This relative increase of eARGs was driven by composting temperature and viral lysis of ARG-carrying bacteria based on metagenome-assembled genome (MAG) analysis. Notably, thermal lysis of mesophilic bacteria carrying ARGs was a key factor in releasing eARGs at the thermophilic phase, while viral lysis played a relatively stronger role during the non-thermal phase of composting. Furthermore, MAG-based tracking of ARGs in combination with direct transformation experiments demonstrated that eARGs released during composting pose a potential transmission risk. Our study provides bioinformatic and experimental evidence of the undiscovered role of temperature and viral lysis in co-driving the spread of ARGs in compost microbiomes via the horizontal transfer of environmentally released DNA.

IMPORTANCE

The spread of antibiotic resistance genes (ARGs) is a critical global health concern. Understanding the factors influencing the release of extracellular ARGs (eARGs) is essential for developing effective strategies. In this study, we investigated the association between viral lysis, heat, and eARG release during composting. Our findings revealed a substantial increase in eARGs despite reduced intracellular ARG abundance. Composting temperature and viral lysis were identified as key drivers, with thermal lysis predominant during the thermophilic phase and viral lysis during non-thermal phases. Moreover, eARGs released during composting posed a transmission risk through horizontal gene transfer. This study highlights the significance of temperature and phage lysis in ARG spread, providing valuable insights for mitigating antibiotic resistance threats.

Bacillus subtilis-derived peptides disrupt quorum sensing and biofilm assembly in multidrug-resistant Staphylococcus aureus

Multidrug-resistant Staphylococcus aureus is one of the most clinically important pathogens in the world, with infections leading to high rates of morbidity and mortality in both humans and animals. The ability of S. aureus to form biofilms protects cells from antibiotics and promotes the transfer of antibiotic resistance genes; therefore, new strategies aimed at inhibiting biofilm growth are urgently needed. Probiotic species, including Bacillus subtilis, are gaining interest as potential therapies against S. aureus for their ability to reduce S. aureus colonization and virulence. Here, we search for strains and microbially derived compounds with strong antibiofilm activity against multidrug-resistant S. aureus by isolating and screening Bacillus strains from a variety of agricultural environments. From a total of 1,123 environmental isolates, we identify a single strain B. subtilis 6D1, with a potent ability to inhibit biofilm growth, disassemble mature biofilm, and improve antibiotic sensitivity of S. aureus biofilms through an Agr quorum sensing interference mechanism. Biochemical and molecular networking analysis of an active organic fraction revealed multiple surfactin isoforms, and an uncharacterized peptide was driving this antibiofilm activity. Compared with commercial high-performance liquid chromatography grade surfactin obtained from B. subtilis, we show these B. subtilis 6D1 peptides are significantly better at inhibiting biofilm formation in all four S. aureus Agr backgrounds and preventing S. aureus-induced cytotoxicity when applied to HT29 human intestinal cells. Our study illustrates the potential of exploring microbial strain diversity to discover novel antibiofilm agents that may help combat multidrug-resistant S. aureus infections and enhance antibiotic efficacy in clinical and veterinary settings.

IMPORTANCE

The formation of biofilms by multidrug-resistant bacterial pathogens, such as Staphylococcus aureus, increases these microorganisms' virulence and decreases the efficacy of common antibiotic regimens. Probiotics possess a variety of strain-specific strategies to reduce biofilm formation in competing organisms; however, the mechanisms and compounds responsible for these phenomena often go uncharacterized. In this study, we identified a mixture of small probiotic-derived peptides capable of Agr quorum sensing interference as one of the mechanisms driving antibiofilm activity against S. aureus. This collection of peptides also improved antibiotic killing and protected human gut epithelial cells from S. aureus-induced toxicity by stimulating an adaptive cytokine response. We conclude that purposeful strain screening and selection efforts can be used to identify unique probiotic strains that possess specially desired mechanisms of action. This information can be used to further improve our understanding of the ways in which probiotic and probiotic-derived compounds can be applied to prevent bacterial infections or improve bacterial sensitivity to antibiotics in clinical and agricultural settings.

Long-read sequencing reveals extensive gut phageome structural variations driven by genetic exchange with bacterial hosts

Genetic variations are instrumental for unraveling phage evolution and deciphering their functional implications. Here, we explore the underlying fine-scale genetic variations in the gut phageome, especially structural variations (SVs). By using virome-enriched long-read metagenomic sequencing across 91 individuals, we identified a total of 14,438 nonredundant phage SVs and revealed their prevalence within the human gut phageome. These SVs are mainly enriched in genes involved in recombination, DNA methylation, and antibiotic resistance. Notably, a substantial fraction of phage SV sequences share close homology with bacterial fragments, with most SVs enriched for horizontal gene transfer (HGT) mechanism. Further investigations showed that these SV sequences were genetic exchanged between specific phage-bacteria pairs, particularly between phages and their respective bacterial hosts. Temperate phages exhibit a higher frequency of genetic exchange with bacterial chromosomes and then virulent phages. Collectively, our findings provide insights into the genetic landscape of the human gut phageome.

Reproducible chemostat cultures to minimize eukaryotic viruses from fecal transplant material

Recent studies indicate an important role of bacteriophages for successful fecal microbiota transplantation (FMT). However, wider clinical applications of FMT are hampered by to donor variability and concerns of infection risks by bacteria and human viruses. To overcome these challenges, mouse cecal and human fecal material were propagated in a chemostat fermentation setup supporting multiplication of bacteria, and phages, while propagation of eukaryotic viruses will be prevented in the absence of eukaryotic host cells. The results showed decrease of the median relative abundance of viral contigs of classified eukaryotic viruses below 0.01%. The corresponding virome profiles showed dilution rate dependency, a reproducibility between biological replicates, and maintained high diversity regarding both the human and mouse inocula. This proof-of-concept cultivation approach may constitute the first step of developing novel therapeutic tools with high reproducibility and with low risk of infection from the donor material to target gut-related diseases.

A natural ANI gap that can define intra-species units of bacteriophages and other viruses

Despite the importance of intra-species variants of viruses for causing disease and/or disrupting ecosystem functioning, there is no universally applicable standard to define these. A (natural) gap in whole-genome average nucleotide identity (ANI) values around 95% is commonly used to define species, especially for bacteriophages, but whether a similar gap exists within species that can be used to define intra-species units has not been evaluated yet. Whole-genome comparisons among members of 1,016 bacteriophage (Caudoviricetes) species revealed a region of low frequency of ANI values around 99.2%-99.8%, showing threefold or fewer pairs than expected for an even distribution. This second gap is prevalent in viruses infecting various cultured or uncultured hosts from a variety of environments, although a few exceptions to this pattern were also observed (3.7% of total species) and are likely attributed to cultivation biases or other factors. Similar results were observed for a limited set of eukaryotic viruses that are adequately sampled, including SARS-CoV-2, whose ANI-based clusters matched well with the WHO-defined variants of concern, indicating that our findings from bacteriophages might be more broadly applicable and the ANI-based clusters may represent functionally and/or ecologically distinct units. These units appear to be predominantly driven by (high) ecological cohesiveness coupled to either frequent recombination for bacteriophages or selection and clonal evolution for other viruses such as SARS-CoV-2, indicating that fundamentally different underlying mechanisms could lead to similar diversity patterns. Accordingly, we propose the ANI gap approach outlined above for defining viral intra-species units, for which we propose the term genomovars.

IMPORTANCE

Viral species are composed of an ensemble of intra-species variants whose individual dynamics may have major implications for human and animal health and/or ecosystem functioning. However, the lack of universally accepted standards to define these intra-species variants has led researchers to use different approaches for this task, creating inconsistent intra-species units across different viral families and confusion in communication. By comparing hundreds of mostly bacteriophage genomes, we show that there is a widely distributed natural gap in whole-genome average nucleotide identity values in most, but not all, of these species that can be used to define intra-species units. Therefore, these results advance the molecular toolbox for tracking viral intra-species units and should facilitate future epidemiological and environmental studies.

Unveiling the unknown viral world in groundwater

Viruses as the prevailing biological entities are poorly understood in underground realms. Here, we establish the first metagenomic Groundwater Virome Catalogue (GWVC) comprising 280,420 viral species ( ≥ 5 kb) detected from 607 monitored wells in seven geo-environmental zones throughout China. In expanding ~10-fold the global portfolio of known groundwater viruses, we uncover over 99% novel viruses and about 95% novel viral clusters. By linking viruses to hosts from 119 prokaryotic phyla, we double the number of microbial phyla known to be virus-infected in groundwater. As keystone ultrasmall symbionts in aquifers, CPR bacteria and DPANN archaea are susceptible to virulent viruses. Certain complete CPR viruses even likely infect non-CPR bacteria, while partial CPR/DPANN viruses harbor cell-surface modification genes that assist symbiont cell adhesion to free-living microbes. This study reveals the unknown viral world and auxiliary metabolism associated with methane, nitrogen, sulfur, and phosphorus cycling in groundwater, and highlights the importance of subsurface virosphere in viral ecology.

Global diversity and ecological functions of viruses inhabiting oil reservoirs

Oil reservoirs, being one of the significant subsurface repositories of energy and carbon, host diverse microbial communities affecting energy production and carbon emissions. Viruses play crucial roles in the ecology of microbiomes, however, their distribution and ecological significance in oil reservoirs remain undetermined. Here, we assemble a catalogue encompassing viral and prokaryotic genomes sourced from oil reservoirs. The catalogue comprises 7229 prokaryotic genomes and 3,886 viral Operational Taxonomic Units (vOTUs) from 182 oil reservoir metagenomes. The results show that viruses are widely distributed in oil reservoirs, and 85% vOTUs in oil reservoir are detected in less than 10% of the samples, highlighting the heterogeneous nature of viral communities within oil reservoirs. Through combined microcosm enrichment experiments and bioinformatics analysis, we validate the ecological roles of viruses in regulating the community structure of sulfate reducing microorganisms, primarily through a virulent lifestyle. Taken together, this study uncovers a rich diversity of viruses and their ecological functions within oil reservoirs, offering a comprehensive understanding of the role of viral communities in the biogeochemical cycles of the deep biosphere.

Longitudinal alterations in the urinary virome of kidney transplant recipients are influenced by BK viremia and patient sex

Little is known about the urinary virome and how it interacts with the host, particularly in renal transplant diseases. Using metagenomic sequencing, we characterized the urinary virome of 23 kidney transplant recipients longitudinally (11 BKV+ patients and 12 BKV- patients). We applied linear mixed effects models, PERMANOVA, k-means clustering, and MaAsLin2 algorithms to determine virome signatures associated with post-transplant time, BK viremia status, and patient sex. We found that the richness and alpha diversity of urinary virome were significantly different in renal transplant recipients with BKV+ over time in comparison to BKV- (richness P = 0.012, alpha P < 0.0001). Female BKV- patients had significantly higher virome richness than males (P = 0.0063). Virome beta diversity was significantly different between patients by BKV status (P < 0.001). Additionally, we identified underlying interactions between patient sex and BKV status, in terms of virome beta diversity (P = 0.008). BK polyomavirus infections were primarily of subtypes IA, IB1, and IB2. The non-BK dominant samples clustered into six urinary virome community states. BKV- samples had more anelloviruses than BKV+ samples though this difference was not statistically significant. Lastly, we identified specific viruses, associated with BKV+ and time in our samples. Our results indicate that dynamic alterations in the urinary virome over the post-transplant period in kidney transplant recipients can be shaped by BK viremia and patient sex. These findings advance our fundamental understanding of the urinary virome and support a new line of investigation in renal disease and transplantation.

IMPORTANCE

The urinary microbiome is increasingly implicated in renal health and disease. While most research focuses on bacteria communities of the microbiome, factors that influence the urinary virome are not understood. Here, we investigated the urinary virome of 23 adult kidney transplant recipients longitudinally over 14 weeks post-transplant. We show that alterations in the urinary virome are associated with kidney transplant recipients with BK polyomavirus viremia that can lead to BK nephropathy and allograft rejection. By modeling the temporal dynamics post-transplant, we delineated specific profiles of the urinary virome associated with patient sex and urinary community states. These findings reveal fundamental aspects of the urinary virome that can be leveraged to better manage kidney diseases.

Viral Communities Suppress the Earthworm Gut Antibiotic Resistome by Lysing Bacteria on a National Scale

Earthworms are critical in regulating soil processes and act as filters for antibiotic resistance genes (ARGs). Yet, the geographic patterns and main drivers of earthworm gut ARGs remain largely unknown. We collected 52 earthworm and soil samples from arable and forest ecosystems along a 3000 km transect across China, analyzing the diversity and abundance of ARGs using shotgun metagenomics. Earthworm guts harbored a lower diversity and abundance of ARGs compared to soil, resulting in a stronger distance-decay rate of ARGs in the gut. Greater deterministic assembly processes of ARGs were found in the gut than in soil. The earthworm gut had a lower frequency of co-occurrence patterns between ARGs and mobile genetic elements (MGEs) in forest than in arable systems. Viral diversity was higher in the gut compared to soil and was negatively correlated with bacterial diversity. Bacteria such as Streptomyces and Pseudomonas were potential hosts of both viruses and ARGs. Viruses had negative effects on the diversity and abundance of ARGs, likely due to the lysis on ARG-bearing bacteria. These findings provide new insights into the variations of ARGs in the earthworm gut and highlight the vital role of viruses in the regulation of ARGs in the soil ecosystem.

Full-genome sequencing of dozens of new DNA viruses found in Spanish bat feces

Bats are natural hosts of multiple viruses, many of which have clear zoonotic potential. The search for emerging viruses has been aided by the implementation of metagenomic tools, which have also enabled the detection of unprecedented viral diversity. Currently, this search is mainly focused on RNA viruses, which are largely over-represented in databases. To compensate for this research bias, we analyzed fecal samples from 189 Spanish bats belonging to 22 different species using viral metagenomics. This allowed us to identify 52 complete or near-complete viral genomes belonging to the families Adenoviridae, Circoviridae, Genomoviridae, Papillomaviridae, Parvoviridae, Polyomaviridae and Smacoviridae. Of these, 30 could constitute new species, doubling the number of viruses currently described in Europe. These findings open the door to a more thorough analysis of bat DNA viruses and their zoonotic potential.

IMPORTANCE

Metagenomics has become a fundamental tool to characterize the global virosphere, allowing us not only to understand the existing viral diversity and its ecological implications but also to identify new and emerging viruses. RNA viruses have a higher zoonotic potential, but this risk is also present for some DNA virus families. In our study, we analyzed the DNA fraction of fecal samples from 22 Spanish bat species, identifying 52 complete or near-complete genomes of different viral families with zoonotic potential. This doubles the number of genomes currently described in Europe. Metagenomic data often produce partial genomes that can be difficult to analyze. Our work, however, has characterized a large number of complete genomes, thus facilitating their taxonomic classification and enabling different analyses to be carried out to evaluate their zoonotic potential. For example, recombination studies are relevant since this phenomenon could play a major role in cross-species transmission.

SegVir: Reconstruction of Complete Segmented RNA Viral Genomes from Metatranscriptomes

Segmented RNA viruses are a complex group of RNA viruses with multisegment genomes. Reconstructing complete segmented viruses is crucial for advancing our understanding of viral diversity, evolution, and public health impact. Using metatranscriptomic data to identify known and novel segmented viruses has sped up the survey of segmented viruses in various ecosystems. However, the high genetic diversity and the difficulty in binning complete segmented genomes present significant challenges in segmented virus reconstruction. Current virus detection tools are primarily used to identify nonsegmented viral genomes. This study presents SegVir, a novel tool designed to identify segmented RNA viruses and reconstruct their complete genomes from complex metatranscriptomes. SegVir leverages both close and remote homology searches to accurately detect conserved and divergent viral segments. Additionally, we introduce a new method that can evaluate the genome completeness and conservation based on gene content. Our evaluations on simulated datasets demonstrate SegVir's superior sensitivity and precision compared to existing tools. Moreover, in experiments using real data, we identified some virus segments missing in the NCBI database, underscoring SegVir's potential to enhance viral metagenome analysis. The source code and supporting data of SegVir are available via https://github.com/HubertTang/SegVir.

Bacterial evolution in Biofiltration of drinking water treatment plant: Different response of phage and plasmid to varied water sources

Biofiltration in drinking water treatment (BDWT) are popular as it holds promise as an alternative to chemical treatments, yet our understanding of the key drivers and trends underlying bacterial evolution within this process remains limited. While plasmids and phages are recognized as the main vectors of horizontal gene transfer (HGT), their roles in shaping bacterial evolution in BDWT remain largely unknown. Here we leverage global metagenomic data to unravel the primary forces driving bacterial evolution in BDWT. Our results revealed that the primary vector of HGT varies depending on the type of source water (groundwater and surface water). Both plasmids and phages accelerated bacterial evolution in BDWT by enhancing genetic diversity within species, but they drove contrasting evolutionary trends in functional redundancy in different source water types. Specifically, trends towards and away from functional redundancy (indicated as gene-protein ratio) were observed in surface-water and groundwater biofilters, respectively. Virulent phages drove bacterial evolution through synergistic interactions with bacterial species capable of natural transformation and with certain natural compounds that disrupt bacterial cytoplasmic membranes. Genes relating to water purification (such as Mn(II)-oxidizing genes), microbial risks (antibiotic resistance genes), and chemical risk (polycyclic aromatic hydrocarbons) were enriched via HGT in BDWT, highlighting the necessity for heighted focus on these useful and risky objects. Overall, these discoveries enhance our understanding of bacterial evolution in BDWT and have implications for the optimization of water treatment strategies.

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 .

Pathogen dynamics and discovery of novel viruses and enzymes by deep nucleic acid sequencing of wastewater

Wastewater contains an extensive reservoir of genetic information, yet largely unexplored. Here, we analyzed by high-throughput sequencing total nucleic acids extracted from wastewater samples collected during a 17 month-period in Berlin, Germany. By integrating global wastewater datasets and applying a novel computational approach to accurately identify viral strains within sewage RNA-sequencing data, we demonstrated the emergence and global dissemination of a specific astrovirus strain. Astrovirus abundance and sequence variation mirrored temporal and spatial patterns of infection, potentially serving as footprints of specific timeframes and geographical locations. Additionally, we revealed more than 100,000 sequence contigs likely originating from novel viral species, exhibiting distinct profiles in total RNA and DNA datasets and including undescribed bunyaviruses and parvoviruses. Finally, we identified thousands of new CRISPR-associated protein sequences, including Transposase B (TnpB), a class of compact, RNA-guided DNA editing enzymes. Collectively, our findings underscore the potential of high-throughput sequencing of total nucleic acids derived from wastewater for a broad range of applications.

The assembly and ecological roles of biofilms attached to plastic debris of Ashmore reef

Plastic pollution in the ocean is a global environmental hazard aggravated by poor management of plastic waste and growth of annual plastic consumption. Microbial communities colonizing the plastic's surface, the plastisphere, has gained global interest resulting in numerous efforts to characterize the plastisphere. However, there are insufficient studies deciphering the underlying metabolic processes governing the function of the plastisphere and the plastic they reside upon. Here, we collected plastic and seawater samples from Ashmore Reef in Australia to examine the planktonic microbes and plastic associated biofilm (PAB) to investigate the ecological impact, pathogenic potential, and plastic degradation capabilities of PAB in Ashmore Reef, as well as the role and impact of bacteriophages on PAB. Using high-throughput metagenomic sequencing, we demonstrated distinct microbial communities between seawater and PAB. Similar numbers of pathogenic bacteria were found in both sample types, yet plastic and seawater select for different pathogen populations. Virulence Factor analysis further illustrated stronger pathogenic potential in PAB, highlighting the pathogenicity of environmental PAB. Furthermore, functional analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathways revealed xenobiotic degradation and fatty acid degradation to be enriched in PABs. In addition, construction of metagenome-assembled genomes (MAG) and functional analysis further demonstrated the presence of a complete Polyethylene (PE) degradation pathway in multiple Proteobacteria MAGs, especially in Rhodobacteriaceae sp. Additionally, we identified viral population presence in PAB, revealing the key role of bacteriophages in shaping these communities within the PAB. Our result provides a comprehensive overview of the various ecological processes shaping microbial community on marine plastic debris.

Hijackers, hitchhikers, or co-drivers? The mysteries of mobilizable genetic elements

Mobile genetic elements shape microbial gene repertoires and populations. Recent results reveal that many, possibly most, microbial mobile genetic elements require helpers to transfer between genomes, which we refer to as Hitcher Genetic Elements (hitchers or HGEs). They may be a large fraction of pathogenicity and resistance genomic islands, whose mechanisms of transfer have remained enigmatic for decades. Together with their helper elements and their bacterial hosts, hitchers form tripartite networks of interactions that evolve rapidly within a parasitism-mutualism continuum. In this emerging view of microbial genomes as communities of mobile genetic elements many questions arise. Which elements are being moved, by whom, and how? How often are hitchers costly hyper-parasites or beneficial mutualists? What is the evolutionary origin of hitchers? Are there key advantages associated with hitchers' lifestyle that justify their unexpected abundance? And why are hitchers systematically smaller than their helpers? In this essay, we start answering these questions and point ways ahead for understanding the principles, origin, mechanisms, and impact of hitchers in bacterial ecology and evolution.

Rheumatoid arthritis patients harbour aberrant enteric bacteriophages with autoimmunity-provoking potential: a paired sibling study

OBJECTIVES

Viruses have been considered as important participants in the development of rheumatoid arthritis (RA). However, the profile of enteric virome and its role in RA remains elusive. This study aimed to investigate the atlas and involvement of virome in RA pathogenesis.

METHODS

Faecal samples from 30 pairs of RA and healthy siblings that minimise genetic interferences were collected for metagenomic sequencing. The α and β diversity of the virome and the virome-bacteriome interaction were analysed. The differential bacteriophages were identified, and their correlations with clinical and immunological features of RA were analysed. The potential involvement of these differential bacteriophages in RA pathogenesis was further investigated by auxiliary metabolic gene annotation and molecular mimicry study. The responses of CD4+ T cells and B cells to the mimotopes derived from the differential bacteriophages were systemically studied.

RESULTS

The composition of the enteric bacteriophageome was distorted in RA. The differentially presented bacteriophages correlated with the immunological features of RA, including anti-CCP autoantibody and HLA-DR shared epitope. Intriguingly, the glycerolipid and purine metabolic genes were highly active in the bacteriophages from RA. Moreover, peptides of RA-enriched phages, in particular Prevotella phage and Oscillibacter phage could provoke the autoimmune responses in CD4+ T cells and plasma cells via molecular mimicry of the disease-associated autoantigen epitopes, especially those of Bip.

CONCLUSIONS

This study provides new insights into enteric bacteriophageome in RA development. In particular, the aberrant bacteriophages demonstrated autoimmunity-provoking potential that would promote the occurrence of the disease.

Genomic insights into redox-driven microbial processes for carbon decomposition in thawing Arctic soils and permafrost

Climate change is rapidly transforming Arctic landscapes where increasing soil temperatures speed up permafrost thaw. This exposes large carbon stocks to microbial decomposition, possibly worsening climate change by releasing more greenhouse gases. Understanding how microbes break down soil carbon, especially under the anaerobic conditions of thawing permafrost, is important to determine future changes. Here, we studied the microbial community dynamics and soil carbon decomposition potential in permafrost and active layer soils under anaerobic laboratory conditions that simulated an Arctic summer thaw. The microbial and viral compositions in the samples were analyzed based on metagenomes, metagenome-assembled genomes, and metagenomic viral contigs (mVCs). Following the thawing of permafrost, there was a notable shift in microbial community structure, with fermentative Firmicutes and Bacteroidota taking over from Actinobacteria and Proteobacteria over the 60-day incubation period. The increase in iron and sulfate-reducing microbes had a significant role in limiting methane production from thawed permafrost, underscoring the competition within microbial communities. We explored the growth strategies of microbial communities and found that slow growth was the major strategy in both the active layer and permafrost. Our findings challenge the assumption that fast-growing microbes mainly respond to environmental changes like permafrost thaw. Instead, they indicate a common strategy of slow growth among microbial communities, likely due to the thermodynamic constraints of soil substrates and electron acceptors, and the need for microbes to adjust to post-thaw conditions. The mVCs harbored a wide range of auxiliary metabolic genes that may support cell protection from ice formation in virus-infected cells.

IMPORTANCE

As the Arctic warms, thawing permafrost unlocks carbon, potentially accelerating climate change by releasing greenhouse gases. Our research delves into the underlying biogeochemical processes likely mediated by the soil microbial community in response to the wet and anaerobic conditions, akin to an Arctic summer thaw. We observed a significant shift in the microbial community post-thaw, with fermentative bacteria like Firmicutes and Bacteroidota taking over and switching to different fermentation pathways. The dominance of iron and sulfate-reducing bacteria likely constrained methane production in the thawing permafrost. Slow-growing microbes outweighed fast-growing ones, even after thaw, upending the expectation that rapid microbial responses to dominate after permafrost thaws. This research highlights the nuanced and complex interactions within Arctic soil microbial communities and underscores the challenges in predicting microbial response to environmental change.

VOGDB-Database of Virus Orthologous Groups

Computational models of homologous protein groups are essential in sequence bioinformatics. Due to the diversity and rapid evolution of viruses, the grouping of protein sequences from virus genomes is particularly challenging. The low sequence similarities of homologous genes in viruses require specific approaches for sequence- and structure-based clustering. Furthermore, the annotation of virus genomes in public databases is not as consistent and up to date as for many cellular genomes. To tackle these problems, we have developed VOGDB, which is a database of virus orthologous groups. VOGDB is a multi-layer database that progressively groups viral genes into groups connected by increasingly remote similarity. The first layer is based on pair-wise sequence similarities, the second layer is based on the sequence profile alignments, and the third layer uses predicted protein structures to find the most remote similarity. VOGDB groups allow for more sensitive homology searches of novel genes and increase the chance of predicting annotations or inferring phylogeny. VOGD B uses all virus genomes from RefSeq and partially reannotates them. VOGDB is updated with every RefSeq release. The unique feature of VOGDB is the inclusion of both prokaryotic and eukaryotic viruses in the same clustering process, which makes it possible to explore old evolutionary relationships of the two groups. VOGDB is freely available at vogdb.org under the CC BY 4.0 license.

Extremely halophilic brine community manipulation shows higher robustness of microbiomes inhabiting human-driven solar saltern than naturally driven lake

Hypersaline ecosystems display taxonomically similar assemblages with low diversities and highly dense accompanying viromes. The ecological implications of viral infection on natural microbial populations remain poorly understood, especially at finer scales of diversity. Here, we sought to investigate the influence of changes in environmental physicochemical conditions and viral predation pressure by autochthonous and allochthonous viruses on host dynamics. For this purpose, we transplanted two microbiomes coming from distant hypersaline systems (solar salterns of Es Trenc in Spain and the thalassohaline lake of Aran-Bidgol lake in Iran), by exchanging the cellular fractions with the sterile-filtered accompanying brines with and without the free extracellular virus fraction. The midterm exposure (1 month) of the microbiomes to the new conditions showed that at the supraspecific taxonomic range, the assemblies from the solar saltern brine more strongly resisted the environmental changes and viral predation than that of the lake. The metagenome-assembled genomes (MAGs) analysis revealed an intraspecific transition at the ecotype level, mainly driven by changes in viral predation pressure, by both autochthonous and allochthonous viruses.

IMPORTANCE

Viruses greatly influence succession and diversification of their hosts, yet the effects of viral infection on the ecological dynamics of natural microbial populations remain poorly understood, especially at finer scales of diversity. By manipulating the viral predation pressure by autochthonous and allochthonous viruses, we uncovered potential phage-host interaction, and their important role in structuring the prokaryote community at an ecotype level.

Viral and Bacterial Community Dynamics in Food Waste and Digestate from Full-Scale Biogas Plants

Anaerobic digestion (AD) is commonly used in food waste treatment. Prokaryotic microbial communities in AD of food waste have been comprehensively studied. The role of viruses, known to affect microbial dynamics and metabolism, remains largely unexplored. This study employed metagenomic analysis and recovered 967 high-quality viral bins within food waste and digestate derived from 8 full-scale biogas plants. The diversity of viral communities was higher in digestate. In silico predictions linked 20.8% of viruses to microbial host populations, highlighting possible virus predators of key functional microbes. Lineage-specific virus-host ratio varied, indicating that viral infection dynamics might differentially affect microbial responses to the varying process parameters. Evidence for virus-mediated gene transfer was identified, emphasizing the potential role of viruses in controlling the microbiome. AD altered the specific process parameters, potentially promoting a shift in viral lifestyle from lysogenic to lytic. Viruses encoding auxiliary metabolic genes (AMGs) were involved in microbial carbon and nutrient cycling, and most AMGs were transcriptionally expressed in digestate, meaning that viruses with active functional states were likely actively involved in AD. These findings provided a comprehensive profile of viral and bacterial communities and expanded knowledge of the interactions between viruses and hosts in food waste and digestate.

Dynamics of the viral community on the surface of a French smear-ripened cheese during maturation and persistence across production years

The surface of smear-ripened cheeses constitutes a dynamic microbial ecosystem resulting from the successive development of different microbial groups such as lactic acid bacteria, fungi, and ripening bacteria. Recent studies indicate that a viral community, mainly composed of bacteriophages, also represents a common and substantial part of the cheese microbiome. However, the composition of this community, its temporal variations, and associations between bacteriophages and their hosts remain poorly characterized. Here, we studied a French smear-ripened cheese by both viral metagenomics and 16S metabarcoding approaches to assess both the succession of phages and bacterial communities on the cheese surface during cheese ripening and their temporal variations in ready-to-eat cheeses over the years of production. We observed a clear transition of the phage community structure during ripening with a decreased relative abundance of viral species (vOTUs) associated with Lactococcus phages, which were replaced by vOTUs associated with phages infecting ripening bacteria such as Brevibacterium, Glutamicibacter, Pseudoalteromonas, and Vibrio. The dynamics of the phage community was strongly associated with bacterial successions observed on the cheese surface. Finally, while some variations in the distribution of phages were observed in ready-to-eat cheeses produced at different dates spanning more than 4 years of production, the most abundant phages were detected throughout. This result revealed the long-term persistence of the dominant phages in the cheese production environment. Together, these findings offer novel perspectives on the ecology of bacteriophages in smear-ripened cheese and emphasize the significance of incorporating bacteriophages in the microbial ecology studies of fermented foods.IMPORTANCEThe succession of diverse microbial populations is critical for ensuring the production of high-quality cheese. We observed a temporal succession of phages on the surface of a smear-ripened cheese, with new phage communities showing up when ripening bacteria start covering this surface. Interestingly, the final phage community of this cheese is also consistent over large periods of time, as the same bacteriophages were found in cheese products from the same manufacturer made over 4 years. This research highlights the importance of considering these bacteriophages when studying the microbial life of fermented foods like cheese.

Distributions, interactions, and dynamics of prokaryotes and phages in a hybrid biological wastewater treatment system

BACKGROUND

Understanding the interactions and dynamics of microbiotas within biological wastewater treatment systems is essential for ensuring their stability and long-term sustainability. In this study, we developed a systematic framework employing multi-omics and Hi-C sequencing to extensively investigate prokaryotic and phage communities within a hybrid biofilm and activated sludge system.

RESULTS

We uncovered distinct distribution patterns, metabolic capabilities, and activities of functional prokaryotes through the analysis of 454 reconstructed prokaryotic genomes. Additionally, we reconstructed a phage catalog comprising 18,645 viral operational taxonomic units (vOTUs) with high length and contiguity using hybrid assembly, and a distinct distribution of phages was depicted between activated sludge (AS) and biofilm. Importantly, 1340 host-phage pairs were established using Hi-C and conventional in silico methods, unveiling the host-determined phage prevalence. The majority of predicted hosts were found to be involved in various crucial metabolic processes, highlighting the potential vital roles of phages in influencing substance metabolism within this system. Moreover, auxiliary metabolic genes (AMGs) related to various categories (e.g., carbohydrate degradation, sulfur metabolism, transporter) were predicted. Subsequent activity analysis emphasized their potential ability to mediate host metabolism during infection. We also profiled the temporal dynamics of phages and their associated hosts using 13-month time-series metagenomic data, further demonstrating their tight interactions. Notably, we observed lineage-specific infection patterns, such as potentially host abundance- or phage/host ratio-driven phage population changes.

CONCLUSIONS

The insights gained from this research contribute to the growing body of knowledge surrounding interactions and dynamics of host-phage and pave the way for further exploration and potential applications in the field of microbial ecology. Video Abstract.