AcetoBase: a functional gene repository and database for formyltetrahydrofolate synthetase sequences

Salinity Impacts the Functional mcrA and dsrA Gene Abundances in Everglades Marshes

Coastal wetlands, such as the Everglades, are increasingly being exposed to stressors that have the potential to modify their existing ecological processes because of global climate change. Their soil microbiomes include a population of organisms important for biogeochemical cycling, but continual stresses can disturb the community's composition, causing functional changes. The Everglades feature wetlands with varied salinity levels, implying that they contain microbial communities with a variety of salt tolerances and microbial functions. Therefore, tracking the effects of stresses on these populations in freshwater and brackish marshes is critical. The study addressed this by utilizing next generation sequencing (NGS) to construct a baseline soil microbial community. The carbon and sulfur cycles were studied by sequencing a microbial functional gene involved in each process, the mcrA and dsrA functional genes, respectively. Saline was introduced over two years to observe the taxonomic alterations that occurred after a long-term disturbance such as seawater intrusion. It was observed that saltwater dosing increased sulfite reduction in freshwater peat soils and decreased methylotrophy in brackish peat soils. These findings add to the understanding of microbiomes by demonstrating how changes in soil qualities impact communities both before and after a disturbance such as saltwater intrusion.

Interpreting tree ensemble machine learning models with endoR

Tree ensemble machine learning models are increasingly used in microbiome science as they are compatible with the compositional, high-dimensional, and sparse structure of sequence-based microbiome data. While such models are often good at predicting phenotypes based on microbiome data, they only yield limited insights into how microbial taxa may be associated. We developed endoR, a method to interpret tree ensemble models. First, endoR simplifies the fitted model into a decision ensemble. Then, it extracts information on the importance of individual features and their pairwise interactions, displaying them as an interpretable network. Both the endoR network and importance scores provide insights into how features, and interactions between them, contribute to the predictive performance of the fitted model. Adjustable regularization and bootstrapping help reduce the complexity and ensure that only essential parts of the model are retained. We assessed endoR on both simulated and real metagenomic data. We found endoR to have comparable accuracy to other common approaches while easing and enhancing model interpretation. Using endoR, we also confirmed published results on gut microbiome differences between cirrhotic and healthy individuals. Finally, we utilized endoR to explore associations between human gut methanogens and microbiome components. Indeed, these hydrogen consumers are expected to interact with fermenting bacteria in a complex syntrophic network. Specifically, we analyzed a global metagenome dataset of 2203 individuals and confirmed the previously reported association between Methanobacteriaceae and Christensenellales. Additionally, we observed that Methanobacteriaceae are associated with a network of hydrogen-producing bacteria. Our method accurately captures how tree ensembles use features and interactions between them to predict a response. As demonstrated by our applications, the resultant visualizations and summary outputs facilitate model interpretation and enable the generation of novel hypotheses about complex systems.

Co-substrate composition is critical for enrichment of functional key species and for process efficiency during biogas production from cattle manure

Cattle manure has a low energy content and high fibre and water content, limiting its value for biogas production. Co-digestion with a more energy-dense material can improve the output, but the co-substrate composition that gives the best results in terms of degree of degradation, gas production and digestate quality has not yet been identified. This study examined the effects of carbohydrate, protein and fat as co-substrates for biogas production from cattle manure. Laboratory-scale semi-continuous mesophilic reactors were operated with manure in mono-digestion or in co-digestion with egg albumin, rapeseed oil, potato starch or a mixture of these, and chemical and microbiological parameters were analysed. The results showed increased gas yield for all co-digestion reactors, but only the reactor supplemented with rapeseed oil showed synergistic effects on methane yield. The reactor receiving potato starch indicated improved fibre degradation, suggesting a priming effect by the easily accessible carbon. Both these reactors showed increased species richness and enrichment of key microbial species, such as fat-degrading Syntrophomonadaceae and families known to include cellulolytic bacteria. The addition of albumin promoted enrichment of known ammonia-tolerant syntrophic acetate- and potential propionate-degrading bacteria, but still caused slight process inhibition and less efficient overall degradation of organic matter in general, and of cellulose in particular.

AcetoBase Version 2: a database update and re-analysis of formyltetrahydrofolate synthetase amplicon sequencing data from anaerobic digesters

AcetoBase is a public repository and database of formyltetrahydrofolate synthetase (FTHFS) sequences. It is the first systematic collection of bacterial FTHFS nucleotide and protein sequences from genomes and metagenome-assembled genomes and of sequences generated by clone library sequencing. At its publication in 2019, AcetoBase (Version 1) was also the first database to establish connections between the FTHFS gene, the Wood–Ljungdahl pathway and 16S ribosomal RNA genes. Since the publication of AcetoBase, there have been significant improvements in the taxonomy of many bacterial lineages and accessibility/availability of public genomics and metagenomics data. The update to the AcetoBase reference database described here (Version 2) provides new sequence data and taxonomy, along with improvements in web functionality and user interface. The evaluation of this latest update by re-analysis of publicly accessible FTHFS amplicon sequencing data previously analysed with AcetoBase Version 1 revealed significant improvements in the taxonomic assignment of FTHFS sequences.

Database URL: https://acetobase.molbio.slu.se

Characterisation of cellulose-degrading organisms in an anaerobic digester

The recalcitrant nature of lignocellulosic biomass hinders efficient exploitation of this fraction for energy production. A better understanding of the microorganisms able to convert plant-based feedstocks is needed to improve anaerobic digestion of lignocellulosic biomass. In this study, active thermophilic cellulose-degrading microorganisms were identified from a full-scale anaerobic digester fed with maize by using metagenome-resolved protein stable isotope probing (protein-SIP). ¹³C-cellulose was converted into ¹³C-methane with a ^13/12C isotope ratio of 0.127 after two days of incubation. Metagenomic analysis revealed 238 different genes coding for carbohydrate-active enzymes (CAZymes), six of which were directly associated with cellulose degradation. The protein-SIP analysis identified twenty heavily labelled peptides deriving from microorganisms actively assimilating labelled carbon from the degradation of ¹³C-cellulose, highlighting several members of the order Clostridiales. Corynebacterium was identified through CAZyme screening, amplicon analysis, and in the metagenome giving a strong identification of being a cellulose degrader.

Syntrophic propionate-oxidizing bacteria in methanogenic systems

The mutual nutritional cooperation underpinning syntrophic propionate degradation provides a scant amount of energy for the microorganisms involved, so propionate degradation often acts as a bottleneck in methanogenic systems. Understanding the ecology, physiology and metabolic capacities of syntrophic propionate-oxidizing bacteria (SPOB) is of interest in both engineered and natural ecosystems, as it offers prospects to guide further development of technologies for biogas production and biomass-derived chemicals, and is important in forecasting contributions by biogenic methane emissions to climate change. SPOB are distributed across different phyla. They can exhibit broad metabolic capabilities in addition to syntrophy (e.g. fermentative, sulfidogenic and acetogenic metabolism) and demonstrate variations in interplay with cooperating partners, indicating nuances in their syntrophic lifestyle. In this review, we discuss distinctions in gene repertoire and organization for the methylmalonyl-CoA pathway, hydrogenases and formate dehydrogenases, and emerging facets of (formate/hydrogen/direct) electron transfer mechanisms. We also use information from cultivations, thermodynamic calculations and omic analyses as the basis for identifying environmental conditions governing propionate oxidation in various ecosystems. Overall, this review improves basic and applied understanding of SPOB and highlights knowledge gaps, hopefully encouraging future research and engineering on propionate metabolism in biotechnological processes.

Miniphocaeibacter halophilus sp. nov., an ammonium-tolerant acetate-producing bacterium isolated from a biogas system

An anaerobic bacterial strain, designated AMB_01T, recovered from mesophilic propionate enrichment of a high-ammonia biogas digester, was characterised using phenotypic and molecular taxonomic methods. Cells of AMB_01T are coccus-shaped and often occur arranged as diplococci or sarcina. Growth occurred at 20-45 °C, initial pH 5.5-8.5 and with up to 0.7 M NH4Cl, with optimum growth at 37-42 °C and pH 8.0. AMB_01T achieved high cell density and highest acetate production when grown on carbohydrates, including monomers, disaccharides and polysaccharides, such as glucose, maltose, cellobiose and starch. The strain was also able to use amino acids and some organic acids and alcoholic compounds for growth. Acetate was formed as the main product and yeast was not required for growth. The major cellular fatty acids were summed feature 4 (iso-C17 : 1I and/or anteiso-C17 : 1B), C18 : 1ω7, C14 : 0, C16 : 0 and summed feature 3 (C16 : 1ω7 and/or iso-C15 : 0 2OH). The highest 16S rRNA gene sequence similarity found was with Miniphocaeibacter massiliensis (96.6 %), within the family Peptoniphilaceae, phylum Bacillota (Firmicutes). The genomic DNA G+C content was 29.0 mol%. An almost complete set of genes for the acetyl-CoA pathway was found. Genome comparisons between AMB_01T and close relatives showed highest digital DNA-DNA hybridisation to Finegoldia magna (23 %), highest average nucleotide identity with genome nucleotide and amino acid sequences to M. massiliensis (72 and 73 %, respectively) and highest average nucleotide identity (87 %) with Schnuerera ultunensis, indicating that AMB_01T represents a novel species. Analysis of genomic, chemotaxonomic, biochemical and physiological data confirmed that strain AMB_01T represents a novel species, for which the name Miniphocaeibacter halophilus sp. nov. is proposed. The type strain is AMB_01T (=DSM 110247T=JCM 39107 T).

Knowing the enemy: homoacetogens in hydrogen production reactors

One of the bottlenecks of the hydrogen production by dark fermentation is the low yields obtained because of the homoacetogenesis persistence, a metabolic pathway where H₂ and CO₂ are consumed to produce acetate. The central reactions of H₂ production and homoacetogenesis are catalyzed by enzyme hydrogenase and the formyltetrahydrofolate synthetase, respectively. In this work, genes encoding for the formyltetrahydrofolate synthetase (fthfs) and hydrogenase (hydA) were used to investigate the diversity of homoacetogens as well as their phylogenetic relationships through quantitative PCR (qPCR) and next-generation amplicon sequencing. A total of 70 samples from 19 different H₂-producing bioreactors with different configurations and operating conditions were analyzed. Quantification through qPCR showed that the abundance of fthfs and hydA was strongly associated with the type of substrate, organic loading rate, and H₂ production performance. In particular, fthfs sequencing revealed that homoacetogens diversity was low with one or two dominant homoacetogens in each sample. Clostridium carboxivorans was detected in the reactors fed with agave hydrolisates; Acetobacterium woodii dominated in systems fed with glucose; Blautia coccoides and unclassified Sporoanaerobacter species were present in reactors fed with cheese whey; finally, Eubacterium limosum and Selenomonas sp. were co-dominant in reactors fed with glycerol. Altogether, quantification and sequencing analysis revealed that the occurrence of homoacetogenesis could take place due to (1) metabolic changes of H₂-producing bacteria towards homoacetogenesis or (2) the displacement of H₂-producing bacteria by homoacetogens. Overall, it was demonstrated that the fthfs gene was a suitable marker to investigate homoacetogens in H₂-producing reactors. KEY POINTS: • qPCR and sequencing analysis revealed two homoacetogenesis phenomena. • fthfs gene was a suitable marker to investigate homoacetogens in H2 reactors.

Microbiological Surveillance of Biogas Plants: Targeting Acetogenic Community

Acetogens play a very important role in anaerobic digestion and are essential in ensuring process stability. Despite this, targeted studies of the acetogenic community in biogas processes remain limited. Some efforts have been made to identify and understand this community, but the lack of a reliable molecular analysis strategy makes the detection of acetogenic bacteria tedious. Recent studies suggest that screening of bacterial genetic material for formyltetrahydrofolate synthetase (FTHFS), a key marker enzyme in the Wood-Ljungdahl pathway, can give a strong indication of the presence of putative acetogens in biogas environments. In this study, we applied an acetogen-targeted analyses strategy developed previously by our research group for microbiological surveillance of commercial biogas plants. The surveillance comprised high-throughput sequencing of FTHFS gene amplicons and unsupervised data analysis with the AcetoScan pipeline. The results showed differences in the acetogenic community structure related to feed substrate and operating parameters. They also indicated that our surveillance method can be helpful in the detection of community changes before observed changes in physico-chemical profiles, and that frequent high-throughput surveillance can assist in management towards stable process operation, thus improving the economic viability of biogas plants. To our knowledge, this is the first study to apply a high-throughput microbiological surveillance approach to visualise the potential acetogenic population in commercial biogas digesters.

Profiling temporal dynamics of acetogenic communities in anaerobic digesters using next-generation sequencing and T-RFLP

Acetogens play a key role in anaerobic degradation of organic material and in maintaining biogas process efficiency. Profiling this community and its temporal changes can help evaluate process stability and function, especially under disturbance/stress conditions, and avoid complete process failure. The formyltetrahydrofolate synthetase (FTHFS) gene can be used as a marker for acetogenic community profiling in diverse environments. In this study, we developed a new high-throughput FTHFS gene sequencing method for acetogenic community profiling and compared it with conventional terminal restriction fragment length polymorphism of the FTHFS gene, 16S rRNA gene-based profiling of the whole bacterial community, and indirect analysis via 16S rRNA profiling of the FTHFS gene-harbouring community. Analyses and method comparisons were made using samples from two laboratory-scale biogas processes, one operated under stable control and one exposed to controlled overloading disturbance. Comparative analysis revealed satisfactory detection of the bacterial community and its changes for all methods, but with some differences in resolution and taxonomic identification. FTHFS gene sequencing was found to be the most suitable and reliable method to study acetogenic communities. These results pave the way for community profiling in various biogas processes and in other environments where the dynamics of acetogenic bacteria have not been well studied.

nZVI Impacts Substrate Conversion and Microbiome Composition in Chain Elongation From D- and L-Lactate Substrates

Medium-chain carboxylates (MCC) derived from biomass biorefining are attractive biochemicals to uncouple the production of a wide array of products from the use of non-renewable sources. Biological conversion of biomass-derived lactate during secondary fermentation can be steered to produce a variety of MCC through chain elongation. We explored the effects of zero-valent iron nanoparticles (nZVI) and lactate enantiomers on substrate consumption, product formation and microbiome composition in batch lactate-based chain elongation. In abiotic tests, nZVI supported chemical hydrolysis of lactate oligomers present in concentrated lactic acid. In fermentation experiments, nZVI created favorable conditions for either chain-elongating or propionate-producing microbiomes in a dose-dependent manner. Improved lactate conversion rates and n-caproate production were promoted at 0.5-2 g nZVI⋅L^-1 while propionate formation became relevant at ≥ 3.5 g nZVI⋅L^-1. Even-chain carboxylates (n-butyrate) were produced when using enantiopure and racemic lactate with lactate conversion rates increased in nZVI presence (1 g⋅L^-1). Consumption of hydrogen and carbon dioxide was observed late in the incubations and correlated with acetate formation or substrate conversion to elongated products in the presence of nZVI. Lactate racemization was observed during chain elongation while isomerization to D-lactate was detected during propionate formation. Clostridium luticellarii, Caproiciproducens, and Ruminococcaceae related species were associated with n-valerate and n-caproate production while propionate was likely produced through the acrylate pathway by Clostridium novyi. The enrichment of different potential n-butyrate producers (Clostridium tyrobutyricum, Lachnospiraceae, Oscillibacter, Sedimentibacter) was affected by nZVI presence and concentrations. Possible theories and mechanisms underlying the effects of nZVI on substrate conversion and microbiome composition are discussed. An outlook is provided to integrate (bio)electrochemical systems to recycle (n)ZVI and provide an alternative reducing power agent as durable control method.

Sugarcane vinasse extreme thermophilic digestion: a glimpse on biogas free management

The high temperature in which sugarcane vinasse (SV) is generated (~ 90 °C) and the positive effect of higher temperatures in biochemical reactions have motivated the evaluation of SV anaerobic digestion (AD) under extreme temperature conditions. Two-stage (acidogenic/methanogenic) and single-stage (methanogenic) AD of SV were evaluated under 70 °C in structured-bed reactors. The extreme temperature was beneficial to the acidogenic step of the two-stage AD process. The methane production, however, was hindered at 70 °C. The VMP of the single and two-stage reactors accounted, respectively, for only 13% and 7% of the production rate reported in sugarcane vinasse AD at 55 °C. At 70 °C, the main genera responsible for methane production was Methanothermobacter and the acetoclastic methanogenesis did not occur, resulting in acetic acid build up (15,800 mg L^-1). These results brought a new perspective for sugarcane vinasse management, with acetic acid production alternatively to methanization. In this perspective, two-stage process would be composed of acidogenic and acetogenic reactors, and beyond acetate, hydrogen and other soluble compounds could be recovered in a complete biorefinery process.

High-Throughput Sequencing and Unsupervised Analysis of Formyltetrahydrofolate Synthetase (FTHFS) Gene Amplicons to Estimate Acetogenic Community Structure

The formyltetrahydrofolate synthetase (FTHFS) gene is a molecular marker of choice to study the diversity of acetogenic communities. However, current analyses are limited due to lack of a high-throughput sequencing approach for FTHFS gene amplicons and a dedicated bioinformatics pipeline for data analysis, including taxonomic annotation and visualization of the sequence data. In the present study, we combined the barcode approach for multiplexed sequencing with unsupervised data analysis to visualize acetogenic community structure. We used samples from a biogas digester to develop proof-of-principle for our combined approach. We successfully generated high-throughput sequence data for the partial FTHFS gene and performed unsupervised data analysis using the novel bioinformatics pipeline “AcetoScan” presented in this study, which resulted in taxonomically annotated OTUs, phylogenetic tree, abundance plots and diversity indices. The results demonstrated that high-throughput sequencing can be used to sequence the FTHFS amplicons from a pool of samples, while the analysis pipeline AcetoScan can be reliably used to process the raw sequence data and visualize acetogenic community structure. The method and analysis pipeline described in this paper can assist in the identification and quantification of known or potentially new acetogens. The AcetoScan pipeline is freely available at https://github.com/abhijeetsingh1704/AcetoScan.