Differential expression of methanogenesis genes of Methanothermobacter thermoautotrophicus (formerly Methanobacterium thermoautotrophicum) in pure culture and in cocultures with fatty acid-oxidizing syntrophs

Core cooperative metabolism in low-complexity CO2-fixing anaerobic microbiota

Biological conversion of carbon dioxide into methane has a crucial role in global carbon cycling and is operated by a specialised set of anaerobic archaea. Although it is known that this conversion is strictly linked with cooperative bacterial activity, such as through syntrophic acetate oxidation, there is also a limited understanding on how this cooperation is regulated and metabolically realised. In this work, we investigate the activity in a microbial community evolved to efficiently convert carbon dioxide into methane and predominantly populated by Methanothermobacter wolfeii. Through multi-omics, biochemical analysis and constraint-based modelling, we identify a potential formate cross-feeding from an uncharacterised Limnochordia species to M. wolfeii, driven by the recently discovered reductive glycine pathway and upregulated when hydrogen and carbon dioxide are limited. The quantitative consistency of this metabolic exchange with experimental data is shown by metagenome-scale metabolic models integrating condition-specific metatranscriptomics, which also indicate a broader three-way interaction involving M. wolfeii, the Limnochordia species, and Sphaerobacter thermophilus. Under limited hydrogen and carbon dioxide, aspartate released by M. wolfeii is fermented by Sphaerobacter thermophilus into acetate, which in turn is convertible into formate by Limnochordia, possibly forming a cooperative loop sustaining hydrogenotrophic methanogenesis. These findings expand our knowledge on the modes of carbon dioxide reduction into methane within natural microbial communities and provide an example of cooperative plasticity surrounding this process.

Energy Availability Determines Strategy of Microbial Amino Acid Synthesis in Volatile Fatty Acid-Fed Anaerobic Methanogenic Chemostats

In natural communities, microbes exchange a variety of metabolites (public goods) with each other, which drives the evolution of auxotroph and shapes interdependent patterns at community-level. However, factors that determine the strategy of public goods synthesis for a given community member still remains to be elucidated. In anaerobic methanogenic communities, energy availability of different community members is largely varied. We hypothesized that this uneven energy availability contributed to the heterogeneity of public goods synthesis ability among the members in these communities. We tested this hypothesis by analyzing the synthetic strategy of amino acids of the bacterial and archaeal members involved in four previously enriched anaerobic methanogenic communities residing in thermophilic chemostats. Our analyses indicate that most of the members in the communities did not possess ability to synthesize all the essential amino acids, suggesting they exchanged these essential public goods to establish interdependent patterns for survival. Importantly, we found that the amino acid synthesis ability of a functional group was largely determined by how much energy it could obtain from its metabolism in the given environmental condition. Moreover, members within a functional group also possessed different amino acid synthesis abilities, which are related to their features of energy metabolism. Our study reveals that energy availability is a key driver of microbial evolution in presence of metabolic specialization at community level and suggests the feasibility of managing anaerobic methanogenic communities for better performance through controlling the metabolic interactions involved.

Co-culturing microbial consortia: approaches for applications in biomanufacturing and bioprocessing

The application of microbial co-cultures is now recognized in the fields of biotechnology, ecology, and medicine. Understanding the biological interactions that govern the association of microorganisms would shape the way in which artificial/synthetic co-cultures or consortia are developed. The ability to accurately predict and control cell-to-cell interactions fully would be a significant enabler in synthetic biology. Co-culturing method development holds the key to strategically engineer environments in which the co-cultured microorganism can be monitored. Various approaches have been employed which aim to emulate the natural environment and gain access to the untapped natural resources emerging from cross-talk between partners. Amongst these methods are the use of a communal liquid medium for growth, use of a solid-liquid interface, membrane separation, spatial separation, and use of microfluidics systems. Maximizing the information content of interactions monitored is one of the major challenges that needs to be addressed by these designs. This review critically evaluates the significance and drawbacks of the co-culturing approaches used to this day in biotechnological applications, relevant to biomanufacturing. It is recommended that experimental results for a co-cultured species should be validated with different co-culture approaches due to variations in interactions that could exist as a result of the culturing method selected.

Analysis of a Methanogen and an Actinobacterium Dominating the Thermophilic Microbial Community of an Electromethanogenic Biocathode

Electromethanogenesis refers to the bioelectrochemical synthesis of methane from CO₂ by biocathodes. In an electromethanogenic system using thermophilic microorganisms, metagenomic analysis along with quantitative real-time polymerase chain reaction and fluorescence in situ hybridization revealed that the biocathode microbiota was dominated by the methanogen Methanothermobacter sp. strain EMTCatA1 and the actinobacterium Coriobacteriaceae sp. strain EMTCatB1. RNA sequencing was used to compare the transcriptome profiles of each strain at the methane-producing biocathodes with those in an open circuit and with the methanogenesis inhibitor 2-bromoethanesulfonate (BrES). For the methanogen, genes related to hydrogenotrophic methanogenesis were highly expressed in a manner similar to those observed under H₂-limited conditions. For the actinobacterium, the expression profiles of genes encoding multiheme c-type cytochromes and membrane-bound oxidoreductases suggested that the actinobacterium directly takes up electrons from the electrode. In both strains, various stress-related genes were commonly induced in the open-circuit biocathodes and biocathodes with BrES. This study provides a molecular inventory of the dominant species of an electromethanogenic biocathode with functional insights and therefore represents the first multiomics characterization of an electromethanogenic biocathode.

Koleobacter methoxysyntrophicus gen. nov., sp. nov., a novel anaerobic bacterium isolated from deep subsurface oil field and proposal of Koleobacteraceae fam. nov. and Koleobacterales ord. nov. within the class Clostridia of the phylum Firmicutes

An anaerobic thermophilic, rod-shaped bacterium possessing a unique non-lipid sheathed-like structure enveloping a single-membraned cell, designated strain NRmbB1^T was isolated from at the deep subsurface oil field located in Yamagata Prefecture, Japan. Growth occurred with 40-60°C (optimum, 55°C), 0-2% (2%), NaCl and pH 6.0-8.5 (8.0). Fermentative growth with various sugars was observed. Glucose-grown cells generated acetate, hydrogen, pyruvate and lactate as the main end products. Syntrophic growth occurred with glucose, pyruvate and 3,4,5-trimethoxybenzoate in the presence of an H₂-scavenging partner, and growth on 3,4,5-trimethoxybenzoate was only observed under syntrophic condition. The predominant cellular fatty acids were C_16:0, iso-C_16:0, anteiso-C_15:0, and iso-C_14:0. Respiratory quinone was not detected. The genomic G+C content was 40.8mol%. Based on 16S rRNA gene phylogeny, strain NRmbB1^T belongs to a distinct order-level clade in the class Clostridia of the phylum Firmicutes, sharing low similarity with other isolated organisms (i.e., 87.5% for top hit Moorella thermoacetica DSM 2955^T). In total, chemotaxonomic, phylogenetic and genomic characterization revealed that strain NRmbB1^T (=KCTC 25035^T, =JCM 39120^T) represents a novel species of a new genus. In addition, we also propose the associated family and order as Koleobacteraceae fam. nov and Koleobacterales ord. nov., respectively.

An iron corrosion-assisted H2-supplying system: a culture method for methanogens and acetogens under low H2 pressures

H₂ is an important fermentation intermediate in anaerobic environments. Although H₂ occurs at very low partial pressures in the environments, the culture and isolation of H₂-utilizing microorganisms is usually carried out under very high H₂ pressures, which might have hampered the discovery and understanding of microorganisms adapting to low H₂ environments. Here we constructed a culture system designated the "iron corrosion-assisted H₂-supplying (iCH) system" by connecting the gas phases of two vials (one for the iron corrosion reaction and the other for culturing microorganisms) to achieve cultures of microorganisms under low H₂ pressures. We conducted enrichment cultures for methanogens and acetogens using rice paddy field soil as the microbial source. In the enrichment culture of methanogens under canonical high H₂ pressures, only Methanobacterium spp. were enriched. By contrast, Methanocella spp. and Methanoculleus spp., methanogens adapting to low H₂ pressures, were specifically enriched in the iCH cultures. We also observed selective enrichment of acetogen species by the iCH system (Acetobacterium spp. and Sporomusa spp.), whereas Clostridium spp. predominated in the high H₂ cultures. These results demonstrate that the iCH system facilitates culture of anaerobic microorganisms under low H₂ pressures, which will enable the selective culture of microorganisms adapting to low H₂ environments.

Integration of absolute multi-omics reveals dynamic protein-to-RNA ratios and metabolic interplay within mixed-domain microbiomes

While the field of microbiology has adapted to the study of complex microbiomes via modern meta-omics techniques, we have not updated our basic knowledge regarding the quantitative levels of DNA, RNA and protein molecules within a microbial cell, which ultimately control cellular function. Here we report the temporal measurements of absolute RNA and protein levels per gene within a mixed bacterial-archaeal consortium. Our analysis of this data reveals an absolute protein-to-RNA ratio of 10²–10⁴ for bacterial populations and 10³–10⁵ for an archaeon, which is more comparable to Eukaryotic representatives’ humans and yeast. Furthermore, we use the linearity between the metaproteome and metatranscriptome over time to identify core functional guilds, hence using a fundamental biological feature (i.e., RNA/protein levels) to highlight phenotypical complementarity. Our findings show that upgrading multi-omic toolkits with traditional absolute measurements unlocks the scaling of core biological questions to dynamic and complex microbiomes, creating a deeper insight into inter-organismal relationships that drive the greater community function.

Here, the authors perform a temporal multi-omic analysis of a minimalistic cellulose-degrading and methane-producing consortium at the strain level and estimate protein-to-RNA ratios and RNA-protein dynamics of the community simultaneously over time.

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.

Energy-Conserving Enzyme Systems Active During Syntrophic Acetate Oxidation in the Thermophilic Bacterium Thermacetogenium phaeum

The thermophilic acetogen Thermacetogenium phaeum uses the Wood-Ljungdahl pathway (WLP) in both directions, either for the production of acetate from various compounds or for the oxidation of acetate in syntrophic cooperation with methanogens. In this study, energy-conserving enzyme systems in T. phaeum were investigated in both metabolic directions. A gene cluster containing a membrane-bound periplasmically oriented formate dehydrogenase directly adjacent to putative menaquinone synthesis genes was identified in the genome. The protein products of these genes were identified by total proteome analysis, and menaquinone MK-7 had been found earlier as the dominant quinone in the membrane. Enzyme assays with membrane preparations and anthraquinone-2,6-disulfonate as electron acceptor verified the presence of a quinone-dependent formate dehydrogenase. A quinone-dependent methylene-THF reductase is active in the soluble fraction and in the membrane fraction. From these results we conclude a reversed electron transport system from methyl-THF oxidation to CO₂ reduction yielding formate as reduced product which is transferred to the methanogenic partner. The redox potential difference between methyl-THF (E_o' = -200 mV) and formate (E_o' = -432 mV) does not allow electron transfer through syntrophic formate removal alone. We postulate that part of the ATP conserved by substrate-level phosphorylation has to be invested into the generation of a transmembrane proton gradient by ATPase. This proton gradient could drive the endergonic oxidation of methyl-THF in an enzyme reaction similar to the membrane-bound reversed electron transport system previously observed in the syntrophically butyrate-oxidizing bacterium Syntrophomonas wolfei. To balance the overall ATP budget in acetate oxidation, we postulate that acetate is activated through an ATP-independent path via aldehyde:ferredoxin oxidoreductase (AOR) and subsequent oxidation of acetaldehyde to acetyl-CoA.

Growth Kinetics, Carbon Isotope Fractionation, and Gene Expression in the Hyperthermophile Methanocaldococcus jannaschii during Hydrogen-Limited Growth and Interspecies Hydrogen Transfer

Hyperthermophilic methanogens are often H2 limited in hot subseafloor environments, and their survival may be due in part to physiological adaptations to low H2 conditions and interspecies H2 transfer. The hyperthermophilic methanogen Methanocaldococcus jannaschii was grown in monoculture at high (80 to 83 μM) and low (15 to 27 μM) aqueous H2 concentrations and in coculture with the hyperthermophilic H2 producer Thermococcus paralvinellae The purpose was to measure changes in growth and CH4 production kinetics, CH4 fractionation, and gene expression in M. jannaschii with changes in H2 flux. Growth and cell-specific CH4 production rates of M. jannaschii decreased with decreasing H2 availability and decreased further in coculture. However, cell yield (cells produced per mole of CH4 produced) increased 6-fold when M. jannaschii was grown in coculture rather than monoculture. Relative to high H2 concentrations, isotopic fractionation of CO2 to CH4 (εCO2-CH4) was 16‰ larger for cultures grown at low H2 concentrations and 45‰ and 56‰ larger for M. jannaschii growth in coculture on maltose and formate, respectively. Gene expression analyses showed H2-dependent methylene-tetrahydromethanopterin (H4MPT) dehydrogenase expression decreased and coenzyme F420-dependent methylene-H4MPT dehydrogenase expression increased with decreasing H2 availability and in coculture growth. In coculture, gene expression decreased for membrane-bound ATP synthase and hydrogenase. The results suggest that H2 availability significantly affects the CH4 and biomass production and CH4 fractionation by hyperthermophilic methanogens in their native habitats.IMPORTANCE Hyperthermophilic methanogens and H2-producing heterotrophs are collocated in high-temperature subseafloor environments, such as petroleum reservoirs, mid-ocean ridge flanks, and hydrothermal vents. Abiotic flux of H2 can be very low in these environments, and there is a gap in our knowledge about the origin of CH4 in these habitats. In the hyperthermophile Methanocaldococcus jannaschii, growth yields increased as H2 flux, growth rates, and CH4 production rates decreased. The same trend was observed increasingly with interspecies H2 transfer between M. jannaschii and the hyperthermophilic H2 producer Thermococcus paralvinellae With decreasing H2 availability, isotopic fractionation of carbon during methanogenesis increased, resulting in isotopically more negative CH4 with a concomitant decrease in H2-dependent methylene-tetrahydromethanopterin dehydrogenase gene expression and increase in F420-dependent methylene-tetrahydromethanopterin dehydrogenase gene expression. The significance of our research is in understanding the nature of hyperthermophilic interspecies H2 transfer and identifying biogeochemical and molecular markers for assessing the physiological state of methanogens and possible source of CH4 in natural environments.

Alternative Pathways of Acetogenic Ethanol and Methanol Degradation in the Thermophilic Anaerobe Thermacetogenium phaeum

Growth of the anaerobic thermophile Thermacetogenium phaeum with methanol, ethanol, ethanolamine, and acetate was investigated in axenic cultures and in syntrophic cultures with Methanothermobacter thermautotrophicus. Microcompartment genes were identified in the T. phaeum genome, and presence of microcompartments was confirmed by transmission electron microscopy and proteome analysis. These genes were expressed only during growth with ethanolamine. Proteome data were compared after growth with all four substrates, and activities of key enzymes of the Wood–Ljungdahl pathway and of enzyme systems leading to production or degradation of acetaldehyde such as alcohol dehydrogenase, aldehyde:ferredoxin oxidoreductase, acetate kinase, and phosphate acetyltransferase were measured in cytoplasmic fractions. Accounting of fermentation stoichiometries and growth yields with all four substrates showed that ethanol and methanol oxidation follow the same stoichiometries as in Acetobacterium woodii. On the other hand, the pathways of ethanol and methanol degradations vary between both organisms. Growth yields of T. phaeum were substantially lower than reported for A. woodii. Since T. phaeum has no Rnf complex encoded in its genome, the mechanisms of ATP synthesis have to be different from those of A. woodii. In addition to the central degradation pathways also found in A. woodii, T. phaeum maintains enzyme systems that compensate for the absence of an Rnf-complex but which on the other hand cause a loss of energy. On the basis of our data, pathways of methanol and ethanol degradation in T. phaeum are discussed.

Comparative proteomic analysis of Methanothermobacter thermautotrophicus reveals methane formation from H2 and CO2 under different temperature conditions

The growth of all methanogens is limited to a specific temperature range. However, Methanothermobacter thermautotrophicus can be found in a variety of natural and artificial environments, the temperatures of which sometimes even exceed the temperature growth ranges of thermophiles. As a result, the extent to which methane production and survival are affected by temperature remains unclear. To investigate the mechanisms of methanogenesis that Archaea have evolved to cope with drastic temperature shifts, the responses of Methanothermobacter thermautotrophicus to temperature were investigated under a high temperature growth (71°C) and cold shock (4°C) using Isobaric tags for relative and absolute quantitation (iTRAQ). The results showed that methane formation is decreased and that protein folding and degradation are increased in both high‐ and low‐temperature treatments. In addition, proteins predicted to be involved in processing environmental information processing and in cell membrane/wall/envelope biogenesis may play key roles in affecting methane formation and enhancing the response of M. thermautotrophicus to temperature stress. Analysis of the genomic locations of the genes corresponding to these temperature‐dependent proteins predicted that 77 of the genes likely to form 32 gene clusters. Here, we assess the response of M. thermautotrophicus to different temperatures and provide a new level of understanding of methane formation and cellular putative adaptive responses.

Effect of ammonia on the active microbiome and metagenome from stable full-scale digesters

Four full-scale anaerobic digesters with a long history of stable operation were characterized in terms of active microbiome and metagenome. Isotopic fractionation of biogas demonstrated that acetotrophy was rather prevalent in reactors operated at <3 gTAN L^-1 while hydrogenotrophy was predominant at >6 gTAN L^-1, suggesting that syntrophic acetate oxidizing bacteria (SAOB) played a significant role in the latter. These results were generally coherent with the observed active bacterial and archaeal communities but no known SAOB were observed. Metagenome descriptions yielded 73 assembled population genomes, of which only 7 could be assigned at the species level. Gene annotation and association to relevant metabolic pathways indicated that the phyla Chloroflexi and Bacteroidales might encompass new, currently undescribed, SAOB/formate producing species that would metabolize acetate via the glycine cleavage system. The predominant hydrogenotrophic counterpart at a high ammonia content belonged to the genus Methanoculleus, which could also grow on acetate to a certain extent.

Non-autotrophic methanogens dominate in anaerobic digesters

Anaerobic digesters are man-made habitats for fermentative and methanogenic microbes, and are characterized by extremely high concentrations of organics. However, little is known about how microbes adapt to such habitats. In the present study, we report phylogenetic, metagenomic, and metatranscriptomic analyses of microbiomes in thermophilic packed-bed digesters fed acetate as the major substrate, and we have shown that acetoclastic and hydrogenotrophic methanogens that utilize acetate as a carbon source dominate there. Deep sequencing and precise binning of the metagenomes reconstructed complete genomes for two dominant methanogens affiliated with the genera Methanosarcina and Methanothermobacter, along with 37 draft genomes. The reconstructed Methanosarcina genome was almost identical to that of a thermophilic acetoclastic methanogen Methanosarcina thermophila TM-1, indicating its cosmopolitan distribution in thermophilic digesters. The reconstructed Methanothermobacter (designated as Met2) was closely related to Methanothermobacter tenebrarum, a non-autotrophic hydrogenotrophic methanogen that grows in the presence of acetate. Met2 lacks the Cdh complex required for CO2 fixation, suggesting that it requires organic molecules, such as acetate, as carbon sources. Although the metagenomic analysis also detected autotrophic methanogens, they were less than 1% in abundance of Met2. These results suggested that non-autotrophic methanogens preferentially grow in anaerobic digesters containing high concentrations of organics.

Complete genome sequence of Methanoculleus bourgensis strain MAB1, the syntrophic partner of mesophilic acetate-oxidising bacteria (SAOB)

Methanoculleus bourgensis strain MAB1 has been identified as the hydrogenotrophic partner of mesophilic acetate-oxidising bacteria, a syntrophic relationship operating close to the thermodynamic equilibrium and of considerable importance in ammonia-rich engineered biogas processes. Methanoculleus bourgensis strain MAB1 belongs to the order Methanomicrobiales, family Methanomicrobiaceae, within the phylum Euryarchaeota. The genome shows a total size of 2,859,299 bp encoding 3450 predicted protein-encoding genes, of which only 1472 (43 %) have been assigned tentative functions. The genome encodes further 44 tRNA genes and three rRNA genes (5S, 16S and 23S rRNA). This study presents assembling and annotation features as well as genomic traits related to ammonia tolerance and methanogenesis.

Bovine Host Genetic Variation Influences Rumen Microbial Methane Production with Best Selection Criterion for Low Methane Emitting and Efficiently Feed Converting Hosts Based on Metagenomic Gene Abundance

Methane produced by methanogenic archaea in ruminants contributes significantly to anthropogenic greenhouse gas emissions. The host genetic link controlling microbial methane production is unknown and appropriate genetic selection strategies are not developed. We used sire progeny group differences to estimate the host genetic influence on rumen microbial methane production in a factorial experiment consisting of crossbred breed types and diets. Rumen metagenomic profiling was undertaken to investigate links between microbial genes and methane emissions or feed conversion efficiency. Sire progeny groups differed significantly in their methane emissions measured in respiration chambers. Ranking of the sire progeny groups based on methane emissions or relative archaeal abundance was consistent overall and within diet, suggesting that archaeal abundance in ruminal digesta is under host genetic control and can be used to genetically select animals without measuring methane directly. In the metagenomic analysis of rumen contents, we identified 3970 microbial genes of which 20 and 49 genes were significantly associated with methane emissions and feed conversion efficiency respectively. These explained 81% and 86% of the respective variation and were clustered in distinct functional gene networks. Methanogenesis genes (e.g. mcrA and fmdB) were associated with methane emissions, whilst host-microbiome cross talk genes (e.g. TSTA3 and FucI) were associated with feed conversion efficiency. These results strengthen the idea that the host animal controls its own microbiota to a significant extent and open up the implementation of effective breeding strategies using rumen microbial gene abundance as a predictor for difficult-to-measure traits on a large number of hosts. Generally, the results provide a proof of principle to use the relative abundance of microbial genes in the gastrointestinal tract of different species to predict their influence on traits e.g. human metabolism, health and behaviour, as well as to understand the genetic link between host and microbiome.

Methane is a highly potent greenhouse gas and ruminants are the major source of methane emissions from anthropogenic activities. Here we show in an experiment with cattle that genetic selection of low-emitting animals is a viable option based on a newly developed selection criterion. The experimental data provided a comprehensive insight into the host additive genetic influence on the microbiome, the impact of nutrition on genetics and the microbiome, and the effect of metagenomic microbial genes on the analysed traits. We developed a selection criterion to change those traits by evaluation of hosts based on the relative abundance of microbial genes. This criterion is shown to be highly informative and it is therefore suggested to be used in studies analysing different traits and species. This study provides a proof of principle that there is an additive genetic influence of the host on its microbiome and that selection for the desired host can be based on the abundance of a suite of genes in the ruminal metagenome associated with the trait. The use of this criterion will allow genetic analysis on a large number of hosts, previously a significant barrier to determination of host genetic effects on such traits.

RUMINANT NUTRITION SYMPOSIUM: Use of genomics and transcriptomics to identify strategies to lower ruminal methanogenesis

Globally, methane (CH4) emissions account for 40% to 45% of greenhouse gas emissions from ruminant livestock, with over 90% of these emissions arising from enteric fermentation. Reduction of carbon dioxide to CH4 is critical for efficient ruminal fermentation because it prevents the accumulation of reducing equivalents in the rumen. Methanogens exist in a symbiotic relationship with rumen protozoa and fungi and within biofilms associated with feed and the rumen wall. Genomics and transcriptomics are playing an increasingly important role in defining the ecology of ruminal methanogenesis and identifying avenues for its mitigation. Metagenomic approaches have provided information on changes in abundances as well as the species composition of the methanogen community among ruminants that vary naturally in their CH4 emissions, their feed efficiency, and their response to CH4 mitigators. Sequencing the genomes of rumen methanogens has provided insight into surface proteins that may prove useful in the development of vaccines and has allowed assembly of biochemical pathways for use in chemogenomic approaches to lowering ruminal CH4 emissions. Metagenomics and metatranscriptomic analysis of entire rumen microbial communities are providing new perspectives on how methanogens interact with other members of this ecosystem and how these relationships may be altered to reduce methanogenesis. Identification of community members that produce antimethanogen agents that either inhibit or kill methanogens could lead to the identification of new mitigation approaches. Discovery of a lytic archaeophage that specifically lyses methanogens is 1 such example. Efforts in using genomic data to alter methanogenesis have been hampered by a lack of sequence information that is specific to the microbial community of the rumen. Programs such as Hungate1000 and the Global Rumen Census are increasing the breadth and depth of our understanding of global ruminal microbial communities, steps that are key to using these tools to further define the science of ruminal methanogenesis.

Response of a rice paddy soil methanogen to syntrophic growth as revealed by transcriptional analyses

Members of Methanocellales are widespread in paddy field soils and play the key role in methane production. These methanogens feature largely in these organisms’ adaptation to low H2 and syntrophic growth with anaerobic fatty acid oxidizers. The adaptive mechanisms, however, remain unknown. In the present study, we determined the transcripts of 21 genes involved in the key steps of methanogenesis and acetate assimilation of Methanocella conradii HZ254, a strain recently isolated from paddy field soil. M. conradii was grown in monoculture and syntrophically with Pelotomaculum thermopropionicum (a propionate syntroph) or Syntrophothermus lipocalidus (a butyrate syntroph). Comparison of the relative transcript abundances showed that three hydrogenase-encoding genes and all methanogenesis-related genes tested were upregulated in cocultures relative to monoculture. The genes encoding formylmethanofuran dehydrogenase (Fwd), heterodisulfide reductase (Hdr), and the membrane-bound energy-converting hydrogenase (Ech) were the most upregulated among the evaluated genes. The expression of the formate dehydrogenase (Fdh)-encoding gene also was significantly upregulated. In contrast, an acetate assimilation gene was downregulated in cocultures. The genes coding for Fwd, Hdr, and the D subunit of F420-nonreducing hydrogenase (Mvh) form a large predicted transcription unit; therefore, the Mvh/Hdr/Fwd complex, capable of mediating the electron bifurcation and connecting the first and last steps of methanogenesis, was predicted to be formed in M. conradii. We propose that Methanocella methanogens cope with low H2 and syntrophic growth by (i) stabilizing the Mvh/Hdr/Fwd complex and (ii) activating formatedependent methanogenesis.

Considerations in the use of fluorescence in situ hybridization (FISH) and confocal laser scanning microscopy to characterize rumen methanogens and define their spatial distributions

In this study, methanogen-specific coenzyme F420autofluorescence and confocal laser scanning microscopy were used to identify rumen methanogens and define their spatial distribution in free-living, biofilm-, and protozoa-associated microenvironments. Fluorescence in situ hybridization (FISH) with temperature-controlled hybridization was used in an attempt to describe methanogen diversity. A heat pretreatment (65 °C, 1 h) was found to be a noninvasive method to increase probe access to methanogen RNA targets. Despite efforts to optimize FISH, 16S rRNA methanogen-specific probes, including Arch915, bound to some cells that lacked F420, possibly identifying uncharacterized Methanomassiliicoccales or reflecting nonspecific binding to other members of the rumen bacterial community. A probe targeting RNA from the methanogenesis-specific methyl coenzyme M reductase (mcr) gene was shown to detect cultured Methanosarcina cells with signal intensities comparable to those of 16S rRNA probes. However, the probe failed to hybridize with the majority of F420-emitting rumen methanogens, possibly because of differences in cell wall permeability among methanogen species. Methanogens were shown to integrate into microbial biofilms and to exist as ecto- and endosymbionts with rumen protozoa. Characterizing rumen methanogens and defining their spatial distribution may provide insight into mitigation strategies for ruminal methanogenesis.

Microbial trophic interactions and mcrA gene expression in monitoring of anaerobic digesters

Anaerobic digestion (AD) is a biological process where different trophic groups of microorganisms break down biodegradable organic materials in the absence of oxygen. A wide range of AD technologies is being used to convert livestock manure, municipal and industrial wastewaters, and solid organic wastes into biogas. AD gains importance not only because of its relevance in waste treatment but also because of the recovery of carbon in the form of methane, which is a renewable energy and is used to generate electricity and heat. Despite the advances on the engineering and design of new bioreactors for AD, the microbiology component always poses challenges. Microbiology of AD processes is complicated as the efficiency of the process depends on the interactions of various trophic groups involved. Due to the complex interdependence of microbial activities for the functionality of the anaerobic bioreactors, the genetic expression of mcrA, which encodes a key enzyme in methane formation, is proposed as a parameter to monitor the process performance in real time. This review evaluates the current knowledge on microbial groups, their interactions, and their relationship to the performance of anaerobic biodigesters with a focus on using mcrA gene expression as a tool to monitor the process.

Application of Anaerobic Digestion Model No. 1 to describe the syntrophic acetate oxidation of poultry litter in thermophilic anaerobic digestion

A molecular analysis found that poultry litter anaerobic digestion was dominated by hydrogenotrophic methanogens which suggests that bacterial acetate oxidation is the primary pathway in the thermophilic digestion of poultry litter. IWA Anaerobic Digestion Model No. 1 (ADM1) was modified to include the bacterial acetate oxidation process in the thermophilic anaerobic digestion (TAD). Two methods for ADM1 parameter estimation were applied: manual calibration with non-linear least squares (MC-NLLS) and an automatic calibration using differential evolution algorithms (DEA). In terms of kinetic parameters for acetate oxidizing bacteria, estimation by MC-NLLS and DEA were, respectively, km 1.12 and 3.25 ± 0.56 kg COD kg COD(-1)d(-1), KS 0.20 and 0.29 ± 0.018 kg COD m(-3) and Yac-st 0.14 and 0.10 ± 0.016 kg COD kg COD(-1). Experimental and predicted volatile fatty acids and biogas composition were in good agreement. Values of BIAS, MSE or INDEX demonstrate that both methods (MC-NLLS and DEA) increased ADM1 accuracy.

Contribution of transcriptomics to systems-level understanding of methanogenic Archaea

Methane-producing Archaea are of interest due to their contribution to atmospheric change and for their roles in technological applications including waste treatment and biofuel production. Although restricted to anaerobic environments, methanogens are found in a wide variety of habitats, where they commonly live in syntrophic relationships with bacterial partners. Owing to tight thermodynamic constraints of methanogenesis alone or in syntrophic metabolism, methanogens must carefully regulate their catabolic pathways including the regulation of RNA transcripts. The transcriptome is a dynamic and important control point in microbial systems. This paper assesses the impact of mRNA (transcriptome) studies on the understanding of methanogenesis with special consideration given to how methanogenesis is regulated to cope with nutrient limitation, environmental variability, and interactions with syntrophic partners. In comparison with traditional microarray-based transcriptome analyses, next-generation high-throughput RNA sequencing is greatly advantageous in assessing transcription start sites, the extent of 5′ untranslated regions, operonic structure, and the presence of small RNAs. We are still in the early stages of understanding RNA regulation but it is already clear that determinants beyond transcript abundance are highly relevant to the lifestyles of methanogens, requiring further study.

Functional responses of methanogenic archaea to syntrophic growth

Methanococcus maripaludis grown syntrophically with Desulfovibrio vulgaris was compared with M. maripaludis monocultures grown under hydrogen limitation using transcriptional, proteomic and metabolite analyses. These measurements indicate a decrease in transcript abundance for energy-consuming biosynthetic functions in syntrophically grown M. maripaludis, with an increase in transcript abundance for genes involved in the energy-generating central pathway for methanogenesis. Compared with growth in monoculture under hydrogen limitation, the response of paralogous genes, such as those coding for hydrogenases, often diverged, with transcripts of one variant increasing in relative abundance, whereas the other was little changed or significantly decreased in abundance. A common theme was an apparent increase in transcripts for functions using H₂ directly as reductant, versus those using the reduced deazaflavin (coenzyme F₄₂₀). The greater importance of direct reduction by H₂ was supported by improved syntrophic growth of a deletion mutant in an F₄₂₀-dependent dehydrogenase of M. maripaludis. These data suggest that paralogous genes enable the methanogen to adapt to changing substrate availability, sustaining it under environmental conditions that are often near the thermodynamic threshold for growth. Additionally, the discovery of interspecies alanine transfer adds another metabolic dimension to this environmentally relevant mutualism.

Rumen microbial (meta)genomics and its application to ruminant production

Meat and milk produced by ruminants are important agricultural products and are major sources of protein for humans. Ruminant production is of considerable economic value and underpins food security in many regions of the world. However, the sector faces major challenges because of diminishing natural resources and ensuing increases in production costs, and also because of the increased awareness of the environmental impact of farming ruminants. The digestion of feed and the production of enteric methane are key functions that could be manipulated by having a thorough understanding of the rumen microbiome. Advances in DNA sequencing technologies and bioinformatics are transforming our understanding of complex microbial ecosystems, including the gastrointestinal tract of mammals. The application of these techniques to the rumen ecosystem has allowed the study of the microbial diversity under different dietary and production conditions. Furthermore, the sequencing of genomes from several cultured rumen bacterial and archaeal species is providing detailed information about their physiology. More recently, metagenomics, mainly aimed at understanding the enzymatic machinery involved in the degradation of plant structural polysaccharides, is starting to produce new insights by allowing access to the total community and sidestepping the limitations imposed by cultivation. These advances highlight the promise of these approaches for characterising the rumen microbial community structure and linking this with the functions of the rumen microbiota. Initial results using high-throughput culture-independent technologies have also shown that the rumen microbiome is far more complex and diverse than the human caecum. Therefore, cataloguing its genes will require a considerable sequencing and bioinformatic effort. Nevertheless, the construction of a rumen microbial gene catalogue through metagenomics and genomic sequencing of key populations is an attainable goal. A rumen microbial gene catalogue is necessary to understand the function of the microbiome and its interaction with the host animal and feeds, and it will provide a basis for integrative microbiome-host models and inform strategies promoting less-polluting, more robust and efficient ruminants.

Comparative proteomic analysis of Methanothermobacter themautotrophicus ΔH in pure culture and in co-culture with a butyrate-oxidizing bacterium

To understand the physiological basis of methanogenic archaea living on interspecies H₂ transfer, the protein expression of a hydrogenotrophic methanogen, Methanothermobacter thermautotrophicus strain ΔH, was investigated in both pure culture and syntrophic coculture with an anaerobic butyrate oxidizer Syntrophothermus lipocalidus strain TGB-C1 as an H₂ supplier. Comparative proteomic analysis showed that global protein expression of methanogen cells in the model coculture was substantially different from that of pure cultured cells. In brief, in syntrophic coculture, although methanogenesis-driven energy generation appeared to be maintained by shifting the pathway to the alternative methyl coenzyme M reductase isozyme I and cofactor F₄₂₀-dependent process, the machinery proteins involved in carbon fixation, amino acid synthesis, and RNA/DNA metabolisms tended to be down-regulated, indicating restrained cell growth rather than vigorous proliferation. In addition, our proteome analysis revealed that α subunits of proteasome were differentially acetylated between the two culture conditions. Since the relevant modification has been suspected to regulate proteolytic activity of the proteasome, the global protein turnover rate could be controlled under syntrophic growth conditions. To our knowledge, the present study is the first report on N-acetylation of proteasome subunits in methanogenic archaea. These results clearly indicated that physiological adaptation of hydrogenotrophic methanogens to syntrophic growth is more complicated than that of hitherto proposed.

Effects of various fatty acid amendments on a microbial digester community in batch culture

Since biogas production is becoming increasingly important the understanding of anaerobic digestion processes is fundamental. However, large-scale digesters often lack online sensor equipment to monitor key parameters. Furthermore the possibility to selectively change fermenting parameter settings in order to investigate methane output or microbial changes is limited. In the present study we examined the possibility to investigate the microbial community of a large scale (750,000 L) digester within a laboratory small-scale approach. We studied the short-term response of the downscaled communities on various fatty acids and its effects on gas production and compared it with data from the original digester sludge. Even high loads of formic acid led to distinct methane formation, whereas high concentrations of other acids (acetic, butyric, propionic acid) caused a marked inhibition of methanogenesis coupled with an increase in hydrogen concentration. Molecular microbial techniques (DGGE/quantitative real-time-PCR) were used to monitor the microbial community changes which were related to data from GC and HPLC analysis. DGGE band patterns showed that the same microorganisms which were already dominant in the original digester re-established again in the lab-scale experiment. Very few microorganisms dominated the whole fermenting process and species diversity was not easily influenced by moderate varying fatty acid amendments--Methanoculleus thermophilus being the most abundant species throughout the variants. MCR-copy number determined via quantitative real-time-PCR--turned out to be a reliable parameter for quantification of methanogens, even in a very complex matrix like fermenter sludge. Generally the downscaled batch approach was shown to be appropriate to investigate microbial communities from large-scale digesters.

Correlation of methane production and functional gene transcriptional activity in a peat soil

The transcription dynamics of subunit A of the key gene in methanogenesis (methyl coenzyme M reductase; mcrA) was studied to evaluate the relationship between process rate (methanogenesis) and gene transcription dynamics in a peat soil ecosystem. Soil methanogen process rates were determined during incubation of peat slurries at temperatures from 4 to 37 degrees C, and real-time quantitative PCR was applied to quantify the abundances of mcrA genes and transcripts; corresponding transcriptional dynamics were calculated from mcrA transcript/gene ratios. Internal standards suggested unbiased recovery of mRNA abundances in comparison to DNA levels. In comparison to those in pure-culture studies, mcrA transcript/gene ratios indicated underestimation by 1 order of magnitude, possibly due to high proportions of inactive or dead methanogens. Methane production rates were temperature dependent, with maxima at 25 degrees C, but changes in abundance and transcription of the mcrA gene showed no correlation with temperature. However, mcrA transcript/gene ratios correlated weakly (regression coefficient = 0.76) with rates of methanogenesis. Methanogen process rates increased over 3 orders of magnitude, while the corresponding maximum transcript/gene ratio increase was only 18-fold. mcrA transcript dynamics suggested steady-state expression in peat soil after incubation for 24 and 48 h, similar to that in stationary-phase cultures. mcrA transcript/gene ratios are therefore potential in situ indicators of methanogen process rate changes in complex soil systems.

Methods for unveiling cryptic microbial partnerships in nature

Syntrophy and mutualism play a central role in carbon and nutrient cycling by microorganisms. Yet our ability to recognize these partnerships in nature or to effectively study their behavior in culture has been hindered by the inherent interdependence of syntrophic associations, their dynamic behavior, and their frequent existence at thermodynamic limits. Now solutions to these challenges are emerging in new methodologies. These include: comparative metagenomics and transcriptomics; discovery-based methods such as Magneto-FISH; and metabolic substrate tracking using stable isotopes coupled either with density gradient separation (SIP) or with FISH-SIMS. These novel approaches are redefining the way we study microbial mutualism and are making intimate microbial associations accessible to both identification and characterization in their native habitats.

Multiple response optimization analysis for pretreatments of Tequila's stillages for VFAs and hydrogen production

The objective of this work was study the effect of three pretreatments (alkalinization, thermical treatment, and sonication) on Tequila's stillages hydrolysis process in acidogenesis stage, through the following response variables: soluble chemical oxygen demand (CODs), total sugar and volatile fatty acids profile and the hydrogen production at the time. The stillages were subject to these pretreatments (according to a 2(3) factorial design); afterward they were transferred to a batch reactor at 35 degrees C and inoculated with an anaerobic digestor sludge. Multiple response optimization (MRO) analysis was done to find the global optimum for the response variables described above. This optimum is able to maximize simultaneously all these variables. It was found adequate to be useful hydrolyzing the organic matter present in Tequila's stillages. Mathematical models were fitted to observe the estimated effects of pretreatments on each response variable, then the MRO was applied.

Specific DNA binding of a potential transcriptional regulator, inosine 5'-monophosphate dehydrogenase-related protein VII, to the promoter region of a methyl coenzyme m reductase I-encoding operon retrieved from Methanothermobacter thermautotrophicus strain DeltaH

Two methyl coenzyme M reductases (MCRs) encoded by the mcr and mrt operons of the hydrogenotrophic methanogen Methanothermobacter thermautotrophicus DeltaH are expressed in response to H(2) availability. In the present study, cis elements and trans-acting factors responsible for the gene expression of MCRs were investigated by using electrophoretic mobility shift assay (EMSA) and affinity particle purification. A survey of their operator regions by EMSA with protein extracts from mrt-expressing cultures restricted them to 46- and 41-bp-long mcr and mrt upstream regions, respectively. Affinity particle purification of DNA-binding proteins conjugated with putative operator regions resulted in the retrieval of a protein attributed to IMP dehydrogenase-related protein VII (IMPDH VII). IMPDH VII is predicted to have a winged helix-turn-helix DNA-binding motif and two cystathionine beta-synthase domains, and it has been suspected to be an energy-sensing module. EMSA with oligonucleotide probes with unusual sequences showed that the binding site of IMPDH VII mostly overlaps the factor B-responsible element-TATA box of the mcr operon. The results presented here suggest that IMPDH VII encoded by MTH126 is a plausible candidate for the transcriptional regulator of the mcr operon in this methanogen.

Evaluation of molecular methods used for establishing the interactions and functions of microorganisms in anaerobic bioreactors

Molecular techniques have unveiled the complexity of the microbial consortium in anaerobic bioreactors and revealed the presence of several uncultivated species. This paper presents a review of the panoply of classical and recent molecular approaches and multivariate analyses that have been, or might be used to establish the interactions and functions of these anaerobic microorganisms. Most of the molecular approaches used so far are based on the analysis of small subunit ribosomal RNA but recent studies also use quantification of functional gene expressions. There are now several studies that have developed quantitative real-time PCR assays to investigate methanogens. With a view to improving the stability and performance of bioreactors, monitoring with molecular methods is also discussed. Advances in metagenomics and proteomics will lead to the development of promising lab-on chip technologies for cost-effective monitoring.

Coupling of Methanothermobacter thermautotrophicus methane formation and growth in fed-batch and continuous cultures under different H2 gassing regimens

In nature, H2- and CO2-utilizing methanogenic archaea have to couple the processes of methanogenesis and autotrophic growth under highly variable conditions with respect to the supply and concentration of their energy source, hydrogen. To study the hydrogen-dependent coupling between methanogenesis and growth, Methanothermobacter thermautotrophicus was cultured in a fed-batch fermentor and in a chemostat under different 80% H(2)-20% CO2 gassing regimens while we continuously monitored the dissolved hydrogen partial pressures (pH2). In the fed-batch system, in which the conditions continuously changed the uptake rates by the growing biomass, the organism displayed a complex and yet defined growth behavior, comprising the consecutive lag, exponential, and linear growth phases. It was found that the in situ hydrogen concentration affected the coupling between methanogenesis and growth in at least two respects. (i) The microorganism could adopt two distinct theoretical maximal growth yields (YCH4 max), notably approximately 3 and 7 g (dry weight) of methane formed mol-1, for growth under low (pH2 < 12 kPa)- and high-hydrogen conditions, respectively. The distinct values can be understood from a theoretical analysis of the process of methanogenesis presented in the supplemental material associated with this study. (ii) The in situ hydrogen concentration affected the "specific maintenance" requirements or, more likely, the degree of proton leakage and proton slippage processes. At low pH2 values, the "specific maintenance" diminished and the specific growth yields approached YCH4 max, indicating that growth and methanogenesis became fully coupled.

Expression of acetylcholine (ACh) and ACh-synthesizing activity in Archaea

Acetylcholine (ACh) is known generally as the neurotransmitter in the mammalian central and peripheral cholinergic nervous systems. However, ACh is also widely expressed in non-neuronal animal tissues and in plants, fungi and bacteria, where it is likely involved in the transport of water, electrolytes and nutrients, and in modulating various other cell functions. We have investigated the expression of ACh and ACh-synthesizing activity in various strains of Archaea, which are situated between Bacteria and Eucarya in the universal phylogenetic tree. Using a sensitive and specific radioimmunoassay, differing levels of ACh were detected in the Hyperthermophiles Thermococcus kodakaraensis KOD1, Sulfolobus tokodaii strain 7 and Pyrobaculum calidifontis VA1; the Methanogens Methanothermobacter thermautotrophicus deltaH and Methanosarcina barkeri; and the Halophiles Halobacterium sp. NRC-1 and Haloferax volcanii. T. kodakaraensis KOD1 expressed the highest levels of ACh among the Archaea tested; moreover, the substance expressed was verified to be ACh using high-performance liquid chromatography with electrochemical detection. Varying degrees of ACh-synthesizing activity were also identified in all of the strains, and the activity of bromoACh-sensitive choline acetyltransferase, an enzyme responsible for ACh synthesis in the nervous system, was detected in T. kodakaraensis KOD1. Our findings demonstrate that ACh and ACh-synthesizing activity are both expressed in evolutionally old Archaea. In the context of the recent discovery of non-neuronal ACh in bacteria, fungi, plants and animals, these findings support the notion that ACh has been expressed in organisms from the origin of life on the earth, functioning as a local mediator as well as a neurotransmitter.

Detecting active methanogenic populations on rice roots using stable isotope probing

Methane is formed on rice roots mainly by CO2 reduction. The present study aimed to identify the active methanogenic populations responsible for this process. Soil-free rice roots were incubated anaerobically under an atmosphere of H2/(13CO2) or N2/(13CO2) with phosphate or carbonate (marble) as buffer medium. Nucleic acids were extracted and fractionated by caesium trifluoroacetate equilibrium density gradient centrifugation after 16-day incubation. Community analyses were performed for gradient fractions using terminal restriction fragment polymorphism analysis (T-RFLP) and sequencing of the 16S rRNA genes. In addition, rRNA was extracted and analysed at different time points to trace the community change during the 16-day incubation. The Methanosarcinaceae and the yet-uncultured archaeal lineage Rice Cluster-I (RC-I) were predominant in the root incubations when carbonate buffer and N2 headspace were used. The analysis of [13C]DNA showed that the relative 16S rRNA gene abundance of RC-I increased whereas that of the Methanosarcinaceae decreased with increasing DNA buoyant density, indicating that members of RC-I were more active than the Methanosarcinaceae. However, an unexpected finding was that RC-I was suppressed in the presence of high H2 concentrations (80%, v/v), which during the early incubation period caused a lower CH4 production compared with that with N2 in the headspace. Eventually, however, CH4 production increased, probably because of the activity of Methanosarcinaceae, which became prevalent. Phosphate buffer appeared to inhibit the activity of the Methanosarcinaceae, resulting in lower CH4 production as compared with carbonate buffer. Under these conditions, Methanobacteriaceae were the prevalent methanogens. Our study suggests that the active methanogenic populations on rice roots change in correspondence to the presence of H2 (80%, v/v) and the type of buffer used in the system.

Global analysis of a "simple" proteome: Methanococcus jannaschii

The completed genome of Methanococcus jannaschii, including the main chromosome and two extra-chromosomal elements, predicts a proteome comprised of 1783 proteins. How many of those proteins are expressed at any given time and the relative abundance of the expressed proteins, however, cannot be predicted solely from the genome sequence. Two-dimensional gel electrophoresis coupled with peptide mass spectrometry is being used to identify the proteins expressed by M. jannaschii cells grown under different conditions as part of an effort to correlate protein expression with regulatory mechanisms. Here we describe the identification of 170 of the most abundant proteins found in total lysates of M. jannaschii grown under optimal fermentation conditions. To optimize the number of proteins detected, two different protein specific stains (Coomassie Blue R250 or silver nitrate) and two different first dimension separation methods (isoelectric focusing or nonequilibrium pH gradient electrophoresis) were used. Thirty-two percent of the proteins identified are annotated as hypothetical (21% conserved hypothetical and 11% hypothetical), 21% are enzymes involved in energy metabolism, 12% are proteins required for protein synthesis, and the remainder include proteins necessary for intermediary metabolism, cell division, and cell structure. Evidence of post-translational modification of numerous M. jannaschii proteins has been found, as well as indications of incomplete dissociation of protein-protein complexes. These results demonstrate the complexity of proteome analysis even when dealing with a relatively simple genome.