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.

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.