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

Significant impact of bleaching treatment on phage-host interaction dynamics in a full-scale wastewater treatment plant

The temporal dynamics of phage-host interactions within full-scale biological wastewater treatment (BWT) plants remain inadequately characterized. Here, we provide an in-depth investigation of viral and bacterial dynamics over a nine-year period in an activated sludge BWT plant, where bleach addition was applied to control sludge foaming. By conducting bioinformatic analyses on 98 metagenomic time-series samples, we reconstructed 3,486 bacterial genomes and 2,435 complete or near-complete viral genomes, which were classified into 361 bacterial and 889 viral clusters, respectively. Our results demonstrate that the primary bleaching event induced significant shifts in both bacterial and viral communities, as well as in virus-host interactions, as evidenced by alterations in bacteria-virus interaction networks and virus-to-host ratio dynamics. Following bleaching, the bacteria-virus network became less interconnected but more compartmentalized. Viral communities mirrored bacterial dynamics, indicating a strong coupling in phage-host interactions. Among the identified virus-host pairs, many exhibited a decelerating rise in viral abundance relative to host abundance, with virus-to-host ratios generally displaying a negative correlation with host abundance. This trend was particularly pronounced in virus-host pairs where viruses harbored integrase genes, indicative of temperate dynamics resembling a "Piggyback-the-Winner" model. Notably, the bleaching intervention appeared to induce a transition from lysogeny to lysis in viruses associated with some foaming-related bacterial species, suggesting a potential virus-involved indirect mechanism by which bleaching mitigates sludge foaming.

Viruses in human-impacted estuarine ecotones: Distribution, metabolic potential, and environmental risks

Estuaries, as dynamic ecological interfaces between marine and terrestrial systems, are characterized by high productivity and intricate microbial communities. Viruses exert critical regulatory effects on microbial processes, influencing ecological functions and contributing to environmental dynamics in estuarine ecosystems. Despite their significance, the diversity and ecological roles of estuarine viruses remain insufficiently understood. This study explored the viral biogeographic patterns, metabolic potential, and influencing factors in 30 subtropical estuaries in China. Few estuarine viruses (< 22 %) exhibited homology with known viruses, and the low overlap of virus clusters with other environments highlights their novelty and habitat specificity. Mantel tests and random forest analysis identified salinity, temperature, nutrients, and pollutants as key factors influencing viral composition and functional profiles. In addition, correlation analysis between virus and host confirmed significant virus-host interactions, while functional analyses highlighted the role of environmental conditions and horizontal gene transfer in shaping auxiliary metabolic genes linked to elemental biogeochemical cycles, particularly phosphorus, sulfur, and nitrogen. The detection of antibiotic resistance genes (ARGs) and virulence factors (VFs) within viral genomes underscores the role of viruses as reservoirs of ARGs and VFs in these ecosystems. These results demonstrate the profound influence of abiotic and host factors on viral community structures in subtropical estuarine ecotones and underscore the ecological significance of metabolic genes in biogeochemical cycling. By clarifying these interactions, this study advances the understanding of viral contributions to ecosystem functioning and biogeochemical dynamics in estuarine environments.

Phages-bacteria interactions underlying the dynamics of polyhydroxyalkanoate-producing mixed microbial cultures via meta-omics study

The dynamics of the structure of polyhydroxyalkanoate-producing mixed microbial cultures (PHA-MMCs) during enrichment and maintenance is an unsolved problem. The effect of phages has been proposed as a cause of dynamic changes in community structure, but evidence is lacking. To address this question, five PHA-MMCs were enriched, and biological samples were sampled temporally to study the interactions between phage and bacterial members by combining metagenomics and metatranscriptomics. A total of 963 metagenome-assembled genomes (MAGs) and 4,294 phage operational taxonomic units (pOTUs) were assembled from bulk metagenomic data. The dynamic changes in the structure of phage and bacterial communities were remarkably consistent. Structural equation modeling analysis showed that phages could infect and lyse dominant species to vacate ecological niches for other species, resulting in a community succession state in which dominant species alternated. Seven key auxiliary metabolic genes (AMGs), phaC, fadJ, acs, ackA, phbB, acdAB, and fadD, potentially contributing to PHA synthesis were identified from phage sequences. Importantly, these AMGs were transcribed, indicating that they were in an active expression state. The meta-analysis provides the first catalog of phages in PHA-MMCs and the AMGs they carry, as well as how they affect the dynamic changes in bacterial communities. This study provides a reference for subsequent studies on understanding and regulating the microbial community structure of open microbial systems.IMPORTANCEThe synthesis of biodegradable plastic PHA from organic waste through mixed microbial cultures (PHA-MMCs), at extremely low cost, has the potential for expanded production. However, the dynamics of dominant species in PHA-MMCs are poorly understood. Our results demonstrate for the first time the impact of phages on the structure of bacterial communities in the PHA-MMCs. There are complex interactions between the PHA producers (e.g., Azomonas, Paracoccus, and Thauera) and phages (e.g., Casadabanvirus and unclassified Hendrixvirinae). Phage communities can regulate the activity and structure of bacterial communities. In addition, the AMGs related to PHA synthesis may hitchhike during phage-host infection cycles, enabling their dissemination across bacterial communities, and phages may act as a critical genetic reservoir for bacterial members, facilitating access to PHA synthesis-related functional traits. This study highlights the impact of phages on bacterial community structure, suggesting that phages have the potential to be used as a tool for better controlling the microbial community structure of PHA-MMCs.

The potential role of viruses in antibiotic resistance gene dissemination in activated sludge viromes

The dissemination of antibiotic resistance genes (ARGs) in activated sludge (AS) systems poses significant environmental and public health challenges. The role of viruses, primarily bacteriophages, in storing and spreading ARGs in AS systems remains largely unexplored. This study characterized the viral community, virus-associated ARGs (vir_ARGs), and mobile genetic elements (MGEs) of aerobic AS viromes from eight wastewater treatment plants (WWTPs) in eastern China. 78,604 viral operational taxonomic units (vOTUs) were identified, including 1685 temperate vOTUs (T-vOTUs). Five ARG types were detected in 37 vOTUs, indicating a low proportion of ARG-carrying viruses. The co-occurrence rate between vir_ARGs and MGEs was 37.83 %, and six ARG-carrying vOTUs contained multiple MGEs, indicating a transfer potential of vir_ARGs. Additionally, ARG and MGE profiles of AS metagenomes were analyzed to evaluate the transfer potential of phage activity on ARGs. The results showed that phage-associated MGEs showed a significant coupling with both the abundance and composition of ARGs, suggesting a potential role of phages in ARG propagation. These findings offer preliminary insights into understanding the viral resistome and its transfer potential in AS systems. Future research necessitates rigorous pure culture and molecular biology experiments to elucidate the precise mechanisms through which viruses contribute to the dissemination and persistence of ARGs.

Electrode functional microorganisms in bioelectrochemical systems and its regulation: A review

Bioelectrochemical systems (BES) as environmental remediation biotechnologies have boomed in the last two decades. Although BESs combined technologies with electro-chemistry, -biology, and -physics, microorganisms and biofilms remain at their core. In this review, various functional microorganisms in BESs for CO2 reduction, dehalogenation, nitrate, phosphate, and sulfate reduction, metal removal, and volatile organic compound oxidation are summarized and compared in detail. Moreover, interrelationship regulation approaches for functional microorganisms and methods for electroactive biofilm development, such as targeted electrode surface modification, chemical treatment, physical revealing, biological optimization, and genetic programming are pointed out. This review provides promising guidance and suggestions for the selection of microbial inoculants and provides an analysis of the role of individual microorganisms in mixed microbial communities and its metabolisms.

Towards a unifying phylogenomic framework for tailed phages

Classifying viruses systematically has remained a key challenge of virology due to the absence of universal genes and vast genetic diversity of viruses. In particular, the most dominant and diverse group of viruses, the tailed double-stranded DNA viruses of prokaryotes belonging to the class Caudoviricetes, lack sufficient similarity in the genetic machinery that unifies them to reconstruct an inclusive, stable phylogeny of these genes. While previous approaches to organize tailed phage diversity have managed to distinguish various taxonomic levels, these methods are limited in scalability, reproducibility, and the inclusion of modes of evolution, like gene gains and losses, remain key challenges. Here, we present a novel, comprehensive, and reproducible framework for examining evolutionary relationships of tailed phages. In this framework, we compare phage genomes based on the presence and absence of a fixed set of gene families which are used as binary trait data that is input into maximum likelihood models. Our resulting phylogeny stably recovers known taxonomic families of tailed phages, with and without the inclusion of metagenome-derived phages. We also quantify the mosaicism of replication and structural genes among known families, and our results suggest that these exchanges likely underpin the emergence of new families. Additionally, we apply this framework to large phages (>100 kilobases) to map emergences of traits associated with genome expansion. Taken together, this evolutionary framework for charting and organizing tailed phage diversity improves the systemization of phage taxonomy, which can unify phage studies and advance our understanding of their evolution.

Virus-prokaryote interactions assist pollutant removal in constructed wetlands

As a vital part of microbial communities, viruses in constructed wetlands (CWs) remain poorly explored, yet they could significantly affect pollutant removal. Here, two pilot-scale CWs were built to investigate the viral community under different hydraulic loading rates (HLRs) using in-depth metagenomic analysis. Gene-sharing networks suggested that the CWs were pools of unexplored viruses. A higher abundance of prokaryotic functional genes related to sulfur cycling and denitrification was observed in the higher HLR condition, which was associated with greater removal of total nitrogen and nitrate nitrogen compared to the lower HLR condition. Viruses also affect nitrogen pollutant removal by potentially infecting functional prokaryotes, such as denitrification bacteria and ammonia-oxidizing bacteria, and by providing auxiliary metabolic genes involved in sulfur and nitrogen cycling. These findings reveal the significance of viruses in pollutant removal in CWs and enhance the understanding of the relationship between engineering design parameters and performance from microbial perspectives.

Unveiling the potential role of virus-encoded polyphosphate kinases in enhancing phosphorus removal in activated sludge systems

While microbial phosphate removal in activated sludge (AS) systems has been extensively studied, the role of viruses in this process remains largely unexplored. In this study, we identified 149 viral auxiliary metabolic genes associated with phosphorus cycling from 2,510 viral contigs (VCs) derived from AS systems. Notably, polyphosphate kinase 1 (ppk1) and polyphosphate kinase 2 (ppk2) genes, which are primarily responsible for phosphate removal, were found in five unclassified VCs. These genes exhibited conserved protein structures and active catalytic sites, indicating a pivotal role of viruses in enhancing phosphorus removal. Phylogenetic analysis demonstrated a close relationship between viral ppk genes and their bacterial counterparts, suggesting the occurrence of horizontal gene transfer. Furthermore, experimental assays validated that viral ppk genes enhanced host phosphate removal capabilities. VCs carrying ppk genes were observed across diverse ecological and geographical contexts, suggesting their potential to bolster host functions in varied environmental and nutrient settings, spanning natural and engineered systems. These findings uncover a previously underappreciated mechanism by which viruses enhance phosphate removal in wastewater treatment plants. Overall, our study highlights the potential for leveraging virus-encoded genes to improve the efficiency of biological phosphorus removal processes, offering new insights into the microbial ecology of AS systems and the role of viruses in biogeochemical cycling.

Characteristics of Phages and Their Interactions With Hosts in Anaerobic Reactors

Anaerobic digestion (AD) of organic wastes relies on the interaction and cooperation of various microorganisms. Phages are crucial components of the microbial community in AD systems, but their diversity and interactions with the prokaryotic populations are still inadequately comprehended. In this study, 2121 viral operational taxonomic units (vOTUs) were recovered from 12 anaerobic fatty acid-fed reactors. Notably, 63.1% of these vOTUs could not be assigned to any known family, revealing a substantial presence of uncharted phages specifically associated with AD environments. Over half of the vOTUs associated with hosts had the capability to infect multiple hosts, ranging from 2 to 49, with a prevalent tendency to infect 2-5 hosts. In silico predictions of phage-host linkages uncovered that only a small fraction of vOTUs were shared across different functional groups, including fermentative bacteria, syntrophic fatty acid-oxidising bacteria (SFOB) and methanogens. Phages linked to hosts in all three groups primarily consisted of generalists and temperate species, especially those linked to SFOB. Additionally, metabolic reconstruction identified auxiliary metabolic genes participating in fatty acid degradation, methanogenesis and energy conservation. The present study provides insights into phage characteristics and their in situ interactions with prokaryotic hosts, highlighting their ecological role in AD systems.

High-throughput single-cell sequencing of activated sludge microbiome

Wastewater treatment plants (WWTPs) represent one of biotechnology's largest and most critical applications, playing a pivotal role in environmental protection and public health. In WWTPs, activated sludge (AS) plays a major role in removing contaminants and pathogens from wastewater. While metagenomics has advanced our understanding of microbial communities, it still faces challenges in revealing the genomic heterogeneity of cells, uncovering the microbial dark matter, and establishing precise links between genetic elements and their host cells as a bulk method. These issues could be largely resolved by single-cell sequencing, which can offer unprecedented resolution to show the unique genetic information. Here we show the high-throughput single-cell sequencing to the AS microbiome. The single-amplified genomes (SAGs) of 15,110 individual cells were clustered into 2,454 SAG bins. We find that 27.5% of the genomes in the AS microbial community represent potential novel species, highlighting the presence of microbial dark matter. Furthermore, we identified 1,137 antibiotic resistance genes (ARGs), 10,450 plasmid fragments, and 1,343 phage contigs, with shared plasmid and phage groups broadly distributed among hosts, indicating a high frequency of horizontal gene transfer (HGT) within the AS microbiome. Complementary analysis using 1,529 metagenome-assembled genomes from the AS samples allowed for the taxonomic classification of 98 SAG bins, which were previously unclassified. Our study establishes the feasibility of single-cell sequencing in characterizing the AS microbiome, providing novel insights into its ecological dynamics, and deepening our understanding of HGT processes, particularly those involving ARGs. Additionally, this valuable tool could monitor the distribution, spread, and pathogenic hosts of ARGs both within AS environments and between AS and other environments, which will ultimately contribute to developing a health risk evaluation system for diverse environments within a One Health framework.

Method evaluation for viruses in activated sludge: Concentration, sequencing, and identification

Activated sludge (AS) in wastewater treatment plants is one of the largest artificial microbial ecosystems on earth and it makes enormous contributions to human societies. Viruses are an important component in AS with a high abundance. However, their communities and functionalities have not been as widely explored as those of other microorganisms, such as bacteria. This gap is mainly due to technical challenges in effective viral concentration, extraction, and sequencing. In this study, we compared four kinds of concentration methods, two sequencing approaches, and four identification bioinformatic tools to evaluate the whole analysis workflow for viruses in AS. Results showed flocculation, filtration, and resuspension (FFR) could get the longest DNA lengths and ultracentrifugation obtained the highest DNA yields for viruses in AS. Based on the results of present study, FFR and tangential flow filtration with the membrane pore size of 100 kDa were most recommended to concentrate viruses in AS samples with huge volumes. Besides, different concentration methods could get different viral catalogs and thus multiple methods should be combined to get the whole picture of viruses in the system. In addition, geNomad was the most recommended identification tool for viruses in the present study and the long-read sequencing could improve the assembly statistics of viruses when compared with the short-read sequencing. For the 8192 viral operational taxonomic units in this study, 95.1 % of them were phages and belonged to the same lineage at the order level of Caudovirales. Virulent phages dominated the AS system and Pseudomonadota were the main host. Taken together, this study provides new insights into methods selection for virus research of AS.

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.

Adaptive expression of phage auxiliary metabolic genes in paddy soils and their contribution toward global carbon sequestration

Habitats with intermittent flooding, such as paddy soils, are crucial reservoirs in the global carbon pool; however, the effect of phage-host interactions on the biogeochemical cycling of carbon in paddy soils remains unclear. Hence, this study applied multiomics and global datasets integrated with validation experiments to investigate phage-host community interactions and the potential of phages to impact carbon sequestration in paddy soils. The results demonstrated that paddy soil phages harbor a diverse and abundant repertoire of auxiliary metabolic genes (AMGs) associated with carbon fixation, comprising 23.7% of the identified AMGs. The successful annotation of protein structures and promoters further suggested an elevated expression potential of these genes within their bacterial hosts. Moreover, environmental stressors, such as heavy metal contamination, cause genetic variation in paddy phages and up-regulate the expression of carbon fixation AMGs, as demonstrated by the significant enrichment of related metabolites (P < 0.05). Notably, the findings indicate that lysogenic phages infecting carbon-fixing hosts increased by 10.7% under heavy metal stress. In addition, in situ isotopic labeling experiments induced by mitomycin-C revealed that by increasing heavy metal concentrations, 13CO2 emissions from the treatment with added lysogenic phage decreased by approximately 17.9%. In contrast, 13C-labeled microbial biomass carbon content increased by an average of 35.4% compared to the control. These results suggest that paddy soil phages prominently influence the global carbon cycle, particularly under global change conditions. This research enhances our understanding of phage-host cooperation in driving carbon sequestration in paddy soils amid evolving environmental conditions.

Strain-level multidrug-resistant pathogenic bacteria in urban wastewater treatment plants: Transmission, source tracking and evolution

Wastewater treatment plants (WWTPs) serve as reservoirs for various pathogens and play a pivotal role in safeguarding environmental safety and public health by mitigating pathogen release. Pathogenic bacteria, known for their potential to cause fatal infections, present a significant and emerging threat to global health and remain poorly understood regarding their origins and transmission in the environment. Using metagenomic approaches, we identified a total of 299 pathogens from three full-scale WWTPs. We comprehensively elucidated the occurrence, dissemination, and source tracking of the pathogens across the WWTPs, addressing deficiencies in traditional detection strategies. While indicator pathogens in current wastewater treatment systems such as Escherichia coli are effectively removed, specific drug-resistant pathogens, including Pseudomonas aeruginosa, Pseudomonas putida, and Aeromonas caviae, persist throughout the treatment process, challenging complete eradication efforts. The anoxic section plays a predominant role in controlling abundance but significantly contributes to downstream pathogen diversity. Additionally, evolution throughout the treatment process enhances pathogen diversity, except for upstream transmission, such as A. caviae str. WP8-S18-ESBL-04 and P. aeruginosa PAO1. Our findings highlight the necessity of expanding current biomonitoring indicators for wastewater treatment to optimize treatment strategies and mitigate the potential health risks posed by emerging pathogens. By addressing these research priorities, we can effectively mitigate risks and safeguard environmental safety and public health.

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.

Seasonal succession, host associations, and biochemical roles of aquatic viruses in a eutrophic lake plagued by cyanobacterial blooms

Viruses are implicated to play key roles as biogeochemical mediators and ecological drivers in freshwater ecosystems. However, the dynamics of viruses and host associations throughout the seasons and during blooming periods in eutrophic freshwater ecosystems remain poorly understood. From the water microbiomes of planktonic biomass from Lake Taihu, a large eutrophic freshwater lake in China that experiences annual Microcystis-dominated harmful algal blooms (HABs), we recovered 41,997 unique viral clusters spanning a wide taxonomic range, including 15,139 Caudovirales clusters targeting bacteria and 1,044 NCLDV clusters targeting eukaryotes. The viral community exhibited clear seasonal succession, driven primarily by microbial communities (particularly Cyanobacteria and Planctomycetes) and environmental factors (mainly nutrients and temperature). Host prediction revealed that viral infection had a more distinct impact on bacteria-driven nitrogen pathways than on phosphate cycling. HAB-induced variations in microbial composition and environmental conditions affected viral strategies including viral lifestyles, host range, and virus-encoded auxiliary metabolic genes (vAMGs) distributions. Viruses infecting Proteobacteria and Actinobacteria showed an enhanced lysogenic lifestyle and a narrower host range during HAB peak in summer, while viruses infecting Bacteroidota adopted an opposite strategy. Notably, vAMGs were most abundant before the HAB outbreak in spring, compensating for bacterial metabolic processes of their hosts such as carbohydrates metabolism, photosynthesis, and phosphate regulation. The findings highlight the intricate relationships between viruses, host microbes, and the bloom-associated environment, underscoring the important biochemical roles viruses play in eutrophic freshwater ecosystems.

Differentiation of the Anammox core microbiome: Unraveling the evolutionary impetus of scalable gene flow

Anaerobic ammonium oxidation bacteria (AAOB), distinguished by their unique autotrophic nitrogen metabolism, hold pivotal positions in the global nitrogen cycle and environmental biotechnologies. However, the ecophysiology and evolution of AAOB remain poorly understood, attributed to the absence of monocultures. Hence, a comprehensive elucidation of the AAOB-dominated core microbiome, anammox core, is imperative to further completing the theory of engineered nitrogen removal and ecological roles of anammox. Performing taxonomic and phylogenetic analyses on collected genome repertoires, we show here that Candidatus Brocadia and Candidatus Kuenenia possesses a more compact core than Candidatus Jettenia, which partly explains why the latter has a less common ecological presence. Evidence of gene flow is particularly striking in functions related to biosynthesis and oxygen detoxification, underscoring the evolutionary forces driving lineage and core differentiation. Furthermore, CRISPR spacer traceback of the AAOB metagenome-assembled genomes (MAGs) reveals a series of genetic traces for the concealed phages. By reconceptualizing the functional divergence of AAOB with the historical role of phages, we ultimately propose a coevolutionary framework to understand the evolutionary trajectory of anammox microecology. The discoveries provided in this study offer new insights into understanding the evolution of AAOB and the ecology of anammox.

Imputing Metagenomic Hi-C Contacts Facilitates the Integrative Contig Binning Through Constrained Random Walk with Restart

Metagenomic Hi-C (metaHi-C) has shown remarkable potential for retrieving high-quality metagenome-assembled genomes from complex microbial communities. Nevertheless, existing metaHi-C-based contig binning methods solely rely on Hi-C interactions between contigs, disregarding crucial biological information such as the presence of single-copy marker genes. To overcome this limitation, we introduce ImputeCC, an integrative contig binning tool optimized for metaHi-C datasets. ImputeCC integrates both Hi-C interactions and the discriminative power of single-copy marker genes to group marker-gene-containing contigs into preliminary bins. It also introduces a novel constrained random walk with restart algorithm to enhance Hi-C connectivity among contigs. Comprehensive assessments using both mock and real metaHi-C datasets from diverse environments demonstrate that ImputeCC consistently outperforms other Hi-C-based contig binning tools. A genus-level analysis of the sheep gut microbiota reconstructed by ImputeCC underlines its capability to recover key species from dominant genera and identify previously unknown genera.

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.

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.

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.

Revealing the Viable Microbial Community of Biofilm in a Sewage Treatment System Using Propidium Monoazide Combined with Real-Time PCR and Metagenomics

Microbial community composition, function, and viability are important for biofilm-based sewage treatment technologies. Most studies of microbial communities mainly rely on the total deoxyribonucleic acid (DNA) extracted from the biofilm. However, nucleotide materials released from dead microorganisms may interfere with the analysis of viable microorganisms and their metabolic potential. In this study, we developed a protocol to assess viability as well as viable community composition and function in biofilm in a sewage treatment system using propidium monoazide (PMA) coupled with real-time quantitative polymerase chain reaction (qPCR) and metagenomic technology. The optimal removal of PMA from non-viable cells was achieved by a PMA concentration of 4 μM, incubation in darkness for 5 min, and exposure for 5 min. Simultaneously, the detection limit can reach a viable bacteria proportion of 1%, within the detection concentration range of 102-108 CFU/mL (colony forming unit/mL), showing its effectiveness in removing interference from dead cells. Under the optimal conditions, the result of PMA-metagenomic sequencing revealed that 6.72% to 8.18% of non-viable microorganisms were influenced and the composition and relative abundance of the dominant genera were changed. Overall, this study established a fast, sensitive, and highly specific biofilm viability detection method, which could provide technical support for accurately deciphering the structural composition and function of viable microbial communities in sewage treatment biofilms.

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.

High-throughput fabrication of antimicrobial phage microgels and example applications in food decontamination

Engineered by nature, biological entities are exceptional building blocks for biomaterials. These entities can impart enhanced functionalities on the final material that are otherwise unattainable. However, preserving the bioactive functionalities of these building blocks during the material fabrication process remains a challenge. We describe a high-throughput protocol for the bottom-up self-assembly of highly concentrated phages into microgels while preserving and amplifying their inherent antimicrobial activity and biofunctionality. Each microgel is comprised of half a million cross-linked phages as the sole structural component, self-organized in aligned bundles. We discuss common pitfalls in the preparation procedure and describe optimization processes to ensure the preservation of the biofunctionality of the phage building blocks. This protocol enables the production of an antimicrobial spray containing the manufactured phage microgels, loaded with potent virulent phages that effectively reduced high loads of multidrug-resistant Escherichia coli O157:H7 on red meat and fresh produce. Compared with other microgel preparation methods, our protocol is particularly well suited to biological materials because it is free of organic solvents and heat. Bench-scale preparation of base materials, namely microporous films (the template for casting microgels) and pure concentrated phage suspension, requires 3.5 h and 5 d, respectively. A single production run, that yields over 1,750,000 microgels, ranges from 2 h to 2 d depending on the rate of cross-linking chemistry. We expect that this platform will address bottlenecks associated with shelf-stability, preservation and delivery of phage for antimicrobial applications, expanding the use of phage for prevention and control of bacterial infections and contaminants.

A panoramic view of the virosphere in three wastewater treatment plants by integrating viral-like particle-concentrated and traditional non-concentrated metagenomic approaches

Wastewater biotreatment systems harbor a rich diversity of microorganisms, and the effectiveness of biotreatment systems largely depends on the activity of these microorganisms. Specifically, viruses play a crucial role in altering microbial behavior and metabolic processes throughout their infection phases, an aspect that has recently attracted considerable interest. Two metagenomic approaches, viral-like particle-concentrated (VPC, representing free viral-like particles) and non-concentrated (NC, representing the cellular fraction), were employed to assess their efficacy in revealing virome characteristics, including taxonomy, diversity, host interactions, lifestyle, dynamics, and functional genes across processing units of three wastewater treatment plants (WWTPs). Our findings indicate that each approach offers unique insights into the viral community and functional composition. Their combined use proved effective in elucidating WWTP viromes. We identified nearly 50,000 viral contigs, with Cressdnaviricota and Uroviricota being the predominant phyla in the VPC and NC fractions, respectively. Notably, two pathogenic viral families, Asfarviridae and Adenoviridae, were commonly found in these WWTPs. We also observed significant differences in the viromes of WWTPs processing different types of wastewater. Additionally, various phage-derived auxiliary metabolic genes (AMGs) were active at the RNA level, contributing to the metabolism of the microbial community, particularly in carbon, sulfur, and phosphorus cycling. Moreover, we identified 29 virus-carried antibiotic resistance genes (ARGs) with potential for host transfer, highlighting the role of viruses in spreading ARGs in the environment. Overall, this study provides a detailed and integrated view of the virosphere in three WWTPs through the application of VPC and NC metagenomic approaches. Our findings enhance the understanding of viral communities, offering valuable insights for optimizing the operation and regulation of wastewater treatment systems.

Discarded diversity: Novel megaphages, auxiliary metabolic genes, and virally encoded CRISPR-Cas systems in landfills

Background

Viruses are the most abundant microbial entity on the planet, impacting microbial community structure and ecosystem services. Despite outnumbering Bacteria and Archaea by an order of magnitude, viruses have been comparatively underrepresented in reference databases. Metagenomic examinations have illustrated that viruses of Bacteria and Archaea have been specifically understudied in engineered environments. Here we employed metagenomic and computational biology methods to examine the diversity, host interactions, and genetic systems of viruses predicted from 27 samples taken from three municipal landfills across North America.

Results

We identified numerous viruses that are not represented in reference databases, including the third largest bacteriophage genome identified to date (~678 kbp), and note a cosmopolitan diversity of viruses in landfills that are distinct from viromes in other systems. Host-virus interactions were examined via host CRISPR spacer to viral protospacer mapping which captured hyper-targeted viral populations and six viral populations predicted to infect across multiple phyla. Virally-encoded auxiliary metabolic genes (AMGs) were identified with the potential to augment hosts' methane, sulfur, and contaminant degradation metabolisms, including AMGs not previously reported in literature. CRISPR arrays and CRISPR-Cas systems were identified from predicted viral genomes, including the two largest bacteriophage genomes to contain these genetic features. Some virally encoded Cas effector proteins appear distinct relative to previously reported Cas systems and are interesting targets for potential genome editing tools.

Conclusions

Our observations indicate landfills, as heterogeneous contaminated sites with unique selective pressures, are key locations for diverse viruses and atypical virus-host dynamics.

Phages in different habitats and their ability to carry antibiotic resistance genes

As the most abundant organisms on Earth, phages play a key role in the evolution of bacterial antibiotic resistance. Although previous studies have demonstrated the molecular mechanisms of horizontal gene transfer mediated by mobile genetic elements, our understanding of the intertwined relationships between antibiotic resistance genes (ARGs) and phages is limited. In this study, we analysed 2781 metagenomic samples to reveal the composition and species interactions of phage communities in different habitats as well as their capacity to carry ARGs with health risks. The composition of phage communities varies in different habitats and mainly depends on environmental conditions. Terrestrial habitats display more complex and robust interactions between phages than aquatic and human-associated habitats, resulting in the highest biodiversity of phages. Several types of phages in certain taxa (4.95-7.67%, mainly belonging to Caudoviricetes) have the capacity to carry specific ARGs and display a high potential risk to human health, especially in human-associated habitats. Overall, our results provide insights into the assembly mechanisms of phage communities and their effects on the dissemination of antibiotic resistance.

Temporal Dynamics and Contribution of Phage Community to the Prevalence of Antibiotic Resistance Genes in a Full-Scale Sludge Anaerobic Digestion Plant

Anaerobic digestion (AD) is an important biological resource recovery process, where microorganisms play key roles for material transformation. There has been some knowledge about the prokaryotic community and antibiotic resistance genes (ARGs) in AD, but there has been very limited knowledge of phages. In this study, samples from a full-scale AD plant were collected over 13 months, sequenced, and analyzed for viral and prokaryotic metagenomes. Totally, 3015 viral operational taxonomic units (vOTUs) were detected, mostly assigned to Caudoviricetes. The phage community had faster temporal variation than the prokaryotic community. Warm seasons harbored a higher abundance of both temperate phages and broad host-range phages. Seven ARGs of 6 subtypes were carried by 20 vOTUs, a representative ermT gene was synthesized and expressed, and the resistance activity in the host was examined, confirming the real activity of virus-carried ARGs in the AD process. Some of the ARGs were horizontally transferred between the phage and prokaryotic genomes. However, phage infection was not found to contribute to ARG transfer. This study provided an insight into the ecological patterns of the phage community, confirmed the antibiotic resistance activity of virus-carried ARGs, evaluated the contribution of phages on the ARG prevalence, and laid the foundation for the control strategies of the community and antibiotic resistance in the AD process.

A Metagenomic Investigation of Potential Health Risks and Element Cycling Functions of Bacteria and Viruses in Wastewater Treatment Plants

The concentration of viruses in sewage sludge is significantly higher (10-1000-fold) than that found in natural environments, posing a potential risk for human and animal health. However, the composition of these viruses and their role in the transfer of pathogenic factors, as well as their role in the carbon, nitrogen, and phosphorus cycles remain poorly understood. In this study, we employed a shotgun metagenomic approach to investigate the pathogenic bacteria and viral composition and function in two wastewater treatment plants located on a campus. Our analysis revealed the presence of 1334 amplicon sequence variants (ASVs) across six sludge samples, with 242 ASVs (41.22% of total reads) identified as pathogenic bacteria. Arcobacter was found to be the most dominant pathogen accounting for 6.79% of total reads. The virome analysis identified 613 viral genera with Aorunvirus being the most abundant genus at 41.85%. Approximately 0.66% of these viruses were associated with human and animal diseases. More than 60% of the virome consisted of lytic phages. Host prediction analysis revealed that the phages primarily infected Lactobacillus (37.11%), Streptococcus (21.11%), and Staphylococcus (7.11%). Furthermore, our investigation revealed an abundance of auxiliary metabolic genes (AMGs) involved in carbon, nitrogen, and phosphorus cycling within the virome. We also detected a total of 113 antibiotic resistance genes (ARGs), covering major classes of antibiotics across all samples analyzed. Additionally, our findings indicated the presence of virulence factors including the clpP gene accounting for approximately 4.78%, along with toxin genes such as the RecT gene representing approximately 73.48% of all detected virulence factors and toxin genes among all samples analyzed. This study expands our understanding regarding both pathogenic bacteria and viruses present within sewage sludge while providing valuable insights into their ecological functions.

Analysis and culturing of the prototypic crAssphage reveals a phage-plasmid lifestyle

The prototypic crAssphage (Carjivirus communis) is one of the most abundant, prevalent, and persistent gut bacteriophages, yet it remains uncultured and its lifestyle uncharacterized. For the last decade, crAssphage has escaped plaque-dependent culturing efforts, leading us to investigate alternative lifestyles that might explain its widespread success. Through genomic analyses and culturing, we find that crAssphage uses a phage-plasmid lifestyle to persist extrachromosomally. Plasmid-related genes are more highly expressed than those implicated in phage maintenance. Leveraging this finding, we use a plaque-free culturing approach to measure crAssphage replication in culture with Phocaeicola vulgatus, Phocaeicola dorei, and Bacteroides stercoris, revealing a broad host range. We demonstrate that crAssphage persists with its hosts in culture without causing major cell lysis events or integrating into host chromosomes. The ability to switch between phage and plasmid lifestyles within a wide range of hosts contributes to the prolific nature of crAssphage in the human gut microbiome.

Investigating the viral ecology and contribution to the microbial ecology in full-scale mesophilic anaerobic digesters

In an attempt to assess the diversity of viruses and their potential to modulate the metabolism of functional microorganisms in anaerobic digesters, we collected digestate from three mesophilic anaerobic digesters in full-scale wastewater treatment plants treating real municipal wastewater. The reads were analyzed using bioinformatics algorithms to elucidate viral diversity, identify their potential role in modulating the metabolism of functional microorganisms, and provide essential genomic information for the potential use of virus-mediated treatment in controlling the anaerobic digester microbiome. We found that Siphoviridae was the dominant family in mesophilic anaerobic digesters, followed by Myoviridae and Podoviridae. Lysogeny was prevalent in mesophilic anaerobic digesters as the majority of metagenome-assembled genomes contained at least one viral genome within them. One virus within the genome of an acetoclastic methanogen (Methanothrix soehngenii) was observed with a gene (fwdE) acquired via lateral transfer from hydrogenotrophic methanogens. The virus-mediated acquisition of fwdE gene enables possibility of mixotrophic methanogenesis in Methanothrix soehngenii. This evidence highlighted that lysogeny provides fitness advantage to methanogens in anaerobic digesters by adding flexibility to changing substrates. Similarly, we found auxiliary metabolic genes, such as cellulase and alpha glucosidase, of bacterial origin responsible for sludge hydrolysis in viruses. Additionally, we discovered novel viral genomes and provided genomic information on viruses infecting acidogenic, acetogenic, and pathogenic bacteria that can potentially be used for virus-mediated treatment to deal with the souring problem in anaerobic digesters and remove pathogens from biosolids before land application. Collectively, our study provides a genome-level understanding of virome in conjunction with the microbiome in anaerobic digesters that can be used to optimize the anaerobic digestion process for efficient biogas generation.

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.

Identification of Huge Phages from Wastewater Metagenomes

Huge phages have genomes larger than 200 kilobases, which are particularly interesting for their genetic inventory and evolution. We screened 165 wastewater metagenomes for the presence of viral sequences. After identifying over 600 potential huge phage genomes, we reduced the dataset using manual curation by excluding viral contigs that did not contain viral protein-coding genes or consisted of concatemers of several small phage genomes. This dataset showed seven fully annotated huge phage genomes. The phages grouped into distinct phylogenetic clades, likely forming new genera and families. A phylogenomic analysis between our huge phages and phages with smaller genomes, i.e., less than 200 kb, supported the hypothesis that huge phages have undergone convergent evolution. The genomes contained typical phage protein-coding genes, sequential gene cassettes for metabolic pathways, and complete inventories of tRNA genes covering all standard and rare amino acids. Our study showed a pipeline for huge phage analyses that may lead to new enzymes for therapeutic or biotechnological applications.

Hi-C metagenome sequencing reveals soil phage-host interactions

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

Application of high-throughput sequencing technologies and analytical tools for pathogen detection in urban water systems: Progress and future perspectives

The ubiquitous presence of pathogenic microorganisms, such as viruses, bacteria, fungi, and protozoa, in urban water systems poses a significant risk to public health. The emergence of infectious waterborne diseases mediated by urban water systems has become one of the leading global causes of mortality. However, the detection and monitoring of these pathogenic microorganisms have been limited by the complexity and diversity in the environmental samples. Conventional methods were restricted by long assay time, high benchmarks of identification, and narrow application sceneries. Novel technologies, such as high-throughput sequencing technologies, enable potentially full-spectrum detection of trace pathogenic microorganisms in complex environmental matrices. This review discusses the current state of high-throughput sequencing technologies for identifying pathogenic microorganisms in urban water systems with a concise summary. Furthermore, future perspectives in pathogen research emphasize the need for detection methods with high accuracy and sensitivity, the establishment of precise detection standards and procedures, and the significance of bioinformatics software and platforms. We have compiled a list of pathogens analysis software/platforms/databases that boast robust engines and high accuracy for preference. We highlight the significance of analyses by combining targeted and non-targeted sequencing technologies, short and long reads technologies, sequencing technologies, and bioinformatic tools in pursuing upgraded biosafety in urban water systems.

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

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

Narrow host range phages infect essential bacteria for water purification reactions in groundwater-fed rapid sand filters

Biofiltration is used worldwide to provide safe potable water due to its low energy demand and excellent treatment performance. For instance, in Denmark, over 95% of drinking water is supplied through groundwater-fed rapid sand filters (RSF). Bacteriophages, viruses that infect bacteria, have been shown to shape the taxonomic and functional composition of microbial communities across a range of natural and engineering systems. However, phages in the biofiltration systems are rarely studied, despite the central role microbes play in water purification. To probe this, metagenomic data from surface water, groundwater and mixed source water biofiltration units (n = 26 from China, Europe and USA) for drinking water production were analysed to characterize prokaryotic viruses and to identify their potential microbial hosts. The source water type and geographical location are found to exert influence on the composition of the phageome in biofilters. Although the viral abundance (71,676 ± 17,841 RPKM) in biofilters is only 14.4% and 17.0% lower than those of the nutrient-rich wastewater treatment plants and fresh surface waters, the richness (1,441 ± 1,046) and diversity (Inverse Simpson: 91 ± 61) in biofiltration units are significantly less by a factor of 2-5 and 3-4, respectively. In depth analysis of data from 24 groundwater-fed RSFs in Denmark revealed a core phageome shared by most RSFs, which was consistently linked to dominant microbial hosts involved in key biological reactions for water purification. Finally, the high number of specific links detected between phages and bacterial species and the large proportion of lytic phages (77%) led to the conjecture that phages regulate bacterial populations through predation, preventing the proliferation of dominant species and contributing to the established functional redundancy among the dominant microbial groups. In conclusion, bacteriophages are likely to play a significant role in water treatment within biofilters, particularly through interactions with key bacterial species.

Metagenomic analysis unveils the underexplored roles of prokaryotic viruses in a full-scale landfill leachate treatment plant

Enormous viral populations have been identified in activated sludge systems, but their ecological and biochemical roles in landfill leachate treatment plants remain poorly understood. To address this knowledge gap, we conducted an in-depth analysis using 36 metagenomic datasets that we collected and sequenced during a half-year time-series sampling campaign at six sites in a full-scale landfill leachate treatment plant (LLTP), elucidating viral distribution, virus‒host dynamics, virus-encoded auxiliary metabolic genes (AMGs), and viral contributions to the spread of virulence and antibiotic resistance genes. Our findings demonstrated that viral and prokaryotic communities differed widely among different treatment units, with stability over time. LLTP viruses were linked to various prokaryotic hosts, spanning 35 bacterial phyla and one archaeal phylum, which included the core microbes involved in biological treatments, as well as some of the less well-characterized microbial dark matter phyla. By encoding 2364 auxiliary metabolic genes (AMGs), viruses harbored the potential to regulate microbial nucleotide metabolism, facilitate the biodegradation of complex organic matter, and enhance flocculation and settling in biological treatment plants. The abundance distribution of AMGs varied considerably across treatment units and showed a lifestyle-dependent pattern with temperate virus-associated AMGs exhibiting a higher average abundance in downstream biological treatment units and effluent water. Meanwhile, temperate viruses tended to carry a higher load of virulence factor genes (VFGs), antibiotic resistance genes (ARGs), and biotic and metal resistance genes (BMRGs), and engaged in more frequent gene exchanges with prokaryotes than lytic viruses, thus acting as a pivotal contributor to the dissemination of pathogenicity and resistance genes in downstream LLTP units. This study provided a comprehensive profile of viral and prokaryotic communities in the LLTP and unveiled the varying roles of different-lifestyle viruses in biochemical processes and water quality safety.

MetaCC allows scalable and integrative analyses of both long-read and short-read metagenomic Hi-C data

Metagenomic Hi-C (metaHi-C) can identify contig-to-contig relationships with respect to their proximity within the same physical cell. Shotgun libraries in metaHi-C experiments can be constructed by next-generation sequencing (short-read metaHi-C) or more recent third-generation sequencing (long-read metaHi-C). However, all existing metaHi-C analysis methods are developed and benchmarked on short-read metaHi-C datasets and there exists much room for improvement in terms of more scalable and stable analyses, especially for long-read metaHi-C data. Here we report MetaCC, an efficient and integrative framework for analyzing both short-read and long-read metaHi-C datasets. MetaCC outperforms existing methods on normalization and binning. In particular, the MetaCC normalization module, named NormCC, is more than 3000 times faster than the current state-of-the-art method HiCzin on a complex wastewater dataset. When applied to one sheep gut long-read metaHi-C dataset, MetaCC binning module can retrieve 709 high-quality genomes with the largest species diversity using one single sample, including an expansion of five uncultured members from the order Erysipelotrichales, and is the only binner that can recover the genome of one important species Bacteroides vulgatus. Further plasmid analyses reveal that MetaCC binning is able to capture multi-copy plasmids.

Presence and role of viruses in anaerobic digestion of food waste under environmental variability

BACKGROUND

The interaction among microorganisms in the anaerobic digestion of food waste (ADFW) reactors lead to the degradation of organics and the recycling of energy. Viruses are an important component of the microorganisms involved in ADFW, but are rarely investigated. Furthermore, little is known about how viruses affect methanogenesis.

RESULTS

Thousands of viral sequences were recovered from five full-scale ADFW reactors. Gene-sharing networks indicated that the ADFW samples contained substantial numbers of unexplored anaerobic-specific viruses. Moreover, the viral communities in five full-scale reactors exhibited both commonalities and heterogeneities. The lab-scale dynamic analysis of typical ADFW scenarios suggested that the viruses had similar kinetic characteristics to their prokaryotic hosts. By associating with putative hosts, a majority of the bacteria and archaea phyla were found to be infected by viruses. Viruses may influence prokaryotic ecological niches, and thus methanogenesis, by infecting key functional microorganisms, such as sulfate-reducing bacteria (SRB), syntrophic acetate-oxidizing bacteria (SAOB), and methanogens. Metabolic predictions for the viruses suggested that they may collaborate with hosts at key steps of sulfur and long-chain fatty acid (LCFA) metabolism and could be involved in typical methanogenesis pathways to participate in methane production.

CONCLUSIONS

Our results expanded the diversity of viruses in ADFW systems and suggested two ways that viral manipulated ADFW biochemical processes. Video Abstract.

Virus impacted community adaptation in oligotrophic groundwater environment revealed by Hi-C coupled metagenomic and viromic study

Viruses play a crucial role in microbial mortality, diversity and biogeochemical cycles. Groundwater is the largest global freshwater and one of the most oligotrophic aquatic systems on Earth, but how microbial and viral communities are shaped in this special habitat is largely unexplored. In this study, we collected groundwater samples from 23 to 60 m aquifers at Yinchuan Plain, China. In total, 1920 non-reductant viral contigs were retrieved from metagenomes and viromes constructed by Illumina and Nanopore hybrid sequencing. Only 3% of them could be clustered with known viruses, most of which were Caudoviricetes. Coupling 1.2 Tb Hi-C sequencing with CRISPR matching and homology search, we connected 469 viruses with their hosts while some viral clusters presented a broad-host-range trait. Meanwhile, a large proportion of biosynthesis related auxiliary metabolism genes were identified. Those characteristics might benefit viruses for a better survival in this special oligotrophic environment. Additionally, the groundwater virome showed genomic features distinct from those of the open ocean and wastewater treatment facilities in GC distribution and unannotated gene compositions. This paper expands the current knowledge of the global viromic records and serves as a foundation for a more thorough understanding of viruses in groundwater.

Nonbacterial Microflora in Wastewater Treatment Plants: an Underappreciated Potential Source of Pathogens

Wastewater treatment plants (WWTPs) receive and treat large volumes of domestic, industrial, and urban wastewater containing pathogenic and nonpathogenic microorganisms, chemical compounds, heavy metals, and other potentially hazardous substances. WWTPs play an essential role in preserving human, animal, and environmental health by removing many of these toxic and infectious agents, particularly biological hazards. Wastewater contains complex consortiums of bacterial, viral, archaeal, and eukaryotic species, and while bacteria in WWTP have been extensively studied, the temporal and spatial distribution of nonbacterial microflora (viruses, archaea, and eukaryotes) is less understood. In this study, we analyzed the viral, archaeal, and eukaryotic microflora in wastewater throughout a treatment plant (raw influent, effluent, oxidation pond water, and oxidation pond sediment) in Aotearoa (New Zealand) using Illumina shotgun metagenomic sequencing. Our results suggest a similar trend across many taxa, with an increase in relative abundance in oxidation pond samples compared to influent and effluent samples, except for archaea, which had the opposite trend. Additionally, some microbial families, such as Podoviridae bacteriophages and Apicomplexa alveolates, appeared largely unaffected by the treatment process, with their relative abundance remaining stable throughout. Several groups encompassing pathogenic species, such as Leishmania, Plasmodium, Toxoplasma, Apicomplexa, Cryptococcus, Botrytis, and Ustilago, were identified. If present, these potentially pathogenic species could be a threat to human and animal health and agricultural productivity; therefore, further investigation is warranted. These nonbacterial pathogens should be considered when assessing the potential for vector transmission, distribution of biosolids to land, and discharge of treated wastewater to waterways or land. IMPORTANCE Nonbacterial microflora in wastewater remain understudied compared to their bacterial counterparts despite their importance in the wastewater treatment process. In this study, we report the temporal and spatial distributions of DNA viruses, archaea, protozoa, and fungi in raw wastewater influent, effluent, oxidation pond water, and oxidation pond sediments by using shotgun metagenomic sequencing. Our study indicated the presence of groups of nonbacterial taxa which encompass pathogenic species that may have potential to cause disease in humans, animals, and agricultural crops. We also observed higher alpha diversity in viruses, archaea, and fungi in effluent samples than in influent samples. This suggests that the resident microflora in the wastewater treatment plant may be making a greater contribution to the diversity of taxa observed in wastewater effluent than previously thought. This study provides important insights to better understand the potential human, animal, and environmental health impacts of discharged treated wastewater.

Mesophilic and thermophilic viruses are associated with nutrient cycling during hyperthermophilic composting

While decomposition of organic matter by bacteria plays a major role in nutrient cycling in terrestrial ecosystems, the significance of viruses remains poorly understood. Here we combined metagenomics and metatranscriptomics with temporal sampling to study the significance of mesophilic and thermophilic bacteria and their viruses on nutrient cycling during industrial-scale hyperthermophilic composting (HTC). Our results show that virus-bacteria density dynamics and activity are tightly coupled, where viruses specific to mesophilic and thermophilic bacteria track their host densities, triggering microbial community succession via top-down control during HTC. Moreover, viruses specific to mesophilic bacteria encoded and expressed several auxiliary metabolic genes (AMGs) linked to carbon cycling, impacting nutrient turnover alongside bacteria. Nutrient turnover correlated positively with virus-host ratio, indicative of a positive relationship between ecosystem functioning, viral abundances, and viral activity. These effects were predominantly driven by DNA viruses as most detected RNA viruses were associated with eukaryotes and not associated with nutrient cycling during the thermophilic phase of composting. Our findings suggest that DNA viruses could drive nutrient cycling during HTC by recycling bacterial biomass through cell lysis and by expressing key AMGs. Viruses could hence potentially be used as indicators of microbial ecosystem functioning to optimize productivity of biotechnological and agricultural systems.

CRISPR-Cas phage defense systems and prophages in Candidatus Accumulibacter

Candidatus Accumulibacter plays a major role in enhanced biological phosphorus removal (EBPR) from wastewater. Although bacteriophages have been shown to represent fatal threats to Ca. Accumulibacter organisms and thus interfere with the stability of the EBPR process, little is known about the ability of different Ca. Accumulibacter strains to resist phage infections. We conducted a systematic analysis of the occurrence and characteristics of clustered regularly interspaced short palindromic repeats and associated proteins (CRISPR-Cas) systems and prophages in Ca. Accumulibacter lineage members (43 in total, including 10 newly recovered genomes). Results indicate that 28 Ca. Accumulibacter genomes encode CRISPR-Cas systems. They were likely acquired via horizontal gene transfer, conveying a distinct adaptivity to phage predation to different Ca. Accumulibacter members. Major differences in the number of spacers show the unique phage resistance of these members. A comparison of the spacers in closely related Ca. Accumulibacter members from distinct geographical locations indicates that habitat isolation may have resulted in the acquisition of resistance to different phages by different Ca. Accumulibacter. Long-term operation of three laboratory-scale EBPR bioreactors revealed high relative abundances of Ca. Accumulibacter with CRISPSR-Cas systems. Their specific resistance to phages in these reactors was indicated by spacer analysis. Metatranscriptomic analyses showed the activation of the CRISPR-Cas system under both anaerobic and aerobic conditions. Additionally, 133 prophage regions were identified in 43 Ca. Accumulibacter genomes. Twenty-seven of them (in 19 genomes) were potentially active. Major differences in the occurrence of CRISPR-Cas systems and prophages in Ca. Accumulibacter will lead to distinct responses to phage predation. This study represents the first systematic analysis of CRISPR-Cas systems and prophages in the Ca. Accumulibacter lineage, providing new perspectives on the potential impacts of phages on Ca. Accumulibacter and EBPR systems.

Deciphering endogenous and exogenous regulations of anammox consortia in responding to lincomycin by multiomics: quorum sensing and CRISPR system

The widespread use of antibiotics has created an antibiotic resistance genes (ARGs)-enriched environment, which causes high risks on human and animal health. Although antibiotics can be partially adsorbed and degraded in wastewater treatment processes, striving for a complete understanding of the microbial adaptive mechanism to antibiotic stress remains urgent. Combined with metagenomics and metabolomics, this study revealed that anammox consortia could adapt to lincomycin by spontaneously changing the preference for metabolite utilization and establishing interactions with eukaryotes, such as Ascomycota and Basidiomycota. Specifically, quorum sensing (QS) based microbial regulation and the ARGs transfer mediated by clustered regularly interspaced short palindromic repeats (CRISPR) system and global regulatory genes were the principal adaptive strategies. Western blotting results validated that Cas9 and TrfA were mainly responsible for the alteration of ARGs transfer pathway. These findings highlight the potential adaptative mechanism of microbes to antibiotic stress and fill gaps in horizontal gene transfer pathways in the anammox process, further facilitating the ARGs control through molecular and synthetic biology techniques.

Potential Auxiliary Metabolic Capabilities and Activities Reveal Biochemical Impacts of Viruses in Municipal Wastewater Treatment Plants

Viruses influence biogeochemical cycles in oceans, freshwater, soil, and human gut through infection and by modulating virocell metabolism through virus-encoded auxiliary metabolic genes (vAMGs). However, the geographical distribution, potential metabolic function, and engineering significance of vAMGs in wastewater treatment plants (WWTPs) remain to be explored. Here, 752 single-contig viral genomes with high confidence, 510 of which belonged to Caudovirales, were recovered from the activated sludge metagenomes of 32 geographically distributed WWTPs. A total of 101 vAMGs involved in various metabolic pathways were identified, the most common of which were the queuosine biosynthesis genes folE, queD, and queE and the sulfur metabolism gene cysH. Phylogenetic analysis and virus-host relationship prediction revealed the probable evolutionary histories of vAMGs involved in carbon (acpP and prsA), nitrogen (amoC), sulfur (cysH), and phosphate (phoH) metabolism, which potentially mediate microbial carbon and nutrient cycling. Notably, 11 of the 38 (28.3%) vAMGs identified in the metagenomes with corresponding metatranscriptomes were transcriptionally expressed, implying an active functional state. This meta-analysis provides the first broad catalog of vAMGs in municipal WWTPs and how they may assist in the basic physiological reactions of their microbial hosts or nutrient cycling in the WWTPs, and therefore, may have important effects on the engineering of wastewater treatment processes.

Microbial density-dependent viral dynamics and low activity of temperate phages in the activated sludge process

The ecological behavior of bacteriophages (phages), the most abundant biological entity in wastewater treatment systems, is poorly understood, especially that of temperate phages. Here, the temporal dynamics of lytic and temperate phages in a laboratory-scale activated sludge reactor with a sludge bulking issue was investigated using coupled sludge metagenomic and viromic analyses. The lysogenic fragments (prophages) identified were widely distributed in the reconstructed metagenome-assembled genomes (61.7%, n = 227). However, only 12.3% of the identified prophages experienced lysogenic-lytic switching, and the abundance contribution of prophages to free virus communities was only 0.02-0.3%, indicating low activity of temperate phages. Although the sludge community changed dramatically during reactor operation, no massive prophage induction events were detected. Statistical analyses showed strong correlations between sludge concentration and free virus and temperate phage communities, suggesting microbial density-dependent virus dynamics in the sludge microbiota.

Enhanced purification of kitchen-oil wastewater driven synergistically by surface microelectric fields and microorganisms

The stable structure and toxic effect of refractory organic pollutants in wastewater lead to the problem of high energy consumption in water treatment technology. Herein, we propose a synergistic purification of refractory wastewater driven by microorganisms and surface microelectric fields (SMEF) over a dual-reaction-center (DRC) catalyst HCLL-S8-M prepared by an in situ growth method of carbon nitride on the Cu-Al₂O₃ surface. Characterization techniques demonstrate the successful construction of SMEF with strong electrostatic force over HCLL-S8-M based on cation-π interactions between metal copper ions and carbon nitride rings. With the catalyst as the core filler, an innovative fixed bed bioreactor is constructed to purify the actual kitchen-oil wastewater. The removal efficiency of the wastewater even with a very low biodegradability (BOD₅/COD = 0.33) can reach 60% after passing through this bioreactor. An innovative reaction mechanism is revealed for the first time that under the condition of a small amount of biodegradable organic matter, the SMEF induces the enrichment of electric active microorganisms (Desulfobulbus and Geobacter) in the wastewater, accelerates the interspecies electron transfer of intertrophic metabolism with the biodegradable bacteria through the extracellular electron transfer mechanism such as cytochrome C and self-secreted electron shuttle. The electrons of the refractory organic pollutants adsorbed on the surface of the catalyst are delocalized by the SMEF, which can be directly utilized by microorganisms through EPS conduction. The SMEF generated by electron polarization can maximize the utilization of pollutants and microorganisms in wastewater and further enhance degradation without adding any external energy, which is of great significance to the development of water self-purification technology.

The microbial dark matter and "wanted list" in worldwide wastewater treatment plants

BACKGROUND

Wastewater treatment plants (WWTPs) are one of the largest biotechnology applications in the world and are of critical importance to modern urban societies. An accurate evaluation of the microbial dark matter (MDM, microorganisms whose genomes remain uncharacterized) proportions in WWTPs is of great value, while there is no such research yet. This study conducted a global meta-analysis of MDM in WWTPs with 317,542 prokaryotic genomes from the Genome Taxonomy Database and proposed a "wanted list" for priority targets in further investigations of activated sludge.

RESULTS

Compared with the Earth Microbiome Project data, WWTPs had relatively lower genome-sequenced proportions of prokaryotes than other ecosystems, such as the animal related environments. Analysis showed that the median proportions of the genome-sequenced cells and taxa (100% identity and 100% coverage in 16S rRNA gene region) in WWTPs reached 56.3% and 34.5% for activated sludge, 48.6% and 28.5% for aerobic biofilm, and 48.3% and 28.5% for anaerobic digestion sludge, respectively. This result meant MDM had high proportions in WWTPs. Besides, all of the samples were occupied by a few predominant taxa, and the majority of the sequenced genomes were from pure cultures. The global-scale "wanted list" for activated sludge contained four phyla that have few representatives and 71 operational taxonomic units with the majority of them having no genome or isolate yet. Finally, several genome mining methods were verified to successfully recover genomes from activated sludge such as hybrid assembly of the second- and third-generation sequencing.

CONCLUSIONS

This work elucidated the proportion of MDM in WWTPs, defined the "wanted list" of activated sludge for future investigations, and certified potential genome recovery methods. The proposed methodology of this study can be applied to other ecosystems and improve understanding of ecosystem structure across diverse habitats. Video Abstract.