A microbial gene catalog of anaerobic digestion from full-scale biogas plants

Metatranscriptomics-guided genome-scale metabolic reconstruction reveals the carbon flux and trophic interaction in methanogenic communities

BACKGROUND

Despite rapid advances in genomic-resolved metagenomics and remarkable explosion of metagenome-assembled genomes (MAGs), the function of uncultivated anaerobic lineages and their interactions in carbon mineralization remain largely uncertain, which has profound implications in biotechnology and biogeochemistry.

RESULTS

In this study, we combined long-read sequencing and metatranscriptomics-guided metabolic reconstruction to provide a genome-wide perspective of carbon mineralization flow from polymers to methane in an anaerobic bioreactor. Our results showed that incorporating long reads resulted in a substantial improvement in the quality of metagenomic assemblies, enabling the effective recovery of 132 high-quality genomes meeting stringent criteria of minimum information about a metagenome-assembled genome (MIMAG). In addition, hybrid assembly obtained 51% more prokaryotic genes in comparison to the short-read-only assembly. Metatranscriptomics-guided metabolic reconstruction unveiled the remarkable metabolic flexibility of several novel Bacteroidales-affiliated bacteria and populations from Mesotoga sp. in scavenging amino acids and sugars. In addition to recovering two circular genomes of previously known but fragmented syntrophic bacteria, two newly identified bacteria within Syntrophales were found to be highly engaged in fatty acid oxidation through syntrophic relationships with dominant methanogens Methanoregulaceae bin.74 and Methanothrix sp. bin.206. The activity of bin.206 preferring acetate as substrate exceeded that of bin.74 with increasing loading, reinforcing the substrate determinantal role.

CONCLUSION

Overall, our study uncovered some key active anaerobic lineages and their metabolic functions in this complex anaerobic ecosystem, offering a framework for understanding carbon transformations in anaerobic digestion. These findings advance the understanding of metabolic activities and trophic interactions between anaerobic guilds, providing foundational insights into carbon flux within both engineered and natural ecosystems. Video Abstract.

MiDAS 5: Global diversity of bacteria and archaea in anaerobic digesters

Anaerobic digestion of organic waste into methane and carbon dioxide (biogas) is carried out by complex microbial communities. Here, we use full-length 16S rRNA gene sequencing of 285 full-scale anaerobic digesters (ADs) to expand our knowledge about diversity and function of the bacteria and archaea in ADs worldwide. The sequences are processed into full-length 16S rRNA amplicon sequence variants (FL-ASVs) and are used to expand the MiDAS 4 database for bacteria and archaea in wastewater treatment systems, creating MiDAS 5. The expansion of the MiDAS database increases the coverage for bacteria and archaea in ADs worldwide, leading to improved genus- and species-level classification. Using MiDAS 5, we carry out an amplicon-based, global-scale microbial community profiling of the sampled ADs using three common sets of primers targeting different regions of the 16S rRNA gene in bacteria and/or archaea. We reveal how environmental conditions and biogeography shape the AD microbiota. We also identify core and conditionally rare or abundant taxa, encompassing 692 genera and 1013 species. These represent 84-99% and 18-61% of the accumulated read abundance, respectively, across samples depending on the amplicon primers used. Finally, we examine the global diversity of functional groups with known importance for the anaerobic digestion process.

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.

Soil microbial community composition and co-occurrence network responses to mild and severe disturbances in volcanic areas

Soil physicochemical properties and vegetation types are the main factors affecting soil microorganisms, but there are few studies on the effects of the disturbance following volcanic eruption. To make up for this lack of knowledge, we used Illumina Miseq high-throughput sequencing to study the characteristics of soil microorganisms on both shores of a volcanically disturbed lake. Soil microorganisms in the two sites were subjected to different degrees of volcanic disturbance and showed significant heterogeneity. Mild volcanic disturbance area had higher enrichment of prokaryotic community. Co-occurrence network analysis showed that a total of 12 keystone taxa (9 prokaryotes and 3 fungi) were identified, suggesting that soil prokaryote may play a more significant role than fungi in overall community structure and function. Compared with severe volcanic disturbance area, the soil microbial community in mild volcanic disturbance area had the higher modular network (0.327 vs 0.291). The competition was stronger (positive/negative link ratio, P/N: 1.422 vs 1.159). Random forest analysis showed that soil superoxide dismutase was the most significant variable associated with soil microbial community. Structural equation model (SEM) results showed that keystone had a directly positive effect on prokaryotic (λ = 0.867, P < 0.001) and fungal (λ = 0.990, P < 0.001) multifunctionality while had also a directly positive effect on fungal diversity (λ = 0.553, P < 0.001), suggesting that keystone taxa played a key role in maintaining ecosystem stability. These results were important for understanding the effects of different levels of volcanic disturbance on soil ecosystems.

Uncovering Microbiome Adaptations in a Full-Scale Biogas Plant: Insights from MAG-Centric Metagenomics and Metaproteomics

The current focus on renewable energy in global policy highlights the importance of methane production from biomass through anaerobic digestion (AD). To improve biomass digestion while ensuring overall process stability, microbiome-based management strategies become more important. In this study, metagenomes and metaproteomes were used for metagenomically assembled genome (MAG)-centric analyses to investigate a full-scale biogas plant consisting of three differentially operated digesters. Microbial communities were analyzed regarding their taxonomic composition, functional potential, as well as functions expressed on the proteome level. Different abundances of genes and enzymes related to the biogas process could be mostly attributed to different process parameters. Individual MAGs exhibiting different abundances in the digesters were studied in detail, and their roles in the hydrolysis, acidogenesis and acetogenesis steps of anaerobic digestion could be assigned. Methanoculleus thermohydrogenotrophicum was an active hydrogenotrophic methanogen in all three digesters, whereas Methanothermobacter wolfeii was more prevalent at higher process temperatures. Further analysis focused on MAGs, which were abundant in all digesters, indicating their potential to ensure biogas process stability. The most prevalent MAG belonged to the class Limnochordia; this MAG was ubiquitous in all three digesters and exhibited activity in numerous pathways related to different steps of AD.

The Potential of Digested Sludge-Assimilating Microflora for Biogas Production from Food Processing Wastes

Food processing wastes (FPWs) are residues generated in food manufacturing, and their composition varies depending on the type of food product being manufactured. Therefore, selecting and acclimatizing seed microflora during the initiation of biogas production is crucial for optimal outcomes. The present study examined the biogas production capabilities of digested sludge-assimilating and biogas-yielding soil (DABYS) and enteric (DABYE) microflorae when used as seed cultures for biogas production from FPWs. After subculturing and feeding these microbial seeds with various FPWs, we assessed their biogas-producing abilities. The subcultures produced biogas from many FPWs, except orange peel, suggesting that the heterogeneity of the bacterial members in the seed microflora facilitates quick adaptation to FPWs. Microflorae fed with animal-derived FPWs contained several methanogenic archaeal families and produced methane. In contrast, microflorae fed with vegetable-, fruit-, and crop-derived FPWs generated hydrogen, and methanogenic archaeal populations were diminished by repeated subculturing. The subcultured microflorae appear to hydrolyze carbohydrates and protein in FPWs using cellulase, pectinase, or protease. Despite needing enhancements in biogas yield for future industrial scale-up, the DABYS and DABYE microflorae demonstrate robust adaptability to various FPWs.

Inter-kingdom interactions and stability of methanogens revealed by machine-learning guided multi-omics analysis of industrial-scale biogas plants

Multi-omics analysis is a powerful tool for the detection and study of inter-kingdom interactions, such as those between bacterial and archaeal members of complex biogas-producing microbial communities. In the present study, the microbiomes of three industrial-scale biogas digesters, each fed with different substrates, were analysed using a machine-learning guided genome-centric metagenomics framework complemented with metatranscriptome data. This data permitted us to elucidate the relationship between abundant core methanogenic communities and their syntrophic bacterial partners. In total, we detected 297 high-quality, non-redundant metagenome-assembled genomes (nrMAGs). Moreover, the assembled 16 S rRNA gene profiles of these nrMAGs showed that the phylum Firmicutes possessed the highest copy number, while the representatives of the archaeal domain had the lowest. Further investigation of the three anaerobic microbial communities showed characteristic alterations over time but remained specific to each industrial-scale biogas plant. The relative abundance of various microorganisms as revealed by metagenome data was independent from corresponding metatranscriptome activity data. Archaea showed considerably higher activity than was expected from their abundance. We detected 51 nrMAGs that were present in all three biogas plant microbiomes with different abundances. The core microbiome correlated with the main chemical fermentation parameters, and no individual parameter emerged as a predominant shaper of community composition. Various interspecies H₂/electron transfer mechanisms were assigned to hydrogenotrophic methanogens in the biogas plants that ran on agricultural biomass and wastewater. Analysis of metatranscriptome data revealed that methanogenesis pathways were the most active of all main metabolic pathways.

Effect of Alkaline and Mechanical Pretreatment of Wheat Straw on Enrichment Cultures from Pachnoda marginata Larva Gut

In order to partially mimic the efficient lignocellulose pretreatment process performed naturally in the gut system of Pachnoda marginata larvae, two wheat straw pretreatments were evaluated: a mechanical pretreatment via cutting the straw into two different sizes and an alkaline pretreatment with calcium hydroxide. After pretreatment, gut enrichment cultures on wheat straw at alkaline pH were inoculated and kept at mesophilic conditions over 45 days. The methanogenic community was composed mainly of the Methanomicrobiaceae and Methanosarcinaceae families. The combined pretreatment, size reduction and alkaline pretreatment, was the best condition for methane production. The positive effect of the straw pretreatment was higher in the midgut cultures, increasing the methane production by 192%, while for hindgut cultures the methane production increased only by 149% when compared to non-pretreated straw. Scanning electron microscopy (SEM) showed that the alkaline pretreatment modified the surface of the wheat straw fibers, which promoted biofilm formation and microbial growth. The enrichment cultures derived from larva gut microbiome were able to degrade larger 1 mm alkaline treated and smaller 250 µm but non-pretreated straw at the same efficiency. The combination of mechanical and alkaline pretreatments resulted in increased, yet not superimposed, methane yield.

Recovery of metagenome-assembled microbial genomes from a full-scale biogas plant of food waste by pacific biosciences high-fidelity sequencing

Background

Anaerobic digestion (AD) is important in treating of food waste, and thousands of metagenome-assembled genomes (MAGs) have been constructed for the microbiome in AD. However, due to the limitations of the short-read sequencing and assembly technologies, most of these MAGs are grouped from hundreds of short contigs by binning algorithms, and the errors are easily introduced.

Results

In this study, we constructed a total of 60 non-redundant microbial genomes from 64.5 Gb of PacBio high-fidelity (HiFi) long reads, generated from the digestate samples of a full-scale biogas plant fed with food waste. Of the 60 microbial genomes, all genomes have at least one copy of rRNA operons (16S, 23S, and 5S rRNA), 54 have ≥18 types of standard tRNA genes, and 39 are circular complete genomes. In comparison with the published short-read derived MAGs for AD, we found 23 genomes with average nucleotide identity less than 95% to any known MAGs. Besides, our HiFi-derived genomes have much higher average contig N50 size, slightly higher average genome size and lower contamination. GTDB-Tk classification of these genomes revealed two genomes belonging to novel genus and four genomes belonging to novel species, since their 16S rRNA genes have identities lower than 95 and 97% to any known 16S rRNA genes, respectively. Microbial community analysis based on the these assembled genomes reveals the most predominant phylum was Thermotogae (70.5%), followed by Euryarchaeota (6.1%), and Bacteroidetes (4.7%), and the most predominant bacterial and archaeal genera were Defluviitoga (69.1%) and Methanothrix (5.4%), respectively. Analysis of the full-length 16S rRNA genes identified from the HiFi reads gave similar microbial compositions to that derived from the 60 assembled genomes.

Conclusion

High-fidelity sequencing not only generated microbial genomes with obviously improved quality but also recovered a substantial portion of novel genomes missed in previous short-read based studies, and the novel genomes will deepen our understanding of the microbial composition in AD of food waste.

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.

Adaptation of gut microbiome and host metabolic systems to lignocellulosic degradation in bamboo rats

Bamboo rats (Rhizomys pruinosus) are among the few mammals that lives on a bamboo-based diet which is mainly composed of lignocellulose. However, the mechanisms of adaptation of their gut microbiome and metabolic systems in the degradation of lignocellulose are largely unknown. Here, we conducted a multi-omics analysis on bamboo rats to investigate the interaction between their gut microbiomes and metabolic systems in the pre- and post-weaning periods, and observed significant relationships between dietary types, gut microbiome, serum metabolome and host gene expression. For comparison, published gut microbial data from the famous bamboo-eating giant panda (Ailuropoda melanoleuca) were also used for analysis. We found that the adaptation of the gut microbiome of the bamboo rat to a lignocellulose diet is related to a member switch in the order Bacteroidales from family Bacteroidaceae to family Muribaculaceae, while for the famous bamboo-eating giant panda, several aerobes and facultative anaerobes increase after weaning. The conversion of bacteria with an increased relative abundance in bamboo rats after weaning enriched diverse carbohydrate-active enzymes (CAZymes) associated with lignocellulose degradation and functionally enhanced the biosynthesis of amino acids and B vitamins. Meanwhile, the circulating concentration of short-chain fatty acids (SCFAs) derived metabolites and the metabolic capacity of linoleic acid in the host were significantly elevated. Our findings suggest that fatty acid metabolism, including linoleic acid and SCFAs, are the main energy sources for bamboo rats in response to the low-nutrient bamboo diet.

Metagenomic and HT-qPCR analysis reveal the microbiome and resistome in pig slurry under storage, composting, and anaerobic digestion

Direct application of pig slurry to agricultural land, as a means of nutrient recycling, introduces pathogens, antibiotic resistant bacteria, or genes, to the environment. With global environmental sustainability policies mandating a reduction in synthetic fertilisation and a commitment to a circular economy it is imperative to find effective on-farm treatments of slurry that maximises its fertilisation value and minimises risk to health and the environment. We assessed and compared the effect of storage, composting, and anaerobic digestion (AD) on pig slurry microbiome, resistome and nutrient content. Shotgun metagenomic sequencing and HT-qPCR arrays were implemented to understand the dynamics across the treatments. Our results identified that each treatment methods have advantages and disadvantages in removal pollutants or increasing nutrients. The data suggests that storage and composting are optimal for the removal of human pathogens and anaerobic digestion for the reduction in antibiotic resistance (AMR) genes and mobile genetic elements. The nitrogen content is increased in storage and AD, while reduced in composting. Thus, depending on the requirement for increased or reduced nitrogen the optimum treatment varies. Combining the results indicates that composting provides the greatest gain by reducing risk to human health and the environment. Network analysis revealed reducing Proteobacteria and Bacteroidetes while increasing Firmicutes will reduce the AMR content. KEGG analysis identified no significant change in the pathways across all treatments. This novel study provides a data driven decision tree to determine the optimal treatment for best practice to minimise pathogen, AMR and excess or increasing nutrient transfer from slurry to environment.

Metaproteomics reveals enzymatic strategies deployed by anaerobic microbiomes to maintain lignocellulose deconstruction at high solids

Economically viable production of cellulosic biofuels requires operation at high solids loadings—on the order of 15 wt%. To this end we characterize Nature’s ability to deconstruct and utilize mid-season switchgrass at increasing solid loadings using an anaerobic methanogenic microbiome. This community exhibits undiminished fractional carbohydrate solubilization at loadings ranging from 30 g/L to 150 g/L. Metaproteomic interrogation reveals marked increases in the abundance of specific carbohydrate-active enzyme classes. Significant enrichment of auxiliary activity family 6 enzymes at higher solids suggests a role for Fenton chemistry. Stress-response proteins accompanying these reactions are similarly upregulated at higher solids, as are β-glucosidases, xylosidases, carbohydrate-debranching, and pectin-acting enzymes—all of which indicate that removal of deconstruction inhibitors is important for observed undiminished solubilization. Our work provides insights into the mechanisms by which natural microbiomes effectively deconstruct and utilize lignocellulose at high solids loadings, informing the future development of defined cultures for efficient bioconversion.

Removal of methylmercury and its potential relationship to microbiota in sludge anaerobic digestion under thermal hydrolysis

Reducing health risk of mercury (Hg)/methylmercury (MeHg) in sewage sludge is vital to its land application. This study revealed that thermal hydrolysis reduced MeHg content both during pretreatment process and subsequent anaerobic digestion (AD), which resulted in decrease of MeHg content from 4.24 ng/g to 0.95 ng/g after thermal hydrolysis (150 ℃) and further decreased to 0.39 ng/g after AD. Notably, thermal hydrolysis at high temperature (120 ℃ and 150 ℃) promoted both Hg methylation and MeHg demethylation rather than the control or at low temperature (100 ℃). Hg methylation dominated in hydrolysis and acidogenesis stage, whereas MeHg demethylation dominated in methanogenesis stage. Though abundance of related genes (HgcA and merA) was dramatically reduced, Ruminococcaceae, Peptococcaceae, and Lachnospiraceae were potentially Hg methylators in hydrolysis and acidogenesis stage. Whereas, MeHg demethylation dominated in the late period of AD due to the improved syntrophic methanogenesis and possibly reduced Hg²⁺ biodegradability by precipitation.

Variation of Metagenome From Feedstock to Digestate in Full-Scale Biogas Plants

Anaerobic digestion (AD) has been widely used to resolve the problem of organic wastes worldwide. Previous studies showed that the types of feedstock have a great influence on the AD microbiome, and a huge number of AD populations are migrated from upstream feedstocks. However, the changes of microbial compositions from feedstock to AD digestate are still less understood. We collected feedstock samples from 56 full-scale biogas plants, generated 1,716 Gb feedstock metagenomic data in total, and constructed the first comprehensive microbial gene catalog of feedstock containing 25.2 million genes. Our result indicated that the predominant phyla in feedstock are Firmicutes, Bacteroidetes, and Proteobacteria, which is similar to that in AD digestate, and the microbial diversity of feedstock samples is higher than that of AD digestate samples. In addition, the relative abundance of most genes involved in methanogenesis increase from feedstock to AD digestate. Besides, the amount of antibiotic resistance genes (ARGs) and pathogenic bacteria in AD are effectively reduced compared to feedstocks. This study provides a comprehensive microbial gene catalog of feedstock, and deepens the understanding of variation of microbial communities from feedstock to AD digestate of full-scale AD. The results also suggest the potential of AD to reduce the level of ARGs and pathogens in animal manure.