Phage-host interactions: The neglected part of biological wastewater treatment

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.

Comparative life cycle assessment of traditional and emerging methods for high concentration organic wastewater treatment

Two traditional methods, including incineration and the sequencing batch reactor activated sludge process (SBR), and two emerging methods, including supercritical water oxidation (SCWO), catalytic wet oxidation (CWO) combined with SBR (CWO-SBR), for high concentration organic wastewater treatment were compared from the view of environmental impact. The printing ink wastewater generated from the surface treatment of electronic products was selected as the typical wastewater and diluted to various concentrations for sensitivity analysis. Simapro 9.3 software was selected to explore the environmental impacts of the four processes with Traci impact assessment method and Ecoinvent 3.8 database. Moreover, two tools including the total environmental impact and comprehensive competitiveness were used for quantitative comparison. The total environmental impacts for the incineration and SCWO processes gradually decrease with increasing chemical oxygen demand (COD) concentration for more steam output and less fuel and electricity consumptions, but opposite trends appear for the CWO-SBR and SBR processes, indicating higher COD concentrations exhibit smaller environmental impacts for the incineration and SCWO processes. The incineration and SBR processes have the slightest and highest total environmental impacts of 0.030 and 1.117 at the base case of COD = 234,000 mg/L, respectively. The trends of the comprehensive indexes for the four processes with increasing COD concentration are similar to those of the total environmental impacts. The SCWO process has the highest comprehensive index of 0.956 at the base case, followed by the incineration, CWO-SBR, and SBR processes, with values of 0.783, 0.460, and 0.158, respectively.

Genome streamlining of Pseudomonas putida B6-2 for bioremediation

Microbial transformation is a favored approach for environmental remediation. However, the effectiveness of microbial remediation has been limited by the lack of chassis cells with satisfactory contaminant degradation performance. Pseudomonas putida B6-2, with a wide substrate spectrum and high solvent tolerance, is a chassis strain with great potential for application in environmental remediation. Here, guided by bioinformatic analyses and genome-scale metabolic model (GEM) predictions, we successfully optimized P. putida B6-2 by rationally reducing its nonessential genetic components and generating a more robust genome-streamlined strain, P. putida BGR4. Several improvements were observed compared with the original P. putida B6-2 strain, including a 1.4 × 105-fold increase in electroporation efficiency, an 8.3-fold increase in conjugation efficiency, improved glycerol utilization capability, and increased phenol utilization after heterologous expression of the phenol monooxygenase encoded by dmpKLMNOP. Additionally, P. putida BGR4 exhibited enhanced tolerance to several stressors, including starvation, oxidative stress, and DNA damage. Transcriptomic analysis revealed that genome streamlining led to the upregulation of genes involved in the "carbon metabolism" and "tricarboxylic acid cycle" pathways in P. putida BGR4, which likely contributed to the superior phenotype of P. putida BGR4 in terms of carbon source utilization and contaminant degradation capabilities. Furthermore, the absence of four prophages was identified as a potential cause of the enhanced stress resistance observed in P. putida BGR4. Overall, we developed a combined genome-streamlining strategy involving bioinformatic analyses and GEM predictions and generated a more robust chassis strain, P. putida BGR4, which expands the repertoire of chassis cells for environmental remediation.IMPORTANCEDespite the development of many chassis cells, there is still a lack of robust chassis cells with satisfactory contaminant degradation performance. Targeted genome streamlining is an effective way to provide powerful chassis cells. However, genome streamlining does not always lead to the improved phenotypes of genome-streamlined chassis cells. In this research, a novel procedure that combined bioinformatic analyses and GEM predictions was proposed to guide genome streamlining and predict the effects of genome streamlining. This genome streamlining procedure was successfully applied to Pseudomonas putida B6-2, which was a chassis cell with great potential for application in environmental remediation and resulted in the generation of a more robust chassis cell, P. putida BGR4, thereby providing a superior chassis cell for efficient and sustainable environmental remediation and a valuable framework for guiding the genome streamlining of strains for other applications.

Genomic mapping of wastewater bacteriophage may predict potential bacterial pathogens infecting the community

Most existing wastewater surveillance studies that focus on viruses have identified a large fraction of bacteriophages. Identifying bacteria by considering bacteriophage-host interactions is a novel method for detecting bacterial pathogens circulating in a community, using wastewater surveillance. This study aims to identify human-related bacterial pathogens in municipal wastewater collected in metro Detroit, using high-throughput sequencing and bioinformatics. Untreated municipal wastewater samples were collected on August 11, 2020, and bacteriophages were concentrated using the VIRus ADsorption-ELution (VIRADEL) method. Bacteriophage-related contigs in samples ranged from 15.53 % to 18.91 %, with 2477 classified and 8853 unclassified contigs. Most identified bacteriophages were from Caudoviricetes and Malgrandaviricetes classes belonging to 19 families. Hosts of bacteriophages were predicted with the PhaBOX (CHERRY) tool. The results indicated that out of the 2477 classified phages, 2373 were associated with known bacterial hosts. Also, out of 8853 unclassified bacteriophages, 8421 were associated with known bacterial hosts, and the remaining 432 were with unknown bacterial hosts. Among all bacteriophage-associated hosts, 399 were identified as pathogenic bacteria at the species level. Approximately, 85 % of the identified pathogenic bacteria are reported to be associated with human diseases. Genome quality assessments showed that 15 bacteriophages had nearly complete genomes, which were further analyzed to understand bacteriophage-bacteria interactions in wastewater. Identified hosts of these complete-genome phages included human pathogens such as Salmonella enterica, Bacillus cereus, Achromobacter xylosoxidans, and Escherichia coli. The S. enterica bacteriophage (k141_1005294) genomic map was annotated, and responsible open reading frames (ORFs) were characterized to illustrate bacteriophage behavior during infection of pathogenic bacteria in untreated wastewater. To the best of our knowledge, this is the first attempt to characterize human bacterial pathogens in wastewater through bacteriophage-pathogen interactions. Novel bioinformatic approaches enhance pathogen detection and improve the understanding of community wastewater microbiomes.

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.

Microbiome regulation for sustainable wastewater treatment

Sustainable wastewater treatment is essential for attaining clean water and sanitation, aligning with UN Sustainable Development Goals. Wastewater treatment plants (WWTPs) have utilized environmental microbiomes in biological treatment processes in this effort for over a century. However, the inherent complexity and redundancy of microbial communities, and emerging chemical and biological contaminants, challenge the biotechnology applications. Over the past decades, understanding and utilization of microbial energy metabolism and interaction relationships have revolutionized the biological system. In this review, we discuss how microbiome regulation strategies are being used to generate actionable performance for low-carbon pollutant removal and resource recovery in WWTPs. The engineering application cases also highlight the real feasibility and promising prospects of the microbiome regulation approaches. In conclusion, we recommend identifying environmental risks associated with chemical and biological contaminants transformation as a prerequisite. We propose the integration of gene editing and enzyme design to precisely regulate microbiomes for the synergistic control of both chemical and biological risks. Additionally, the development of integrated technologies and engineering equipment is crucial in addressing the ongoing water crisis. This review advocates for the innovation of conventional wastewater treatment biotechnology to ensure sustainable wastewater treatment.

Supporting data-enhanced hybrid ordinary differential equation model for phosphate dynamics in municipal wastewater treatment

A parallel hybrid ordinary differential equation (ODE) integrating the Activated Sludge Model No. 2d (ASM2d) and an artificial neural network (ANN) was developed to simulate biological phosphorus removal (BPR) with high accuracy and interpretability. Two novelties were introduced; first, the involved supporting data (i.e., phosphate-release activity) were incorporated as an input in the ANN. Second, the outputs of the ANN were selective. Three models were implemented using different ANN outputs, and all three outperformed ASM2d in phosphate estimation for anaerobic/aerobic sequencing batch reactor operation. In particular, the incorporation of four variables responsible for BPR into the ANN enabled the highest performance (R2 = 0.93) owing to the capture of increasing phosphate-accumulating organisms (PAOs). The ANN with the supporting data worked satisfactorily to compensate for ASM2d by adding proper PAOs, resulting in improvement in phosphate estimation. The novel parallel hybrid ODE can simulate BPR while maintaining physical meaning.

Application of machine learning models to improve the prediction of pesticide photodegradation in water by ZnO-based photocatalysts

Pesticide pollution has been posing a significant risk to human and ecosystems, and photocatalysis is widely applied for the degradation of pesticides. Machine learning (ML) emerges as a powerful method for modeling complex water treatment processes. For the first time, this study developed novel ML models that improved the estimation of the photocatalytic degradation of various pesticides using ZnO-based photocatalysts. The input parameters encompassed the source of light, mass proportion of dopants to Zn, initial pesticide concentration (C0), pH of the solution, catalyst dosage and irradiation time. Additionally, physicochemical properties such as the molecular weight of the dopants and pesticides, as well as the water solubility of both dopants and pesticides, were considered. Notably, the numerical data were extracted from the literature via relevant tables (directly) or graphs (indirectly) using the web-based tool WebPlotDigitizer. Four ML models including multi-layer perceptron artificial neural network (MLP-ANN), particle swarm optimization-adaptive neuro fuzzy inference system (PSO-ANFIS), radial basis function (RBF), and coupled simulated annealing-least squares support vector machine (CSA-LSSVM) were developed. In comparison, RBF showed the best accuracy of modeling among all models, with the highest determination coefficient (R2) of 0.978 and average absolute relative deviation (AARD) of 4.80%. RBF model was effective in estimating the photocatalytic degradation of pesticides except for 2-chlorophenol, triclopyr and lambda-cyhalothrin, where CSA-LSSVM model demonstrated superior performance. Dichlorvos was completely degraded by ZnO photocatalyst under visible light. The sensitivity analysis by relevancy factor exhibited that light irradiation time and initial pesticide concentration were the most important parameters influencing photocatalytic degradation of pesticides positively and negatively, respectively. The new ML models provide a powerful tool for predicting pesticide degradation in wastewater treatment, which will reduce photochemical experiments and promote sustainable development.

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.

Correlation between the development of phage resistance and the original antibiotic resistance of host bacteria under the co-exposure of antibiotic and bacteriophage

Bacteriophages (phages) are viruses capable of regulating the proliferation of antibiotic resistant bacteria (ARB). However, phages that directly cause host lethality may quickly select for phage resistant bacteria, and the co-evolutionary trade-offs under varying environmental conditions, including the presence of antibiotics, remains unclear as to their impact on phage and antibiotic resistance. Here, we report the emergence of phage resistance in three distinct E. coli strains with varying resistance to β-lactam antibiotics, treated with different ampicillin (AMP) concentrations. Hosts exhibiting stronger antibiotic resistance demonstrated a higher propensity to develop and maintain stable phage resistance. When exposed to polyvalent phage KNT-1, the growth of AMP-sensitive E. coli K12 was nearly suppressed within 18 h, while the exponential growth of AMP-resistant E. coli TEM and super-resistant E. coli NDM-1 was delayed by 12 h and 8 h, respectively. The mutation frequency and mutated colony count of E. coli NDM-1 were almost unaffected by co-existing AMP, whereas for E. coli TEM and K12, these metrics significantly decreased with increasing AMP concentration from 8 to 50 μg/mL, becoming unquantifiable at 100 μg/mL. Furthermore, the fitness costs of phage resistance mutation and its impact on initial antibiotic resistance in bacteria were further examined, through analyzing AMP susceptibility, biofilm formation and EPS secretion of the isolated phage resistant mutants. The results indicated that acquiring phage resistance could decrease antibiotic resistance, particularly for hosts lacking strong antibiotic resistance. The ability of mutants to form biofilm contributes to antibiotic resistance, but the correlation is not entirely positive, while the secretion of extracellular polymeric substance (EPS), especially the protein content, plays a crucial role in protecting the bacteria from both antibiotic and phage exposure. This study explores phage resistance development in hosts with different antibiotic resistance and helps to understand the limitations and possible solutions of phage-based technologies.

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.

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

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

Short-, long-read metagenome and virome reveal the profile of phage-mediated ARGs in anoxic-oxic processes for swine wastewater treatment

Phages are among the most widely spread viruses, but their profiles and the antibiotic resistance genes (ARGs) they carry in swine wastewater remain underexplored. The present study investigated the distribution characteristics of phages and their ARG risk in anoxic/oxic (A/O) wastewater treatment processes of swine farms using short- and long-read metagenome and virome. The results demonstrated that the virome could extract more phage sequences than the total metagenome; thus, it was more suited for studying phages in wastewater settings. Intriguingly, phages had significantly lower abundance of ARG than ARGs harbored by total microorganisms (P < 0.01). Eleven ARGs co-occurred with phages and bacteria (R > 0.6 and P < 0.05), with Siphoviridae being the phage co-occurring with the most ARGs (5). Horizontal gene transfer (HGT) events were observed between Proteobacteria and the major phyla except for Bacteroidota. Furthermore, there were prophage sequences and ARGs on the same contig in bacterial MAGs. These data strongly demonstrate that phages promote horizontal transfer of ARG between bacterial hosts in A/O processes for swine wastewater treatment. Therefore, the risk of phage-mediated horizontal transfer of ARGs cannot be overlooked despite the low abundance of phage ARGs (pARG).

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.

Recent advances and future perspective on lignocellulose-based materials as adsorbents in diverse water treatment applications

The growing shortage of non-renewable resources and the burden of toxic pollutants in water have gradually become stumbling blocks in the path of sustainable human development. To this end, there has been great interest in finding renewable and environmentally friendly materials to promote environmental sustainability and combat harmful pollutants in wastewater. Of the many options, lignocellulose, as an abundant, biocompatible and renewable material, is the most attractive candidate for water remediation due to the unique physical and chemical properties of its constituents. Herein, we review the latest research advances in lignocellulose-based adsorbents, focusing on lignocellulosic composition, material modification, application of adsorbents. The modification and preparation methods of lignin, cellulose and hemicellulose and their applications in the treatment of diverse contaminated water are systematically and comprehensively presented. Also, the detailed description of the adsorption model, the adsorption mechanism and the adsorbent regeneration technique provides an excellent reference for understanding the underlying adsorption mechanism and the adsorbent recycling. Finally, the challenges and limitations of lignocellulosic adsorbents are evaluated from a practical application perspective, and future developments in the related field are discussed. In summary, this review offers rational insights to develop lignocellulose-based environmentally-friendly reactive materials for the removal of hazardous aquatic contaminants.

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.

An Eco-evolutionary Model on Surviving Lysogeny Through Grounding and Accumulation of Prophages

Temperate phages integrate into the bacterial genomes propagating along with the bacterial genomes. Multiple phage elements, representing diverse prophages, are present in most bacterial genomes. The evolutionary events and the ecological dynamics underlying the accumulation of prophage elements in bacterial genomes have yet to be understood. Here, we show that the local wastewater had 7% of lysogens (hosting mitomycin C-inducible prophages), and they showed resistance to superinfection by their corresponding lysates. Genomic analysis of four lysogens and four non-lysogens revealed the presence of multiple prophages (belonging to Myoviridae and Siphoviridae) in both lysogens and non-lysogens. For large-scale comparison, 2180 Escherichia coli genomes isolated from various sources across the globe and 523 genomes specifically isolated from diverse wastewaters were analyzed. A total of 15,279 prophages were predicted among 2180 E. coli genomes and 2802 prophages among 523 global wastewater isolates, with a mean of ~ 5 prophages per genome. These observations indicate that most putative prophages are relics of past bacteria-phage conflicts; they are "grounded" prophages that cannot excise from the bacterial genome. Prophage distribution analysis based on the sequence homology suggested the random distribution of E. coli prophages within and between E. coli clades. The independent occurrence pattern of these prophages indicates extensive horizontal transfers across the genomes. We modeled the eco-evolutionary dynamics to reconstruct the events that could have resulted in the prophage accumulation accounting for infection, superinfection immunity, and grounding. In bacteria-phage conflicts, the bacteria win by grounding the prophage, which could confer superinfection immunity.

Lysogenic bacteriophages encoding arsenic resistance determinants promote bacterial community adaptation to arsenic toxicity

Emerging evidence from genomics gives us a glimpse into the potential contribution of lysogenic bacteriophages (phages) to the environmental adaptability of their hosts. However, it is challenging to quantify this kind of contribution due to the lack of appropriate genetic markers and the associated controllable environmental factors. Here, based on the unique transformable nature of arsenic (the controllable environmental factor), a series of flooding microcosms was established to investigate the contribution of arsM-bearing lysogenic phages to their hosts' adaptation to trivalent arsenic [As(III)] toxicity, where arsM is the marker gene associated with microbial As(III) detoxification. In the 15-day flooding period, the concentration of As(III) was significantly increased, and this elevated As(III) toxicity visibly inhibited the bacterial population, but the latter quickly adapted to As(III) toxicity. During the flooding period, some lysogenic phages re-infected new hosts after an early burst, while others persistently followed the productive cycle (i.e., lytic cycle). The unique phage-host interplay contributed to the rapid spread of arsM among soil microbiota, enabling the quick recovery of the bacterial community. Moreover, the higher abundance of arsM imparted a greater arsenic methylation capability to soil microbiota. Collectively, this study provides experimental evidence for lysogenic phages assisting their hosts in adapting to an extreme environment, which highlights the ecological perspectives on lysogenic phage-host mutualism.

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.