Complete genome sequence of Methanothermobacter marburgensis, a methanoarchaeon model organism

Complete genome sequence of Methanothermobacter sp. DP, a hydrogenotrophic and thermophilic methanogen

Here, we report the complete genome sequence of the thermophilic hydrogenotrophic methanogen Methanothermobacter sp. DP isolated in South Korea from an anaerobic digester fed with cigarette waste. The genome consists of 1,693,285 bp, with 1,772 protein-coding genes and a GC content of 48.8%.

An integrated systems biology approach reveals differences in formate metabolism in the genus Methanothermobacter

Methanogenesis allows methanogenic archaea to generate cellular energy for their growth while producing methane. Thermophilic hydrogenotrophic species of the genus Methanothermobacter have been recognized as robust biocatalysts for a circular carbon economy and are already applied in power-to-gas technology with biomethanation, which is a platform to store renewable energy and utilize captured carbon dioxide. Here, we generated curated genome-scale metabolic reconstructions for three Methanothermobacter strains and investigated differences in the growth performance of these same strains in chemostat bioreactor experiments with hydrogen and carbon dioxide or formate as substrates. Using an integrated systems biology approach, we identified differences in formate anabolism between the strains and revealed that formate anabolism influences the diversion of carbon between biomass and methane. This finding, together with the omics datasets and the metabolic models we generated, can be implemented for biotechnological applications of Methanothermobacter in power-to-gas technology, and as a perspective, for value-added chemical production.

Molecular characterization of microorganisms with industrial potential for methane production in sludge from Kangemi sewage treatment plant, Nyeri county-Kenya

Microbial consortia under anaerobic conditions are involved in oxidizing organic matter in the sludge to produce methane gas. However, in developing countries like Kenya, these microbes have not been fully identified to target them for the efficient harnessing of biofuel. This study collected wet sludge from two anaerobic digestion lagoons 1 and 2 that were operational during sampling at Kangemi Sewage Treatment Plant, in Nyeri County, Kenya. DNA was extracted from samples using commercially available ZymoBIOMICS™ DNA Miniprep Kit and sequenced using Shotgun metagenomics. Samples were analyzed using MG-RAST software (Project ID: mgp100988), which allowed for identifying microorganisms directly involved in various stages of methanogenesis pathways. The study found hydrogenotrophic methanogens, such as Methanospirillum (32%), Methanobacterium (27%), Methanobrevibacter (27%), and Methanosarcina (32%), being predominant in the lagoon communities, whereas acetoclastic microorganisms such as the Methanoregula (22%) and the acetate oxidazing bacteria such as Clostridia (68%) were the key microbes for that pathway in the sewage digester sludge. Furthermore, Methanothermobacter (18%), Methanosarcina (21%), Methanosaeta (15%), and Methanospirillum (13%) carried out the methylotrophic pathway. In contrast, Methanosarcina (23%),Methanoregula (14%), methanosaeta (13%), and methnanoprevibacter (13%) seemed to play an important role in the final step of methane release. This study concluded that the sludge produced from the Nyeri-Kangemi WWTP harbors microbes with significant potential for biogas production. The study recommends a pilot study to investigate the efficiency of the identified microbes for biogas production.

Random transposon mutagenesis identifies genes essential for transformation in Methanococcus maripaludis

Natural transformation, the process whereby a cell acquires DNA directly from the environment, is an important driver of evolution in microbial populations, yet the mechanism of DNA uptake is only characterized in bacteria. To expand our understanding of natural transformation in archaea, we undertook a genetic approach to identify a catalog of genes necessary for transformation in Methanococcus maripaludis. Using an optimized method to generate random transposon mutants, we screened 6144 mutant strains for defects in natural transformation and identified 25 transformation-associated candidate genes. Among these are genes encoding components of the type IV-like pilus, transcription/translation associated genes, genes encoding putative membrane bound transport proteins, and genes of unknown function. Interestingly, similar genes were identified regardless of whether replicating or integrating plasmids were provided as a substrate for transformation. Using allelic replacement mutagenesis, we confirmed that several genes identified in these screens are essential for transformation. Finally, we identified a homolog of a membrane bound substrate transporter in Methanoculleus thermophilus and verified its importance for transformation using allelic replacement mutagenesis, suggesting a conserved mechanism for DNA transfer in multiple archaea. These data represent an initial characterization of the genes important for transformation which will inform efforts to understand gene flow in natural populations. Additionally, knowledge of the genes necessary for natural transformation may assist in identifying signatures of transformation machinery in archaeal genomes and aid the establishment of new model genetic systems for studying archaea.

Quantitative Analysis of Core Lipid Production in Methanothermobacter marburgensis at Different Scales

Archaeal lipids have a high biotechnological potential, caused by their high resistance to oxidative stress, extreme pH values and temperatures, as well as their ability to withstand phospholipases. Further, methanogens, a specific group of archaea, are already well-established in the field of biotechnology because of their ability to use carbon dioxide and molecular hydrogen or organic substrates. In this study, we show the potential of the model organism Methanothermobacter marburgensis to act both as a carbon dioxide based biological methane producer and as a potential supplier of archaeal lipids. Different cultivation settings were tested to gain an insight into the optimal conditions to produce specific core lipids. The study shows that up-scaling at a constant particle number (n/n = const.) seems to be a promising approach. Further optimizations regarding the length and number of the incubation periods and the ratio of the interaction area to the total liquid volume are necessary for scaling these settings for industrial purposes.

Comparative genomic analysis reveals starvation survival systems in Methanothermobacter thermautotrophicus ΔH

Methanothermobacter thermautotrophicus ΔH (MTH) is a thermophilic hydrogenotrophic methanogenic archaeon capable of reducing CO₂ with H₂ to produce methane gas. It is the potential candidate in the biomethanation of CO₂ and CO in anaerobic reactors and biogas upgrading process. However, systematic studies addressing its genome conservation and function remain scant in this genome. In this study, we have evaluated its evolutionary resemblance and metabolic discrepancy, particularly in starvation survival systems by comparing the genomic contexts with Methanothermobacter marburgensis str. Marburg (MMG) and Methanobacterium formicicum DSM 1535 (MFO). The phylogenomic analysis of this study indicated that there was a strong phylogenomic signal among MTH, MMG, and MFO in the whole-genome tree. DNA replication machinery was conserved in the MTH genome and might have evolved at different evolution rates. Genome synteny analysis observed collinearity of either gene orders or gene families has to be maintained with syntenic blocks located in the syntenic out-paralogs. A genome-wide metabolic analysis identified some unique putative metabolic subsystems in MTH, which are proposed to determine its growth characteristics in diverse environments. MTH genome comprised of 93 unique genes-coding for starvation survival and stress-response proteins. These proteins confer its adaptation to nutritional deprivation and other abiotic stresses. MTH has a typical system to withstand its growth and cell viability during stable operation and recovery after prolonged starvation. Thus, the present work will provide an insight to improve the genome refinement and metabolic reconstruction in parallel to other closely related species.

Complete Genomic Sequence of the Thermophilic and Hydrogenotrophic Methanogen Methanothermobacter sp. Strain KEPCO-1

Here, we describe the complete genome of Methanothermobacter sp. strain KEPCO-1, a thermophilic and hydrogenotrophic methanogen that was isolated from an anaerobic digester in Seoul, Republic of Korea. The genome of KEPCO-1 shares 96.98% of its sequence with Methanothermobacter marburgensis strain DSM 2133 and consists of 1,741,029 bp, with 1,822 protein-coding genes, 44 noncoding RNAs, and a GC content of 48.47%. The development of this genome will facilitate future genomic studies of KEPCO-1.

Here, we describe the complete genome of Methanothermobacter sp. strain KEPCO-1, a thermophilic and hydrogenotrophic methanogen that was isolated from an anaerobic digester in Seoul, Republic of Korea. The genome of KEPCO-1 shares 96.98% of its sequence with Methanothermobacter marburgensis strain DSM 2133 and consists of 1,741,029 bp, with 1,822 protein-coding genes, 44 noncoding RNAs, and a GC content of 48.47%. The development of this genome will facilitate future genomic studies of KEPCO-1.

Genome Analyses and Genome-Centered Metatranscriptomics of Methanothermobacter wolfeii Strain SIV6, Isolated from a Thermophilic Production-Scale Biogas Fermenter

In the thermophilic biogas-producing microbial community, the genus Methanothermobacter was previously described to be frequently abundant. The aim of this study was to establish and analyze the genome sequence of the archaeal strain Methanothermobacter wolfeii SIV6 originating from a thermophilic industrial-scale biogas fermenter and compare it to related reference genomes. The circular chromosome has a size of 1,686,891 bases, featuring a GC content of 48.89%. Comparative analyses considering three completely sequenced Methanothermobacter strains revealed a core genome of 1494 coding sequences and 16 strain specific genes for M. wolfeii SIV6, which include glycosyltransferases and CRISPR/cas associated genes. Moreover, M. wolfeii SIV6 harbors all genes for the hydrogenotrophic methanogenesis pathway and genome-centered metatranscriptomics indicates the high metabolic activity of this strain, with 25.18% of all transcripts per million (TPM) belong to the hydrogenotrophic methanogenesis pathway and 18.02% of these TPM exclusively belonging to the mcr operon. This operon encodes the different subunits of the enzyme methyl-coenzyme M reductase (EC: 2.8.4.1), which catalyzes the final and rate-limiting step during methanogenesis. Finally, fragment recruitment of metagenomic reads from the thermophilic biogas fermenter on the SIV6 genome showed that the strain is abundant (1.2%) within the indigenous microbial community. Detailed analysis of the archaeal isolate M. wolfeii SIV6 indicates its role and function within the microbial community of the thermophilic biogas fermenter, towards a better understanding of the biogas production process and a microbial-based management of this complex process.

Effect of N2 on Biological Methanation in a Continuous Stirred-Tank Reactor with Methanothermobacter marburgensis

In this contribution, the effect of the presence of a presumed inert gas like N2 in the feed gas on the biological methanation of hydrogen and carbon dioxide with Methanothermobacter marburgensis was investigated. N2 can be found as a component besides CO2 in possible feed gases like mine gas, weak gas, or steel mill gas. To determine whether there is an effect on the biological methanation of CO2 and H2 from renewable sources or not, the process was investigated using feed gases containing CO2, H2, and N2 in different ratios, depending on the CO2 content. A possible effect can be a lowered conversion rate of CO2 and H2 to CH4. Feed gases containing up to 47N2 were investigated. The conversion of hydrogen and carbon dioxide was possible with a conversion rate of up to 91 but was limited by the amount of H2 when feeding a stoichiometric ratio of 4:1 and not by adding N2 to the feed gas.

Comparative Genomic Analysis of Members of the Genera Methanosphaera and Methanobrevibacter Reveals Distinct Clades with Specific Potential Metabolic Functions

Methanobrevibacter and Methanosphaera species represent some of the most prevalent methanogenic archaea in the gastrointestinal tract of animals and humans and play an important role in this environment. The aim of this study was to identify genomic features that are shared or specific for members of each genus with a special emphasis of the analysis on the assimilation of nitrogen and acetate and the utilization of methanol and ethanol for methanogenesis. Here, draft genome sequences of Methanobrevibacter thaueri strain DSM 11995^T, Methanobrevibacter woesei strain DSM 11979^T, and Methanosphaera cuniculi strain 4103^T are reported and compared to those of 16 other Methanobrevibacter and Methanosphaera genomes, including genomes of the 13 currently available types of strains of the two genera. The comparative genome analyses indicate that among other genes, the absence of molybdopterin cofactor biosynthesis is conserved in Methanosphaera species but reveals also that the three species share a core set of more than 300 genes that distinguishes the genus Methanosphaera from the genus Methanobrevibacter. Multilocus sequence analysis shows that the genus Methanobrevibacter can be subdivided into clades, potentially new genera, which may display characteristic specific metabolic features. These features include not only the potential ability of nitrogen fixation and acetate assimilation in a clade comprised of Methanobrevibacter species from the termite gut and Methanobrevibacter arboriphilus strains but also the potential capability to utilize ethanol and methanol in a clade comprising Methanobrevibacter wolinii strain DSM 11976^T, Mbb. sp. AbM4, and Mbb. boviskoreani strain DSM 25824^T.

Cofactor Tail Length Modulates Catalysis of Bacterial F420-Dependent Oxidoreductases

F₄₂₀ is a microbial cofactor that mediates a wide range of physiologically important and industrially relevant redox reactions, including in methanogenesis and tetracycline biosynthesis. This deazaflavin comprises a redox-active isoalloxazine headgroup conjugated to a lactyloligoglutamyl tail. Here we studied the catalytic significance of the oligoglutamate chain, which differs in length between bacteria and archaea. We purified short-chain F₄₂₀ (two glutamates) from a methanogen isolate and long-chain F₄₂₀ (five to eight glutamates) from a recombinant mycobacterium, confirming their different chain lengths by HPLC and LC/MS analysis. F₄₂₀ purified from both sources was catalytically compatible with purified enzymes from the three major bacterial families of F₄₂₀-dependent oxidoreductases. However, long-chain F₄₂₀ bound to these enzymes with a six- to ten-fold higher affinity than short-chain F₄₂₀. The cofactor side chain also significantly modulated the kinetics of the enzymes, with long-chain F₄₂₀ increasing the substrate affinity (lower K_m) but reducing the turnover rate (lower k_cat) of the enzymes. Molecular dynamics simulations and comparative structural analysis suggest that the oligoglutamate chain of F₄₂₀ makes dynamic electrostatic interactions with conserved surface residues of the oxidoreductases while the headgroup binds the catalytic site. In conjunction with the kinetic data, this suggests that electrostatic interactions made by the oligoglutamate tail result in higher-affinity, lower-turnover catalysis. Physiologically, we propose that bacteria have selected for long-chain F₄₂₀ to better control cellular redox reactions despite tradeoffs in catalytic rate. Conversely, this suggests that industrial use of shorter-length F₄₂₀ will greatly increase the rates of bioremediation and biocatalysis processes relying on purified F₄₂₀-dependent oxidoreductases.

Draft Genome Sequence of Methanothermobacter sp. Strain EMTCatA1, Reconstructed from the Metagenome of a Thermophilic Electromethanogenesis-Catalyzing Biocathode

A draft genome of Methanothermobacter sp. strain EMTCatA1 was reconstructed from a metagenome of a thermophilic electromethanogenic biocathode. This genome will provide information about methanogens catalyzing methanogenesis at the biocathodes.

The Draft Genome of the Non-Host-Associated Methanobrevibacter arboriphilus Strain DH1 Encodes a Large Repertoire of Adhesin-Like Proteins

Methanobrevibacter arboriphilus strain DH1 is an autotrophic methanogen that was isolated from the wetwood of methane-emitting trees. This species has been of considerable interest for its unusual oxygen tolerance and has been studied as a model organism for more than four decades. Strain DH1 is closely related to other host-associated Methanobrevibacter species from intestinal tracts of animals and the rumen, making this strain an interesting candidate for comparative analysis to identify factors important for colonizing intestinal environments. Here, the genome sequence of M. arboriphilus strain DH1 is reported. The draft genome is composed of 2.445.031 bp with an average GC content of 25.44% and predicted to harbour 1964 protein-encoding genes. Among the predicted genes, there are also more than 50 putative genes for the so-called adhesin-like proteins (ALPs). The presence of ALP-encoding genes in the genome of this non-host-associated methanogen strongly suggests that target surfaces for ALPs other than host tissues also need to be considered as potential interaction partners. The high abundance of ALPs may also indicate that these types of proteins are more characteristic for specific phylogenetic groups of methanogens rather than being indicative for a particular environment the methanogens thrives in.

Proteomic Analysis of the Hydrogen and Carbon Monoxide Metabolism of Methanothermobacter marburgensis

Hydrogenotrophic methanogenic archaea are efficient H₂ utilizers, but only a few are known to be able to utilize CO. Methanothermobacter thermoautotrophicus is one of the hydrogenotrophic methanogens able to grow on CO, albeit about 100 times slower than on H₂ + CO₂. In this study, we show that the hydrogenotrophic methanogen Methanothermobacter marburgensis, is able to perform methanogenic growth on H₂/CO₂/CO and on CO as a sole substrate. To gain further insight in its carboxydotrophic metabolism, the proteome of M. marburgensis, grown on H₂/CO₂ and H₂/CO₂/CO, was analyzed. Cultures grown with H₂/CO₂/CO showed relative higher abundance of enzymes involved in the reductive acetyl-CoA pathway and proteins involved in redox metabolism. The data suggest that the strong reducing capacity of CO negatively affects hydrogenotrophic methanogenesis, making growth on CO as a sole substrate difficult for this type of methanogens. M. marburgensis appears to partly deal with this by up-regulating co-factor regenerating reactions and activating additional pathways allowing for formation of other products, like acetate.

Identifying novel protein phenotype annotations by hybridizing protein-protein interactions and protein sequence similarities

Studies of protein phenotypes represent a central challenge of modern genetics in the post-genome era because effective and accurate investigation of protein phenotypes is one of the most critical procedures to identify functional biological processes in microscale, which involves the analysis of multifactorial traits and has greatly contributed to the development of modern biology in the post genome era. Therefore, we have developed a novel computational method that identifies novel proteins associated with certain phenotypes in yeast based on the protein-protein interaction network. Unlike some existing network-based computational methods that identify the phenotype of a query protein based on its direct neighbors in the local network, the proposed method identifies novel candidate proteins for a certain phenotype by considering all annotated proteins with this phenotype on the global network using a shortest path (SP) algorithm. The identified proteins are further filtered using both a permutation test and their interactions and sequence similarities to annotated proteins. We compared our method with another widely used method called random walk with restart (RWR). The biological functions of proteins for each phenotype identified by our SP method and the RWR method were analyzed and compared. The results confirmed a large proportion of our novel protein phenotype annotation, and the RWR method showed a higher false positive rate than the SP method. Our method is equally effective for the prediction of proteins involving in all the eleven clustered yeast phenotypes with a quite low false positive rate. Considering the universality and generalizability of our supporting materials and computing strategies, our method can further be applied to study other organisms and the new functions we predicted can provide pertinent instructions for the further experimental verifications.

A novel cytosolic NADH:quinone oxidoreductase from Methanothermobacter marburgensis

Methanothermobacter marburgensis is a strictly anaerobic, thermophilic methanogenic archaeon that uses methanogenesis to convert H₂ and CO₂ to energy. M. marburgensis is one of the best-studied methanogens, and all genes required for methanogenic metabolism have been identified. Nonetheless, the present study describes a gene (Gene ID 9704440) coding for a putative NAD(P)H:quinone oxidoreductase that has not yet been identified as part of the metabolic machinery. The gene product, MmNQO, was successfully expressed, purified and characterized biochemically, as well as structurally. MmNQO was identified as a flavin-dependent NADH:quinone oxidoreductase with the capacity to oxidize NADH in the presence of a wide range of electron acceptors, whereas NADPH was oxidized with only three acceptors. The 1.50 Å crystal structure of MmNQO features a homodimeric enzyme where each monomer comprises 196 residues folding into flavodoxin-like α/β domains with non-covalently bound FMN (flavin mononucleotide). The closest structural homologue is the modulator of drug activity B from Streptococcus mutans with 1.6 Å root-mean-square deviation on 161 Cα atoms and 28% amino-acid sequence identity. The low similarity at sequence and structural level suggests that MmNQO is unique among NADH:quinone oxidoreductases characterized to date. Based on preliminary bioreactor experiments, MmNQO could provide a useful tool to prevent overflow metabolism in applications that require cells with high energy demand.

A novel NADH:quinone oxidoreductase, MmNQO, from Methanothermobacter marburgensis was identified. MmNQO oxidizes NADH with several electron acceptors and is structurally similar to bacterial MdaB. It is localized in the cytosol and may provide a useful tool to prevent overflow metabolism.

Annotating functional RNAs in genomes using Infernal

Many different types of functional non-coding RNAs participate in a wide range of important cellular functions but the large majority of these RNAs are not routinely annotated in published genomes. Several programs have been developed for identifying RNAs, including specific tools tailored to a particular RNA family as well as more general ones designed to work for any family. Many of these tools utilize covariance models (CMs), statistical models of the conserved sequence, and structure of an RNA family. In this chapter, as an illustrative example, the Infernal software package and CMs from the Rfam database are used to identify RNAs in the genome of the archaeon Methanobrevibacter ruminantium, uncovering some additional RNAs not present in the genome's initial annotation. Analysis of the results and comparison with family-specific methods demonstrate some important strengths and weaknesses of this general approach.

In vivo activation of methyl-coenzyme M reductase by carbon monoxide

Methyl-coenzyme M reductase (MCR) from methanogenic archaea catalyzes the rate-limiting and final step in methane biosynthesis. Using coenzyme B as the two-electron donor, MCR reduces methyl-coenzyme M (CH₃-SCoM) to methane and the mixed disulfide, CoBS-SCoM. MCR contains an essential redox-active nickel tetrahydrocorphinoid cofactor, Coenzyme F₄₃₀, at its active site. The active form of the enzyme (MCR_red1) contains Ni(I)-F₄₃₀. Rapid and efficient conversion of MCR to MCR_red1 is important for elucidating the enzymatic mechanism, yet this reduction is difficult because the Ni(I) state is subject to oxidative inactivation. Furthermore, no in vitro methods have yet been described to convert Ni(II) forms into MCR_red1. Since 1991, it has been known that MCR_red1 from Methanothermobacter marburgensis can be generated in vivo when cells are purged with 100% H₂. Here we show that purging cells or cell extracts with CO can also activate MCR. The rate of in vivo activation by CO is about 15 times faster than by H₂ (130 and 8 min^-1, respectively) and CO leads to twofold higher MCR_red1 than H₂. Unlike H₂-dependent activation, which exhibits a 10-h lag time, there is no lag for CO-dependent activation. Based on cyanide inhibition experiments, carbon monoxide dehydrogenase is required for the CO-dependent activation. Formate, which also is a strong reductant, cannot activate MCR in M. marburgensis in vivo.

De novo modeling of the F(420)-reducing [NiFe]-hydrogenase from a methanogenic archaeon by cryo-electron microscopy

Methanogenic archaea use a [NiFe]-hydrogenase, Frh, for oxidation/reduction of F420, an important hydride carrier in the methanogenesis pathway from H2 and CO2. Frh accounts for about 1% of the cytoplasmic protein and forms a huge complex consisting of FrhABG heterotrimers with each a [NiFe] center, four Fe-S clusters and an FAD. Here, we report the structure determined by near-atomic resolution cryo-EM of Frh with and without bound substrate F420. The polypeptide chains of FrhB, for which there was no homolog, was traced de novo from the EM map. The 1.2-MDa complex contains 12 copies of the heterotrimer, which unexpectedly form a spherical protein shell with a hollow core. The cryo-EM map reveals strong electron density of the chains of metal clusters running parallel to the protein shell, and the F420-binding site is located at the end of the chain near the outside of the spherical structure. DOI:http://dx.doi.org/10.7554/eLife.00218.001.

Methanothermobacter tenebrarum sp. nov., a hydrogenotrophic, thermophilic methanogen isolated from gas-associated formation water of a natural gas field

A thermophilic and hydrogenotrophic methanogen, strain RMAST, was isolated from gas-associated formation water of a gas-producing well in a natural gas field in Japan. Strain RMAST grew solely on H2/CO2 but required Casamino acids, tryptone, yeast extract or vitamins for growth. Growth of strain RMAST was stimulated by acetate. Cells were non-motile, straight rods (0.5×3.5–10.5 µm) and occurred singly or in pairs. Bundles of fimbriae occurred at both poles of cells and the cell wall was thick (approximately 21 nm, as revealed by ultrathin section electron microscopy). Strain RMAST grew at 45–80 °C (optimum, 70 °C), at pH 5.8–8.7 (optimum, pH 6.9–7.7) and with 0.001–20 g NaCl l−1 (optimum, 2.5 g NaCl l−1). Phylogenetic analysis revealed that Methanothermobacter thermautotrophicus ΔHT was most closely related to the isolate (95.7 % 16S rRNA gene sequence similarity). On the basis of morphological, phenotypic and phylogenetic characteristics, it is clear that strain RMAST represents a novel species of the genus Methanothermobacter , for which we propose the name Methanothermobacter tenebrarum sp. nov. The type strain is RMAST ( = DSM 23052T = JCM 16532T = NBRC 106236T).

More than 200 genes required for methane formation from H₂ and CO₂ and energy conservation are present in Methanothermobacter marburgensis and Methanothermobacter thermautotrophicus

The hydrogenotrophic methanogens Methanothermobacter marburgensis and Methanothermobacter thermautotrophicus can easily be mass cultured. They have therefore been used almost exclusively to study the biochemistry of methanogenesis from H₂ and CO₂, and the genomes of these two model organisms have been sequenced. The close relationship of the two organisms is reflected in their genomic architecture and coding potential. Within the 1,607 protein coding sequences (CDS) in common, we identified approximately 200 CDS required for the synthesis of the enzymes, coenzymes, and prosthetic groups involved in CO₂ reduction to methane and in coupling this process with the phosphorylation of ADP. Approximately 20 additional genes, such as those for the biosynthesis of F₄₃₀ and methanofuran and for the posttranslational modifications of the two methyl-coenzyme M reductases, remain to be identified.