Exploring Hydrogenotrophic Methanogenesis: a Genome Scale Metabolic Reconstruction of Methanococcus maripaludis.

Hydrogenotrophic methanogenesis occurs in multiple environments, ranging from the intestinal tracts of animals to anaerobic sediments and hot springs. Energy conservation in hydrogenotrophic methanogens was long a mystery; only within the last decade was it reported that net energy conservation for growth depends on electron bifurcation. In this work, we focus on Methanococcus maripaludis, a well-studied hydrogenotrophic marine methanogen. To better understand hydrogenotrophic methanogenesis and compare it with methylotrophic methanogenesis that utilizes oxidative phosphorylation rather than electron bifurcation, we have built iMR539, a genome scale metabolic reconstruction that accounts for 539 of the 1,722 protein-coding genes of M. maripaludis strain S2. Our reconstructed metabolic network uses recent literature to not only represent the central electron bifurcation reaction but also incorporate vital biosynthesis and assimilation pathways, including unique cofactor and coenzyme syntheses. We show that our model accurately predicts experimental growth and gene knockout data, with 93% accuracy and a Matthews correlation coefficient of 0.78. Furthermore, we use our metabolic network reconstruction to probe the implications of electron bifurcation by showing its essentiality, as well as investigating the infeasibility of aceticlastic methanogenesis in the network. Additionally, we demonstrate a method of applying thermodynamic constraints to a metabolic model to quickly estimate overall free-energy changes between what comes in and out of the cell. Finally, we describe a novel reconstruction-specific computational toolbox we created to improve usability. Together, our results provide a computational network for exploring hydrogenotrophic methanogenesis and confirm the importance of electron bifurcation in this process. IMPORTANCE: Understanding and applying hydrogenotrophic methanogenesis is a promising avenue for developing new bioenergy technologies around methane gas. Although a significant portion of biological methane is generated through this environmentally ubiquitous pathway, existing methanogen models portray the more traditional energy conservation mechanisms that are found in other methanogens. We have constructed a genome scale metabolic network of Methanococcus maripaludis that explicitly accounts for all major reactions involved in hydrogenotrophic methanogenesis. Our reconstruction demonstrates the importance of electron bifurcation in central metabolism, providing both a window into hydrogenotrophic methanogenesis and a hypothesis-generating platform to fuel metabolic engineering efforts.

A systems level predictive model for global gene regulation of methanogenesis in a hydrogenotrophic methanogen.

Methanogens catalyze the critical methane-producing step (called methanogenesis) in the anaerobic decomposition of organic matter. Here, we present the first predictive model of global gene regulation of methanogenesis in a hydrogenotrophic methanogen, Methanococcus maripaludis. We generated a comprehensive list of genes (protein-coding and noncoding) for M. maripaludis through integrated analysis of the transcriptome structure and a newly constructed Peptide Atlas. The environment and gene-regulatory influence network (EGRIN) model of the strain was constructed from a compendium of transcriptome data that was collected over 58 different steady-state and time-course experiments that were performed in chemostats or batch cultures under a spectrum of environmental perturbations that modulated methanogenesis. Analyses of the EGRIN model have revealed novel components of methanogenesis that included at least three additional protein-coding genes of previously unknown function as well as one noncoding RNA. We discovered that at least five regulatory mechanisms act in a combinatorial scheme to intercoordinate key steps of methanogenesis with different processes such as motility, ATP biosynthesis, and carbon assimilation. Through a combination of genetic and environmental perturbation experiments we have validated the EGRIN-predicted role of two novel transcription factors in the regulation of phosphate-dependent repression of formate dehydrogenase-a key enzyme in the methanogenesis pathway. The EGRIN model demonstrates regulatory affiliations within methanogenesis as well as between methanogenesis and other cellular functions.

Essential anaplerotic role for the energy-converting hydrogenase Eha in hydrogenotrophic methanogenesis.

Despite decades of study, electron flow and energy conservation in methanogenic Archaea are still not thoroughly understood. For methanogens without cytochromes, flavin-based electron bifurcation has been proposed as an essential energy-conserving mechanism that couples exergonic and endergonic reactions of methanogenesis. However, an alternative hypothesis posits that the energy-converting hydrogenase Eha provides a chemiosmosis-driven electron input to the endergonic reaction. In vivo evidence for both hypotheses is incomplete. By genetically eliminating all nonessential pathways of H(2) metabolism in the model methanogen Methanococcus maripaludis and using formate as an additional electron donor, we isolate electron flow for methanogenesis from flux through Eha. We find that Eha does not function stoichiometrically for methanogenesis, implying that electron bifurcation must operate in vivo. We show that Eha is nevertheless essential, and a substoichiometric requirement for H(2) suggests that its role is anaplerotic. Indeed, H(2) via Eha stimulates methanogenesis from formate when intermediates are not otherwise replenished. These results fit the model for electron bifurcation, which renders the methanogenic pathway cyclic, and as such requires the replenishment of intermediates. Defining a role for Eha and verifying electron bifurcation provide a complete model of methanogenesis where all necessary electron inputs are accounted for.

Parallel evolution of transcriptome architecture during genome reorganization.

Assembly of genes into operons is generally viewed as an important process during the continual adaptation of microbes to changing environmental challenges. However, the genome reorganization events that drive this process are also the roots of instability for existing operons. We have determined that there exists a statistically significant trend that correlates the proportion of genes encoded in operons in archaea to their phylogenetic lineage. We have further characterized how microbes deal with operon instability by mapping and comparing transcriptome architectures of four phylogenetically diverse extremophiles that span the range of operon stabilities observed across archaeal lineages: a photoheterotrophic halophile (Halobacterium salinarum NRC-1), a hydrogenotrophic methanogen (Methanococcus maripaludis S2), an acidophilic and aerobic thermophile (Sulfolobus solfataricus P2), and an anaerobic hyperthermophile (Pyrococcus furiosus DSM 3638). We demonstrate how the evolution of transcriptional elements (promoters and terminators) generates new operons, restores the coordinated regulation of translocated, inverted, and newly acquired genes, and introduces completely novel regulation for even some of the most conserved operonic genes such as those encoding subunits of the ribosome. The inverse correlation (r=-0.92) between the proportion of operons with such internally located transcriptional elements and the fraction of conserved operons in each of the four archaea reveals an unprecedented view into varying stages of operon evolution. Importantly, our integrated analysis has revealed that organisms adapted to higher growth temperatures have lower tolerance for genome reorganization events that disrupt operon structures.

VhuD facilitates electron flow from H2 or formate to heterodisulfide reductase in Methanococcus maripaludis.

Flavin-based electron bifurcation has recently been characterized as an essential energy conservation mechanism that is utilized by hydrogenotrophic methanogenic Archaea to generate low-potential electrons in an ATP-independent manner. Electron bifurcation likely takes place at the flavin associated with the α subunit of heterodisulfide reductase (HdrA). In Methanococcus maripaludis the electrons for this reaction come from either formate or H2 via formate dehydrogenase (Fdh) or Hdr-associated hydrogenase (Vhu). However, how these enzymes bind to HdrA to deliver electrons is unknown. Here, we present evidence that the δ subunit of hydrogenase (VhuD) is central to the interaction of both enzymes with HdrA. When M. maripaludis is grown under conditions where both Fdh and Vhu are expressed, these enzymes compete for binding to VhuD, which in turn binds to HdrA. Under these conditions, both enzymes are fully functional and are bound to VhuD in substoichiometric quantities. We also show that Fdh copurifies specifically with VhuD in the absence of other hydrogenase subunits. Surprisingly, in the absence of Vhu, growth on hydrogen still occurs; we show that this involves F420-reducing hydrogenase. The data presented here represent an initial characterization of specific protein interactions centered on Hdr in a hydrogenotrophic methanogen that utilizes multiple electron donors for growth.

Effects of H2 and formate on growth yield and regulation of methanogenesis in Methanococcus maripaludis.

Hydrogenotrophic methanogenic Archaea are defined by an H2 requirement for growth. Despite this requirement, many hydrogenotrophs are also capable of growth with formate as an electron donor for methanogenesis. While certain responses of these organisms to hydrogen availability have been characterized, responses to formate starvation have not been reported. Here we report that during continuous culture of Methanococcus maripaludis under defined nutrient conditions, growth yields relative to methane production decreased markedly with either H2 excess or formate excess. Analysis of the growth yields of several mutants suggests that this phenomenon occurs independently of the storage of intracellular carbon or a transcriptional response to methanogenesis. Using microarray analysis, we found that the expression of genes encoding coenzyme F420-dependent steps of methanogenesis, including one of two formate dehydrogenases, increased with H2 starvation but with formate occurred at high levels regardless of limitation or excess. One gene, encoding H2-dependent methylene-tetrahydromethanopterin dehydrogenase, decreased in expression with either H2 limitation or formate limitation. Expression of genes for the second formate dehydrogenase, molybdenum-dependent formylmethanofuran dehydrogenase, and molybdenum transport increased specifically with formate limitation. Of the two formate dehydrogenases, only the first could support growth on formate in batch culture where formate was in excess.

Protein complexing in a methanogen suggests electron bifurcation and electron delivery from formate to heterodisulfide reductase.

In methanogenic Archaea, the final step of methanogenesis generates methane and a heterodisulfide of coenzyme M and coenzyme B (CoM-S-S-CoB). Reduction of this heterodisulfide by heterodisulfide reductase to regenerate HS-CoM and HS-CoB is an exergonic process. Thauer et al. [Thauer, et al. 2008 Nat Rev Microbiol 6:579-591] recently suggested that in hydrogenotrophic methanogens the energy of heterodisulfide reduction powers the most endergonic reaction in the pathway, catalyzed by the formylmethanofuran dehydrogenase, via flavin-based electron bifurcation. Here we present evidence that these two steps in methanogenesis are physically linked. We identify a protein complex from the hydrogenotrophic methanogen, Methanococcus maripaludis, that contains heterodisulfide reductase, formylmethanofuran dehydrogenase, F(420)-nonreducing hydrogenase, and formate dehydrogenase. In addition to establishing a physical basis for the electron-bifurcation model of energy conservation, the composition of the complex also suggests that either H(2) or formate (two alternative electron donors for methanogenesis) can donate electrons to the heterodisulfide-H(2) via F(420)-nonreducing hydrogenase or formate via formate dehydrogenase. Electron flow from formate to the heterodisulfide rather than the use of H(2) as an intermediate represents a previously unknown path of electron flow in methanogenesis. We further tested whether this path occurs by constructing a mutant lacking F(420)-nonreducing hydrogenase. The mutant displayed growth equal to wild-type with formate but markedly slower growth with hydrogen. The results support the model of electron bifurcation and suggest that formate, like H(2), is closely integrated into the methanogenic pathway.

A Genetic Study of Nif-Associated Genes in a Hyperthermophilic Methanogen.

Methanocaldococcus sp. strain FS406-22, a hyperthermophilic methanogen, fixes nitrogen with a minimal set of known nif genes. Only four structural nif genes, nifH, nifD, nifK, and nifE, are present in a cluster, and a nifB homolog is present elsewhere in the genome. nifN, essential for the final synthesis of the iron-molybdenum cofactor of nitrogenase in well-characterized diazotrophs, is absent from FS406-22. In addition, FS406-22 encodes four novel hypothetical proteins, and a ferredoxin, in the nif cluster. Here, we develop a set of genetic tools for FS406-22 and test the functionality of genes in the nif cluster by making markerless in-frame deletion mutations. Deletion of the gene for one hypothetical protein, designated Hp4, delayed the initiation of diazotrophic growth and decreased the growth rate, an effect we confirmed by genetic complementation. NifE also appeared to play a role in diazotrophic growth, and the encoding of Hp4 and NifE in a single operon suggested they may work together in some way in the synthesis of the nitrogenase cofactor. No role could be discerned for any of the other hypothetical proteins, nor for the ferredoxin, despite the presence of these genes in a variety of related organisms. Possible pathways and evolutionary scenarios for the synthesis of the nitrogenase cofactor in an organism that lacks nifN are discussed. IMPORTANCEMethanocaldococcus has been considered a model genus, but genetic tools have not been forthcoming until recently. Here, we develop and illustrate the utility of positive selection with either of two selective agents (simvastatin and neomycin), negative selection, generation of markerless in-frame deletion mutations, and genetic complementation. These genetic tools should be useful for a variety of related species. We address the question of the minimal set of nif genes, which has implications for how nitrogen fixation evolved.

Overlapping repressor binding sites regulate expression of the Methanococcus maripaludis glnK(1) operon.

The euryarchaeal transcriptional repressor NrpR regulates a variety of nitrogen assimilation genes by 2-oxoglutarate-reversible binding to conserved palindromic operators. The number and positioning of these operators varies among promoter regions of regulated genes, suggesting NrpR can bind in different patterns. Particularly intriguing is the contrast between the nif and glnK(1) promoter regions of Methanococcus maripaludis, where two operators are present but with different configurations. Here we study NrpR binding and regulation at the glnK(1) promoter, where the two operator sequences overlap and occur on opposite faces of the double helix. We find that both operators function in binding, with a dimer of NrpR binding simultaneously to each overlapping operator. We show in vivo that the first operator plays a primary role in regulation and the second operator plays an enhancing role. This is the first demonstration of overlapping operators functioning in Archaea.

Interdomain conjugal transfer of DNA from bacteria to archaea.

Escherichia coli transforms the methanogenic archaeon Methanococcus maripaludis at frequencies ranging from 0.2 x 10(-6) to 2 x 10(-6) per recipient cell. Transformation requires cell-to-cell contact, oriT, and tra functions, is insensitive to DNase I, and otherwise displays hallmarks of conjugation.

Comprehensive computational analysis of Hmd enzymes and paralogs in methanogenic Archaea.

BACKGROUND: Methanogenesis is the sole means of energy production in methanogenic Archaea. H2-forming methylenetetrahydromethanopterin dehydrogenase (Hmd) catalyzes a step in the hydrogenotrophic methanogenesis pathway in class I methanogens. At least one hmd paralog has been identified in nine of the eleven complete genome sequences of class I hydrogenotrophic methanogens. The products of these paralog genes have thus far eluded any detailed functional characterization. RESULTS: Here we present a thorough computational analysis of Hmd enzymes and paralogs that includes state of the art phylogenetic inference, structure prediction, and functional site prediction techniques. We determine that the Hmd enzymes are phylogenetically distinct from Hmd paralogs but share a common overall structure. We predict that the active site of the Hmd enzyme is conserved as a functional site in Hmd paralogs and use this observation to propose possible molecular functions of the paralog that are consistent with previous experimental evidence. We also identify an uncharacterized site in the N-terminal domains of both proteins that is predicted by our methods to directly impart function. CONCLUSION: This study contributes to our understanding of the evolutionary history, structural conservation, and functional roles, of the Hmd enzymes and paralogs. The results of our phylogenetic and structural analysis constitute datasets that will aid in the future study of the Hmd protein family. Our functional site predictions generate several testable hypotheses that will guide further experimental characterization of the Hmd paralog. This work also represents a novel approach to protein function prediction in which multiple computational methods are integrated to achieve a detailed characterization of proteins that are not well understood.

Quantitative proteomics of nutrient limitation in the hydrogenotrophic methanogen Methanococcus maripaludis.

BACKGROUND: Methanogenic Archaea play key metabolic roles in anaerobic ecosystems, where they use H2 and other substrates to produce methane. Methanococcus maripaludis is a model for studies of the global response to nutrient limitations. RESULTS: We used high-coverage quantitative proteomics to determine the response of M. maripaludis to growth-limiting levels of H2, nitrogen, and phosphate. Six to ten percent of the proteome changed significantly with each nutrient limitation. H2 limitation increased the abundance of a wide variety of proteins involved in methanogenesis. However, one protein involved in methanogenesis decreased: a low-affinity [Fe] hydrogenase, which may dominate over a higher-affinity mechanism when H2 is abundant. Nitrogen limitation increased known nitrogen assimilation proteins. In addition, the increased abundance of molybdate transport proteins suggested they function for nitrogen fixation. An apparent regulon governed by the euryarchaeal nitrogen regulator NrpR is discussed. Phosphate limitation increased the abundance of three different sets of proteins, suggesting that all three function in phosphate transport. CONCLUSION: The global proteomic response of M. maripaludis to each nutrient limitation suggests a wider response than previously appreciated. The results give new insight into the function of several proteins, as well as providing information that should contribute to the formulation of a regulatory network model.

Roles of coenzyme F420-reducing hydrogenases and hydrogen- and F420-dependent methylenetetrahydromethanopterin dehydrogenases in reduction of F420 and production of hydrogen during methanogenesis.

Reduced coenzyme F420 (F420H2) is an essential intermediate in methanogenesis from CO2. During methanogenesis from H2 and CO2, F420H2 is provided by the action of F420-reducing hydrogenases. However, an alternative pathway has been proposed, where H2-dependent methylenetetrahydromethanopterin dehydrogenase (Hmd) and F420H2-dependent methylenetetrahydromethanopterin dehydrogenase (Mtd) together reduce F420 with H2. Here we report the construction of mutants of Methanococcus maripaludis that are defective in each putative pathway. Their analysis demonstrates that either pathway supports growth on H2 and CO2. Furthermore, we show that during growth on formate instead of H2, where F420H2 is a direct product of formate oxidation, H2 production occurs. H2 presumably arises from the oxidation of F420H2, and the analysis of the mutants during growth on formate suggests that this too can occur by either pathway. We designate the alternative pathway for the interconversion of H2 and F420H2 the Hmd-Mtd cycle.

Global responses of Methanococcus maripaludis to specific nutrient limitations and growth rate.

Continuous culture, transcriptome arrays, and measurements of cellular amino acid pools and tRNA charging levels were used to determine the response of Methanococcus maripaludis to leucine limitation. For comparison, the responses to phosphate and H2 limitations were measured as well. In addition, the effect of growth rate was determined. Leucine limitation resulted in a broad response. tRNA(Leu) charging decreased, but only small increases in mRNA were seen for amino acid biosynthesis genes. However, the cellular levels of free isoleucine and valine showed significant increases, indicating a coordinate regulation of branched-chain amino acids at a post-mRNA level. Leucine limitation also resulted in increased mRNA abundance for ribosomal protein genes, increased rRNA abundance, and decreased mRNA abundance for genes of methanogenesis. In contrast, phosphate limitation induced a specific response, a marked increase in mRNA levels for a phosphate transporter. Some mRNA levels responded to more than one factor; for example, transcripts for flagellum synthesis genes decreased under conditions of leucine limitation and increased under H2 limitation. Increased growth rate resulted in increased mRNA levels for ribosomal protein genes, increased rRNA abundance, and increased mRNA for a gene encoding an S-layer protein.

Formate-dependent H2 production by the mesophilic methanogen Methanococcus maripaludis.

Methanococcus maripaludis, an H(2)- and formate-utilizing methanogen, produced H(2) at high rates from formate. The rates and kinetics of H(2) production depended upon the growth conditions, and H(2) availability during growth was a major factor. Specific activities of resting cells grown with formate or H(2) were 0.4 to 1.4 U mg(-1) (dry weight). H(2) production in formate-grown cells followed Michaelis-Menten kinetics, and the concentration of formate required for half-maximal activity (K(f)) was 3.6 mM. In contrast, in H(2)-grown cells this process followed sigmoidal kinetics, and the K(f) was 9 mM. A key enzyme for formate-dependent H(2) production was formate dehydrogenase, Fdh. H(2) production and growth were severely reduced in a mutant containing a deletion of the gene encoding the Fdh1 isozyme, indicating that it was the primary Fdh. In contrast, a mutant containing a deletion of the gene encoding the Fdh2 isozyme possessed near-wild-type activities, indicating that this isozyme did not play a major role. H(2) production by a mutant containing a deletion of the coenzyme F(420)-reducing hydrogenase Fru was also severely reduced, suggesting that the major pathway of H(2) production comprised Fdh1 and Fru. Because a Deltafru-Deltafrc mutant retained 10% of the wild-type activity, an additional pathway is present. Mutants possessing deletions of the gene encoding the F(420)-dependent methylene-H(4)MTP dehydrogenase (Mtd) or the H(2)-forming methylene-H(4)MTP dehydrogenase (Hmd) also possessed reduced activity, which suggested that this second pathway was comprised of Fdh1-Mtd-Hmd. In contrast to H(2) production, the cellular rates of methanogenesis were unaffected in these mutants, which suggested that the observed H(2) production was not a direct intermediate of methanogenesis. In conclusion, high rates of formate-dependent H(2) production demonstrated the potential of M. maripaludis for the microbial production of H(2) from formate.

Metabolic modeling of a mutualistic microbial community.

The rate of production of methane in many environments depends upon mutualistic interactions between sulfate-reducing bacteria and methanogens. To enhance our understanding of these relationships, we took advantage of the fully sequenced genomes of Desulfovibrio vulgaris and Methanococcus maripaludis to produce and analyze the first multispecies stoichiometric metabolic model. Model results were compared to data on growth of the co-culture on lactate in the absence of sulfate. The model accurately predicted several ecologically relevant characteristics, including the flux of metabolites and the ratio of D. vulgaris to M. maripaludis cells during growth. In addition, the model and our data suggested that it was possible to eliminate formate as an interspecies electron shuttle, but hydrogen transfer was essential for syntrophic growth. Our work demonstrated that reconstructed metabolic networks and stoichiometric models can serve not only to predict metabolic fluxes and growth phenotypes of single organisms, but also to capture growth parameters and community composition of simple bacterial communities.

NifI inhibits nitrogenase by competing with Fe protein for binding to the MoFe protein.

Reduction of substrate by nitrogenase requires direct electron transfer from the Fe protein to the MoFe protein. Inhibition of nitrogenase activity in Methanococcus maripaludis occurs when the regulatory protein NifI(1,2) binds the MoFe protein. This inhibition is relieved by 2-oxoglutarate. Here we present evidence that NifI(1,2) binding prevents association of the two nitrogenase components. Increasing amounts of Fe protein competed with NifI(1,2), decreasing its inhibitory effect. NifI(1,2) prevented the co-purification of MoFe protein with a mutant form of the Fe protein that forms a stable complex with the MoFe protein, and NifI(1,2) was unable to bind to an AlF4--stabilized Fe protein:MoFe protein complex. NifI(1,2) inhibited ATP- and MoFe protein-dependent oxidation of the Fe protein, and 2OG relieved this inhibition. These results support a model where NifI(1,2) competes with the Fe protein for binding to MoFe protein and prevents electron transfer.

Diverse homologues of the archaeal repressor NrpR function similarly in nitrogen regulation.

NrpR is a transcriptional repressor of nitrogen assimilation genes that was recently discovered and characterized in the methanogenic archaeon Methanococcus maripaludis. NrpR homologues are widely distributed in Euryarchaeota and present in a few bacterial species. They exist in three different domain configurations: a single ORF encoding one NrpR domain following an N-terminal helix-turn-helix (HTH); a single ORF encoding two NrpR domains fused in tandem following an N-terminal HTH; and two separate ORFs, one with a single domain following an N-terminal HTH and one with a single domain without a HTH. Phylogenetic analysis indicated that the NrpR family forms five distinct groups: the single domain HTH type, the two domains of the double domain HTH type and the two separately encoded domains. To determine the function of diverse NrpR homologues, representative genes in were expressed an Methanococcus maripaludis nrpR deletion mutant. Homologues from species that possess a single gene restored regulated repression, regardless of domain structure. In the case of Methanosarcina acetivorans that contains two genes, both were required. The results show that distantly related NrpR homologues that are present in widely dispersed phyla regulate the expression of nitrogen assimilation genes in a similar fashion.

Continuous culture of Methanococcus maripaludis under defined nutrient conditions.

To study global regulation in the methanogenic archaeon Methanococcus maripaludis, we devised a system for steady-state growth in chemostats. New Brunswick Bioflo 110 bioreactors were equipped with controlled delivery of hydrogen, nitrogen, carbon dioxide, hydrogen sulfide, and anaerobic medium. We determined conditions and media compositions for growth with three different limiting nutrients, hydrogen, phosphate, and leucine. To investigate leucine limitation we constructed and characterized a mutant in the leuA gene for 2-isopropylmalate synthase, demonstrating for the first time the function of this gene in the Archaea. Steady state specific growth rates in these studies ranged from 0.042 to 0.24 h(-1). Plots of culture density vs. growth rate for each condition showed the behavior predicted by growth modeling. The results show that growth behavior is normal and reproducible and validate the use of the chemostat system for metabolic and global regulation studies in M. maripaludis.

Metabolic versatility in methanogens.

Methanogenesis is an anaerobic metabolism responsible for the generation of >90% of the methane formed on Earth today, with important implications for fuels production and global warming. Although methanogenic Archaea have been cultured for over 70 years, key insights regarding electron flow and energy conservation in methanogenesis have only recently emerged. Fundamental differences between two metabolic types of methanogenesis, hydrogenotrophic and methylotrophic, are now understood, with implications for metabolic versatility and the potential for engineering of methanogens to utilize new substrates. The development of model species with genetic and bioinformatic tools has advanced the field and holds potential for further characterizing and engineering of methanogenesis. Our understanding of a related pathway, anaerobic methane oxidation, is in its infancy.