Methanogenesis byMethanosarcina acetivoransinvolves two structurally and functionally distinct classes of heterodisulfide reductase

Physiological and transcriptomic response to methyl-coenzyme M reductase limitation in Methanosarcina acetivorans

Methyl-coenzyme M reductase (MCR) catalyzes the final step of methanogenesis, the microbial metabolism responsible for nearly all biological methane emissions to the atmosphere. Decades of biochemical and structural research studies have generated detailed insights into MCR function in vitro, yet very little is known about the interplay between MCR and methanogen physiology. For instance, while it is routinely stated that MCR catalyzes the rate-limiting step of methanogenesis, this has not been categorically tested. In this study, to gain a more direct understanding of MCR's control on the growth of Methanosarcina acetivorans, we generate a strain with an inducible mcr operon on the chromosome, allowing for careful control of MCR expression. We show that MCR is not growth rate-limiting in substrate-replete batch cultures. However, through careful titration of MCR expression, growth-limiting state(s) can be obtained. Transcriptomic analysis of M. acetivorans experiencing MCR limitation reveals a global response with hundreds of differentially expressed genes across diverse functional categories. Notably, MCR limitation leads to strong induction of methylsulfide methyltransferases, likely due to insufficient recycling of metabolic intermediates. In addition, the mcr operon is not transcriptionally regulated, i.e., it is constitutively expressed, suggesting that the overabundance of MCR might be beneficial when cells experience nutrient limitation or stressful conditions. Altogether, we show that there is a wide range of cellular MCR concentrations that can sustain optimal growth, suggesting that other factors such as anabolic reactions might be rate-limiting for methanogenic growth.

IMPORTANCE

Methane is a potent greenhouse gas that has contributed to ca. 25% of global warming in the post-industrial era. Atmospheric methane is primarily of biogenic origin, mostly produced by microorganisms called methanogens. Methyl-coenzyme M reductase (MCR) catalyzes methane formatio in methanogens. Even though MCR comprises ca. 10% of the cellular proteome, it is hypothesized to be growth-limiting during methanogenesis. In this study, we show that Methanosarcina acetivorans cells grown in substrate-replicate batch cultures produce more MCR than its cellular demand for optimal growth. The tools outlined in this study can be used to refine metabolic models of methanogenesis and assay lesions in MCR in a higher-throughput manner than isolation and biochemical characterization of pure protein.

Fusion/fission protein family identification in Archaea

The majority of newly discovered archaeal lineages remain without a cultivated representative, but scarce experimental data from the cultivated organisms show that they harbor distinct functional repertoires. To unveil the ecological as well as evolutionary impact of Archaea from metagenomics, new computational methods need to be developed, followed by in-depth analysis. Among them is the genome-wide protein fusion screening performed here. Natural fusions and fissions of genes not only contribute to microbial evolution but also complicate the correct identification and functional annotation of sequences. The products of these processes can be defined as fusion (or composite) proteins, the ones consisting of two or more domains originally encoded by different genes and split proteins, and the ones originating from the separation of a gene in two (fission). Fusion identifications are required for proper phylogenetic reconstructions and metabolic pathway completeness assessments, while mappings between fused and unfused proteins can fill some of the existing gaps in metabolic models. In the archaeal genome-wide screening, more than 1,900 fusion/fission protein clusters were identified, belonging to both newly sequenced and well-studied lineages. These protein families are mainly associated with different types of metabolism, genetic, and cellular processes. Moreover, 162 of the identified fusion/fission protein families are archaeal specific, having no identified fused homolog within the bacterial domain. Our approach was validated by the identification of experimentally characterized fusion/fission cases. However, around 25% of the identified fusion/fission families lack functional annotations for both composite and split states, showing the need for experimental characterization in Archaea.IMPORTANCEGenome-wide fusion screening has never been performed in Archaea on a broad taxonomic scale. The overlay of multiple computational techniques allows the detection of a fine-grained set of predicted fusion/fission families, instead of rough estimations based on conserved domain annotations only. The exhaustive mapping of fused proteins to bacterial organisms allows us to capture fusion/fission families that are specific to archaeal biology, as well as to identify links between bacterial and archaeal lineages based on cooccurrence of taxonomically restricted proteins and their sequence features. Furthermore, the identification of poorly characterized lineage-specific fusion proteins opens up possibilities for future experimental and computational investigations. This approach enhances our understanding of Archaea in general and provides potential candidates for in-depth studies in the future.

The novel regulator HdrR controls the transcription of the heterodisulfide reductase operon hdrBCA in Methanosarcina barkeri

Methanogenic archaea play a key role in the global carbon cycle because these microorganisms remineralize organic compounds in various anaerobic environments. The microorganism Methanosarcina barkeri is a metabolically versatile methanogen, which can utilize acetate, methanol, and H2/CO2 to synthesize methane. However, the regulatory mechanisms underlying methanogenesis for different substrates remain unknown. In this study, RNA-seq analysis was used to investigate M. barkeri growth and gene transcription under different substrate regimes. According to the results, M. barkeri showed the best growth under methanol, followed by H2/CO2 and acetate, and these findings corresponded well with the observed variations in genes transcription abundance for different substrates. In addition, we identified a novel regulator, MSBRM_RS03855 (designated as HdrR), which specifically activates the transcription of the heterodisulfide reductase hdrBCA operon in M. barkeri. HdrR was able to bind to the hdrBCA operon promoter to regulate transcription. Furthermore, the structural model analyses revealed a helix-turn-helix domain, which is likely involved in DNA binding. Taken together, HdrR serves as a model to reveal how certain regulatory factors control the expression of key enzymes in the methanogenic pathway.IMPORTANCEThe microorganism Methanosarcina barkeri has a pivotal role in the global carbon cycle and contributes to global temperature homeostasis. The consequences of biological methanogenesis are far-reaching, including impacts on atmospheric methane and CO2 concentrations, agriculture, energy production, waste treatment, and human health. As such, reducing methane emissions is crucial to meeting set climate goals. The methanogenic activity of certain microorganisms can be drastically reduced by inhibiting the transcription of the hdrBCA operon, which encodes heterodisulfide reductases. Here, we provide novel insight into the mechanisms regulating hdrBCA operon transcription in the model methanogen M. barkeri. The results clarified that HdrR serves as a regulator of heterodisulfide reductase hdrBCA operon transcription during methanogenesis, which expands our understanding of the unique regulatory mechanisms that govern methanogenesis. The findings presented in this study can further our understanding of how genetic regulation can effectively reduce the methane emissions caused by methanogens.

Unveiling the unique role of iron in the metabolism of methanogens: A review

Methanogens are the main participants in the carbon cycle, catalyzing five methanogenic pathways. Methanogens utilize different iron-containing functional enzymes in different methanogenic processes. Iron is a vital element in methanogens, which can serve as a carrier or reactant in electron transfer. Therefore, iron plays an important role in the growth and metabolism of methanogens. In this paper, we cast light on the types and functions of iron-containing functional enzymes involved in different methanogenic pathways, and the roles iron play in energy/substance metabolism of methanogenesis. Furthermore, this review provides certain guiding significance for lowering CH4 emissions, boosting the carbon sink capacity of ecosystems and promoting green and low-carbon development in the future.

Microbiome assembly and stability during start-up of a full-scale, two-phase anaerobic digester fed cow manure and mixed organic feedstocks

Carbon transformations during anaerobic digestion are mediated by complex microbiomes, but their assembly is poorly understood, especially in full-scale digesters. Gene-centric metagenomics combining functional and taxonomic classification was performed for an on-farm digester during start-up. Cow manure and organic waste pre-treated in a hydrolysis tank were fed to the methane-producing digester and the volatile solids loading rate was slowly increased from 0 to 3.5 kg volatile solids m-3 d-1 over one year. The microbial community in the anaerobic digester exhibited a high ratio of archaea, which were dominated by hydrogenotrophic methanogens. Bacteria in the anaerobic digester had a high abundance of genes for ferredoxin cycling, H2 generation, and more metabolically complex fermentations than in the hydrolysis tank. In total, the results show that a functionally stable microbiome was achieved quickly during start-up and that the microbiome created in the low-pH hydrolysis tank did not persist in the downstream anaerobic digester.

Thiol-Disulfide Exchange Kinetics and Redox Potential of the Coenzyme M and Coenzyme B Heterodisulfide, an Electron Acceptor Coupled to Energy Conservation in Methanogenic Archaea

Methanogenic and methanotrophic archaea play important roles in the global carbon cycle by interconverting CO2 and methane. To conserve energy from these metabolic pathways that happen close to the thermodynamic equilibrium, specific electron carriers have evolved to balance the redox potentials between key steps. Reduced ferredoxins required to activate CO2 are provided by energetical coupling to the reduction of the high-potential heterodisulfide (HDS) of coenzyme M (2-mercaptoethanesulfonate) and coenzyme B (7-mercaptoheptanoylthreonine phosphate). While the standard redox potential of this important HDS has been determined previously to be -143 mV (Tietze et al. 2003 DOI: 10.1002/cbic.200390053), we have measured thiol disulfide exchange kinetics and reassessed this value by equilibrating thiol-disulfide mixtures of coenzyme M, coenzyme B, and mercaptoethanol. We determined the redox potential of the HDS of coenzyme M and coenzyme B to be -16.4±1.7 mV relative to the reference thiol mercaptoethanol (E0 '=-264 mV). The resulting E0 ' values are -281 mV for the HDS, -271 mV for the homodisulfide of coenzyme M, and -270 mV for the homodisulfide of coenzyme B. We discuss the importance of these updated values for the physiology of methanogenic and methanotrophic archaea and their implications in terms of energy conservation.

Evolving understanding of rumen methanogen ecophysiology

Production of methane by methanogenic archaea, or methanogens, in the rumen of ruminants is a thermodynamic necessity for microbial conversion of feed to volatile fatty acids, which are essential nutrients for the animals. On the other hand, methane is a greenhouse gas and its production causes energy loss for the animal. Accordingly, there are ongoing efforts toward developing effective strategies for mitigating methane emissions from ruminant livestock that require a detailed understanding of the diversity and ecophysiology of rumen methanogens. Rumen methanogens evolved from free-living autotrophic ancestors through genome streamlining involving gene loss and acquisition. The process yielded an oligotrophic lifestyle, and metabolically efficient and ecologically adapted descendants. This specialization poses serious challenges to the efforts of obtaining axenic cultures of rumen methanogens, and consequently, the information on their physiological properties remains in most part inferred from those of their non-rumen representatives. This review presents the current knowledge of rumen methanogens and their metabolic contributions to enteric methane production. It also identifies the respective critical gaps that need to be filled for aiding the efforts to mitigate methane emission from livestock operations and at the same time increasing the productivity in this critical agriculture sector.

Model Organisms To Study Methanogenesis, a Uniquely Archaeal Metabolism

Methanogenic archaea are the only organisms that produce CH4 as part of their energy-generating metabolism. They are ubiquitous in oxidant-depleted, anoxic environments such as aquatic sediments, anaerobic digesters, inundated agricultural fields, the rumen of cattle, and the hindgut of termites, where they catalyze the terminal reactions in the degradation of organic matter. Methanogenesis is the only metabolism that is restricted to members of the domain Archaea. Here, we discuss the importance of model organisms in the history of methanogen research, including their role in the discovery of the archaea and in the biochemical and genetic characterization of methanogenesis. We also discuss outstanding questions in the field and newly emerging model systems that will expand our understanding of this uniquely archaeal metabolism.

The dual role of a multi-heme cytochrome in methanogenesis: MmcA is important for energy conservation and carbon metabolism in Methanosarcina acetivorans

Methanogenic archaea belonging to the Order Methanosarcinales conserve energy using an electron transport chain (ETC). In the genetically tractable strain Methanosarcina acetivorans, ferredoxin donates electrons to the ETC via the Rnf (Rhodobacter nitrogen fixation) complex. The Rnf complex in M. acetivorans, unlike its counterpart in Bacteria, contains a multiheme c-type cytochrome (MHC) subunit called MmcA. Early studies hypothesized MmcA is a critical component of Rnf, however recent work posits that the primary role of MmcA is facilitating extracellular electron transport. To explore the physiological role of MmcA, we characterized M. acetivorans mutants lacking either the entire Rnf complex (∆mmcA-rnf) or just the MmcA subunit (∆mmcA). Our data show that MmcA is essential for growth during acetoclastic methanogenesis but neither Rnf nor MmcA is required for methanogenic growth on methylated compounds. On methylated compounds, the absence of MmcA alone leads to a more severe growth defect compared to a Rnf deletion likely due to different strategies for ferredoxin oxidation that arise in each strain. Transcriptomic data suggest that the ∆mmcA mutant might oxidize ferredoxin by upregulating the cytosolic Wood-Ljundahl pathway for acetyl-CoA synthesis, whereas the ∆mmcA-rnf mutant may repurpose the F₄₂₀ dehydrogenase complex (Fpo) to oxidize ferredoxin coupled to proton translocation. Beyond energy conservation, the deletion of rnf or mmcA leads to global transcriptional changes of genes involved in methanogenesis, carbon assimilation and regulation. Overall, our study provides systems-level insights into the non-overlapping roles of the Rnf bioenergetic complex and the associated MHC, MmcA.

Supplementation of Sulfide or Acetate and 2-Mercaptoethane Sulfonate Restores Growth of the Methanosarcina acetivorans ΔhdrABC Deletion Mutant during Methylotrophic Methanogenesis

Methanogenic archaea are important organisms in the global carbon cycle that grow by producing methane gas. Methanosarcina acetivorans is a methanogenic archaeum that can grow using methylated compounds, carbon monoxide, or acetate and produces renewable methane as a byproduct. However, there is limited knowledge of how combinations of substrates may affect metabolic fluxes in methanogens. Previous studies have shown that heterodisulfide reductase, the terminal oxidase in the electron transport system, is an essential enzyme in all methanogens. Deletion of genes encoding the nonessential methylotrophic heterodisulfide reductase enzyme (HdrABC) results in slower growth rate but increased metabolic efficiency. We hypothesized that increased sulfide, supplementation of mercaptoethanesulfonate (coenzyme M, CoM-SH), or acetate would metabolically alleviate the effect of the ΔhdrABC mutation. Increased sulfide improved growth of the mutant as expected; however, supplementation of both CoM-SH and acetate together were necessary to reduce the effect of the ΔhdrABC mutation. Supplementation of CoM-SH or acetate alone did not improve growth. These results support our model for the role of HdrABC in methanogenesis and suggest M.acetivorans is more efficient at conserving energy when supplemented with acetate. Our study suggests decreased Hdr enzyme activity can be overcome by nutritional supplementation with sulfide or coenzyme M and acetate, which are abundant in anaerobic environments.

Modular structure of complex II: An evolutionary perspective

Succinate dehydrogenases (SDHs) and fumarate reductases (FRDs) catalyse the interconversion of succinate and fumarate, a reaction highly conserved in all domains of life. The current classification of SDH/FRDs is based on the structure of the membrane anchor subunits and their cofactors. It is, however, unknown whether this classification would hold in the context of evolution. In this work, a large-scale comparative genomic analysis of complex II addresses the questions of its taxonomic distribution and phylogeny. Our findings report that for types C, D, and F, structural classification and phylogeny go hand in hand, while for types A, B and E the situation is more complex, highlighting the possibility for their classification into subgroups. Based on these findings, we proposed a revised version of the evolutionary scenario for these enzymes in which a primordial soluble module, corresponding to the cytoplasmatic subunits, would give rise to the current diversity via several independent membrane anchor attachment events.

Insights into the biotechnology potential of Methanosarcina

Methanogens are anaerobic archaea which conserve energy by producing methane. Found in nearly every anaerobic environment on earth, methanogens serve important roles in ecology as key organisms of the global carbon cycle, and in industry as a source of renewable biofuels. Environmentally, methanogenic archaea play an essential role in the reintroducing unavailable carbon to the carbon cycle by anaerobically converting low-energy, terminal metabolic degradation products such as one and two-carbon molecules into methane which then returns to the aerobic portion of the carbon cycle. In industry, methanogens are commonly used as an inexpensive source of renewable biofuels as well as serving as a vital component in the treatment of wastewater though this is only the tip of the iceberg with respect to their metabolic potential. In this review we will discuss how the efficient central metabolism of methanoarchaea could be harnessed for future biotechnology applications.

Evidence for Assimilatory Nitrate Reduction as a Previously Overlooked Pathway of Reactive Nitrogen Transformation in Estuarine Suspended Particulate Matter

Suspended particulate matter (SPM) contributes to the loss of reactive nitrogen (Nr) in estuarine ecosystems. Although denitrification and anaerobic ammonium oxidation in SPM compensate for the current imbalance of global nitrogen (N) inputs and sinks, it is largely unclear whether other pathways for Nr transformation exist in SPM. Here, we combined stable isotope measurements with metagenomics and metatranscriptomics to verify the occurrence of dissimilatory nitrate reduction to ammonium (DNRA) in the SPM of the Pearl River Estuary (PRE). Surprisingly, the conventional functional genes of DNRA (nirBD) were abundant and highly expressed in SPM, which was inconsistent with a low potential rate. Through taxonomic and comparative genomic analyses, we demonstrated that nitrite reductase (NirBD) in conjunction with assimilatory nitrate reductase (NasA) performed assimilatory nitrate reduction (ANR) in SPM, and diverse alpha- and gamma-proteobacterial lineages were identified as key active heterotrophic ANR bacteria. Moreover, ANR was predicted to have a relative higher occurrence than denitrification and DNRA in a survey of Nr transformation pathways in SPM across the PRE spanning 65 km. Collectively, this study characterizes a previously overlooked pathway of Nr transformation mediated by heterotrophic ANR bacteria in SPM and has important implications for our understanding of N cycling in estuaries.

Multivariate insights into enhanced biogas production in thermophilic dry anaerobic co-digestion of food waste with kitchen waste or garden waste: Process properties, microbial communities and metagenomic analyses

Multisubstrate synergetic anaerobic co-digestion can effectively overcome low efficiency of food waste (FW) mono-digestion. This study investigated the effect of supplementing FW with kitchen waste (KW) or garden waste (GW) on thermophilic dry anaerobic co-digestion. FW-KW and FW-GW co-digestion enhanced biogas production by 24.69 % and 44.96 % at organic loading rate (OLR) of 3 g VS L^-1 d^-1, and increased OLR tolerance from 3 to 4 g VS L^-1 d^-1 through mitigating ammonia nitrogen inhibition and volatile fatty acids accumulation. Co-digestion enriched the dominant hydrolytic bacteria Defluviitoga, resulting in an acceleration of substrate hydrolysis. FW-KW co-digestion improved biogas production by increasing gene abundance related to key enzymes in methanogenesis pathways and promoting the conversion of intermediate products into methane. FW-GW co-digestion enhanced biogas production by enriching ABC transporters-associated genes, leading to efficient substrate utilization. This study provides a promising approach for FW treatment with multivariate insights into thermophilic dry anaerobic co-digestion.

Deconstructing Methanosarcina acetivorans into an acetogenic archaeon

The reductive acetyl-coenzyme A (acetyl-CoA) pathway is the only carbon fixation pathway that can also be used for energy conservation like it is known for acetogenic bacteria. In methanogenic archaea, this pathway is extended with one route toward acetyl-CoA formation for anabolism and another route toward methane formation for catabolism. Which of these traits is ancestral in evolution has not been resolved. By diverging virtually all substrate carbon from methanogenesis to flow through acetyl-CoA, Methanosarcina acetivorans can be converted to an acetogenic organism. Being able to deconstruct methanogenic into the seemingly simpler acetogenic energy metabolism provides compelling evidence that methanogens are not nearly as metabolically limited as previously thought and suggests that methanogenesis might have evolved from the acetyl-CoA pathway.

The reductive acetyl-coenzyme A (acetyl-CoA) pathway, whereby carbon dioxide is sequentially reduced to acetyl-CoA via coenzyme-bound C1 intermediates, is the only autotrophic pathway that can at the same time be the means for energy conservation. A conceptually similar metabolism and a key process in the global carbon cycle is methanogenesis, the biogenic formation of methane. All known methanogenic archaea depend on methanogenesis to sustain growth and use the reductive acetyl-CoA pathway for autotrophic carbon fixation. Here, we converted a methanogen into an acetogen and show that Methanosarcina acetivorans can dispense with methanogenesis for energy conservation completely. By targeted disruption of the methanogenic pathway, followed by adaptive evolution, a strain was created that sustained growth via carbon monoxide–dependent acetogenesis. A minute flux (less than 0.2% of the carbon monoxide consumed) through the methane-liberating reaction remained essential, indicating that currently living methanogens utilize metabolites of this reaction also for anabolic purposes. These results suggest that the metabolic flexibility of methanogenic archaea might be much greater than currently known. Also, our ability to deconstruct a methanogen into an acetogen by merely removing cellular functions provides experimental support for the notion that methanogenesis could have evolved from the reductive acetyl-coenzyme A pathway.

Comparative genomics reveals electron transfer and syntrophic mechanisms differentiating methanotrophic and methanogenic archaea

The anaerobic oxidation of methane coupled to sulfate reduction is a microbially mediated process requiring a syntrophic partnership between anaerobic methanotrophic (ANME) archaea and sulfate-reducing bacteria (SRB). Based on genome taxonomy, ANME lineages are polyphyletic within the phylum Halobacterota, none of which have been isolated in pure culture. Here, we reconstruct 28 ANME genomes from environmental metagenomes and flow sorted syntrophic consortia. Together with a reanalysis of previously published datasets, these genomes enable a comparative analysis of all marine ANME clades. We review the genomic features that separate ANME from their methanogenic relatives and identify what differentiates ANME clades. Large multiheme cytochromes and bioenergetic complexes predicted to be involved in novel electron bifurcation reactions are well distributed and conserved in the ANME archaea, while significant variations in the anabolic C1 pathways exists between clades. Our analysis raises the possibility that methylotrophic methanogenesis may have evolved from a methanotrophic ancestor.

A comparative genomics study of anaerobic methanotrophic (ANME) archaea reveals the genetic "parts list" associated with the repeated evolutionary transition between methanogenic and methanotrophic metabolism in the archaeal domain of life.

Anaerobic Production of Isoprene by Engineered Methanosarcina Species Archaea

Isoprene is a valuable petrochemical used for a wide variety of consumer goods, such as adhesives and synthetic rubber. We were able to achieve a high yield of renewable isoprene by taking advantage of the naturally high-flux mevalonate lipid synthesis pathway in anaerobic methane-producing archaea (methanogens). Our study illustrates that by genetically manipulating Methanosarcina species methanogens, it is possible to create organisms that grow by producing the hemiterpene isoprene. Mass balance measurements show that engineered methanogens direct up to 4% of total carbon flux to isoprene, demonstrating that methanogens produce higher isoprene yields than engineered yeast, bacteria, or cyanobacteria, and from inexpensive feedstocks. Expression of isoprene synthase resulted in increased biomass and changes in gene expression that indicate that isoprene synthesis depletes membrane precursors and redirects electron flux, enabling isoprene to be a major metabolic product. Our results demonstrate that methanogens are a promising engineering chassis for renewable isoprene synthesis.IMPORTANCE A significant barrier to implementing renewable chemical technologies is high production costs relative to those for petroleum-derived products. Existing technologies using engineered organisms have difficulty competing with petroleum-derived chemicals due to the cost of feedstocks (such as glucose), product extraction, and purification. The hemiterpene monomer isoprene is one such chemical that cannot currently be produced using cost-competitive renewable biotechnologies. To reduce the cost of renewable isoprene, we have engineered methanogens to synthesize it from inexpensive feedstocks such as methane, methanol, acetate, and carbon dioxide. The "isoprenogen" strains we developed have potential to be used for industrial production of inexpensive renewable isoprene.

Functional diversity of prokaryotic HdrA(BC) modules: Role in flavin-based electron bifurcation processes and beyond

In methanogenic archaea, the archetypical complex of heterodisulfide reductase (HdrABC) and hydrogenase (MvhAGD) couples the endergonic reduction of CO₂ by H₂ to the exergonic reduction of the CoB-S-S-CoM heterodisulfide by H₂ via flavin-based electron bifurcation. Presently known enzymes containing HdrA(BC)-like components play key roles in methanogenesis, acetogenesis, respiratory sulfate reduction, lithotrophic reduced sulfur compound oxidation, aromatic compound degradation, fermentations, and probably many further processes. This functional diversity is achieved by a modular architecture of HdrA(BC) enzymes, where a big variety of electron input/output modules may be connected either directly or via adaptor modules to the HdrA(BC) components. Many, but not all HdrA(BC) complexes are proposed to catalyse a flavin-based electron bifurcation/confurcation. Despite the availability of HdrA(BC) crystal structures, fundamental questions of electron transfer and energy coupling processes remain. Here, we address the common properties and functional diversity of HdrA(BC) core modules integrated into electron-transfer machineries of outstanding complexity.

Newly discovered Asgard archaea Hermodarchaeota potentially degrade alkanes and aromatics via alkyl/benzyl-succinate synthase and benzoyl-CoA pathway

Asgard archaea are widely distributed in anaerobic environments. Previous studies revealed the potential capability of Asgard archaea to utilize various organic substrates including proteins, carbohydrates, fatty acids, amino acids and hydrocarbons, suggesting that Asgard archaea play an important role in sediment carbon cycling. Here, we describe a previously unrecognized archaeal phylum, Hermodarchaeota, affiliated with the Asgard superphylum. The genomes of these archaea were recovered from metagenomes generated from mangrove sediments, and were found to encode alkyl/benzyl-succinate synthases and their activating enzymes that are similar to those identified in alkane-degrading sulfate-reducing bacteria. Hermodarchaeota also encode enzymes potentially involved in alkyl-coenzyme A and benzoyl-coenzyme A oxidation, the Wood–Ljungdahl pathway and nitrate reduction. These results indicate that members of this phylum have the potential to strictly anaerobically degrade alkanes and aromatic compounds, coupling the reduction of nitrate. By screening Sequence Read Archive, additional genes encoding 16S rRNA and alkyl/benzyl-succinate synthases analogous to those in Hermodarchaeota were identified in metagenomic datasets from a wide range of marine and freshwater sediments. These findings suggest that Asgard archaea capable of degrading alkanes and aromatics via formation of alkyl/benzyl-substituted succinates are ubiquitous in sediments.

Evidence for non-methanogenic metabolisms in globally distributed archaeal clades basal to the Methanomassiliicoccales

Recent discoveries of mcr and mcr-like genes in genomes from diverse archaeal lineages suggest that methane metabolism is an ancient pathway with a complicated evolutionary history. One conventional view is that methanogenesis is an ancestral metabolism of the class Thermoplasmata. Through comparative genomic analysis of 12 Thermoplasmata metagenome-assembled genomes (MAGs) basal to the Methanomassiliicoccales, we show that these microorganisms do not encode the genes required for methanogenesis. Further analysis of 770 Ca. Thermoplasmatota genomes/MAGs found no evidence of mcrA homologues outside of the Methanomassiliicoccales. Together, these results suggest that methanogenesis was laterally acquired by an ancestor of the Methanomassiliicoccales. The 12 analysed MAGs include representatives from four orders basal to the Methanomassiliicoccales, including a high-quality MAG that likely represents a new order, Ca. Lunaplasma lacustris ord. nov. sp. nov. These MAGs are predicted to use diverse energy conservation pathways, including heterotrophy, sulfur and hydrogen metabolism, denitrification, and fermentation. Two lineages are widespread among anoxic, sedimentary environments, whereas Ca. Lunaplasma lacustris has thus far only been detected in alpine caves and subarctic lake sediments. These findings advance our understanding of the metabolic potential, ecology, and global distribution of the Thermoplasmata and provide insight into the evolutionary history of methanogenesis within the Ca. Thermoplasmatota.

Methanosarcina acetivorans: A Model for Mechanistic Understanding of Aceticlastic and Reverse Methanogenesis

Acetate-utilizing methanogens are responsible for approximately two-thirds of the one billion metric tons of methane produced annually in Earth's anaerobic environments. Methanosarcina acetivorans has emerged as a model organism for the mechanistic understanding of aceticlastic methanogenesis and reverse methanogenesis applicable to understanding the methane and carbon cycles in nature. It has the largest genome in the Archaea, supporting a metabolic complexity that enables a remarkable ability for adapting to environmental opportunities and challenges. Biochemical investigations have revealed an aceticlastic pathway capable of fermentative and respiratory energy conservation that explains how Ms. acetivorans is able to grow and compete in the environment. The mechanism of respiratory energy conservation also plays a role in overcoming endothermic reactions that are key to reversing methanogenesis.

Energy Conservation and Hydrogenase Function in Methanogenic Archaea, in Particular the Genus Methanosarcina

The biological production of methane is vital to the global carbon cycle and accounts for ca. 74% of total methane emissions. The organisms that facilitate this process, methanogenic archaea, belong to a large and phylogenetically diverse group that thrives in a wide range of anaerobic environments. Two main subgroups exist within methanogenic archaea: those with and those without cytochromes. Although a variety of metabolisms exist within this group, the reduction of growth substrates to methane using electrons from molecular hydrogen is, in a phylogenetic sense, the most widespread methanogenic pathway. Methanogens without cytochromes typically generate methane by the reduction of CO₂ with electrons derived from H₂, formate, or secondary alcohols, generating a transmembrane ion gradient for ATP production via an Na⁺-translocating methyltransferase (Mtr). These organisms also conserve energy with a novel flavin-based electron bifurcation mechanism, wherein the endergonic reduction of ferredoxin is facilitated by the exergonic reduction of a disulfide terminal electron acceptor coupled to either H₂ or formate oxidation. Methanogens that utilize cytochromes have a broader substrate range, and can convert acetate and methylated compounds to methane, in addition to the ability to reduce CO₂ Cytochrome-containing methanogens are able to supplement the ion motive force generated by Mtr with an H⁺-translocating electron transport system. In both groups, enzymes known as hydrogenases, which reversibly interconvert protons and electrons to molecular hydrogen, play a central role in the methanogenic process. This review discusses recent insight into methanogen metabolism and energy conservation mechanisms with a particular focus on the genus Methanosarcina.

Immediate Effects of Ammonia Shock on Transcription and Composition of a Biogas Reactor Microbiome

The biotechnological process of biogas production from organic material is carried out by a diverse microbial community under anaerobic conditions. However, the complex and sensitive microbial network present in anaerobic degradation of organic material can be disturbed by increased ammonia concentration introduced into the system by protein-rich substrates and imbalanced feeding. Here, we report on a simulated increase of ammonia concentration in a fed batch lab-scale biogas reactor experiment. Two treatment conditions were used simulating total ammonia nitrogen concentrations of 4.9 and 8.0 g/L with four replicate reactors. Each reactor was monitored concerning methane generation and microbial composition using 16S rRNA gene amplicon sequencing, while the transcriptional activity of the overall process was investigated by metatranscriptomic analysis. This allowed investigating the response of the microbial community in terms of species composition and transcriptional activity to a rapid upshift to high ammonia conditions. Clostridia and Methanomicrobiales dominated the microbial community throughout the entire experiment under both experimental conditions, while Methanosarcinales were only present in minor abundance. Transcription analysis demonstrated clostridial dominance with respect to genes encoding for enzymes of the hydrolysis step (cellulase, EC 3.2.1.4) as well as dominance of key genes for enzymes of the methanogenic pathway (methyl-CoM reductase, EC 2.8.4.1; heterodisulfide reductase, EC 1.8.98.1). Upon ammonia shock, the selected marker genes showed significant changes in transcriptional activity. Cellulose hydrolysis as well as methanogenesis were significantly reduced at high ammonia concentrations as indicated by reduced transcription levels of the corresponding genes. Based on these experiments we concluded that, apart from the methanogenic archaea, hydrolytic cellulose-degrading microorganisms are negatively affected by high ammonia concentrations. Further, Acholeplasma and Erysipelotrichia showed lower abundance under increased ammonia concentrations and thus might serve as indicator species for an earlier detection in order to counteract against ammonia crises.

A Membrane-Bound Cytochrome Enables Methanosarcina acetivorans To Conserve Energy from Extracellular Electron Transfer

The discovery of a methanogen that can conserve energy to support growth solely from the oxidation of organic carbon coupled to the reduction of an extracellular electron acceptor expands the possible environments in which methanogens might thrive. The potential importance of c-type cytochromes for extracellular electron transfer to syntrophic bacterial partners and/or Fe(III) minerals in some Archaea was previously proposed, but these studies with Methanosarcina acetivorans provide the first genetic evidence for cytochrome-based extracellular electron transfer in Archaea. The results suggest parallels with Gram-negative bacteria, such as Shewanella and Geobacter species, in which multiheme outer-surface c-type cytochromes are an essential component for electrical communication with the extracellular environment. M. acetivorans offers an unprecedented opportunity to study mechanisms for energy conservation from the anaerobic oxidation of one-carbon organic compounds coupled to extracellular electron transfer in Archaea with implications not only for methanogens but possibly also for Archaea that anaerobically oxidize methane.

Extracellular electron exchange in Methanosarcina species and closely related Archaea plays an important role in the global carbon cycle and enhances the speed and stability of anaerobic digestion by facilitating efficient syntrophic interactions. Here, we grew Methanosarcina acetivorans with methanol provided as the electron donor and the humic analogue, anthraquione-2,6-disulfonate (AQDS), provided as the electron acceptor when methane production was inhibited with bromoethanesulfonate. AQDS was reduced with simultaneous methane production in the absence of bromoethanesulfonate. Transcriptomics revealed that expression of the gene for the transmembrane, multiheme, c-type cytochrome MmcA was higher in AQDS-respiring cells than in cells performing methylotrophic methanogenesis. A strain in which the gene for MmcA was deleted failed to grow via AQDS reduction but grew with the conversion of methanol or acetate to methane, suggesting that MmcA has a specialized role as a conduit for extracellular electron transfer. Enhanced expression of genes for methanol conversion to methyl-coenzyme M and the Rnf complex suggested that methanol is oxidized to carbon dioxide in AQDS-respiring cells through a pathway that is similar to methyl-coenzyme M oxidation in methanogenic cells. However, during AQDS respiration the Rnf complex and reduced methanophenazine probably transfer electrons to MmcA, which functions as the terminal reductase for AQDS reduction. Extracellular electron transfer may enable the survival of methanogens in dynamic environments in which oxidized humic substances and Fe(III) oxides are intermittently available. The availability of tools for genetic manipulation of M. acetivorans makes it an excellent model microbe for evaluating c-type cytochrome-dependent extracellular electron transfer in Archaea.

Methanogens: pushing the boundaries of biology

Methanogens are anaerobic archaea that grow by producing methane gas. These microbes and their exotic metabolism have inspired decades of microbial physiology research that continues to push the boundary of what we know about how microbes conserve energy to grow. The study of methanogens has helped to elucidate the thermodynamic and bioenergetics basis of life, contributed our understanding of evolution and biodiversity, and has garnered an appreciation for the societal utility of studying trophic interactions between environmental microbes, as methanogens are important in microbial conversion of biogenic carbon into methane, a high-energy fuel. This review discusses the theoretical basis for energy conservation by methanogens and identifies gaps in methanogen biology that may be filled by undiscovered or yet-to-be engineered organisms.

Electron and Proton Flux for Carbon Dioxide Reduction in Methanosarcina barkeri During Direct Interspecies Electron Transfer

Direct interspecies electron transfer (DIET) is important in diverse methanogenic environments, but how methanogens participate in DIET is poorly understood. Therefore, the transcriptome of Methanosarcina barkeri grown via DIET in co-culture with Geobacter metallireducens was compared with its transcriptome when grown via H2 interspecies transfer (HIT) with Pelobacter carbinolicus. Notably, transcripts for the F420H2 dehydrogenase, Fpo, and the heterodisulfide reductase, HdrABC, were more abundant during growth on DIET. A model for CO2 reduction was developed from these results in which electrons delivered to methanophenazine in the cell membrane are transferred to Fpo. The external proton gradient necessary to drive the otherwise thermodynamically unfavorable reverse electron transport for Fpo-catalyzed F420 reduction is derived from protons released from G. metallireducens metabolism. Reduced F420 is a direct electron donor in the carbon dioxide reduction pathway and also serves as the electron donor for the proposed HdrABC-catalyzed electron bifurcation reaction in which reduced ferredoxin (also required for carbon dioxide reduction) is generated with simultaneous reduction of CoM-S-S-CoB. Expression of genes for putative redox-active proteins predicted to be localized on the outer cell surface was higher during growth on DIET, but further analysis will be required to identify the electron transfer route to methanophenazine. The results indicate that the pathways for electron and proton flux for CO2 reduction during DIET are substantially different than for HIT and suggest that gene expression patterns may also be useful for determining whether Methanosarcina are directly accepting electrons from other extracellular electron donors, such as corroding metals or electrodes.

Genetic, Biochemical, and Molecular Characterization of Methanosarcina barkeri Mutants Lacking Three Distinct Classes of Hydrogenase

The methanogenic archaeon Methanosarcina barkeri encodes three distinct types of hydrogenase, whose functions vary depending on the growth substrate. These include the F₄₂₀-dependent (Frh), methanophenazine-dependent (Vht), and ferredoxin-dependent (Ech) hydrogenases. To investigate their physiological roles, we characterized a series of mutants lacking each hydrogenase in various combinations. Mutants lacking Frh, Vht, or Ech in any combination failed to grow on H₂-CO₂, whereas only Vht and Ech were essential for growth on acetate. In contrast, a mutant lacking all three grew on methanol with a final growth yield similar to that of the wild type and produced methane and CO₂ in the expected 3:1 ratio but had a ca. 33% lower growth rate. Thus, hydrogenases play a significant, but nonessential, role during growth on this substrate. As previously observed, mutants lacking Ech failed to grow on methanol-H₂ unless they were supplemented with biosynthetic precursors. Interestingly, this phenotype was abolished in the Δech Δfrh and Δech Δfrh Δvht mutants, consistent with the idea that hydrogenases inhibit methanol oxidation in the presence of H₂, which prevents production of the reducing equivalents needed for biosynthesis. Quantification of the methane and CO₂ produced from methanol by resting cell suspensions of various mutants supported this conclusion. On the basis of the global transcriptional profiles, none of the hydrogenases were upregulated to compensate for the loss of the others. However, the transcript levels of the F₄₂₀ dehydrogenase operon were significantly higher in all strains lacking frh, suggesting a mechanism to sense the redox state of F₄₂₀ The roles of the hydrogenases in energy conservation during growth with each methanogenic pathway are discussed.IMPORTANCE Methanogenic archaea are key players in the global carbon cycle due to their ability to facilitate the remineralization of organic substrates in many anaerobic environments. The consequences of biological methanogenesis are far-reaching, with impacts on atmospheric methane and CO₂ concentrations, agriculture, energy production, waste treatment, and human health. The data presented here clarify the in vivo function of hydrogenases during methanogenesis, which in turn deepens our understanding of this unique form of metabolism. This knowledge is critical for a variety of important issues ranging from atmospheric composition to human health.

Origin and Evolution of Flavin-Based Electron Bifurcating Enzymes

Twelve evolutionarily unrelated oxidoreductases form enzyme complexes that catalyze the simultaneous coupling of exergonic and endergonic oxidation–reduction reactions to circumvent thermodynamic barriers and minimize free energy loss in a process known as flavin-based electron bifurcation. Common to these 12 bifurcating (Bf) enzymes are protein-bound flavin, the proposed site of bifurcation, and the electron carrier ferredoxin. Despite the documented role of Bf enzymes in balancing the redox state of intracellular electron carriers and in improving the efficiency of cellular metabolism, a comprehensive description of the diversity and evolutionary history of Bf enzymes is lacking. Here, we report the taxonomic distribution, functional diversity, and evolutionary history of Bf enzyme homologs in 4,588 archaeal, bacterial, and eukaryal genomes and 3,136 community metagenomes. Bf homologs were primarily detected in the genomes of anaerobes, including those of sulfate-reducers, acetogens, fermenters, and methanogens. Phylogenetic analyses of Bf enzyme catalytic subunits (oxidoreductases) suggest they were not a property of the Last Universal Common Ancestor of Archaea and Bacteria, which is consistent with the limited and unique taxonomic distributions of enzyme homologs among genomes. Further, phylogenetic analyses of oxidoreductase subunits reveal that non-Bf homologs predate Bf homologs. These observations indicate that multiple independent recruitments of flavoproteins to existing oxidoreductases enabled coupling of numerous new electron Bf reactions. Consistent with the role of these enzymes in the energy metabolism of anaerobes, homologs of Bf enzymes were enriched in metagenomes from subsurface environments relative to those from surface environments. Phylogenetic analyses of homologs from metagenomes reveal that the earliest evolving homologs of most Bf enzymes are from subsurface environments, including fluids from subsurface rock fractures and hydrothermal systems. Collectively, these data suggest strong selective pressures drove the emergence of Bf enzyme complexes via recruitment of flavoproteins that allowed for an increase in the efficiency of cellular metabolism and improvement in energy capture in anaerobes inhabiting a variety of subsurface anoxic habitats where the energy yield of oxidation-reduction reactions is generally low.

Electron Bifurcation and Confurcation in Methanogenesis and Reverse Methanogenesis

Reduction of the disulfide of coenzyme M and coenzyme B (CoMS–SCoB) by heterodisulfide reductases (HdrED and HdrABC) is the final step in all methanogenic pathways. Flavin-based electron bifurcation (FBEB) by soluble HdrABC homologs play additional roles in driving essential endergonic reactions at the expense of the exergonic reduction of CoMS–SCoM. In the first step of the CO₂ reduction pathway, HdrABC complexed with hydrogenase or formate dehydrogenase generates reduced ferredoxin (Fdx^2-) for the endergonic reduction of CO₂ coupled to the exergonic reduction of CoMS–SCoB dependent on FBEB of electrons from H₂ or formate. Roles for HdrABC:hydrogenase complexes are also proposed for pathways wherein the methyl group of methanol is reduced to methane with electrons from H₂. The HdrABC complexes catalyze FBEB-dependent oxidation of H₂ for the endergonic reduction of Fdx driven by the exergonic reduction of CoMS–SCoB. The Fdx^2- supplies electrons for reduction of the methyl group to methane. In H₂^- independent pathways, three-fourths of the methyl groups are oxidized producing Fdx^2- and reduced coenzyme F₄₂₀ (F₄₂₀H₂). The F₄₂₀H₂ donates electrons for reduction of the remaining methyl groups to methane requiring transfer of electrons from Fdx^2- to F₄₂₀. HdrA1B1C1 is proposed to catalyze FBEB-dependent oxidation of Fdx^2- for the endergonic reduction of F₄₂₀ driven by the exergonic reduction of CoMS–SCoB. In H₂^- independent acetotrophic pathways, the methyl group of acetate is reduced to methane with electrons derived from oxidation of the carbonyl group mediated by Fdx. Electron transport involves a membrane-bound complex (Rnf) that oxidizes Fdx^2- and generates a Na⁺ gradient driving ATP synthesis. It is postulated that F₄₂₀ is reduced by Rnf requiring HdrA2B2C2 catalyzing FBEB-dependent oxidation of F₄₂₀H₂ for the endergonic reduction of Fdx driven by the exergonic reduction of CoMS–SCoB. The Fdx^2- is recycled by Rnf and HdrA2B2C2 thereby conserving energy. The HdrA2B2C2 is also proposed to play a role in Fe(III)-dependent reverse methanogenesis. A flavin-based electron confurcating (FBEC) HdrABC complex is proposed for nitrate-dependent reverse methanogenesis in which the oxidation of CoM-SH/CoB-SH and Fdx^2- is coupled to reduction of F₄₂₀. The F₄₂₀H₂ donates electrons to a membrane complex that generates a proton gradient driving ATP synthesis.

A biochemical framework for anaerobic oxidation of methane driven by Fe(III)-dependent respiration

Consumption of methane by aerobic and anaerobic microbes governs the atmospheric level of this powerful greenhouse gas. Whereas a biochemical understanding of aerobic methanotrophy is well developed, a mechanistic understanding of anaerobic methanotrophy has been prevented by the unavailability of pure cultures. Here we report a biochemical investigation of Methanosarcina acetivorans, a methane-producing species capable of anaerobic methanotrophic growth dependent on reduction of Fe(III). Our findings support a pathway anchored by Fe(III)-dependent mechanisms for energy conservation driving endergonic reactions that are key to methanotrophic growth. The pathway is remarkably similar to pathways hypothesized for uncultured anaerobic methanotrophic archaea. The results contribute to an improved understanding of the methane cycle that is paramount to understanding human interventions influencing Earth’s climate. Finally, the pathway enables advanced development and optimization of biotechnologies converting methane to value-added products through metabolic engineering of M. acetivorans.

The unavailability of pure cultures has prevented a mechanistic understanding of anaerobic methanotrophy. Here the authors report a biochemical investigation of Methanosarcina acetivorans that supports a pathway anchored by Fe(III)-dependent mechanisms for energy conservation and driving endergonic reactions.

Identification, Comparison, and Validation of Robust Rumen Microbial Biomarkers for Methane Emissions Using Diverse Bos Taurus Breeds and Basal Diets

Previous shotgun metagenomic analyses of ruminal digesta identified some microbial information that might be useful as biomarkers to select cattle that emit less methane (CH4), which is a potent greenhouse gas. It is known that methane production (g/kgDMI) and to an extent the microbial community is heritable and therefore biomarkers can offer a method of selecting cattle for low methane emitting phenotypes. In this study a wider range of Bos Taurus cattle, varying in breed and diet, was investigated to determine microbial communities and genetic markers associated with high/low CH4 emissions. Digesta samples were taken from 50 beef cattle, comprising four cattle breeds, receiving two basal diets containing different proportions of concentrate and also including feed additives (nitrate or lipid), that may influence methane emissions. A combination of partial least square analysis and network analysis enabled the identification of the most significant and robust biomarkers of CH4 emissions (VIP > 0.8) across diets and breeds when comparing all potential biomarkers together. Genes associated with the hydrogenotrophic methanogenesis pathway converting carbon dioxide to methane, provided the dominant biomarkers of CH4 emissions and methanogens were the microbial populations most closely correlated with CH4 emissions and identified by metagenomics. Moreover, these genes grouped together as confirmed by network analysis for each independent experiment and when combined. Finally, the genes involved in the methane synthesis pathway explained a higher proportion of variation in CH4 emissions by PLS analysis compared to phylogenetic parameters or functional genes. These results confirmed the reproducibility of the analysis and the advantage to use these genes as robust biomarkers of CH4 emissions. Volatile fatty acid concentrations and ratios were significantly correlated with CH4, but these factors were not identified as robust enough for predictive purposes. Moreover, the methanotrophic Methylomonas genus was found to be negatively correlated with CH4. Finally, this study confirmed the importance of using robust and applicable biomarkers from the microbiome as a proxy of CH4 emissions across diverse production systems and environments.

Physiological Evidence for Isopotential Tunneling in the Electron Transport Chain of Methane-Producing Archaea

Many, but not all, organisms use quinones to conserve energy in their electron transport chains. Fermentative bacteria and methane-producing archaea (methanogens) do not produce quinones but have devised other ways to generate ATP. Methanophenazine (MPh) is a unique membrane electron carrier found in Methanosarcina species that plays the same role as quinones in the electron transport chain. To extend the analogy between quinones and MPh, we compared the MPh pool sizes between two well-studied Methanosarcina species, Methanosarcina acetivorans C2A and Methanosarcina barkeri Fusaro, to the quinone pool size in the bacterium Escherichia coli We found the quantity of MPh per cell increases as cultures transition from exponential growth to stationary phase, and absolute quantities of MPh were 3-fold higher in M. acetivorans than in M. barkeri The concentration of MPh suggests the cell membrane of M. acetivorans, but not of M. barkeri, is electrically quantized as if it were a single conductive metal sheet and near optimal for rate of electron transport. Similarly, stationary (but not exponentially growing) E. coli cells also have electrically quantized membranes on the basis of quinone content. Consistent with our hypothesis, we demonstrated that the exogenous addition of phenazine increases the growth rate of M. barkeri three times that of M. acetivorans Our work suggests electron flux through MPh is naturally higher in M. acetivorans than in M. barkeri and that hydrogen cycling is less efficient at conserving energy than scalar proton translocation using MPh.IMPORTANCE Can we grow more from less? The ability to optimize and manipulate metabolic efficiency in cells is the difference between commercially viable and nonviable renewable technologies. Much can be learned from methane-producing archaea (methanogens) which evolved a successful metabolic lifestyle under extreme thermodynamic constraints. Methanogens use highly efficient electron transport systems and supramolecular complexes to optimize electron and carbon flow to control biomass synthesis and the production of methane. Worldwide, methanogens are used to generate renewable methane for heat, electricity, and transportation. Our observations suggest Methanosarcina acetivorans, but not Methanosarcina barkeri, has electrically quantized membranes. Escherichia coli, a model facultative anaerobe, has optimal electron transport at the stationary phase but not during exponential growth. This study also suggests the metabolic efficiency of bacteria and archaea can be improved using exogenously supplied lipophilic electron carriers. The enhancement of methanogen electron transport through methanophenazine has the potential to increase renewable methane production at an industrial scale.

Cas9-mediated genome editing in the methanogenic archaeon Methanosarcina acetivorans

Although Cas9-mediated genome editing has proven to be a powerful genetic tool in eukaryotes, its application in Bacteria has been limited because of inefficient targeting or repair; and its application to Archaea has yet to be reported. Here we describe the development of a Cas9-mediated genome-editing tool that allows facile genetic manipulation of the slow-growing methanogenic archaeon Methanosarcina acetivorans Introduction of both insertions and deletions by homology-directed repair was remarkably efficient and precise, occurring at a frequency of approximately 20% relative to the transformation efficiency, with the desired mutation being found in essentially all transformants examined. Off-target activity was not observed. We also observed that multiple single-guide RNAs could be expressed in the same transcript, reducing the size of mutagenic plasmids and simultaneously simplifying their design. Cas9-mediated genome editing reduces the time needed to construct mutants by more than half (3 vs. 8 wk) and allows simultaneous construction of double mutants with high efficiency, exponentially decreasing the time needed for complex strain constructions. Furthermore, coexpression the nonhomologous end-joining (NHEJ) machinery from the closely related archaeon, Methanocella paludicola, allowed efficient Cas9-mediated genome editing without the need for a repair template. The NHEJ-dependent mutations included deletions ranging from 75 to 2.7 kb in length, most of which appear to have occurred at regions of naturally occurring microhomology. The combination of homology-directed repair-dependent and NHEJ-dependent genome-editing tools comprises a powerful genetic system that enables facile insertion and deletion of genes, rational modification of gene expression, and testing of gene essentiality.

A Ferredoxin- and F420H2-Dependent, Electron-Bifurcating, Heterodisulfide Reductase with Homologs in the Domains Bacteria and Archaea

Heterodisulfide reductases (Hdr) of the HdrABC class are ancient enzymes and a component of the anaerobic core belonging to the prokaryotic common ancestor. The ancient origin is consistent with the widespread occurrence of genes encoding putative HdrABC homologs in metabolically diverse prokaryotes predicting diverse physiological functions; however, only one HdrABC has been characterized and that was from a narrow metabolic group of obligate CO₂-reducing methanogenic anaerobes (methanogens) from the domain Archaea. Here we report the biochemical characterization of an HdrABC homolog (HdrA2B2C2) from the acetate-utilizing methanogen Methanosarcina acetivorans with unusual properties structurally and functionally distinct from the only other HdrABC characterized. Homologs of the HdrA2B2C2 archetype are present in phylogenetically and metabolically diverse species from the domains Bacteria and Archaea. The expression of the individual HdrA2, HdrB2, and HdrB2C2 enzymes in Escherichia coli, and reconstitution of an active HdrA2B2C2 complex, revealed an intersubunit electron transport pathway dependent on ferredoxin or coenzyme F₄₂₀ (F₄₂₀H₂) as an electron donor. Remarkably, HdrA2B2C2 couples the previously unknown endergonic oxidation of F₄₂₀H₂ and reduction of ferredoxin with the exergonic oxidation of F₄₂₀H₂ and reduction of the heterodisulfide of coenzyme M and coenzyme B (CoMS-SCoB). The unique electron bifurcation predicts a role for HdrA2B2C2 in Fe(III)-dependent anaerobic methane oxidation (ANME) by M. acetivorans and uncultured species from ANME environments. HdrA2B2C2, ubiquitous in acetotrophic methanogens, was shown to participate in electron transfer during acetotrophic growth of M. acetivorans and proposed to be essential for growth in the environment when acetate is limiting.

Discovery of the archetype HdrA2B2C2 heterodisulfide reductase with categorically unique properties extends the understanding of this ancient family beyond CO₂-reducing methanogens to include diverse prokaryotes from the domains Bacteria and Archaea. The unprecedented coenzyme F₄₂₀-dependent electron bifurcation, an emerging fundamental principle of energy conservation, predicts a role for HdrA2B2C2 in diverse metabolisms, including anaerobic CH₄-oxidizing pathways. The results document an electron transport role for HdrA2B2C2 in acetate-utilizing methanogens responsible for at least two-thirds of the methane produced in Earth’s biosphere. The previously unavailable heterologous production of individual subunits and the reconstitution of HdrA2B2C2 with activity have provided an understanding of intersubunit electron transfer in the HdrABC class and a platform for investigating the principles of electron bifurcation.

Genome-Scale Metabolic Modeling of Archaea Lends Insight into Diversity of Metabolic Function

Decades of biochemical, bioinformatic, and sequencing data are currently being systematically compiled into genome-scale metabolic reconstructions (GEMs). Such reconstructions are knowledge-bases useful for engineering, modeling, and comparative analysis. Here we review the fifteen GEMs of archaeal species that have been constructed to date. They represent primarily members of the Euryarchaeota with three-quarters comprising representative of methanogens. Unlike other reviews on GEMs, we specially focus on archaea. We briefly review the GEM construction process and the genealogy of the archaeal models. The major insights gained during the construction of these models are then reviewed with specific focus on novel metabolic pathway predictions and growth characteristics. Metabolic pathway usage is discussed in the context of the composition of each organism's biomass and their specific energy and growth requirements. We show how the metabolic models can be used to study the evolution of metabolism in archaea. Conservation of particular metabolic pathways can be studied by comparing reactions using the genes associated with their enzymes. This demonstrates the utility of GEMs to evolutionary studies, far beyond their original purpose of metabolic modeling; however, much needs to be done before archaeal models are as extensively complete as those for bacteria.

Genome-wide gene expression and RNA half-life measurements allow predictions of regulation and metabolic behavior in Methanosarcina acetivorans

BACKGROUND

While a few studies on the variations in mRNA expression and half-lives measured under different growth conditions have been used to predict patterns of regulation in bacterial organisms, the extent to which this information can also play a role in defining metabolic phenotypes has yet to be examined systematically. Here we present the first comprehensive study for a model methanogen.

RESULTS

We use expression and half-life data for the methanogen Methanosarcina acetivorans growing on fast- and slow-growth substrates to examine the regulation of its genes. Unlike Escherichia coli where only small shifts in half-lives were observed, we found that most mRNA have significantly longer half-lives for slow growth on acetate compared to fast growth on methanol or trimethylamine. Interestingly, half-life shifts are not uniform across functional classes of enzymes, suggesting the existence of a selective stabilization mechanism for mRNAs. Using the transcriptomics data we determined whether transcription or degradation rate controls the change in transcript abundance. Degradation was found to control abundance for about half of the metabolic genes underscoring its role in regulating metabolism. Genes involved in half of the metabolic reactions were found to be differentially expressed among the substrates suggesting the existence of drastically different metabolic phenotypes that extend beyond just the methanogenesis pathways. By integrating expression data with an updated metabolic model of the organism (iST807) significant differences in pathway flux and production of metabolites were predicted for the three growth substrates.

CONCLUSIONS

This study provides the first global picture of differential expression and half-lives for a class II methanogen, as well as provides the first evidence in a single organism that drastic genome-wide shifts in RNA half-lives can be modulated by growth substrate. We determined which genes in each metabolic pathway control the flux and classified them as regulated by transcription (e.g. transcription factor) or degradation (e.g. post-transcriptional modification). We found that more than half of genes in metabolism were controlled by degradation. Our results suggest that M. acetivorans employs extensive post-transcriptional regulation to optimize key metabolic steps, and more generally that degradation could play a much greater role in optimizing an organism's metabolism than previously thought.

Metatranscriptome analysis of active microbial communities in produced water samples from the Marcellus Shale

Controlling microbial activity is a primary concern during the management of the large volumes of wastewater (produced water) generated during high-volume hydraulic fracturing. In this study we analyzed the transcriptional activity (metatranscriptomes) of three produced water samples from the Marcellus Shale. The goal of this study was to describe active metabolic pathways of industrial concern for produced water management and reuse, and to improve understanding of produced water microbial activity. Metatranscriptome analysis revealed active biofilm formation, sulfide production, and stress management mechanisms of the produced water microbial communities. Biofilm-formation and sulfate-reduction pathways were identified in all samples. Genes related to a diverse array of stress response mechanisms were also identified with implications for biocide efficacy. Additionally, active expression of a methanogenesis pathway was identified in a sample of produced water collected prior to holding pond storage. The active microbial community identified by metatranscriptome analysis was markedly different than the community composition as identified by 16S rRNA sequencing, highlighting the value of evaluating the active microbial fraction during assessments of produced water biofouling potential and evaluation of biocide application strategies. These results indicate biofouling and corrosive microbial processes are active in produced water and should be taken into consideration while designing produced water reuse strategies.

Genomic composition and dynamics among Methanomicrobiales predict adaptation to contrasting environments

Members of the order Methanomicrobiales are abundant, and sometimes dominant, hydrogenotrophic (H₂-CO₂ utilizing) methanoarchaea in a broad range of anoxic habitats. Despite their key roles in greenhouse gas emissions and waste conversion to methane, little is known about the physiological and genomic bases for their widespread distribution and abundance. In this study, we compared the genomes of nine diverse Methanomicrobiales strains, examined their pangenomes, reconstructed gene flow and identified genes putatively mediating their success across different habitats. Most strains slowly increased gene content whereas one, Methanocorpusculum labreanum, evidenced genome downsizing. Peat-dwelling Methanomicrobiales showed adaptations centered on improved transport of scarce inorganic nutrients and likely use H⁺ rather than Na⁺ transmembrane chemiosmotic gradients during energy conservation. In contrast, other Methanomicrobiales show the potential to concurrently use Na⁺ and H⁺ chemiosmotic gradients. Analyses also revealed that the Methanomicrobiales lack a canonical electron bifurcation system (MvhABGD) known to produce low potential electrons in other orders of hydrogenotrophic methanogens. Additional putative differences in anabolic metabolism suggest that the dynamics of interspecies electron transfer from Methanomicrobiales syntrophic partners can also differ considerably. Altogether, these findings suggest profound differences in electron trafficking in the Methanomicrobiales compared with other hydrogenotrophs, and warrant further functional evaluations.

Assessing methanotrophy and carbon fixation for biofuel production by Methanosarcina acetivorans

Background

Methanosarcina acetivorans is a model archaeon with renewed interest due to its unique reversible methane production pathways. However, the mechanism and relevant pathways implicated in (co)utilizing novel carbon substrates in this organism are still not fully understood. This paper provides a comprehensive inventory of thermodynamically feasible routes for anaerobic methane oxidation, co-reactant utilization, and maximum carbon yields of major biofuel candidates by M. acetivorans.

Results

Here, an updated genome-scale metabolic model of M. acetivorans is introduced (iMAC868 containing 868 genes, 845 reactions, and 718 metabolites) by integrating information from two previously reconstructed metabolic models (i.e., iVS941 and iMB745), modifying 17 reactions, adding 24 new reactions, and revising 64 gene-protein-reaction associations based on newly available information. The new model establishes improved predictions of growth yields on native substrates and is capable of correctly predicting the knockout outcomes for 27 out of 28 gene deletion mutants. By tracing a bifurcated electron flow mechanism, the iMAC868 model predicts thermodynamically feasible (co)utilization pathway of methane and bicarbonate using various terminal electron acceptors through the reversal of the aceticlastic pathway.

Conclusions

This effort paves the way in informing the search for thermodynamically feasible ways of (co)utilizing novel carbon substrates in the domain Archaea.

Electronic supplementary material

The online version of this article (doi:10.1186/s12934-015-0404-4) contains supplementary material, which is available to authorized users.

Rerouting Cellular Electron Flux To Increase the Rate of Biological Methane Production

Methanogens are anaerobic archaea that grow by producing methane, a gas that is both an efficient renewable fuel and a potent greenhouse gas. We observed that overexpression of the cytoplasmic heterodisulfide reductase enzyme HdrABC increased the rate of methane production from methanol by 30% without affecting the growth rate relative to the parent strain. Hdr enzymes are essential in all known methane-producing archaea. They function as the terminal oxidases in the methanogen electron transport system by reducing the coenzyme M (2-mercaptoethane sulfonate) and coenzyme B (7-mercaptoheptanoylthreonine sulfonate) heterodisulfide, CoM-S-S-CoB, to regenerate the thiol-coenzymes for reuse. In Methanosarcina acetivorans, HdrABC expression caused an increased rate of methanogenesis and a decrease in metabolic efficiency on methylotrophic substrates. When acetate was the sole carbon and energy source, neither deletion nor overexpression of HdrABC had an effect on growth or methane production rates. These results suggest that in cells grown on methylated substrates, the cell compensates for energy losses due to expression of HdrABC with an increased rate of substrate turnover and that HdrABC lacks the appropriate electron donor in acetate-grown cells.

Genetic basis for metabolism of methylated sulfur compounds in Methanosarcina species

Methanosarcina acetivorans uses a variety of methylated sulfur compounds as carbon and energy sources. Previous studies implicated the mtsD, mtsF, and mtsH genes in catabolism of dimethylsulfide, but the genes required for use of other methylsulfides have yet to be established. Here, we show that a four-gene locus, designated mtpCAP-msrH, is specifically required for growth on methylmercaptopropionate (MMPA). The mtpC, mtpA, and mtpP genes encode a putative corrinoid protein, a coenzyme M (CoM) methyltransferase, and a major facilitator superfamily (MFS) transporter, respectively, while msrH encodes a putative transcriptional regulator. Mutants lacking mtpC or mtpA display a severe growth defect in MMPA medium but are unimpaired during growth on other substrates. The mtpCAP genes comprise a transcriptional unit that is highly and specifically upregulated during growth on MMPA, whereas msrH is monocistronic and constitutively expressed. Mutants lacking msrH fail to transcribe mtpCAP and grow poorly in MMPA medium, consistent with the assignment of its product as a transcriptional activator. The mtpCAP-msrH locus is conserved in numerous marine methanogens, including eight Methanosarcina species that we showed are capable of growth on MMPA. Mutants lacking the mtsD, mtsF, and mtsH genes display a 30% reduction in growth yield when grown on MMPA, suggesting that these genes play an auxiliary role in MMPA catabolism. A quadruple ΔmtpCAP ΔmtsD ΔmtsF ΔmtsH mutant strain was incapable of growth on MMPA. Reanalysis of mtsD, mtsF, and mtsH mutants suggests that the preferred substrate for MtsD is dimethylsulfide, while the preferred substrate for MtsF is methanethiol.

IMPORTANCE

Methylated sulfur compounds play pivotal roles in the global sulfur and carbon cycles and contribute to global temperature homeostasis. Although the degradation of these molecules by aerobic bacteria has been well studied, relatively little is known regarding their fate in anaerobic ecosystems. In this study, we identify the genetic basis for metabolism of methylmercaptopropionate, dimethylsulfide, and methanethiol by strictly anaerobic methanogens of the genus Methanosarcina. These data will aid the development of predictive sulfur cycle models and enable molecular ecological approaches for the study of methylated sulfur metabolism in anaerobic ecosystems.

New mode of energy metabolism in the seventh order of methanogens as revealed by comparative genome analysis of “Candidatus methanoplasma termitum”

The recently discovered seventh order of methanogens, the Methanomassiliicoccales (previously referred to as “Methanoplasmatales”), so far consists exclusively of obligately hydrogen-dependent methylotrophs. We sequenced the complete genome of “Candidatus Methanoplasma termitum” from a highly enriched culture obtained from the intestinal tract of termites and compared it with the previously published genomes of three other strains from the human gut, including the first isolate of the order. Like all other strains, “Ca. Methanoplasma termitum” lacks the entire pathway for CO2 reduction to methyl coenzyme Mand produces methane by hydrogen-dependent reduction of methanol or methylamines, which is consistent with additional physiological data. However, the shared absence of cytochromes and an energy-converting hydrogenase for the reoxidation of the ferredoxin produced by the soluble heterodisulfide reductase indicates that Methanomassiliicoccales employ a new mode of energy metabolism, which differs from that proposed for the obligately methylotrophic Methanosphaera stadtmanae. Instead, all strains possess a novel complex that is related to the F420:methanophenazine oxidoreductase (Fpo) of Methanosarcinales butlacks an F420-oxidizing module, resembling the apparently ferredoxin-dependent Fpo-like homolog in Methanosaeta thermophila. Since all Methanomassiliicoccales also lack the subunit E of the membrane-bound heterodisulfide reductase (HdrDE), wepropose that the Fpo-like complex interacts directly with subunit D, forming an energy-converting ferredoxin: heterodisulfideoxidoreductase. The dual function of heterodisulfide in Methanomassiliicoccales, which serves both in electron bifurcation and as terminal acceptor in a membrane-associated redox process, may be a unique characteristic of the novel order.

A multienzyme complex channels substrates and electrons through acetyl-CoA and methane biosynthesis pathways in Methanosarcina

Multienzyme complexes catalyze important metabolic reactions in many organisms, but little is known about the complexes involved in biological methane production (methanogenesis). A crosslinking-mass spectrometry (XL-MS) strategy was employed to identify proteins associated with coenzyme M-coenzyme B heterodisulfide reductase (Hdr), an essential enzyme in all methane-producing archaea (methanogens). In Methanosarcina acetivorans, Hdr forms a multienzyme complex with acetyl-CoA decarbonylase synthase (ACDS), and F420-dependent methylene-H4MPT reductase (Mer). ACDS is essential for production of acetyl-CoA during growth on methanol, or for methanogenesis from acetate, whereas Mer is essential for methanogenesis from all substrates. Existence of a Hdr:ACDS:Mer complex is consistent with growth phenotypes of ACDS and Mer mutant strains in which the complex samples the redox status of electron carriers and directs carbon flux to acetyl-CoA or methanogenesis. We propose the Hdr:ACDS:Mer complex comprises a special class of multienzyme redox complex which functions as a "biological router" that physically links methanogenesis and acetyl-CoA biosynthesis pathways.

Specific inhibitors for identifying pathways for methane production from carbon monoxide by a nonadapted anaerobic mixed culture

Specific inhibitors such as 2-bromoethanesulfonate (BES) and vancomycin were employed in activity batch tests to decipher metabolic pathways that are preferentially used by a mixed anaerobic consortium (sludge from an anaerobic digester) to transform carbon monoxide (CO) into methane (CH4). We first evaluated the inhibitory effect of both BES and vancomycin on the microbial community, as well as the efficiency and stability of vancomycin at 35 °C, over time. The activity tests with CO2-H2, CO, glucose, acetate, formate, propionate, butyrate, methanol, and ethanol showed that vancomycin does not inhibit some Gram-negative bacteria, and 50 mmol/L BES effectively blocks CH4 production in the sludge. However, when sludge was incubated with propionate, butyrate, methanol, or ethanol as the sole energy and carbon source, methanogenesis was only partially inhibited by BES. Separate tests showed that 0.07 mmol/L vancomycin is enough to maintain its inhibitory efficiency and stability in the population for at least 32 days at 35 °C. Using the inhibitors above, it was demonstrated that CO conversion to CH4 is an indirect, 2-step process, in which the CO is converted first to acetate and subsequently to CH4.

Characterization of the RnfB and RnfG subunits of the Rnf complex from the archaeon Methanosarcina acetivorans

Rnf complexes are redox-driven ion pumps identified in diverse species from the domains Bacteria and Archaea, biochemical characterizations of which are reported for two species from the domain Bacteria. Here, we present characterizations of the redox-active subunits RnfG and RnfB from the Rnf complex of Methanosarcina acetivorans, an acetate-utilizing methane-producing species from the domain Archaea. The purified RnfG subunit produced in Escherichia coli fluoresced in SDS-PAGE gels under UV illumination and showed a UV-visible spectrum typical of flavoproteins. The Thr166Gly variant of RnfG was colorless and failed to fluoresce under UV illumination confirming a role for Thr166 in binding FMN. Redox titration of holo-RnfG revealed a midpoint potential of −129 mV for FMN with n = 2. The overproduced RnfG was primarily localized to the membrane of E. coli and the sequence contained a transmembrane helix. A topological analysis combining reporter protein fusion and computer predictions indicated that the C-terminal domain containing FMN is located on the outer aspect of the cytoplasmic membrane. The purified RnfB subunit produced in E. coli showed a UV-visible spectrum typical of iron-sulfur proteins. The EPR spectra of reduced RnfB featured a broad spectral shape with g values (2.06, 1.94, 1.90, 1.88) characteristic of magnetically coupled 3Fe-4S and 4Fe-4S clusters in close agreement with the iron and acid-labile sulfur content. The ferredoxin specific to the aceticlastic pathway served as an electron donor to RnfB suggesting this subunit is the entry point of electrons to the Rnf complex. The results advance an understanding of the organization and biochemical properties of the Rnf complex and lay a foundation for further understanding the overall mechanism in the pathway of methane formation from acetate.

Southern Appalachian peatlands support high archaeal diversity

Mid-latitude peatlands with a temperate climate are sparsely studied and as such represent a gap in the current knowledge base regarding archaeal populations present and their roles in these environments. Phylogenetic analysis of the archaeal populations among three peatlands in the Southern Appalachians reveal not only methanogenic species but also significant populations of thaumarchaeal and crenarchaeal-related organisms of the uncultured miscellaneous crenarchaeotal group (MCG) and the terrestrial group 1.1c, as well as deep-branching Euryarchaeota primarily within the Lake Dagow sediment and rice cluster V lineages. The Thaum/Crenarchaea and deep-branching Euryarchaea represented approximately 24-83% and 2-18%, respectively, of the total SSU rRNA clones retrieved in each library, and methanogens represented approximately 14-72% of the clones retrieved. Several taxa that are either rare or novel to acidic peatlands were detected including the euryarchaeal SM1K20 cluster and thaumarchaeal/crenarchaeal-related clusters 1.1a, C3, SAGMCG-1, pSL12, and AK59. All three major groups (methanogens, Thaumarchaea/Crenarchaea, and deep-branching Euryarchaea) were detected in the RNA library, suggesting at least a minimum level of maintenance activity. Compared to their northern counterparts, Southern Appalachian peatlands appear to harbor a relatively high diversity of Archaea and exhibit a high level of intra-site heterogeneity.

Towards a computational model of a methane producing archaeum

Progress towards a complete model of the methanogenic archaeum Methanosarcina acetivorans is reported. We characterized size distribution of the cells using differential interference contrast microscopy, finding them to be ellipsoidal with mean length and width of 2.9 μm and 2.3 μm, respectively, when grown on methanol and 30% smaller when grown on acetate. We used the single molecule pull down (SiMPull) technique to measure average copy number of the Mcr complex and ribosomes. A kinetic model for the methanogenesis pathways based on biochemical studies and recent metabolic reconstructions for several related methanogens is presented. In this model, 26 reactions in the methanogenesis pathways are coupled to a cell mass production reaction that updates enzyme concentrations. RNA expression data (RNA-seq) measured for cell cultures grown on acetate and methanol is used to estimate relative protein production per mole of ATP consumed. The model captures the experimentally observed methane production rates for cells growing on methanol and is most sensitive to the number of methyl-coenzyme-M reductase (Mcr) and methyl-tetrahydromethanopterin:coenzyme-M methyltransferase (Mtr) proteins. A draft transcriptional regulation network based on known interactions is proposed which we intend to integrate with the kinetic model to allow dynamic regulation.