Construction of multiple metagenome assembled genomes containing carbon monoxide dehydrogenases from anaerobic carbon monoxide enrichment cultures

Despite its toxicity to many organisms, including most prokaryotes, carbon monoxide (CO) is utilized by some aerobic and anaerobic prokaryotes. Hydrogenogenic CO utilizers employ carbon monoxide dehydrogenase (CODH) and energy-converting hydrogenase (ECH) to oxidize CO and reduce protons to produce H2. Those prokaryotes constitute a rare biosphere and are difficult to detect even with PCR amplification and with metagenomic analyses. In this study, anaerobic CO-enrichment cultures followed by construction of metagenome assembled genomes (MAGs) detected high-quality MAGs from potential hydrogenogenic CO utilizers. Of 32 MAGs constructed, 5 were potential CO utilizer harboring CODH genes. Of the five MAGs, two were classified into the genus Thermolithobacter on the basis of 16S rRNA sequence identity, related to Carboxydocella tharmautotrophica 41, with an average nucleotide identity (ANI) of approximately 72%. Additionally, two were related to Geoglobus acetivorans with ANI values ranging from 75 to 77% to G. acetivorans SBH6, and one MAG was identified as Desulfotomaculum kuznetsovii with an ANI > 96% to D. kuznetsovii DSM 6115. The two Thermolithobacter MAGs identified in this study contained CODH-ECH gene clusters, and were therefore identified as potential hydrogenogenic CO utilizers. However, these MAGs harbored three CODH gene clusters that showed distinct physiological functions in addition to CODH-ECH gene clusters. In total, the five potential CO utilizer MAGs contained sixteen CODH genes. Among those CODHs, four sets did not cluster with any known CODH protein sequences (with an identity of > 90%), and the CODH database was expanded.

Linking DOM characteristics to microbial community: The potential role of DOM mineralization for arsenic release in shallow groundwater

Dissolved organic matter (DOM) play critical roles in arsenic (As) biotransformation in groundwater, but its compositional characteristics and interactions with indigenous microbial communities remain unclear. In this study, DOM signatures coupled with taxonomy and functions of microbial community were characterized in As-enriched groundwater by excitation-emission matrix, Fourier transform ion cyclotron resonance mass spectrometry and metagenomic sequencing. Results showed that As concentrations were significantly positively correlated with DOM humification (r = 0.707, p < 0.01) and the most dominant humic acid-like DOM components (r = 0.789, p < 0.01). Molecular characterization further demonstrated high DOM oxidation degree, with the prevalence of unsaturated oxygen-low aromatics, nitrogen (N₁/N₂)-containing compounds and unique CHO molecules in high As groundwater. These DOM properties were consistent with microbial composition and functional potentials. Both taxonomy and binning analyses demonstrated the dominance of Pseudomonas stutzeri, Microbacterium and Sphingobium xenophagum in As-enriched groundwater which possessed abundant As-reducing gene, with organic carbon degrading genes capable of labile to recalcitrant compounds degradation and high potentials of organic nitrogen mineralization to generate ammonium. Besides, most assembled bins in high As groundwater presented strong fermentation potentials which could facilitate carbon utilization by heterotrophic microbes. This study provides better insight into the potential role of DOM mineralization for As release in groundwater system.

Genome-centric metagenomic insights into the role of Chloroflexi in anammox, activated sludge and methanogenic reactors

BACKGROUND

The phylum Chloroflexi is highly abundant in a wide variety of wastewater treatment bioreactors. It has been suggested that they play relevant roles in these ecosystems, particularly in degrading carbon compounds and on structuring flocs or granules. Nevertheless, their function is not yet well understood as most species have not been isolated in axenic cultures. Here we used a metagenomic approach to investigate Chloroflexi diversity and their metabolic potential in three environmentally different bioreactors: a methanogenic full-scale reactor, a full-scale activated sludge reactor and a lab scale anammox reactor.

RESULTS

Differential coverage binning approach was used to assemble the genomes of 17 new Chloroflexi species, two of which are proposed as new Candidatus genus. In addition, we recovered the first representative genome belonging to the genus 'Ca. Villigracilis'. Even though samples analyzed were collected from bioreactors operating under different environmental conditions, the assembled genomes share several metabolic features: anaerobic metabolism, fermentative pathways and several genes coding for hydrolytic enzymes. Interestingly, genome analysis from the anammox reactor indicated a putative role of Chloroflexi in nitrogen conversion. Genes related to adhesiveness and exopolysaccharides production were also detected. Complementing sequencing analysis, filamentous morphology was detected by Fluorescent in situ hybridization.

CONCLUSION

Our results suggest that Chloroflexi participate in organic matter degradation, nitrogen removal and biofilm aggregation, playing different roles according to the environmental conditions.

Insights into carbon-fixation pathways through metagonomics in the sediments of deep-sea cold seeps

Carbon fixation by chemoautotrophic microorganisms in the dark ocean has a major impact on global carbon cycling and ecological relationships in the ocean's interior. At present, six pathways of autotrophic carbon fixation have been found: the Calvin cycle, the reductive Acetyl-CoA or Wood-Ljungdahl pathway (rAcCoA), the reductive tricarboxylic acid cycle (rTCA), the 3-hydroxypropionate bicycle (3HP), the 3-hydroxypropionate/4-hydroxybutyrate cycle (3HP/4HB), and the dicarboxylate/4-hydroxybutyrate cycle (DC/4HB). Although our knowledge about carbon fixation pathways in the ocean has increased significantly, carbon fixation pathways in the cold seeps are still unknown. In this study, we collected sediment samples from two cold seeps and one trough in the south China sea (SCS), and investigated with metagenomic and metagenome assembled genomes (MAGs). We found that six autotrophic carbon fixation pathways present in the cold seeps and trough with rTCA cycle was the most common pathway, whose genes were particularly high in the cold seeps and increased with sediment depths; the rAcCoA cycle mainly occurred in the cold seep regions, and the abundance of module genes increased with sediment depths. We also elucidated members of chemoautotrophic microorganisms involved in these six carbon-fixation pathways. The rAcCoA, rTCA and DC/4-HB cycles required significantly less energy probably play an important role in the deep-sea environments, especially in the cold seeps. This study provided metabolic insights into the carbon fixation pathways in the cold seeps, and laid the foundation for future detailed study on processes and rates of carbon fixation in the deep-sea ecosystems.

Ecogenomics and metabolic potential of the South Atlantic Ocean microbiome

The unique combination of depth, salinity, and water masses make the South Atlantic Ocean an ecosystem of special relevance within the global ocean. Yet, the microbiome of this ecosystem has received less attention than other regions of the global Ocean. This has hampered our understanding of the diversity and metabolic potential of the microorganisms that dwell in this habitat. To fill this knowledge gap, we analyzed a collection of 31 metagenomes from the Atlantic Ocean that spanned the epipelagic, mesopelagic and bathypelagic zones (surface to 4000 m). Read-centric and gene-centric analysis revealed the unique taxonomic and functional composition of metagenomes from each depth zone, which was driven by differences in physical and chemical parameters. In parallel, a total of 40 metagenome-assembled genomes were obtained, which recovered one third of the total community. Phylogenomic reconstruction revealed that many of these genomes are derived from poorly characterized taxa of Bacteria and Archaea. Genomes derived from heterotrophic bacteria of the aphotic zone displayed a large apparatus of genes suited for the utilization of recalcitrant organic compounds such as cellulose, chitin and alkanes. In addition, we found genomic evidence suggesting that mixotrophic bacteria from the bathypelagic zone could perform carbon fixation through the Calvin-Benson-Bassham cycle, fueled by sulfur oxidation. Finally, we found that the viral communities shifted throughout the water column regarding their targeted hosts and virus-to-microbe ratio, in response to shifts in the composition and functioning their microbial counterparts. Our findings shed light on the microbial and viral drivers of important biogeochemical processes that take place in the South Atlantic Ocean.

Estimating the quality of eukaryotic genomes recovered from metagenomic analysis with EukCC

Microbial eukaryotes constitute a significant fraction of biodiversity and have recently gained more attention, but the recovery of high-quality metagenomic assembled eukaryotic genomes is limited by the current availability of tools. To help address this, we have developed EukCC, a tool for estimating the quality of eukaryotic genomes based on the automated dynamic selection of single copy marker gene sets. We demonstrate that our method outperforms current genome quality estimators, particularly for estimating contamination, and have applied EukCC to datasets derived from two different environments to enable the identification of novel eukaryote genomes, including one from the human skin.

Genome-centric microbiome analysis reveals solid retention time (SRT)-shaped species interactions and niche differentiation in food waste and sludge co-digesters

Co-digestion of food waste with sewage sludge is widely applied for waste stabilization and energy recovery around the world. However, the effect of solid retention time (SRT) on the microbial population dynamics, metabolism and interspecies interaction have not been fully elucidated. Here, the influence of SRTs (5-25 days) on the performance of the co-digestion system was investigated and state-of-the-art genome-centric metagenomic analysis was employed to uncover the dynamics and metabolic network of the key players underlying the well-functioned and poorly-functioned co-digestion microbial communities. The results of the microbial analyses indicated that SRT largely shaped microbial community structure by enriching the syntrophic specialist Syntrophomonas and CO₂/H₂ ( formate)-using methanogen Methanocorpusculum in the well-functioned co-digester operated at SRT of 25 days, while selecting acid-tolerant populations Lactobacillus at SRT of 5 days. The metagenome assembled genomes (MAGs) of key players, such as Syntrophomonadaceae, Methanocorpusculum, and Mesotoga, were retrieved, additionally, the syntrophic acetate oxidation plus hydrogenotrophic methanogenesis (SAO-HM) were proposed as the dominant pathway for methane production. The metabolic interaction in the co-digestion microbial consortia was profiled by assigning MAGs into functional guilds. Functional redundancy was found in the bacterial groups in hydrolysis step, and the members in these groups reduced the direct competition by niche differentiation.

Impact of process temperature and organic loading rate on cellulolytic / hydrolytic biofilm microbiomes during biomethanation of ryegrass silage revealed by genome-centered metagenomics and metatranscriptomics

BACKGROUND

Anaerobic digestion (AD) of protein-rich grass silage was performed in experimental two-stage two-phase biogas reactor systems at low vs. increased organic loading rates (OLRs) under mesophilic (37 °C) and thermophilic (55 °C) temperatures. To follow the adaptive response of the biomass-attached cellulolytic/hydrolytic biofilms at increasing ammonium/ammonia contents, genome-centered metagenomics and transcriptional profiling based on metagenome assembled genomes (MAGs) were conducted.

RESULTS

In total, 78 bacterial and archaeal MAGs representing the most abundant members of the communities, and featuring defined quality criteria were selected and characterized in detail. Determination of MAG abundances under the tested conditions by mapping of the obtained metagenome sequence reads to the MAGs revealed that MAG abundance profiles were mainly shaped by the temperature but also by the OLR. However, the OLR effect was more pronounced for the mesophilic systems as compared to the thermophilic ones. In contrast, metatranscriptome mapping to MAGs subsequently normalized to MAG abundances showed that under thermophilic conditions, MAGs respond to increased OLRs by shifting their transcriptional activities mainly without adjusting their proliferation rates. This is a clear difference compared to the behavior of the microbiome under mesophilic conditions. Here, the response to increased OLRs involved adjusting of proliferation rates and corresponding transcriptional activities. The analysis led to the identification of MAGs positively responding to increased OLRs. The most outstanding MAGs in this regard, obviously well adapted to higher OLRs and/or associated conditions, were assigned to the order Clostridiales (Acetivibrio sp.) for the mesophilic biofilm and the orders Bacteroidales (Prevotella sp. and an unknown species), Lachnospirales (Herbinix sp. and Kineothrix sp.) and Clostridiales (Clostridium sp.) for the thermophilic biofilm. Genome-based metabolic reconstruction and transcriptional profiling revealed that positively responding MAGs mainly are involved in hydrolysis of grass silage, acidogenesis and / or acetogenesis.

CONCLUSIONS

An integrated -omics approach enabled the identification of new AD biofilm keystone species featuring outstanding performance under stress conditions such as increased OLRs. Genome-based knowledge on the metabolic potential and transcriptional activity of responsive microbiome members will contribute to the development of improved microbiological AD management strategies for biomethanation of renewable biomass.

Ecogenomic characterization of widespread, closely-related SAR11 clades of the freshwater genus "Candidatus Fonsibacter" and proposal of Ca. Fonsibacter lacus sp. nov

The ubiquitous alpha-proteobacteria of the order "Candidatus Pelagibacterales" (SAR11) are highly abundant in aquatic environments, and among them, members of the monophyletic lineage LD12 (also known as SAR11 clade IIIb) are specifically found in lacustrine ecosystems. Clade IIIb bacteria are some of the most prominent members of freshwater environments, but little is known about their biology due to the lack of genome representatives. Only recently, the first non-marine isolate was cultured and described as "Candidatus Fonsibacter ubiquis". Here, we expand the collection of freshwater IIIb representatives and describe a new IIIb species of the genus "Ca. Fonsibacter". Specifically, we assembled a collection of 67 freshwater metagenomic datasets from the interconnected lakes of the Chattahoochee River basin (GA, USA) and obtained nearly complete metagenome-assembled genomes (MAGs) representing 5 distinct IIIb subclades, roughly equivalent to species based on genomic standards, including the previously described "Ca. F. ubiquis". Genomic comparisons between members of the IIIb species revealed high similarity in gene content. However, when comparing their abundance profiles in the Chattahoochee basin and various aquatic environments, differences in temporal and spatial distributions among the distinct species were observed implying niche differentiation might be underlying the coexistence of the highly functionally similar representatives. The name Ca. Fonsibacter lacus sp. nov. is proposed for the most abundant and widespread species in the Chattahoochee River basin and various freshwater ecosystems.

Microbial activity response to hydrogen injection in thermophilic anaerobic digesters revealed by genome-centric metatranscriptomics

BACKGROUND

The expansion of renewable energy produced by windmills and photovoltaic panels has generated a considerable electricity surplus, which can be utilized in water electrolysis systems for hydrogen production. The resulting hydrogen can then be funneled to anaerobic digesters for biogas upgrading (biomethanation) purposes (power-to-methane) or to produce high value-added compounds such as short-chain fatty acids (power-to-chemicals). Genome-centric metagenomics and metatranscriptomic analyses were performed to better understand the metabolic dynamics associated with H₂ injection in two different configurations of anaerobic digesters treating acidic wastes, specifically cheese manufacturing byproducts. These approaches revealed the key-genes involved in methanation and carbon fixation pathways at species level.

RESULTS

The biogas upgrading process in the single-stage configuration increased the CH₄ content by 7%. The dominant methanogenic species responsible for the upregulation of the hydrogenotrophic pathway in this reactor was Methanothermobacter wolfeii UC0008. In the two-stage configuration, H₂ injection induced an upregulation of CO₂ fixation pathways producing short-chain fatty acids, mainly acetate and butyrate. In this configuration, the abundant species Anaerobaculum hydrogeniformans UC0046 and Defluviitoga tunisiensis UC0050 primarily upregulated genes related to electron transport chains, suggesting putative syntrophisms with hydrogen scavenger microbes. Interestingly, Tepidanaerobacter acetatoxydans UC0018 did not act as an acetate-oxidizer in either reactor configurations, and instead regulated pathways involved in acetate production and uptake. A putative syntrophic association between Coprothermobacter proteolyticus UC0011 and M. wolfeii UC0008 was proposed in the two-stage reactor. In order to support the transcriptomic findings regarding the hydrogen utilization routes, an advanced bioconversion model was adapted for the simulation of the single- and two-stage reactor setups.

CONCLUSIONS

This is the first study investigating biogas reactor metatranscriptome dynamics following hydrogen injection for biomethanation and carbon fixation to short-chain fatty acids purposes. The same microbes showed different patterns of metabolic regulation in the two reactor configurations. It was observed an effect of the specialized acidogenic reactor on the overall microbial consortium composition and activity in the two-stage digester. There were also suggested the main species responsible for methanation, short-chain fatty acids production, and electron transport chain mechanisms, in both reactor configurations.