Phage Genome Diversity in a Biogas-Producing Microbiome Analyzed by Illumina and Nanopore GridION Sequencing

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.

A unified compendium of prokaryotic and viral genomes from over 300 anaerobic digestion microbiomes

BACKGROUND

The anaerobic digestion process degrades organic matter into simpler compounds and occurs in strictly anaerobic and microaerophilic environments. The process is carried out by a diverse community of microorganisms where each species has a unique role and it has relevant biotechnological applications since it is used for biogas production. Some aspects of the microbiome, including its interaction with phages, remains still unclear: a better comprehension of the community composition and role of each species is crucial for a cured understanding of the carbon cycle in anaerobic systems and improving biogas production.

RESULTS

The primary objective of this study was to expand our understanding on the anaerobic digestion microbiome by jointly analyzing its prokaryotic and viral components. By integrating 192 additional datasets into a previous metagenomic database, the binning process generated 11,831 metagenome-assembled genomes from 314 metagenome samples published between 2014 and 2022, belonging to 4,568 non-redundant species based on ANI calculation and quality verification. CRISPR analysis on these genomes identified 76 archaeal genomes with active phage interactions. Moreover, single-nucleotide variants further pointed to archaea as the most critical members of the community. Among the MAGs, two methanogenic archaea, Methanothrix sp. 43zhSC_152 and Methanoculleus sp. 52maCN_3230, had the highest number of SNVs, with the latter having almost double the density of most other MAGs.

CONCLUSIONS

This study offers a more comprehensive understanding of microbial community structures that thrive at different temperatures. The findings revealed that the fraction of archaeal species characterized at the genome level and reported in public databases is higher than that of bacteria, although still quite limited. The identification of shared spacers between phages and microbes implies a history of phage-bacterial interactions, and specifically lysogenic infections. A significant number of SNVs were identified, primarily comprising synonymous and nonsynonymous variants. Together, the findings indicate that methanogenic archaea are subject to intense selective pressure and suggest that genomic variants play a critical role in the anaerobic digestion process. Overall, this study provides a more balanced and diverse representation of the anaerobic digestion microbiota in terms of geographic location, temperature range and feedstock utilization.

Viral Diversity in Benthic Abyssal Ecosystems: Ecological and Methodological Considerations

Viruses are the most abundant 'biological entities' in the world's oceans. However, technical and methodological constraints limit our understanding of their diversity, particularly in benthic abyssal ecosystems (>4000 m depth). To verify advantages and limitations of analyzing virome DNA subjected either to random amplification or unamplified, we applied shotgun sequencing-by-synthesis to two sample pairs obtained from benthic abyssal sites located in the North-eastern Atlantic Ocean at ca. 4700 m depth. One amplified DNA sample was also subjected to single-molecule long-read sequencing for comparative purposes. Overall, we identified 24,828 viral Operational Taxonomic Units (vOTUs), belonging to 22 viral families. Viral reads were more abundant in the amplified DNA samples (38.5-49.9%) compared to the unamplified ones (4.4-5.8%), with the latter showing a greater viral diversity and 11-16% of dsDNA viruses almost undetectable in the amplified samples. From a procedural point of view, the viromes obtained by direct sequencing (without amplification step) provided a broader overview of both ss and dsDNA viral diversity. Nevertheless, our results suggest that the contextual use of random amplification of the same sample and long-read technology can improve the assessment of viral assemblages by reducing off-target reads.

Characterization of phage vB_EcoS-EE09 infecting E. coli DSM613 Isolated from Wastewater Treatment Plant Effluent and Comparative Proteomics of the Infected and Non-Infected Host

Phages influence microbial communities, can be applied in phage therapy, or may serve as bioindicators, e.g., in (waste)water management. We here characterized the Escherichia phage vB_EcoS-EE09 isolated from an urban wastewater treatment plant effluent. Phage vB_EcoS-EE09 belongs to the genus Dhillonvirus, class Caudoviricetes. It has an icosahedral capsid with a long non-contractile tail and a dsDNA genome with an approximate size of 44 kb and a 54.6% GC content. Phage vB_EcoS-EE09 infected 12 out of the 17 E. coli strains tested. We identified 16 structural phage proteins, including the major capsid protein, in cell-free lysates by protein mass spectrometry. Comparative proteomics of protein extracts of infected E. coli cells revealed that proteins involved in amino acid and protein metabolism were more abundant in infected compared to non-infected cells. Among the proteins involved in the stress response, 74% were less abundant in the infected cultures compared to the non-infected controls, with six proteins showing significant less abundance. Repressing the expression of these proteins may be a phage strategy to evade host defense mechanisms. Our results contribute to diversifying phage collections, identifying structural proteins to enable better reliability in annotating taxonomically related phage genomes, and understanding phage-host interactions at the protein level.

Special Issue “Bacteriophage Genomics”: Editorial

Virus genomics as a separate branch of biology has emerged relatively recently [...]

Analysis of the anaerobic digestion metagenome under environmental stresses stimulating prophage induction

BACKGROUND

The viral community has the potential to influence the structure of the microbiome and thus the yield of the anaerobic digestion process. However, the virome composition in anaerobic digestion is still under-investigated. A viral induction experiment was conducted on separate batches undergoing a series of DNA-damaging stresses, in order to coerce temperate viruses to enter the lytic cycle.

RESULTS

The sequencing of the metagenome revealed a viral community almost entirely composed of tailed bacteriophages of the order Caudovirales. Following a binning procedure 1,092 viral and 120 prokaryotic genomes were reconstructed, 64 of which included an integrated prophage in their sequence. Clustering of coverage profiles revealed the presence of species, both viral and microbial, sharing similar reactions to shocks. A group of viral genomes, which increase under organic overload and decrease under basic pH, uniquely encode the yopX gene, which is involved in the induction of temperate prophages. Moreover, the in-silico functional analysis revealed an enrichment of sialidases in viral genomes. These genes are associated with tail proteins and, as such, are hypothesised to be involved in the interaction with the host. Archaea registered the most pronounced changes in relation to shocks and featured behaviours not shared with other species. Subsequently, data from 123 different samples of the global anaerobic digestion database was used to determine coverage profiles of host and viral genomes on a broader scale.

CONCLUSIONS

Viruses are key components in anaerobic digestion environments, shaping the microbial guilds which drive the methanogenesis process. In turn, environmental conditions are pivotal in shaping the viral community and the rate of induction of temperate viruses. This study provides an initial insight into the complexity of the anaerobic digestion virome and its relation with the microbial community and the diverse environmental parameters. Video Abstract.