Across-Experiment Transcriptomics of Sheep Rumen Identifies Expression of Lipid/Oxo-Acid Metabolism and Muscle Cell Junction Genes Associated With Variation in Methane-Related Phenotypes

Abstract

<p>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, <i>n</i> = 62) and New Zealand (NZ, <i>n</i> = 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 (<i>p</i> < 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 (<i>p</i> < 0.05) and yield (<i>p</i> < 0.01). These genes encode key enzymes in the ketone body synthesis pathway and included members of the poorly characterized aldo-keto reductase 1C (<i>AKR1C</i> ) family. Several <i>AKR1C</i> 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, <i>p</i> < 0.01) and <i>NEXN</i> (nexilin, <i>p</i> < 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.</p>