Methane production from rumen methanogenesis contributes approximately 71% of greenhouse gas emissions from the agricultural sector. This study has performed genomic predictions for methane production from 99 sheep across 3 yr using a residual methane phenotype that is log methane yield corrected for live weight, rumen volume, and feed intake. Using genomic relationships, the prediction accuracies (as determined by the correlation between predicted and observed residual methane production) ranged from 0.058 to 0.220 depending on the time point being predicted. The best linear unbiased prediction algorithm was then applied to relationships between animals that were built on the rumen metabolome and microbiome. Prediction accuracies for the metabolome-based relationships for the two available time points were 0.254 and 0.132; the prediction accuracy for the first microbiome time point was 0.142. The second microbiome time point could not successfully predict residual methane production. When the metabolomic relationships were added to the genomic relationships, the accuracy of predictions increased to 0.274 (from 0.201 when only the genomic relationship was used) and 0.158 (from 0.081 when only the genomic relationship was used) for the two time points, respectively. When the microbiome relationships from the first time point were added to the genomic relationships, the maximum prediction accuracy increased to 0.247 (from 0.216 when only the genomic relationship was used), which was achieved by giving the genomic relationships 10 times more weighting than the microbiome relationships. These accuracies were higher than the genomic, metabolomic, and microbiome relationship matrixes achieved alone when identical sets of animals were used.

Ruminants are significant contributors to the livestock generated component of the greenhouse gas, methane (CH₄). The CH₄ is primarily produced by the rumen microbes. Although the composition of the diet and animal intake amount have the largest effect on CH₄ production and yield (CH₄ production/dry matter intake, DMI), the host also influences CH₄ yield. Shorter rumen feed mean retention time (MRT) is associated with higher dry matter intake and lower CH₄ yield, but the molecular mechanism(s) by which the host affects CH₄ production remain unclear. We integrated rumen wall transcriptome data and CH₄ phenotypes from two independent experiments conducted with sheep in Australia (AUS, n = 62) and New Zealand (NZ, n = 24). The inclusion of the AUS data validated the previously identified clusters and gene sets representing rumen epithelial, metabolic and muscular functions. In addition, the expression of the cell cycle genes as a group was consistently positively correlated with acetate and butyrate concentrations (p < 0.05, based on AUS and NZ data together). The expression of a group of metabolic genes showed positive correlations in both AUS and NZ datasets with CH₄ production (p < 0.05) and yield (p < 0.01). These genes encode key enzymes in the ketone body synthesis pathway and included members of the poorly characterized aldo-keto reductase 1C (AKR1C ) family. Several AKR1C family genes appear to have ruminant specific evolution patterns, supporting their specialized roles in the ruminants. Combining differential gene expression in the rumen wall muscle of the shortest and longest MRT AUS animals (no data available for the NZ animals) with correlation and network analysis, we identified a set of rumen muscle genes involved in cell junctions as potential regulators of MRT, presumably by influencing contraction rates of the smooth muscle component of the rumen wall. Higher rumen expression of these genes, including SYNPO (synaptopodin, p < 0.01) and NEXN (nexilin, p < 0.05), was associated with lower CH₄ yield in both AUS and NZ datasets. Unlike the metabolic genes, the variations in the expression of which may reflect the availability of rumen metabolites, the muscle genes are currently our best candidates for causal genes that influence CH₄ yield.