Application of meta-omics techniques to understand greenhouse gas emissions originating from ruminal metabolism
Genetics Selection Evolution volume 49, Article number: 9 (2017)
Methane emissions from ruminal fermentation contribute significantly to total anthropological greenhouse gas (GHG) emissions. New meta-omics technologies are beginning to revolutionise our understanding of the rumen microbial community structure, metabolic potential and metabolic activity. Here we explore these developments in relation to GHG emissions. Microbial rumen community analyses based on small subunit ribosomal RNA sequence analysis are not yet predictive of methane emissions from individual animals or treatments. Few metagenomics studies have been directly related to GHG emissions. In these studies, the main genes that differed in abundance between high and low methane emitters included archaeal genes involved in methanogenesis, with others that were not apparently related to methane metabolism. Unlike the taxonomic analysis up to now, the gene sets from metagenomes may have predictive value. Furthermore, metagenomic analysis predicts metabolic function better than only a taxonomic description, because different taxa share genes with the same function. Metatranscriptomics, the study of mRNA transcript abundance, should help to understand the dynamic of microbial activity rather than the gene abundance; to date, only one study has related the expression levels of methanogenic genes to methane emissions, where gene abundance failed to do so. Metaproteomics describes the proteins present in the ecosystem, and is therefore arguably a better indication of microbial metabolism. Both two-dimensional polyacrylamide gel electrophoresis and shotgun peptide sequencing methods have been used for ruminal analysis. In our unpublished studies, both methods showed an abundance of archaeal methanogenic enzymes, but neither was able to discriminate high and low emitters. Metabolomics can take several forms that appear to have predictive value for methane emissions; ruminal metabolites, milk fatty acid profiles, faecal long-chain alcohols and urinary metabolites have all shown promising results. Rumen microbial amino acid metabolism lies at the root of excessive nitrogen emissions from ruminants, yet only indirect inferences for nitrogen emissions can be drawn from meta-omics studies published so far. Annotation of meta-omics data depends on databases that are generally weak in rumen microbial entries. The Hungate 1000 project and Global Rumen Census initiatives are therefore essential to improve the interpretation of sequence/metabolic information.
Many terms employ the ‘meta-’ prefix and ‘-omics’ or ‘-ome’ suffixes. Arguably, among all these, the four most relevant to the rumen microbial community and ruminal metabolism are metagenomics, metatranscriptomics, metaproteomics and metabolomics. All four take advantage of technologies that have only recently become generally available. Metagenomics, the study of all the genes present in the ecosystem, and metatranscriptomics, the study of transcribed genes, employ high-throughput DNA-sequencing, which has become incredibly fast and inexpensive over the last decade. Metaproteomics, which catalogues the total protein complement of the community—the translated genes—now uses high-resolution mass spectrometry to identify peptides derived from these proteins by shotgun hydrolysis. Metabolomics uses a variety of spectroscopic and mass spectrometric methods and separation techniques to quantify the metabolites that are present. Each of the meta-omics technologies tells us something different about the microbial community and its activities. Here we assess how they may help to provide effective strategies to mitigate the pressing environmental problems associated with greenhouse gas (GHG) emissions from ruminant livestock production.
Concerns about methane and nitrogen emissions from ruminants
The 2006 publication  by the Food and Agriculture Organisation (FAO) and the Livestock, Environment and Development Initiative, ‘Livestock’s Long Shadow, Environmental Issues and Options’, marked a watershed in public and political views on livestock and the environment. The following highly emotive paragraph in the Executive Summary encapsulates its message—“Livestock’s contribution to environmental problems is on a massive scale and its potential contribution to their solution is equally large. The impact is so significant that it needs to be addressed with urgency. Major reductions in impact could be achieved at reasonable cost.” Land degradation, water shortage and biodiversity are important, and also the atmosphere and climate change. Ruminants loom large in the last concern, because they, and their excreta, produce large amounts of methane and nitrous oxide emitted to the atmosphere. The report concluded that the livestock sector is responsible for 18% of total greenhouse gas (GHG) emissions and 37% of total anthropogenic methane, which is largely responsible for the total amount. While the exact numbers have varied in the interim, and more aspects of the whole system have been factored into the models, it is clear that ruminant methane and nitrogen (N) emissions, which originate largely from rumen microbial activity, must be addressed in our efforts to limit climate change.
The rumen microbial community and methane
The rumen is home to a vast array of microbes from the three great domains of life. Their abundance per g of digesta ranges from 104 to 106 ciliate protozoa (although sometimes there are none), 103 to 105 anaerobic fungi, 1010 to 1011 anaerobic bacteria and 108 to 109 archaea. The protozoa can comprise up to half the rumen microbial biomass, the fungi about 7%, the archaea 1 to 4% and the bacteria form the remainder. In a recent publication , we reviewed the composition of the ruminal community relating to methanogenesis. Briefly, the abundance of archaea has only a weak correlation with methane emissions from individual cattle and sheep. The composition of the archaeal community appears to have a stronger effect, with animals that harbour the Methanobrevibacterium gottschalkii clade tending to be associated with greater methane emissions. Although ciliate protozoa are well known to produce H2 and harbour abundant archaea, their numbers do not have a strong relation to methane emissions. A meta-analysis of defaunation revealed methane emissions to be on average 11% lower than in faunated animals . Methane emissions are greater from ruminants that have high abundance of H2-producing bacteria, and lower when non-H2-producers, such as Succinovibrionaceae, are more numerous. Individual taxa correlate with methane emissions, but not necessarily in the manner expected. Fundamental questions regarding the physiology and metabolism of individual species, both cultivated and those not yet cultivated, need to be addressed in order to understand how methane emissions are affected by the microbiome.
Methane is a GHG that is 28 times more potent than CO2 . Around 90% of the methane produced by ruminants is derived from the rumen , where methanogenic archaea convert the H2 and CO2 produced by the protozoa, bacteria and fungi to methane . Worldwide research efforts have investigated various mitigation strategies, particularly feed additives that might inhibit H2 production, provide an alternative H sink or inhibit the archaea themselves [7–10]. Other strategies include chemogenomics and immunization [11–13]. A strategy that could be most sustainable, because of its persistence and ease of implementation, is genetic selection for low methane-emitting animals [14–16]. If it can be demonstrated that the different volumes of methane emissions from different animals can be explained by their differing ruminal microbiomes, and that the property is persistent and heritable, it should be possible to select future generations of cattle and sheep that have genetically determined lower methane emissions. Thus far, it has been demonstrated that methane emissions in sheep [14, 15, 17], dairy cows  and beef steers [19, 20] are significantly heritable. Indeed, the prediction of methane emissions via milk fatty acid composition, as described below, is heritable . It had been expected that lower methane emissions would improve the efficiency of energy retention and thereby increase feed efficiency. However, unfortunately that largely intuitive prediction does not seem to hold in practice [22, 23], thus weakening the incentive to farmers to adopt measures that would lower methane emissions. However, the reverse is undoubtedly true, i.e. that more efficient cattle will produce less methane per unit product (meat, milk), thus a focus on feed efficiency may be more fruitful, rather than simply methane or N emissions alone.
Nitrous oxide is about ten times as potent a GHG as methane . It is formed by microbial denitrification in soil and in anaerobic slurries, both of which are exacerbated by the oversupply of dietary protein to cattle. The quantity of protein flowing from the rumen is a major factor that limits the productivity of ruminant livestock production [24, 25]. The protein reaching the abomasum consists of a mixture of dietary and microbial protein and, following digestion and absorption, it provides the amino acids upon which ruminants depend for their amino acid requirements. Rumen wall tissue protein turnover also contributes to the protein drain imposed by ruminal microorganisms, because ruminal bacteria tend to invade and digest ruminal epithelial tissues [26, 27]. In order to compensate for these inefficiencies, ruminant livestock producers tend to oversupply the animals with relatively cheap protein sources such as soybean meal. The excess N is excreted in urine and faeces, which then present a disposal problem.
Nitrous oxide emissions are equivalent to methane emissions in Scotland in terms of GHG from agriculture . Nitrogenous excretion from ruminants is therefore another area that needs to be addressed. Part of the inefficiency stems from the animal itself, with inefficient amino acid metabolism, but the main inefficiency arises from the proteolytic and bacteriolytic activities of ruminal microorganisms [24, 25].
Rumen microbial metagenomics and GHG emissions
The first application of the metagenome concept to the rumen microbiota was gene mining, whereby gene libraries that were sequenced from the total DNA of ruminal digesta were screened for target activities. This approach proved successful in the discovery and characterisation of many key microbial enzymes such as glycosyl hydrolases [29–33], polyphenol oxidases , and lipases [35, 36]. During annotation of whole metagenomes in rumen studies, it was apparent that the majority of the open reading frames (ORF) encoded genes that were unknown or not yet included in reference databases. Furthermore, with the vast majority of ruminal species yet to be cultivated in vitro [37, 38], the potential of metagenome mining in the rumen is vast.
Pioneering papers to explore the wider potential of metagenomics applied to the rumen were those by Brulc et al.  and Hess et al. . Brulc et al.  were the first to report the results of deep sequencing of the ruminal metagenome. They focussed mainly on glycosyl hydrolase sequence analysis, in a comparative metagenomics exercise that was the first of its kind in the rumen. A comparison of the glycosyl hydrolase and cellulosome functional genes in digesta from three steers revealed that, in the rumen microbiome, initial colonization of fibre appears to be by organisms that possess enzymes that attack the easily available side chains of complex plant polysaccharides rather than the more recalcitrant main chains, especially cellulose. In an interesting cross-species comparison, Brulc et al.  compared their rumen data with that of the termite hindgut microbiome. Fundamental differences in the glycosyl hydrolase content appeared to be diet-dependent, with cattle consuming forages and legumes compared to the consumption of wood by termites.
Hess et al.  were also driven largely by the potential discovery of new glycosyl hydrolases that might be of value in the biofuels industry, but they demonstrated also the depth of new information that could be extracted from metagenomic deep sequencing. Only one cow was used in this experiment, yet the wealth of new discoveries was immense. At least five operational taxonomic units (OTU) were enriched on the switchgrass. None of these was identified to be a cultivated species, indicating a major opportunity to isolate the enriched species that by implication could be involved in switchgrass degradation and therefore be useful in the biofuels industry. Only 12% of the 27,755 carbohydrate-active genes that were assembled from the ruminal metagenome of switchgrass-adherent microorganisms were more than 75% identical to genes deposited in the NCBI non-redundant database, whereas 43% of the genes had less than 50% identity to any known protein. Ninety of the candidate proteins were expressed in vitro, of which 57% were enzymatically active against cellulosic substrates. It might be argued that, since glycosyl hydrolases are by far the best characterised enzymes from the ruminal ecosystem, even more novelty would be seen when mining enzymes with different functions that are important to ruminal microorganisms, such as protein or lipid metabolism. The gene mining so far accomplished has barely scratched the surface of such a complex enzymatic ecosystem.
Perhaps the most remarkable demonstration of the Hess et al.  analysis was the assembly of 15 bacterial genomes from uncultured species at completeness that ranged from 60 to 93%. The assemblies were validated by complementary methods including single-cell genome sequencing. This kind of genome assembly, by analysing the genes that are present, can help us to understand the metabolic role and ecological niche of bacteria that have yet to be cultivated.
Things are now moving rapidly in relating metagenomics to methane emissions. Denman and McSweeney  and McAllister et al.  published extensive reviews less than two years ago, to which the reader is referred. Since then, several fundamental research papers have been published using metagenomics to understand GHG emissions. The methanogenic archaeal community and its gene complement were characterized by metagenomics analysis in the buffalo rumen . Genes encoding all the key steps of methanogenesis were found. Moreover, and a potentially significant finding was the discovery of genes involved in the acetogenesis pathway, a possible alternative to methanogenesis in the rumen. However, in goats, the contribution of reductive acetogenesis in redirecting H2 away from methanogenesis was minimal, even when methanogenesis was inhibited by bromochloromethane . Instead, genes involved in propionate formation via the randomizing pathway, and numbers of corresponding bacteria among Prevotella and Selenomonas spp., increased in the presence of bromochloromethane, while the genes involved in methanogenesis decreased.
Another example study in beef demonstrated significant differences (P < 0.05) in the abundance of 21 of the most numerous (>0.1%) genes when the rumen microbial metagenomes from high and low methane-emitting beef steers were compared . Eight of the nine most significantly differing genes were associated with methane metabolism, but the others were not. Indeed, their link with methanogenesis was not obvious. The abundance of the 21 genes in total explained 88% of the variation in methane production, thus possibly forming the basis for genetic selection of animals with a low-methane genotype. The same experiments showed that sire-progeny groups differed in their methane emissions. Further analysis  demonstrated that the abundance of 49 genes explained 86% of the variation in feed efficiency. Once again, the reasons that underlie these correlations were not obvious, although it was noted that host-microbiota crosstalk gene expression (TSTA3 and FucI) were significantly associated with feed efficiency. These results suggest, as proposed by Taxis et al. , that future studies of the whole animal-gut microbiome networks hold high promise for understanding the ‘superorganism’ .
A significant recent paper on ruminal metagenomics explored feed efficiency in dairy cows . Methane emissions were also measured ex vivo, while metagenomes were studied from deep sequencing . Species diversity was lower in the more efficient animals, as was gene diversity (Fig. 1). Moreover, methane emissions were also significantly lower in the efficient animals, as was found previously in cattle . Ruminal digesta contained more propionate, butyrate and isovalerate in efficient animals. Most striking of all, metabolic pathway analysis showed that genes of the non-randomizing acrylate pathway of propionate production were much more prevalent in the efficient cattle. The acrylate pathway is found principally in the distinctive, large Gram-negative coccus, Megasphaera elsdenii, which has been identified with a stabilising effect on ruminal fermentation because of its rapid conversion of lactate to propionate and butyrate [48, 49]; M. elsdenii also produces isovalerate and valerate as end-products of its amino acid-fermenting ability [50, 51]. rRNA gene amplicon analysis showed that M. elsdenii abundance was much greater in efficient animals, corresponding to the acrylate gene abundance. In the study by Wallace et al.  in beef cattle, although not reported in the paper itself, the abundance of M. elsdenii was 13-fold higher in the low-methane steers, thus entirely consistent with the results of Shabat et al. . M. elsdenii has been trialled with some success as a probiotic for ruminants on the grounds of its pH-stabilizing properties [52, 53]. Thus, thanks to these studies, a picture is emerging whereby it can be seen that differences in the abundance of H2-producing bacteria, non-H2-producing bacteria and H2 utilisers, together with the abundance of pH-stabilizing bacteria, affect the quantity of methane that a ruminant animal produces and its feed efficiency.
The Hungate 1000 project and Global Rumen Census
Thus far, combined understanding of function and phylogenetic identity in metagenomics data has been limited by the relatively few completed rumen bacterial genomes and in turn by the number of annotated genes and protein sequences of ruminal species. This issue is being addressed by the Hungate 1000 project (www.hungate1000.org.nz). The project title refers to the pioneering work in culturing strictly anaerobic ruminal bacteria carried out by Robert E. Hungate . The aim of the project is to produce a reference set of 1000 rumen microbial genome sequences from cultivated rumen bacteria and methanogenic archaea, together with representative cultures of rumen anaerobic fungi and ciliate protozoa. The project is funded by the New Zealand Government in support of the Livestock Research Group of the Global Research Alliance on Agricultural Greenhouse Gases. The sequencing effort obtained support from the US Department of Energy Joint Genome Institute Community Sequencing Program, and the overall project is a global collaboration between members of the Rumen Microbial Genomics Network, established to accelerate knowledge development and mitigation solutions in the rumen microbial genomics research area. The reference genome information gathered will be used to facilitate genome-enabled research aimed at understanding rumen function in order to find a balance between food production and GHG emissions, and to support international efforts to develop methane mitigation and rumen adaptation technologies. Once the Hungate 1000 project is completed, genes discovered from deep metagenome sequencing will be able to be pinned with much greater certainty to known species.
The most extensive exploration of the ruminal microbiome that was recently published was the Global Rumen Census, an international effort that analysed the microbial community in 742 samples from 32 animal species from 35 countries . The results revealed a common core microbiome in all samples, and prompted the conclusion that significant new taxonomic groups were unlikely to be discovered. The authors also commented on likely functional redundancy, with different taxa performing essentially the same function using related genes, a topic that has also been reviewed recently . Indeed future understanding of rumen function, and its relationship with the host genome, is likely to be expanded most significantly by exploring and linking gene networks . An important overall conclusion of the Global Rumen Census was that diet, rather than genetics or geographical location, had the greatest influence on the ruminal microbiome.
Similar to metagenomics, metatranscriptomics was used first as a tool for gene mining by Qi et al. , again with the principal objective of identifying novel lignocellulolytic and glycosyl hydrolase genes in the muskoxen rumen, with the interesting hypothesis that new genes, particularly from the eukaryotic community, might be found. The investigation was highly successful, achieving an 8.7× higher rate of total carbohydrate active enzyme discovery than that found in previous metagenomics analyses. The metatranscriptomic approach offers the unique possibility to restrict the analysis only to transcribed genes, thus removing the often very high noise of non-transcribed portions of the genome, which are by contrast always present in metagenomic experiments.
Shi et al.  investigated methane production in a cohort of New Zealand sheep using metagenomics and metatranscriptomic techniques that aimed at understanding microbiological differences between animals that produced low and high amounts of methane. The paper illustrated the power of deep sequencing in understanding the microbial community and its activity. Four rams with a high-methane phenotype, identified from a pool of 22 animals, were compared with four rams with a low-methane phenotype and four rams with an intermediate phenotype. The difference in methane production between the high and low phenotypes was about 1.7-fold, similar to the beef cattle study  discussed above. Microbial community structures were compared by extracting rRNA gene sequence information from deep sequencing and also by qPCR of rRNA and mcrA/mcrT genes. No differences were detected in the different microbial groups. Further detailed analysis of the archaeal community found higher abundances of Methanobrevibacter gottschalkii in high producers, an observation that has been reported in other studies . Methanogenic gene abundances also did not differ between the animal groups. It was only the metatranscriptome that differed, where the abundance of mRNA sequences was compared. Three of the ten most increased transcripts in the high producers coded for enzymes in the methanogenesis pathway. The idea that the transcriptome is more responsive as a measurement of methane emissions has gained currency. This argument was challenged  because ATP production in methanogens is entirely dependent on methane formation and the growth yield, molar growth yield (g biomass/mol ATP utilised) (YATP), is proportional to ATP production. However, molar growth yield [g biomass/mol CH4 produced (Ymethane)] varies significantly according to growth conditions, with excess H2 apparently leading to uncoupling  analogous to that observed in bacteria where growth is limited by a nutrient other than a sugar as energy source . Given the extensive nature of electron-transport-linked metabolism in methanogens , ruminal archaea may well use similar mechanisms to maintain cellular metabolites during periods of stress, and their abundance may therefore not be proportional to the quantity of methane formed. Further metatranscriptomic studies, linked possibly to metabolomics analysis, might be useful in investigating this point.
The proteome differs from the previous -omes in that, while the others predict what genes are present and how they are transcribed, the proteome reflects the end-product, the proteins that are actually expressed. There has not been a concerted effort to characterise the proteomes of different pure cultures of ruminal microorganisms. Indeed, it appears that the technology and interest have jumped that particular step to study the metaproteome, i.e. the entire complement of proteins that is expressed by the ruminal microbiome. Metaproteomic analysis aims at characterising the entire protein content of an environmental sample at a given point in time . At first, it may seem improbable that such a complex community, comprising hundreds of species each with thousands of genes, would present a proteome that would be sufficiently discriminated to enable the identification of individual proteins. Nonetheless, earlier examples of metaproteomic analyses from the human gut  and soil  have shown that it is in fact technically feasible.
Two main technical approaches are available in proteomics. The first is the long established two-dimensional SDS polyacrylamide gel electrophoresis (2D SDS-PAGE) technology that was originated by O’Farrell . Separation of the total protein is accomplished by isoelectric point in the first dimension and molecular size in the second. The proteome is visualised using a stain, revealing individual spots that can be identified by mass spectroscopic analysis following trypsinisation of spots cut from the gel. Protein identification depends heavily on searches of reference databases that contain relatively few rumen microbial proteomes. A recent development in proteome technology uses state-of-the-art mass spectrometers that are capable of analysing complex mixtures of peptides derived by partial hydrolysis of total protein mixtures. Raw data are generated as a massive set of mass spectra which are converted into a long list of short peptide sequences (the metapeptidome). These are assembled into proteins by mapping to a reference database in a similar way to shotgun DNA sequencing, hence the name shotgun metaproteomics. Many believe that the shotgun method, with the much larger volume of data generated, will supplant the gel-based method. There are still a number of technical issues that need to be addressed before shotgun metaproteomics can be used for comparative analysis. The first is a reliable method to quantitate data. This has been carried out previously using spectrum counting but can also be achieved by labelling samples with stable isotopes. Moreover, there is a lack of bioinformatics analysis support and, similar to 2D SDS-PAGE, the identification of proteins relies on mapping data to amino acid sequence databases in which the great majority of ruminal species are not represented.
The rumen ecosystem shares some characteristics with microbial communities in the environment and human gut that have previously been characterised using metaproteomics, such as microbial diversity and relative abundance of microorganisms in some studies [64, 66–68] and the abundance of nutrients in others , but it provides a unique challenge in the combination of these properties. The metaproteome will provide a different insight of the function of the rumen microbial community compared to the nucleic acid meta-omes, arguably one that might prove more useful as part of the campaign to lower methane emissions and to better understand the role of key enzymes involved in feed utilisation efficiency in ruminants.
The RuminOmics project (www.ruminomics.eu) investigated SDS-PAGE methods for generating metaproteomic information from ruminal digesta . Results were variable according to the sample. In some gels, distinct spots were observed, while in others interference by humic substances that are derived from the plant materials consumed by the animal, resulted in no distinct protein spot pattern. In the gels where spots were resolved, tandem mass spectrum analysis indicated that structural proteins from protozoa were most abundant, an expected result considering the high proportion of their biomass in the rumen. A surprising discovery was the strong resolution of key enzymes associated with methanogenesis from the archaea that form a relatively small proportion of the rumen microbial community. In a comparison of the 2D PAGE metaproteomes of high- and low-methane emitting dairy cows, no significant difference was evident although this was possibly due to the lack of precision using this technique.
The first analysis of the ruminal metaproteome using shotgun peptide methodology was published in 2015 . Remarkably, taxonomic information assigned to the predicted proteins enabled a community analysis to be carried out, in which the relative abundance of different bacterial and archaeal families and eukaryote phyla were calculated. The composition of the microbial community was different from those most commonly seen in the rumen, but no comparative DNA-based analysis was presented. It would be very interesting to examine the correspondence of the microbiome deduced from the metagenome to that predicted from the corresponding metaproteome. Another feature of the analysis was the high abundance of plant-derived peptides detected by the shotgun method (Fig. 2) compared to none being detected by the 2-D method . The plant-derived peptides would be the proteolytic products of microbial digestion of plant protein in the feed. Methanogenesis-associated proteins were only mentioned in passing, but presumably the metaproteome may be as useful in predicting metabolic pathways as it was in describing pathways of starch metabolism . Once again, no differences in the microbial community based on the metaproteome between high- and low-methane emitters were evident in dairy cows in the RuminOmics project (Fig. 2).
Metabolomics provides a detailed array of information about ruminal metabolic activity, which is complementary to the DNA- and protein-based methods described above. The reader is directed to the ground-breaking paper by Saleem et al.  on this topic. Here, four metabolomic analyses will be discussed that relate specifically to methanogenesis.
Shabat et al.  analysed the ruminal metabolome in cows that varied in feed efficiency, i.e. the conversion of feed to product. Higher concentrations of short-chain fatty acids were observed in more efficient cows, accompanied by lower methane emissions. Other metabolites were not significantly different, except for putrescine, which was present at higher concentrations in efficient cows. Whether this reflects a key aspect of ruminal metabolism that affects efficiency is unclear. A large number of metabolites was detected by Zhao et al. , related to both nitrogen and protein metabolism. The metabolome was highly dependent on the dietary composition, and no attempt was made to correlate metabolites with emissions, but significant differences were seen in amino acid metabolites and in methylamines that are substrates for methylotrophic methanogenesis, suggesting possible future value in these measurements to studies of ruminant GHG emissions.
The faecal metabolome associated with methane includes the distinctive membrane lipids of the archaea, namely dialkyl glycerol diethers (DGDG) and glycerol dialkyl glycerol tetraether (GDGT). The most common forms of DGDG and GDGT are archaeol and caldarchaeol, respectively (Fig. 3). Archaeol has received the most attention, and has had its relationship with methane production analysed across a range of diets in studies on beef and dairy cattle [73–75]. These studies concluded that there is considerable between-animal variation in the relationship, although the relationship is significant when comparing the treatment means. Between-animal variation could be attributed to differences in the location and kinetics of methanogens in the ruminant digestive tract, and a lack of relationship between archaeol measurements in the rumen and the faeces . Another potential cause of the variation could be the oversight of the presence of caldarchaeol in the methanogen membrane. In comparison to archaeol, which forms a bilayer, caldarchaeol forms a monolayer and is less permeable to protons. McCartney et al.  found that the proportion of caldarchaeol in the faeces increased markedly when the animal was fed a diet high in starch, and thus perhaps protecting the methanogens from the resultant drop in ruminal pH. Furthermore, concentrations of caldarchaeol and total ether lipids were found to be more proportional to measured methane production than archaeol concentrations. In summary, archaeol is potentially a useful alternative marker for determining methanogen abundance, however, as a methane proxy, more work is needed to further investigate both archaeol and caldarchaeol.
The urinary metabolome has provided much useful information on N retention and fluxes in the animal for many years . The purines provide a useful proxy measure of microbial protein flow from the rumen [79, 80] and urea itself is of course the most important metabolite associated with the efficiency of N retention. However, a recent study in protozoa-depleted lambs revealed that the urinary metabolomes of faunated and the protozoa-depleted animals were almost completely polarised in terms of protein-derived metabolites following discriminant analysis . In spite of the complexity of the data, the clear separation of the metabolome according to the different treatments gives an indication of the value of further investigation into the urinary metabolome. Correlation with the composition of the microbiota also suggests the possibility of using the urinary metabolome to predict rumen microbial metabolism. For instance, metabolites of tryptophan were linked not only to the abundance of protozoa but also to bacterial taxa mostly distantly related to known species. Intriguingly, the urinary metabolome study also revealed a possible link to methanogenesis. Methane emissions were not measured, but a negative relationship was found between urinary trimethylamine-N-oxide and the ruminal abundance of the methylotrophic methanogenesis order Methanomassiliicoccales.
The fatty acid composition of milk, sometimes called the milk lipidome, can also be useful in predicting ruminal metabolism, including methanogenesis, in dairy cows. A meta-analysis  concluded that milk fatty acid concentrations of C10:0, C12:0, C14:0-iso, C14:0, cis-9 C14:1, C15:0, and C16:0 were positively related to methane yield per unit of milk, while C4:0, C18:0, trans-10 + 11 C18:1, cis-9 C18:1, cis-11 C18:1, and cis-9,12 C18:2 in milk fat were negatively related. Mathematical analysis enabled prediction equations to be formulated that had moderate potential for predicting methane yield per unit of feed and a slightly lower potential for predicting methane yield per unit of milk. Subsequent experiments suggested that mid-infrared spectroscopy was a useful tool in predicting methane emissions from milk fatty acid composition . In spite of these observed correlations, in the RuminOmics project, the predictive value of individual fatty acid concentrations in more than 200 fatty acids measured was weak for methane emissions. The link between milk fatty acids and methane is the ruminal microbiota. Key species differ in their fatty acid composition , thus different minor fatty acids derived from these species appear in milk depending on the abundance of different members of the microbial community. Since the community of H2-producing bacteria, for example, has an influence on methanogenesis, the quantities of their fatty acids in milk can indicate their abundance in the rumen and therefore indirectly their effect on methanogenesis. In the RuminOmics project, the predictive value of individual fatty acid concentrations for methane emissions was weak. Multiple correlations were found, however, few of them had been observed previously.
Current -omics technologies can provide detailed information about the animal genome, the ruminal metagenome and their respective functional activities from the metatranscriptome and metaproteome. Comparative analysis using these technologies allows us to characterise the interaction between the animal and its rumen microbiota. At the present time, it is mainly the power and potential of metagenomics, metatranscriptomics and metaproteomics that are being investigated, with fewer studies investigating their application to problems associated with animal production. Furthermore, integrating the results of various meta-omics analyses remains a challenge. Improving the data present in public databases to include progressively more information on rumen microbial species is a priority. Indeed, research groups around the world are joining forces to meet these challenges. A much larger knowledge base for rumen microbial genomics will allow these methods to become more robust for the detection of relevant species as well as for a correct identification and quantification of microbial genes and proteins directly related to rumen metabolic pathways, which could have an important role in the improvement of livestock productions and breeding programmes. Some progress has been made with methane emissions. However, an arguably more acceptable strategy, particularly to the livestock producer, would be to focus on the efficiency of feed utilisation rather than methane itself. The equally important issue of N emissions has received too little attention.
Food and Agriculture Organisation of the United Nations. Livestock’s long shadow: environmental issues and options. Rome: FAO; 2006.
Tapio I, Snelling TJ, Strozzi F, Wallace RJ. The ruminal microbiome associated with methane emissions from ruminant livestock. J Anim Sci Biotechnol. 2017 (accepted)
Newbold CJ, de la Fuente G, Belanche A, Ramos-Morales E, McEwan NR. The role of ciliate protozoa in the rumen. Front Microbiol. 2015;6:1313.
IPCC. Climate change 2014: synthesis report. Contribution of working groups I, II and III to the fifth assessment report of the intergovernmental panel on climate change. In: Pachauri RK, Meyer LA, editors. Core writing team. Geneva: IPCC; 2014. p. 151.
Murray RM, Bryant AM, Leng RA. Rates of production of methane in the rumen and large intestine of sheep. Br J Nutr. 1976;36:1–14.
Martin C, Morgavi DP, Doreau M. Methane mitigation in ruminants: from microbe to the farm scale. Animal. 2010;4:351–65.
McAllister TA, Meale SJ, Valle E, Guan LL, Zhou M, Kelly WJ, et al. Ruminant nutrition symposium: use of genomics and transcriptomics to identify strategies to lower ruminal methanogenesis. J Anim Sci. 2015;93:1431–49.
Morgavi DP, Kelly WJ, Janssen PH, Attwood GT. Rumen microbial (meta)genomics and its application to ruminant production. Animal. 2013;7:184–201.
Kumar S, Choudhury PK, Carro MD, Griffith GW, Dagar SS, Puniya M, et al. New aspects and strategies for methane mitigation from ruminants. Appl Microbiol Biotechnol. 2014;98:31–44.
Yang C, Rooke JA, Cabeza I, Wallace RJ. Nitrate and inhibition of ruminal methanogenesis: microbial ecology, obstacles, and opportunities for lowering methane emissions from ruminant livestock. Front Microbiol. 2016;7:132.
Attwood GT, Altermann E, Kelly WJ, Leahy SC, Zhang L, Morrison M. Exploring rumen methanogen genomes to identify targets for methane mitigation strategies. Anim Feed Sci Technol. 2011;166–67:65–75.
Leahy SC, Kelly WJ, Altermann E, Ronimus RS, Yeoman CJ, Pacheco DM, et al. The genome sequence of the rumen methanogen Methanobrevibacter ruminantium reveals new possibilities for controlling ruminant methane emissions. PLoS One. 2010;5:e8926.
Wright ADG, Kennedy P, O’Neill CJ, Toovey AF, Popovski S, Rea SM, et al. Reducing methane emissions in sheep by immunization against rumen methanogens. Vaccine. 2004;22:3976–85.
Pinares-Patino CS, Hickey SM, Young EA, Dodds KG, MacLean S, Molano G, et al. Heritability estimates of methane emissions from sheep. Animal. 2013;7:316–21.
Goopy JP, Robinson DL, Woodgate RT, Donaldson AJ, Oddy VH, Vercoe PE, et al. Estimates of repeatability and heritability of methane production in sheep using portable accumulation chambers. Anim Prod Sci. 2015;56:116–22.
de Haas Y, Windig JJ, Calus MPL, Dijkstra J, de Haan M, Bannink A, et al. Genetic parameters for predicted methane production and potential for reducing enteric emissions through genomic selection. J Dairy Sci. 2011;94:6122–34.
Robinson DL, Goopy JP, Hegarty RS, Oddy VH, Thompson AN, Toovey AF, et al. Genetic and environmental variation in methane emissions of sheep at pasture. J Anim Sci. 2014;92:4349–63.
Pickering NK, Chagunda MG, Banos G, Mrode R, McEwan JC, Wall E. Genetic parameters for predicted methane production and laser methane detector measurements. J Anim Sci. 2015;93:11–20.
Herd RM, Arthur PF, Donoghue KA, Bird SH, Bird-Gardiner T, Hegarty RS. Measures of methane production and their phenotypic relationships with dry matter intake, growth, and body composition traits in beef cattle. J Anim Sci. 2014;92:5267–74.
Roehe R, Dewhurst RJ, Duthie CA, Rooke JA, McKain N, Ross DW, et al. 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. PLoS Genet. 2016;12:e1005846.
van Engelen S, Bovenhuis H, Dijkstra J, van Arendonk JA, Visker MH. Short communication: genetic study of methane production predicted from milk fat composition in dairy cows. J Dairy Sci. 2015;98:8223–6.
Knapp JR, Laur GL, Vadas PA, Weiss WP, Tricarico JM. Invited review: enteric methane in dairy cattle production—quantifying the opportunities and impact of reducing emissions. J Dairy Sci. 2014;97:3231–61.
Hristov AN, Oh J, Giallongo F, Frederick TW, Harper MT, Weeks HL, et al. An inhibitor persistently decreased enteric methane emission from dairy cows with no negative effect on milk production. Proc Natl Acad Sci USA. 2015;112:10663–8.
Leng RA, Nolan JV. Nitrogen metabolism in the rumen. J Dairy Sci. 1984;67:1072–89.
Broderick GA, Wallace RJ, Orskov ER. Control of rate and extent of protein degradation. In: Tsuda T, Sasaki Y, Kawashima R, editors. Physiological aspects of digestion and metabolism in ruminants. London: Academic Press Ltd; 1991. p. 541–92.
Cheng KJ, Costerton JW. Adherent rumen bacteria—their role in the digestion of plant material, urea and epithelial cells. In: Ruckebusch Y, Thivend P, editors. Digestive physiology and metabolism in ruminants. Lancaster: MTP Press; 1980. p. 227–50.
Wallace RJ, Newbold CJ, Bequette BJ, MacRae JC, Lobley GE. Increasing the flow of protein from ruminal fermentation. Asian Australas J Anim Sci. 2001;14:885–93.
Patton M, Moss J, Zhang L, Kim IS, Binfield J, Westhoff P. FAPRI-UK greenhouse gas emission modelling system for England, Wales, Scotland and Northern Ireland. http://randd.defra.gov.uk/Document.aspx?Document=9859_GHGPaper(Dec2010).pdf (2010). Accessed on 30 Dec 2016.
Ferrer M, Beloqui A, Golyshina OV, Plou FJ, Neef A, Chernikova TN, et al. Biochemical and structural features of a novel cyclodextrinase from cow rumen metagenome. Biotechnol J. 2007;2:207–13.
Ferrer M, Golyshina OV, Chernikova TN, Khachane AN, Reyes-Duarte D, Dos Santos VAPM, et al. Novel hydrolase diversity retrieved from a metagenome library of bovine rumen microflora. Environ Microbiol. 2005;7:1996–2010.
Bao L, Huang Q, Chang L, Sun Q, Zhou J, Lu H. Cloning and characterization of two beta-glucosidase/xylosidase enzymes from yak rumen metagenome. Appl Biochem Biotechnol. 2012;166:72–86.
Ko KC, Lee JH, Han Y, Choi JH, Song JJ. A novel multifunctional cellulolytic enzyme screened from metagenomic resources representing ruminal bacteria. Biochem Biophys Res Commun. 2013;441:567–72.
Rashamuse KJ, Visser DF, Hennessy F, Kemp J, Roux-van der Merwe MP, Badenhorst J, et al. Characterisation of two bifunctional cellulase-xylanase enzymes isolated from a bovine rumen metagenome library. Curr Microbiol. 2013;66:145–51.
Beloqui A, Pita M, Polaina J, Martinez-Arias A, Golyshina OV, Zumarraga M, et al. Novel polyphenol oxidase mined from a metagenome expression library of bovine rumen—biochemical properties, structural analysis, and phylogenetic relationships. J Biol Chem. 2006;281:22933–42.
Liu K, Wang J, Bu D, Zhao S, McSweeney C, Yu P, et al. Isolation and biochemical characterization of two lipases from a metagenomic library of China Holstein cow rumen. Biochem Biophys Res Commun. 2009;385:605–11.
Privé F, Newbold CJ, Kaderbhai NN, Girdwood SG, Golyshina OV, Golyshin PN, et al. Isolation and characterization of novel lipases/esterases from a bovine rumen metagenome. Appl Microbiol Biotechnol. 2015;99:5475–85.
Nelson KE, Zinder SH, Hance I, Burr P, Odongo D, Wasawo D, et al. Phylogenetic analysis of the microbial populations in the wild herbivore gastrointestinal tract: insights into an unexplored niche. Environ Microbiol. 2003;5:1212–20.
Kenters N, Henderson G, Jeyanathan J, Kittelmann S, Janssen PH. Isolation of previously uncultured rumen bacteria by dilution to extinction using a new liquid culture medium. J Microbiol Methods. 2011;84:52–60.
Brulc JM, Antonopoulos DA, Miller MEB, Wilson MK, Yannarell AC, Dinsdale EA, et al. Gene-centric metagenomics of the fiber-adherent bovine rumen microbiome reveals forage specific glycoside hydrolases. Proc Natl Acad Sci USA. 2009;106:1948–53.
Hess M, Sczyrba A, Egan R, Kim TW, Chokhawala H, Schroth G, et al. Metagenomic discovery of biomass-degrading genes and genomes from cow rumen. Science. 2011;331:463–7.
Denman SE, McSweeney CS. The early impact of genomics and metagenomics on ruminal microbiology. Annu Rev Anim Biosci. 2015;3:447–65.
Singh KM, Patel AK, Shah RK, Reddy B, Joshi CG. Potential functional gene diversity involved in methanogenesis and methanogenic community structure in Indian buffalo (Bubalus bubalis) rumen. J Appl Genet. 2015;56:411–26.
Denman SE, Martinez Fernandez G, Shinkai T, Mitsumori M, McSweeney CS. Metagenomic analysis of the rumen microbial community following inhibition of methane formation by a halogenated methane analog. Front Microbiol. 2015;6:1087.
Wallace RJ, Rooke JA, McKain N, Duthie CA, Hyslop JJ, Ross DW, et al. The rumen microbial metagenome associated with high methane production in cattle. BMC Genomics. 2015;16:839.
Taxis TM, Wolff S, Gregg SJ, Minton NO, Zhang C, Dai J, et al. The players may change but the game remains: network analyses of ruminal microbiomes suggest taxonomic differences mask functional similarity. Nucleic Acids Res. 2015;43:9600–12.
Shabat SK, Sasson G, Doron-Faigenboim A, Durman T, Yaacoby S, Berg Miller ME, et al. Specific microbiome-dependent mechanisms underlie the energy harvest efficiency of ruminants. ISME J. 2016;10:2958–72.
Hegarty RS, Goopy JP, Herd RM, McCorkell B. Cattle selected for lower residual feed intake have reduced daily methane production. J Anim Sci. 2007;85:1479–86.
Counotte GHM, Prins RA, Janssen RHAM, Debie MJA. Role of Megasphaera elsdenii in the fermentation of DL-[2-13C] lactate in the rumen of dairy cattle. Appl Environ Microbiol. 1981;42:649–55.
Stewart CS, Flint HJ, Bryant MP. The rumen bacteria. In: Hobson PN, Stewart CS, editors. The rumen microbial ecosystem. London: Chapman & Hall; 1997. p. 10–72.
Rychlik JL, LaVera R, Russell JB. Amino acid deamination by ruminal Megasphaera elsdenii strains. Curr Microbiol. 2002;45:340–5.
Wallace RJ. Catabolism of amino acids by Megasphaera elsdenii LC1. Appl Environ Microbiol. 1986;51:1141–3.
Hagg FM, Erasmus LJ, Henning PH, Coertze RJ. The effect of a direct fed microbial (Megasphaera elsdenii) on the productivity and health of Holstein cows. S Afr J Anim Sci. 2010;40:101–12.
Henning PH, Horn CH, Steyn DG, Meissner HH, Hagg FM. The potential of Megasphaera elsdenii isolates to control ruminal acidosis. Anim Feed Sci Technol. 2010;157:13–9.
Ling JR. The golden anniversary of rumen and ruminant studies. Lett Appl Microbiol. 1994;19:291–3.
Henderson G, Cox F, Ganesh S, Jonker A, Young W, Janssen PH. Rumen microbial community composition varies with diet and host, but a core microbiome is found across a wide geographical range. Sci Rep. 2015;5:14567.
Weimer PJ. Redundancy, resilience, and host specificity of the ruminal microbiota: implications for engineering improved ruminal fermentations. Front Microbiol. 2015;6:296.
Qi M, Wang P, O’Toole N, Barboza PS, Ungerfeld E, Leigh MBA, et al. Snapshot of the eukaryotic gene expression in muskoxen rumen-a metatranscriptomic approach. PLoS One. 2011;6:e20521.
Shi WB, Moon CD, Leahy SC, Kang DW, Froula J, Kittelmann S, et al. Methane yield phenotypes linked to differential gene expression in the sheep rumen microbiome. Genome Res. 2014;24:1517–52.
Costa KC, Yoon SH, Pan M, Burn JA, Baliga NS, Leigh JA. Effects of H2 and formate on growth yield and regulation of methanogenesis in Methanococcus maripaludis. J Bacteriol. 2013;195:1456–62.
Neijssel OM, Teixeira de Mattos MJ. The energetics of bacterial growth: a reassessment. Mol Microbiol. 1994;13:179–82.
Thauer RK, Kaster AK, Seedorf H, Buckel W, Hedderich R. Methanogenic archaea: ecologically relevant differences in energy conservation. Nat Rev Microbiol. 2008;6:579–91.
Wilmes P, Bond PL. Metaproteomics: studying functional gene expression in microbial ecosystems. Trends Microbiol. 2006;14:92–7.
Klaassens ES, De Vos WM, Vaughan EE. Metaproteomics approach to study the functionality of the microbiota in the human infant gastrointestinal tract. Appl Environ Microbiol. 2007;73:1388–92.
Chourey K, Jansson J, VerBerkmoes N, Shah M, Chavarria KL, Tom LM, et al. Direct cellular lysis/protein extraction protocol for soil metaproteomics. J Proteome Res. 2010;9:6615–22.
O’Farrell PH. High resolution two-dimensional electrophoresis of proteins. J Biol Chem. 1975;250:4007–21.
Abram F, Gunnigle E, O’Flaherty V. Optimisation of protein extraction and 2-DE for metaproteomics of microbial communities from anaerobic wastewater treatment biofilms. Electrophoresis. 2009;30:4149–51.
Benndorf D, Vogt C, Jehmlich N, Schmidt Y, Thomas H, Woffendin G, et al. Improving protein extraction and separation methods for investigating the metaproteome of anaerobic benzene communities within sediments. Biodegradation. 2009;20:737–50.
Wu L, Wang H, Zhang Z, Lin R, Zhang Z, Lin W. Comparative metaproteomic analysis on consecutively Rehmannia glutinosa-monocultured rhizosphere soil. PLoS One. 2011;6:e20611.
Snelling TJ, Wallace RJ. The rumen microbial metaproteome as revealed by SDS-PAGE. BMC Microbiol. 2017;17:9. doi:10.1186/s12866-016-0917-y
Deusch S, Seifert J. Catching the tip of the iceberg: evaluation of sample preparation protocols for metaproteomic studies of the rumen microbiota. Proteomics. 2015;15:3590–5.
Saleem F, Bouatra S, Guo AC, Psychogios N, Mandal R, Dunn SM, et al. The bovine ruminal fluid metabolome. Metabolomics. 2013;9:360–78.
Zhao S, Zhao J, Bu D, Sun P, Wang J, Dong Z. Metabolomics analysis reveals large effect of roughage types on rumen microbial metabolic profile in dairy cows. Lett Appl Microbiol. 2014;59:79–85.
McCartney CA, Bull ID, Yan T, Dewhurst RJ. Assessment of archaeol as a molecular proxy for methane production in cattle. J Dairy Sci. 2013;96:1211–7.
McCartney CA, Bull ID, Dewhurst RJ. Chemical markers for rumen methanogens and methanogenesis. Animal. 2013;7:409–17.
Schwarm A, Schweigel-Rontgen M, Kreuzer M, Ortmann S, Gill F, Kuhla B, et al. Methane emission, digestive characteristics and faecal archaeol in heifers fed diets based on silage from brown midrib maize as compared to conventional maize. Arch Anim Nutr. 2015;69:159–76.
McCartney CA, Bull ID, Dewhurst RJ. Using archaeol to investigate the location of methanogens in the ruminant digestive tract. Livest Sci. 2014;164:39–45.
McCartney CA, Dewhurst RJ, Bull ID. Changes in the ratio of tetraether to diether lipids in cattle feces in response to altered dietary ratio of grass silage and concentrates. J Anim Sci. 2014;92:4095–8.
Nolan JV. Quantitative models of nitrogen metabolism in sheep. In: McDonald IW, Warner ACI, editors. Digestion and metabolism in the ruminant. Armidale: University of New England Publishing Unit; 1975. p. 416–31.
Chen XB, Mathieson J, Hovell FDD, Reeds PJ. Measurement of purine derivatives in urine of ruminants using automated methods. J Sci Food Agric. 1990;53:23–33.
Chen XB, Chen YK, Franklin MF, Orskov ER, Shand WJ. The effect of feed intake and body weight on purine derivative excretion and microbial protein supply in sheep. J Anim Sci. 1992;70:1534–42.
Morgavi DP, Rathahao-Paris E, Popova M, Boccard J, Nielsen KF, Boudra H. Rumen microbial communities influence metabolic phenotypes in lambs. Front Microbiol. 2015;6:1060.
van Lingen HJ, Crompton LA, Hendriks WH, Reynolds CK, Dijkstra J. Meta-analysis of relationships between enteric methane yield and milk fatty acid profile in dairy cattle. J Dairy Sci. 2014;97:7115–32.
van Gastelen S, Dijkstra J. Prediction of methane emission from lactating dairy cows using milk fatty acids and mid-infrared spectroscopy. J Sci Food Agric. 2016;96:3963–8.
Vlaeminck B, Fievez V, Tamminga S, Dewhurst RJ, van Vuuren A, De Brabander D, et al. Milk odd- and branched-chain fatty acids in relation to the rumen fermentation pattern. J Dairy Sci. 2006;89:3954–64.
The authors wrote the manuscript together, RJW having initiated the project. All authors read and approved the final manuscript.
The Rowett Institute of Nutrition and Health is funded by the Rural and Environment Science and Analytical Services Division (RESAS) of the Scottish Government. This study was financially supported by RuminOmics (Project No. 289319 of EC 7th Framework Programme: Food, Agriculture, Fisheries and Biotechnology).
The authors declare that they have no competing interests.
An erratum to this article is available at http://dx.doi.org/10.1186/s12711-017-0304-7.
About this article
Cite this article
Wallace, R.J., Snelling, T.J., McCartney, C.A. et al. Application of meta-omics techniques to understand greenhouse gas emissions originating from ruminal metabolism. Genet Sel Evol 49, 9 (2017). https://doi.org/10.1186/s12711-017-0285-6
- Methane Emission
- Glycosyl Hydrolase
- Archaeal Community
- Milk Fatty Acid
- Urinary Metabolome