Contribution of transcriptomics to systems-level understanding of methanogenic Archaea

Comparative Analysis of rRNA Removal Methods for RNA-Seq Differential Expression in Halophilic Archaea

Despite intense recent research interest in archaea, the scientific community has experienced a bottleneck in the study of genome-scale gene expression experiments by RNA-seq due to the lack of commercial and specifically designed rRNA depletion kits. The high rRNA:mRNA ratio (80–90%: ~10%) in prokaryotes hampers global transcriptomic analysis. Insufficient ribodepletion results in low sequence coverage of mRNA, and therefore, requires a substantially higher number of replicate samples and/or sequencing reads to achieve statistically reliable conclusions regarding the significance of differential gene expression between case and control samples. Here, we show that after the discontinuation of the previous version of RiboZero (Illumina, San Diego, CA, USA) that was useful in partially or completely depleting rRNA from archaea, archaeal transcriptomics studies have experienced a slowdown. To overcome this limitation, here, we analyze the efficiency for four different hybridization-based kits from three different commercial suppliers, each with two sets of sequence-specific probes to remove rRNA from four different species of halophilic archaea. We conclude that the key for transcriptomic success with the currently available tools is the probe-specificity for the rRNA sequence hybridization. With this paper, we provide insights into the archaeal community for selecting certain reagents and strategies over others depending on the archaeal species of interest. These methods yield improved RNA-seq sensitivity and enhanced detection of low abundance transcripts.

Proteomic Sample Preparation and Data Analysis in Line with the Archaeal Proteome Project

Despite the ecological, evolutionary and economical significance of archaea, key aspects of their cell biology, metabolic pathways, and adaptations to a wide spectrum of environmental conditions, remain to be elucidated. Proteomics allows for the system-wide analysis of proteins, their changes in abundance between different conditions, as well as their post-translational modifications, providing detailed insights into the function of proteins and archaeal cell biology. In this chapter, we describe a sample preparation and mass spectrometric analysis workflow that has been designed for Haloferax volcanii but can be applied to a broad range of archaeal species. Furthermore, proteomics experiments provide a wealth of data that is invaluable to various disciplines. Therefore, we previously initiated the Archaeal Proteome Project (ArcPP), a community project that combines the analysis of multiple datasets with expert knowledge in various fields of archaeal research. The corresponding bioinformatic analysis, allowing for the integration of new proteomics data into the ArcPP, as well as the interactive exploration of ArcPP results is also presented here. In combination, these protocols facilitate an optimized, detailed and collaborative approach to archaeal proteomics.

Review: Fifty years of research on rumen methanogenesis: lessons learned and future challenges for mitigation

Meat and milk from ruminants provide an important source of protein and other nutrients for human consumption. Although ruminants have a unique advantage of being able to consume forages and graze lands not suitable for arable cropping, 2% to 12% of the gross energy consumed is converted to enteric CH4 during ruminal digestion, which contributes approximately 6% of global anthropogenic greenhouse gas emissions. Thus, ruminant producers need to find cost-effective ways to reduce emissions while meeting consumer demand for food. This paper provides a critical review of the substantial amount of ruminant CH4-related research published in past decades, highlighting hydrogen flow in the rumen, the microbiome associated with methanogenesis, current and future prospects for CH4 mitigation and insights into future challenges for science, governments, farmers and associated industries. Methane emission intensity, measured as emissions per unit of meat and milk, has continuously declined over the past decades due to improvements in production efficiency and animal performance, and this trend is expected to continue. However, continued decline in emission intensity will likely be insufficient to offset the rising emissions from increasing demand for animal protein. Thus, decreases in both emission intensity (g CH4/animal product) and absolute emissions (g CH4/day) are needed if the ruminant industries continue to grow. Providing producers with cost-effective options for decreasing CH4 emissions is therefore imperative, yet few cost-effective approaches are currently available. Future abatement may be achieved through animal genetics, vaccine development, early life programming, diet formulation, use of alternative hydrogen sinks, chemical inhibitors and fermentation modifiers. Individually, these strategies are expected to have moderate effects (<20% decrease), with the exception of the experimental inhibitor 3-nitrooxypropanol for which decreases in CH4 have consistently been greater (20% to 40% decrease). Therefore, it will be necessary to combine strategies to attain the sizable reduction in CH4 needed, but further research is required to determine whether combining anti-methanogenic strategies will have consistent additive effects. It is also not clear whether a decrease in CH4 production leads to consistent improved animal performance, information that will be necessary for adoption by producers. Major constraints for decreasing global enteric CH4 emissions from ruminants are continued expansion of the industry, the cost of mitigation, the difficulty of applying mitigation strategies to grazing ruminants, the inconsistent effects on animal performance and the paucity of information on animal health, reproduction, product quality, cost-benefit, safety and consumer acceptance.

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.

Model Microbial Consortia as Tools for Understanding Complex Microbial Communities

A major biological challenge in the postgenomic era has been untangling the composition and functions of microbes that inhabit complex communities or microbiomes. Multi-omics and modern bioinformatics have provided the tools to assay molecules across different cellular and community scales; however, mechanistic knowledge over microbial interactions often remains elusive. This is due to the immense diversity and the essentially undiminished volume of not-yet-cultured microbes. Simplified model communities hold some promise in enabling researchers to manage complexity so that they can mechanistically understand the emergent properties of microbial community interactions. In this review, we surveyed several approaches that have effectively used tractable model consortia to elucidate the complex behavior of microbial communities. We go further to provide some perspectives on the limitations and new opportunities with these approaches and highlight where these efforts are likely to lead as advances are made in molecular ecology and systems biology.

Unique Archaeal Small RNAs

Advances in genome-wide sequence technologies allow for detailed insights into the complexity of RNA landscapes of organisms from all three domains of life. Recent analyses of archaeal transcriptomes identified interaction and regulation networks of noncoding RNAs in this understudied domain. Here, we review current knowledge of small, noncoding RNAs with important functions for the archaeal lifestyle, which often requires adaptation to extreme environments. One focus is RNA metabolism at elevated temperatures in hyperthermophilic archaea, which reveals elevated amounts of RNA-guided RNA modification and virus defense strategies. Genome rearrangement events result in unique fragmentation patterns of noncoding RNA genes that require elaborate maturation pathways to yield functional transcripts. RNA-binding proteins, e.g., L7Ae and LSm, are important for many posttranscriptional control functions of RNA molecules in archaeal cells. We also discuss recent insights into the regulatory potential of their noncoding RNA partners.

Methanol metabolism and archaeal community changes in a bioelectrochemical anaerobic digestion sequencing batch reactor with copper-coated graphite cathode

In this study, the metabolism of methanol and changes in an archaeal community were examined in a bioelectrochemical anaerobic digestion sequencing batch reactor with a copper-coated graphite cathode (BEAD-SBR_Cu). Copper-coated graphite cathode produced methanol from food waste. The BEAD-SBR_Cu showed higher methanol removal and methane production than those of the anaerobic digestion (AD)-SBR. The methane production and pH of the BEAD-SBR_Cu were stable even under a high organic loading rate (OLR). The hydrogenotrophic methanogens increased from 32.2 to 60.0%, and the hydrogen-dependent methylotrophic methanogens increased from 19.5 to 37.7% in the bulk of BEAD-SBR_Cu at high OLR. Where methanol was directly injected as a single substrate into the BEAD-SBR_Cu, the main metabolism of methane production was hydrogenotrophic methanogenesis using carbon dioxide and hydrogen released by the oxidation of methanol on the anode through bioelectrochemical reactions.

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.

Horizontal gene transfer drives the evolution of Rh50 permeases in prokaryotes

Background

Rh50 proteins belong to the family of ammonia permeases together with their Amt/MEP homologs. Ammonia permeases increase the permeability of NH₃/NH₄⁺ across cell membranes and are believed to be involved in excretion of toxic ammonia and in the maintenance of pH homeostasis. RH50 genes are widespread in eukaryotes but absent in land plants and fungi, and remarkably rare in prokaryotes. The evolutionary history of RH50 genes in prokaryotes is just beginning to be unveiled.

Results

Here, a molecular phylogenetic approach suggests horizontal gene transfer (HGT) as a primary force driving the evolution and spread of RH50 among prokaryotes. In addition, the taxonomic distribution of the RH50 gene among prokaryotes turned out to be very narrow; a single-copy RH50 is present in the genome of only a small proportion of Bacteria, and, first evidence to date, in only three methanogens among Euryarchaea. The coexistence of RH50 and AMT in prokaryotes seems also a rare event. Finally, phylogenetic analyses were used to reconstruct the HGT network along which prokaryotic RH50 evolution has taken place.

Conclusions

The eukaryotic or bacterial “origin” of the RH50 gene remains unsolved. The RH50 prokaryotic HGT network suggests a preferential directionality of transfer from aerobic to anaerobic organisms. The observed HGT events between archaeal methanogens, anaerobic and aerobic ammonia-oxidizing bacteria suggest that syntrophic relationships play a major role in the structuring of the network, and point to oxygen minimum zones as an ecological niche that might be of crucial importance for HGT-driven evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0850-6) contains supplementary material, which is available to authorized users.

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.

Cattle Manure Enhances Methanogens Diversity and Methane Emissions Compared to Swine Manure under Rice Paddy

Livestock manures are broadly used in agriculture to improve soil quality. However, manure application can increase the availability of organic carbon, thereby facilitating methane (CH4) production. Cattle and swine manures are expected to have different CH4 emission characteristics in rice paddy soil due to the inherent differences in composition as a result of contrasting diets and digestive physiology between the two livestock types. To compare the effect of ruminant and non-ruminant animal manure applications on CH4 emissions and methanogenic archaeal diversity during rice cultivation (June to September, 2009), fresh cattle and swine manures were applied into experimental pots at 0, 20 and 40 Mg fresh weight (FW) ha-1 in a greenhouse. Applications of manures significantly enhanced total CH4 emissions as compared to chemical fertilization, with cattle manure leading to higher emissions than swine manure. Total organic C contents in cattle (466 g kg-1) and swine (460 g kg-1) manures were of comparable results. Soil organic C (SOC) contents were also similar between the two manure treatments, but dissolved organic C (DOC) was significantly higher in cattle than swine manure. The mcrA gene copy numbers were significantly higher in cattle than swine manure. Diverse groups of methanogens which belong to Methanomicrobiaceae were detected only in cattle-manured but not in swine-manured soil. Methanogens were transferred from cattle manure to rice paddy soils through fresh excrement. In conclusion, cattle manure application can significantly increase CH4 emissions in rice paddy soil during cultivation, and its pretreatment to suppress methanogenic activity without decreasing rice productivity should be considered.

Methane yield phenotypes linked to differential gene expression in the sheep rumen microbiome

Ruminant livestock represent the single largest anthropogenic source of the potent greenhouse gas methane, which is generated by methanogenic archaea residing in ruminant digestive tracts. While differences between individual animals of the same breed in the amount of methane produced have been observed, the basis for this variation remains to be elucidated. To explore the mechanistic basis of this methane production, we measured methane yields from 22 sheep, which revealed that methane yields are a reproducible, quantitative trait. Deep metagenomic and metatranscriptomic sequencing demonstrated a similar abundance of methanogens and methanogenesis pathway genes in high and low methane emitters. However, transcription of methanogenesis pathway genes was substantially increased in sheep with high methane yields. These results identify a discrete set of rumen methanogens whose methanogenesis pathway transcription profiles correlate with methane yields and provide new targets for CH₄ mitigation at the levels of microbiota composition and transcriptional regulation.

Towards a computational model of a methane producing archaeum

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