Large-scale analysis of sheep rumen metagenome profiles captured by reduced representation sequencing reveals individual profiles are influenced by the environment and genetics of the host.

BACKGROUND: Producing animal protein while reducing the animal's impact on the environment, e.g., through improved feed efficiency and lowered methane emissions, has gained interest in recent years. Genetic selection is one possible path to reduce the environmental impact of livestock production, but these traits are difficult and expensive to measure on many animals. The rumen microbiome may serve as a proxy for these traits due to its role in feed digestion. Restriction enzyme-reduced representation sequencing (RE-RRS) is a high-throughput and cost-effective approach to rumen metagenome profiling, but the systematic (e.g., sequencing) and biological factors influencing the resulting reference based (RB) and reference free (RF) profiles need to be explored before widespread industry adoption is possible. RESULTS: Metagenome profiles were generated by RE-RRS of 4,479 rumen samples collected from 1,708 sheep, and assigned to eight groups based on diet, age, time off feed, and country (New Zealand or Australia) at the time of sample collection. Systematic effects were found to have minimal influence on metagenome profiles. Diet was a major driver of differences between samples, followed by time off feed, then age of the sheep. The RF approach resulted in more reads being assigned per sample and afforded greater resolution when distinguishing between groups than the RB approach. Normalizing relative abundances within the sampling Cohort abolished structures related to age, diet, and time off feed, allowing a clear signal based on methane emissions to be elucidated. Genus-level abundances of rumen microbes showed low-to-moderate heritability and repeatability and were consistent between diets. CONCLUSIONS: Variation in rumen metagenomic profiles was influenced by diet, age, time off feed and genetics. Not accounting for environmental factors may limit the ability to associate the profile with traits of interest. However, these differences can be accounted for by adjusting for Cohort effects, revealing robust biological signals. The abundances of some genera were consistently heritable and repeatable across different environments, suggesting that metagenomic profiles could be used to predict an individual's future performance, or performance of its offspring, in a range of environments. These results highlight the potential of using rumen metagenomic profiles for selection purposes in a practical, agricultural setting.

Combining host and rumen metagenome profiling for selection in sheep: prediction of methane, feed efficiency, production, and health traits.

BACKGROUND: Rumen microbes break down complex dietary carbohydrates into energy sources for the host and are increasingly shown to be a key aspect of animal performance. Host genotypes can be combined with microbial DNA sequencing to predict performance traits or traits related to environmental impact, such as enteric methane emissions. Metagenome profiles were generated from 3139 rumen samples, collected from 1200 dual purpose ewes, using restriction enzyme-reduced representation sequencing (RE-RRS). Phenotypes were available for methane (CH4) and carbon dioxide (CO2) emissions, the ratio of CH4 to CH4 plus CO2 (CH4Ratio), feed efficiency (residual feed intake: RFI), liveweight at the time of methane collection (LW), liveweight at 8 months (LW8), fleece weight at 12 months (FW12) and parasite resistance measured by faecal egg count (FEC1). We estimated the proportion of phenotypic variance explained by host genetics and the rumen microbiome, as well as prediction accuracies for each of these traits. RESULTS: Incorporating metagenome profiles increased the variance explained and prediction accuracy compared to fitting only genomics for all traits except for CO2 emissions when animals were on a grass diet. Combining the metagenome profile with host genotype from lambs explained more than 70% of the variation in methane emissions and residual feed intake. Predictions were generally more accurate when incorporating metagenome profiles compared to genetics alone, even when considering profiles collected at different ages (lamb vs adult), or on different feeds (grass vs lucerne pellet). A reference-free approach to metagenome profiling performed better than metagenome profiles that were restricted to capturing genera from a reference database. We hypothesise that our reference-free approach is likely to outperform other reference-based approaches such as 16S rRNA gene sequencing for use in prediction of individual animal performance. CONCLUSIONS: This paper shows the potential of using RE-RRS as a low-cost, high-throughput approach for generating metagenome profiles on thousands of animals for improved prediction of economically and environmentally important traits. A reference-free approach using a microbial relationship matrix from log10 proportions of each tag normalized within cohort (i.e., the group of animals sampled at the same time) is recommended for future predictions using RE-RRS metagenome profiles.

Impact of breeding for reduced methane emissions in New Zealand sheep on maternal and health traits.

Enteric methane emissions from ruminants account for ∼35% of New Zealand's greenhouse gas emissions. This poses a significant threat to the pastoral sector. Breeding has been shown to successfully lower methane emissions, and genomic prediction for lowered methane emissions has been introduced at the national level. The long-term genetic impacts of including low methane in ruminant breeding programs, however, are unknown. The success of the New Zealand sheep industry is currently heavily reliant on the prolificacy, fecundity and survival of adult ewes. The objective of this study was to determine genetic and phenotypic correlations between adult maternal ewe traits (live weight, body condition score, number of lambs born, litter survival to weaning, pregnancy scanning and fleece weight), faecal and Nematodirus egg counts and measures of methane in respiration chambers. More than 9,000 records for methane from over 2,200 sheep measured in respiration chambers were collected over 10 years. Sheep were fed on a restricted diet calculated as approximately twice the maintenance. Methane measures were converted to absolute daily emissions of methane measured in g per day (CH4/day). Two measures of methane yield were recorded: the ratio of CH4 to dry matter intake (g CH4/kg DMI; CH4/DMI) and the ratio of CH4 to total gas emissions (CH4/(CH4 + CO2)). Ewes were maintained in the flocks for at least two parities. Non-methane trait data from over 8,000 female relatives were collated to estimate genetic correlations. Results suggest that breeding for low CH4/DMI is unlikely to negatively affect faecal egg counts, adult ewe fertility and litter survival traits, with no evidence for significant genetic correlations. Fleece weight was unfavourably (favourably) correlated with CH4/DMI (rg = -0.21 ± 0.09). Live weight (rg = 0.3 ± 0.1) and body condition score (rg = 0.2 ± 0.1) were positively correlated with methane yield. Comparing the two estimates of methane yield, CH4/DMI had lower heritability and repeatability. However, correlations of both measures with adult ewe traits were similar. This suggests that breeding is a suitable mitigation strategy for lowering methane yield, but wool, live weight and fat deposition traits may be affected over time and should be monitored.

A restriction enzyme reduced representation sequencing approach for low-cost, high-throughput metagenome profiling.

Microbial community profiles have been associated with a variety of traits, including methane emissions in livestock. These profiles can be difficult and expensive to obtain for thousands of samples (e.g. for accurate association of microbial profiles with traits), therefore the objective of this work was to develop a low-cost, high-throughput approach to capture the diversity of the rumen microbiome. Restriction enzyme reduced representation sequencing (RE-RRS) using ApeKI or PstI, and two bioinformatic pipelines (reference-based and reference-free) were compared to bacterial 16S rRNA gene sequencing using repeated samples collected two weeks apart from 118 sheep that were phenotypically extreme (60 high and 58 low) for methane emitted per kg dry matter intake (n = 236). DNA was extracted from freeze-dried rumen samples using a phenol chloroform and bead-beating protocol prior to RE-RRS. The resulting sequences were used to investigate the repeatability of the rumen microbial community profiles, the effect of laboratory and analytical method, and the relationship with methane production. The results suggested that the best method was PstI RE-RRS analyzed with the reference-free approach, which accounted for 53.3±5.9% of reads, and had repeatabilities of 0.49±0.07 and 0.50±0.07 for the first two principal components (PC1 and PC2), phenotypic correlations with methane yield of 0.43±0.06 and 0.46±0.06 for PC1 and PC2, and explained 41±8% of the variation in methane yield. These results were significantly better than for bacterial 16S rRNA gene sequencing of the same samples (p<0.05) except for the correlation between PC2 and methane yield. A Sensitivity study suggested approximately 2000 samples could be sequenced in a single lane on an Illumina HiSeq 2500, meaning the current work using 118 samples/lane and future proposed 384 samples/lane are well within that threshold. With minor adaptations, our approach could be used to obtain microbial profiles from other metagenomic samples.

Genetic parameters of plasma and ruminal volatile fatty acids in sheep fed alfalfa pellets and genetic correlations with enteric methane emissions1.

Animal-to-animal variation in methane (CH4) emissions determined in respiration chambers has a genetic basis, but rapid phenotyping methods that can be applied on-farm are required to enable increased genetic progress by the farming industry. Fermentation of carbohydrates in the rumen results in the formation of VFA with hydrogen (H2) as a byproduct that is used for CH4 formation. Generally, fermentation pathways leading to acetate are associated with the most H2 production, less H2 formation is associated with butyrate production, and propionate and valerate production are associated with reduced H2 production. Therefore, VFA may constitute a potential correlated proxy for CH4 emissions to enable high-throughput animal screening. The objective of the present study was to determine the genetic parameters for ruminal and plasma VFA concentrations in sheep fed alfalfa (Medicago sativa L.) pellets and their genetic (rg) and phenotypic (rp) correlations with CH4 emissions. Measurements of CH4 emissions in respiration chambers and ruminal (stomach tubing 18 h from last meal) and blood plasma (3 h post-feeding) VFA concentrations were made on 1,538 lambs from 5 birth years (2007 and 2009 to 2012) aged between 5 and 10 mo, while the animals were fed alfalfa pellets at 2.0 times maintenance requirements in 2 equal size meals (0900 and 1500 h). These measurements were repeated twice (rounds) 14 d apart. Mean (± SD) CH4 production was 24.4 ± 3.08 g/d, and the mean CH4 yield was 15.8 ± 1.51 g/kg DMI. Mean concentration of total ruminal VFA was 52.2 mM, with concentrations of acetate, propionate and butyrate of 35.97, 8.83, and 4.02 mM, respectively. Ruminal total VFA concentration had heritability (h2) and repeatability estimates (± SE) of 0.24 ± 0.05 and 0.35 ± 0.03, respectively, and similar estimates were found for acetate, propionate, and butyrate. Blood plasma concentrations of VFA had much lower estimates of h2 and repeatability than ruminal VFA. Genetic correlations with CH4 yield were greatest for total concentrations of ruminal VFA and acetate, with 0.54 ± 0.12 and 0.56 ± 0.12, respectively, which were much greater than their corresponding rp. The rp and rg of ruminal VFA proportions and blood VFAs with CH4 emissions were in general lower than for ruminal VFA concentrations. However, minor ruminal VFA proportions had also moderate rg with CH4 yield. Pre-feeding concentrations of total VFA and acetate were the strongest correlated proxies to select sheep that are genetically low CH4 emitters.

Genetic parameters of methane emissions determined using portable accumulation chambers in lambs and ewes grazing pasture and genetic correlations with emissions determined in respiration chambers.

Methane (CH4) emission traits were previously found to be heritable and repeatable in sheep fed alfalfa pellets in respiration chambers (RC). More rapid screening methods are, however, required to increase genetic progress and to provide a cost-effective method to the farming industry for maintaining the generation of breeding values in the future. The objective of the current study was to determine CH4 and carbon dioxide (CO2) emissions using several 1-h portable accumulation chamber (PAC) measurements from lambs and again as ewes while grazing ryegrass-based pasture. Many animals with PAC measurements were also measured in RC while fed alfalfa pellets at 2.0 × maintenance metabolizable energy requirements (MEm). Heritability estimates from mixed models for CH4 and CO2 production (g/d) were 0.19 and 0.16, respectively, when measured using PAC with lambs; 0.20 and 0.27, respectively, when measured using PAC with ewes; and 0.23 and 0.34, respectively, when measured using RC with lambs. For measured gas traits, repeatabilities of measurements collected 14 d apart ranged from 0.33 to 0.55 for PAC (combined lambs and ewes) and were greater at 0.65 to 0.76 for the same traits measured using RC. Genetic correlations (rg) between PAC in lambs and ewes were 0.99 for CH4, 0.93 for CH4 + CO2, and 0.85 for CH4/(CH4 + CO2), suggesting that CH4 emissions in lambs and ewes are the same trait. Genetic correlations between PAC and RC measurements were lower, at 0.62 to 0.67 for CH4 and 0.41 to 0.42 for CH4 + CO2, likely reflecting different environmental conditions associated with the protocols used with the 2 measurement methods. The CH4/(CH4 + CO2) ratio was the most similar genetic trait measured using PAC (both lambs and ewes, 63% and 66% selection efficiency, respectively) compared with CH4 yield (g/kg DMI) measured using RC. These results suggest that PAC measurements have considerable value as a rapid low-cost method to estimate breeding values for CH4 emissions in sheep.

Fatty acid synthase effects on bovine adipose fat and milk fat.

A quantitative trait locus (QTL) was identified by linkage analysis on bovine Chromosome 19 that affects the fatty acid, myristic acid (C14:0), in subcutaneous adipose tissue of pasture-fed beef cattle (99% level: experiment-wise significance). The QTL was also shown to have significant effects on ten fatty acids in the milk fat of pasture-fed dairy cattle. A positional candidate gene for this QTL was identified as fatty acid synthase (FASN), which is a multifunctional enzyme with a central role in the metabolism of lipids. Five single nucleotide polymorphisms (SNPs) were identified in the bovine FASN gene, and animals were genotyped for FASN SNPs in three different cattle resource populations. Linkage and association mapping results using these SNPs were consistent with FASN being the gene underlying the QTL. SNP substitution effects for C14:0 percentage were found to have an effect in the opposite direction in adipose fat to that in milk fat. It is concluded that SNPs in the bovine FASN gene are associated with variation in the fatty acid composition of adipose fat and milk fat.