Integrating metagenomic binning with flux balance analysis to unravel syntrophies in anaerobic CO2 methanation

Background Carbon fixation through biological methanation has emerged as a promising technology to produce renewable energy in the context of the circular economy. The anaerobic digestion microbiome is the fundamental biological system operating biogas upgrading and is paramount in power-to-gas conversion. Carbon dioxide (CO2) methanation is frequently performed by microbiota attached to solid supports generating biofilms. Despite the apparent simplicity of the microbial community involved in biogas upgrading, the dynamics behind most of the interspecies interaction remain obscure. To understand the role of the microbial species in CO2 fixation, the biofilm generated during the biogas upgrading process has been selected as a case study. The present work investigates via genome-centric metagenomics, based on a hybrid Nanopore-Illumina approach the biofilm developed on the diffusion devices of four ex situ biogas upgrading reactors. Moreover, genome-guided metabolic reconstruction and flux balance analysis were used to propose a biological role for the dominant microbes. Results The combined microbiome was composed of 59 species, with five being dominant (> 70% of total abundance); the metagenome-assembled genomes representing these species were refined to reach a high level of completeness. Genome-guided metabolic analysis appointed Firmicutes sp. GSMM966 as the main responsible for biofilm formation. Additionally, species interactions were investigated considering their co-occurrence in 134 samples, and in terms of metabolic exchanges through flux balance simulation in a simplified medium. Some of the most abundant species (e.g., Limnochordia sp. GSMM975) were widespread (~ 67% of tested experiments), while others (e.g., Methanothermobacter wolfeii GSMM957) had a scattered distribution. Genome-scale metabolic models of the microbial community were built with boundary conditions taken from the biochemical data and showed the presence of a flexible interaction network mainly based on hydrogen and carbon dioxide uptake and formate exchange. Conclusions Our work investigated the interplay between five dominant species within the biofilm and showed their importance in a large spectrum of anaerobic biogas reactor samples. Flux balance analysis provided a deeper insight into the potential syntrophic interaction between species, especially Limnochordia sp. GSMM975 and Methanothermobacter wolfeii GSMM957. Finally, it suggested species interactions to be based on formate and amino acids exchanges.