digestome: a licence-clean marker-gene panel for functional profiling of anaerobic digestion microbiomes
Abstract
Functional profiling of anaerobic digestion microbiomes is routinely performed against KEGG or MetaCyc, whose paid commercial licences block fee-for-service analysis for biogas operators, the setting where the results have the most immediate operational value. We present digestome, a curated marker-gene panel for anaerobic digestion built exclusively on sources that are free for commercial use (NCBIfam, public domain; Pfam, CC0; Rhea, CC BY), together with a scorer that reports pathway completeness, branch capability and, for hydrolysis, whether the enzyme is built for export. Benchmarked against 1,401 metagenome-assembled genomes from 134 anaerobic digesters, the panel assigned no methanogenesis route to any of 1,361 non-methanogens, and every acetoclastic call fell within Methanosarcina or Methanothrix. The result held for 3,043 species representatives from the Genome Taxonomy Database, spanning 203 phyla, which are built differently from binned metagenomes: the only background genomes given a route were two archaea carrying mcrA, and none of the 90 anaerobic methane and alkane oxidisers was given one. Specificity against non-methanogens does not show that a methanogen gets the right route: with a family-level check, 59 of 156 acetoclastic calls in that sample fell on methylotrophic genera that do not use acetate, and a lineage policy removed them, with 14 more left unconfirmed in unnamed genera. Testing whether a catalytic domain shares a polypeptide with an export module reduced the genomes credited with cellulolytic capacity from 597 to 50. A causal graph of the measurement shows that a profile reports the community together with three measurement choices (sampling, DNA extraction and sequencing depth, and reference genomes), which fixes what one profile can attribute and what a prospective validation would have to record.
Related
- Viromic and metagenomic analyses reveal distribution and dissemination of antibiotic resistance genes in anaerobic digestion-composting systems: roles of quorum sensing and viruses
- Response of the bacterial consortium under supersaturation conditions in an anaerobic co-digestion environment between Typha domingensis and bovine manure
- Biochar-enhanced direct interspecies electron transfer in anaerobic digestion: A review
- Coupling of stable isotope probing and metagenomics reveals the diversity of viruses infecting formatotrophs during anaerobic digestion process
- Invasive Chinese tallow is associated with distinct rhizosphere microbial communities, core microbiomes, and functional potential relative to native Loblolly pine
BibTeX
@misc{pall2026digestome,
title = {{digestome: a licence-clean marker-gene panel for functional profiling of anaerobic digestion microbiomes}},
author = {T. Päll},
howpublished = {bioRxiv},
year = {2026},
month = sep,
doi = {10.64898/2026.09.28.755039},
url = {https://www.biorxiv.org/content/10.64898/2026.09.28.755039v1},
note = {CC BY},
}
RIS
TY - UNPB
TI - digestome: a licence-clean marker-gene panel for functional profiling of anaerobic digestion microbiomes
AU - T. Päll
PY - 2026
DA - 2026/09/28
PB - bioRxiv
DO - 10.64898/2026.09.28.755039
UR - https://www.biorxiv.org/content/10.64898/2026.09.28.755039v1
AB - Functional profiling of anaerobic digestion microbiomes is routinely performed against KEGG or MetaCyc, whose paid commercial licences block fee-for-service analysis for biogas operators, the setting where the results have the most immediate operational value. We present digestome, a curated marker-gene panel for anaerobic digestion built exclusively on sources that are free for commercial use (NCBIfam, public domain; Pfam, CC0; Rhea, CC BY), together with a scorer that reports pathway completeness, branch capability and, for hydrolysis, whether the enzyme is built for export. Benchmarked against 1,401 metagenome-assembled genomes from 134 anaerobic digesters, the panel assigned no methanogenesis route to any of 1,361 non-methanogens, and every acetoclastic call fell within Methanosarcina or Methanothrix. The result held for 3,043 species representatives from the Genome Taxonomy Database, spanning 203 phyla, which are built differently from binned metagenomes: the only background genomes given a route were two archaea carrying mcrA, and none of the 90 anaerobic methane and alkane oxidisers was given one. Specificity against non-methanogens does not show that a methanogen gets the right route: with a family-level check, 59 of 156 acetoclastic calls in that sample fell on methylotrophic genera that do not use acetate, and a lineage policy removed them, with 14 more left unconfirmed in unnamed genera. Testing whether a catalytic domain shares a polypeptide with an export module reduced the genomes credited with cellulolytic capacity from 597 to 50. A causal graph of the measurement shows that a profile reports the community together with three measurement choices (sampling, DNA extraction and sequencing depth, and reference genomes), which fixes what one profile can attribute and what a prospective validation would have to record.
ER -