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.

Stringent Response of Cyanobacteria and Other Bacterioplankton during Different Stages of a Harmful Cyanobacterial Bloom.

We conducted a field study to investigate the role of stringent response in cyanobacteria and coexisting bacterioplankton during nutrient-deprived periods at various stages of bloom in a freshwater lake (Utah Lake) for the first time. Using metagenomics and metatranscriptomics analyses, we examined the cyanobacterial ecology and expression of important functional genes related to stringent response, N and P metabolism, and regulation. Our findings mark a significant advancement in understanding the mechanisms by which toxic cyanobacteria survive and proliferate during nitrogen (N) and phosphorus (P) limitations. We successfully identified and analyzed the metagenome-assembled genomes (MAGs) of the dominant bloom-forming cyanobacteria, namely, Dolichospermum circinale, Aphanizomenon flos-aquae UKL13-PB, Planktothrix agardhii, and Microcystis aeruginosa. By mapping RNA-seq data to the coding sequences of the MAGs, we observed that these four prevalent cyanobacteria species activated multiple functions to adapt to the depletion of inorganic nutrients. During and after the blooms, the four dominant cyanobacteria species expressed high levels of transcripts related to toxin production, such as microcystins (mcy), anatoxins (ana), and cylindrospermopsins (cyr). Additionally, genes associated with polyphosphate (poly-P) storage and the stringent response alarmone (p)ppGpp synthesis/hydrolysis, including ppk, relA, and spoT, were highly activated in both cyanobacteria and bacterioplankton. Under N deficiency, the main N pathways shifted from denitrification and dissimilatory nitrate reduction in bacterioplankton toward N2-fixing and assimilatory nitrate reduction in certain cyanobacteria with a corresponding shift in the community composition. P deprivation triggered a stringent response mediated by spoT-dependent (p)ppGpp accumulation and activation of the Pho regulon in both cyanobacteria and bacterioplankton, facilitating inorganic and organic P uptake. The dominant cyanobacterial MAGs exhibited the presence of multiple alkaline phosphatase (APase) transcripts (e.g., phoA in Dolichospermum, phoX in Planktothrix, and Microcystis), suggesting their ability to synthesize and release APase enzymes to convert ambient organic P into bioavailable forms. Conversely, transcripts associated with bacterioplankton-dominated pathways like denitrification were low and did not align with the occurrence of intense cyanoHABs. The strong correlations observed among N, P, stringent response metabolisms and the succession of blooms caused by dominant cyanobacterial species provide evidence that the stringent response, induced by nutrient limitation, may activate unique N and P functions in toxin-producing cyanobacteria, thereby sustaining cyanoHABs.

Viruses and Their Interactions With Bacteria and Archaea of Hypersaline Great Salt Lake.

Viruses play vital biogeochemical and ecological roles by (a) expressing auxiliary metabolic genes during infection, (b) enhancing the lateral transfer of host genes, and (c) inducing host mortality. Even in harsh and extreme environments, viruses are major players in carbon and nutrient recycling from organic matter. However, there is much that we do not yet understand about viruses and the processes mediated by them in the extreme environments such as hypersaline habitats. The Great Salt Lake (GSL) in Utah, United States is a hypersaline ecosystem where the biogeochemical role of viruses is poorly understood. This study elucidates the diversity of viruses and describes virus-host interactions in GSL sediments along a salinity gradient. The GSL sediment virosphere consisted of Haloviruses (32.07 ± 19.33%) and members of families Siphoviridae (39.12 ± 19.8%), Myoviridae (13.7 ± 6.6%), and Podoviridae (5.43 ± 0.64%). Our results demonstrate that salinity alongside the concentration of organic carbon and inorganic nutrients (nitrogen and phosphorus) governs the viral, bacteria, and archaeal diversity in this habitat. Computational host predictions for the GSL viruses revealed a wide host range with a dominance of viruses that infect Proteobacteria, Actinobacteria, and Firmicutes. Identification of auxiliary metabolic genes for photosynthesis (psbA), carbon fixation (rbcL, cbbL), formaldehyde assimilation (SHMT), and nitric oxide reduction (NorQ) shed light on the roles played by GSL viruses in biogeochemical cycles of global relevance.

Differences in chlorine and peracetic acid disinfection kinetics of Enterococcus faecalis and Escherichia fergusonii and their susceptible strains based on gene expressions and genomics.

This study was conducted to investigate mechanisms of cross-resistance to chlorine and peracetic acid (PAA) disinfectants by antibiotic-resistant bacteria. Our study evaluated chlorine and PAA based disinfection kinetics of erythromycin-resistant Enterococcus faecalis, meropenem-resistant Escherichia fergusonii, and susceptible strains of these species. Using the integrated second-order disinfectant decay model and first-order Chick-Watson's Law, it was found that the meropenem-resistant Escherichia fergusonii strain showed significantly less log inactivation compared to the susceptible E. fergusonii strain in response to both chlorine and PAA disinfection (p-value = 0.059, 3.5 × 10-6). On the other hand, the susceptible Enterococcus faecalis strain showed similar log inactivation compared to the erythromycin-resistant strain in response to either treatment (p-value = 0.075, 0.28). Meropenem-resistant E. fergusonii showed an increase in gene expression of New Delhi metallo-ß-lactamase (blaNDM-1) gene to chlorine, but there was no increase in expression to PAA. Whole genome sequencing (WGS) was then conducted to elucidate the differences in genes among both resistant and susceptible table E. fergusonii strains. The average nucleotide identity (ANI) analysis of the draft genomes (>97% similarity) suggests that meropenem-resistant E. fergusonii (S1) and meropenem-susceptible E. fergusonii (S2) are the same species but different strains. Both strains have the same genes for oxidative stress pathways, oxidative scavenger genes, and nearly 40 different antibiotic efflux pump genes. The chromosomal and plasmid draft genomes of meropenem-resistant and susceptible E. fergusonii strains each have 65 and 52 antibiotic resistance genes, respectively. Of these, the resistant E. fergusonii strain harbored the extended-spectrum beta-lactamases blaCTX-M-15 and blaTEM-1 genes located on the chromosome, and a blaTEM-1 gene on the plasmid. The overall findings of this study are significant, as they reveal that antibiotic-resistant and susceptible strains of E. fergusonii exhibit different responses towards chlorine and PAA disinfection.

Prevalence of SARS-CoV-2 genes in water reclamation facilities: From influent to anaerobic digester.

Several treatment plants were sampled for influent, primary clarifier sludge, return activated sludge (RAS), and anaerobically digested sludge throughout nine weeks during the summer of the COVID-19 pandemic. Primary clarifier sludge had a significantly higher number of SARS-CoV-2 gene copy number per liter (GC/L) than other sludge samples, within a range from 1.0 × 105 to 1.0 × 106 GC/L. Gene copy numbers in raw influent significantly correlated with gene copy numbers in RAS in Silver Creek (p-value = 0.007, R2 = 0.681) and East Canyon (p-value = 0.009, R2 = 0.775) WRFs; both of which lack primary clarifiers or industrial pretreatment processes. This data indicates that SARS-CoV-2 gene copies tend to partition into primary clarifier sludges, at which point a significant portion of them are removed through sedimentation. Furthermore, it was found that East Canyon WRF gene copy numbers in influent were a significant predictor of daily cases (p-value = 0.0322, R2 = 0.561), and gene copy numbers in RAS were a significant predictor of weekly cases (p-value = 0.0597, R2 = 0.449). However, gene copy numbers found in primary sludge samples from other plants significantly predicted the number of COVID-19 cases for the following week (t = 2.279) and the week after that (t = 2.122) respectively. These data indicate that SARS-CoV-2 extracted from WRF biosolids may better suit epidemiological monitoring that exhibits a time lag. It also supports the observation that primary sludge removes a significant portion of SARS-CoV-2 marker genes. In its absence, RAS can also be used to predict the number of COVID-19 cases due to direct flow through from influent. This research represents the first of its kind to thoroughly examine SARS-CoV-2 gene copy numbers in biosolids throughout the wastewater treatment process and the relationship between primary, return activated, and anaerobically digested sludge and reported positive COVID-19 cases.