Including microbiome information in a multi-trait genomic evaluation: a case study on longitudinal growth performance in beef cattle

BACKGROUND

Growth rate is an important component of feed conversion efficiency in cattle and varies across the different stages of the finishing period. The metabolic effect of the rumen microbiome is essential for cattle growth, and investigating the genomic and microbial factors that underlie this temporal variation can help maximize feed conversion efficiency at each growth stage.

RESULTS

By analysing longitudinal body weights during the finishing period and genomic and metagenomic data from 359 beef cattle, our study demonstrates that the influence of the host genome on the functional rumen microbiome contributes to the temporal variation in average daily gain (ADG) in different months (ADG1, ADG2, ADG3, ADG4). Five hundred and thirty-three additive log-ratio transformed microbial genes (alr-MG) had non-zero genomic correlations (rg) with at least one ADG-trait (ranging from |0.21| to |0.42|). Only a few alr-MG correlated with more than one ADG-trait, which suggests that a differential host-microbiome determinism underlies ADG at different stages. These alr-MG were involved in ribosomal biosynthesis, energy processes, sulphur and aminoacid metabolism and transport, or lipopolysaccharide signalling, among others. We selected two alternative subsets of 32 alr-MG that had a non-uniform or a uniform rg sign with all the ADG-traits, regardless of the rg magnitude, and used them to develop a microbiome-driven breeding strategy based on alr-MG only, or combined with ADG-traits, which was aimed at shaping the rumen microbiome towards increased ADG at all finishing stages. Combining alr-MG information with ADG records increased prediction accuracy of genomic estimated breeding values (GEBV) by 11 to 22% relative to the direct breeding strategy (using ADG-traits only), whereas using microbiome information, only, achieved lower accuracies (from 7 to 41%). Predicted selection responses varied consistently with accuracies. Restricting alr-MG based on their rg sign (uniform subset) did not yield a gain in the predicted response compared to the non-uniform subset, which is explained by the absence of alr-MG showing non-zero rg at least with more than one of the ADG-traits.

CONCLUSIONS

Our work sheds light on the role of the microbial metabolism in the growth trajectory of beef cattle at the genomic level and provides insights into the potential benefits of using microbiome information in future genomic breeding programs to accurately estimate GEBV and increase ADG at each finishing stage in beef cattle.

Using in-abattoir 3-dimensional measurements from images of beef carcasses for the prediction of EUROP classification grade and carcass weight

Imaging technology can aid the automatic extraction of measurements from beef carcasses, which can be used for objective grading. Many abattoirs, however, rely on manual grading due to the required infrastructure and cost, making technology unfeasible. This study explores 3-dimensional (3D) imaging technology, requiring limited infrastructure, and its ability to predict carcass weight, conformation class and fat class for non-invasive, objective classification. Time-of-flight near-infrared cameras captured 3-dimensional point clouds of beef carcasses, on-line in one commercial abattoir in Scotland, over a 6-month period. Thirty-five 3D images were captured per carcass and processed using machine vison software. Seventy-four measurements were extracted from each point cloud. Removal of extreme outliers resulted in 285,109 datapoints for 17,250 carcasses. Coefficients of variation (CV) for each measurement on a per-animal basis were low and consistent, and measurements were averaged across images. Using a training and validation dataset (70:30), multiple linear regression models predicted EUROP conformation class, fat class, and carcass weight. Stepwise models included fixed effects (sex, breed type, kill date (and cold carcass weight for conformation and fat class)), and 3D image measurements. Including 3D measurements resulted in prediction accuracies of 70%, 50% and 23% for cold carcass weight, conformation, and fat class respectively. Mapping predictions on the traditional EUROP grid used in the UK showed that 99% of conformation classes and 93% of fat classes were classified within the correct or neighbouring grade. The results of this study indicate the potential for non-invasive, in-abattoir technology requiring limited infrastructure to predict carcass traits objectively.

Intrinsic calf factors associated with the behavior of healthy pre-weaned group-housed dairy-bred calves

Technology-derived behaviors are researched for disease detection in artificially-reared calves. Whilst existing studies demonstrate differences in behaviors between healthy and diseased calves, intrinsic calf factors (e.g., sex and birthweight) that may affect these behaviors have received little systematic study. This study aimed to understand the impact of a range of calf factors on milk feeding and activity variables of dairy-bred calves. Calves were group-housed from ~7 days to 39 days of age. Seven liters of milk replacer was available daily from an automatic milk feeder, which recorded feeding behaviors and live-weight. Calves were health scored daily and a tri-axial accelerometer used to record activity variables. Healthy calves were selected by excluding data collected 3 days either side of a poor health score or a treatment event. Thirty-one calves with 10 days each were analyzed. Mixed models were used to identify which of live-weight, age, sex, season of birth, age of inclusion into the group, dam parity, birthweight, and sire breed type (beef or dairy), had a significant influence on milk feeding and activity variables. Heavier calves visited the milk machine more frequently for shorter visits, drank faster and were more likely to drink their daily milk allowance than lighter calves. Older calves had a shorter mean standing bout length and were less active than younger calves. Calves born in summer had a longer daily lying time, performed more lying and standing bouts/day and had shorter mean standing bouts than those born in autumn or winter. Male calves had a longer mean lying bout length, drank more slowly and were less likely to consume their daily milk allowance than their female counterparts. Calves that were born heavier had fewer lying and standing bouts each day, a longer mean standing bout length and drank less milk per visit. Beef-sired calves had a longer mean lying bout length and drank more slowly than their dairy sired counterparts. Intrinsic calf factors influence different healthy calf behaviors in different ways. These factors must be considered in the design of research studies and the field application of behavior-based disease detection tools in artificially reared calves.

Enteric Methane Emissions from Dairy-Beef Steers Supplemented with the Essential Oil Blend Agolin Ruminant

Agriculture is the largest source of methane globally, and enteric methane accounts for 32% of methane emissions globally. Dairy-beef is an increasingly important contributor to the beef industry. The objective of this study was to investigate if supplementation with a blend of essential oils (Agolin Ruminant) reduced enteric methane emissions from dairy-bred steers. Methane was measured from thirty-six Holstein Friesian steers (18 control and 18 treatment) in open-circuit respiration chambers, at three time-points relative to the introduction of Agolin Ruminant: (i) -3 (pre-additive introduction co-variate), (ii) 46 days after introduction, and (iii) 116 days after introduction. A significantly lower methane yield was observed in treated animals compared to control animals at both 46 days (p < 0.05) and 116 days (p < 0.01) after the introduction of Agolin Ruminant, although there was no difference in methane production (g/day). Control animals appeared to be more affected by isolation in respiration chambers than animals receiving Agolin Ruminant, as indicated by a significant reduction in dry matter intake by control animals in respiration chambers.

Feed Conversion Ratio (FCR) and Performance Group Estimation Based on Predicted Feed Intake for the Optimisation of Beef Production

This paper reports on the use of estimates of individual animal feed intake (made using time spent feeding measurements) to predict the Feed Conversion Ratio (FCR), a measure of the amount of feed consumed to produce 1 kg of body mass, for an individual animal. Reported research to date has evaluated the ability of statistical methods to predict daily feed intake based on measurements of time spent feeding measured using electronic feeding systems. The study collated data of the time spent eating for 80 beef animals over a 56-day period as the basis for the prediction of feed intake. A Support Vector Regression (SVR) model was trained to predict feed intake and the performance of the approach was quantified. Here, feed intake predictions are used to estimate individual FCR and use this information to categorise animals into three groups based on the estimated Feed Conversion Ratio value. Results provide evidence of the feasibility of utilising the 'time spent eating' data to estimate feed intake and in turn Feed Conversion Ratio (FCR), the latter providing insights that guide farmer decisions on the optimisation of production costs.

Microbiome-driven breeding strategy potentially improves beef fatty acid profile benefiting human health and reduces methane emissions

BACKGROUND

Healthier ruminant products can be achieved by adequate manipulation of the rumen microbiota to increase the flux of beneficial fatty acids reaching host tissues. Genomic selection to modify the microbiome function provides a permanent and accumulative solution, which may have also favourable consequences in other traits of interest (e.g. methane emissions). Possibly due to a lack of data, this strategy has never been explored.

RESULTS

This study provides a comprehensive identification of ruminal microbial mechanisms under host genomic influence that directly or indirectly affect the content of unsaturated fatty acids in beef associated with human dietary health benefits C18:3n-3, C20:5n-3, C22:5n-3, C22:6n-3 or cis-9, trans-11 C18:2 and trans-11 C18:1 in relation to hypercholesterolemic saturated fatty acids C12:0, C14:0 and C16:0, referred to as N3 and CLA indices. We first identified that ~27.6% (1002/3633) of the functional core additive log-ratio transformed microbial gene abundances (alr-MG) in the rumen were at least moderately host-genomically influenced (HGFC). Of these, 372 alr-MG were host-genomically correlated with the N3 index (n=290), CLA index (n=66) or with both (n=16), indicating that the HGFC influence on beef fatty acid composition is much more complex than the direct regulation of microbial lipolysis and biohydrogenation of dietary lipids and that N3 index variation is more strongly subjected to variations in the HGFC than CLA. Of these 372 alr-MG, 110 were correlated with the N3 and/or CLA index in the same direction, suggesting the opportunity for enhancement of both indices simultaneously through a microbiome-driven breeding strategy. These microbial genes were involved in microbial protein synthesis (aroF and serA), carbohydrate metabolism and transport (galT, msmX), lipopolysaccharide biosynthesis (kdsA, lpxD, lpxB), or flagellar synthesis (flgB, fliN) in certain genera within the Proteobacteria phyla (e.g. Serratia, Aeromonas). A microbiome-driven breeding strategy based on these microbial mechanisms as sole information criteria resulted in a positive selection response for both indices (1.36±0.24 and 0.79±0.21 sd of N3 and CLA indices, at 2.06 selection intensity). When evaluating the impact of our microbiome-driven breeding strategy to increase N3 and CLA indices on the environmental trait methane emissions (g/kg of dry matter intake), we obtained a correlated mitigation response of -0.41±0.12 sd.

CONCLUSION

This research provides insight on the possibility of using the ruminal functional microbiome as information for host genomic selection, which could simultaneously improve several microbiome-driven traits of interest, in this study exemplified with meat quality traits and methane emissions. Video Abstract.

Bovine host genome acts on rumen microbiome function linked to methane emissions

Our study provides substantial evidence that the host genome affects the comprehensive function of the microbiome in the rumen of bovines. Of 1,107/225/1,141 rumen microbial genera/metagenome assembled uncultured genomes (RUGs)/genes identified from whole metagenomics sequencing, 194/14/337 had significant host genomic effects (heritabilities ranging from 0.13 to 0.61), revealing that substantial variation of the microbiome is under host genomic control. We found 29/22/115 microbial genera/RUGs/genes host-genomically correlated (|0.59| to |0.93|) with emissions of the potent greenhouse gas methane (CH4), highlighting the strength of a common host genomic control of specific microbial processes and CH4. Only one of these microbial genes was directly involved in methanogenesis (cofG), whereas others were involved in providing substrates for archaea (e.g. bcd and pccB), important microbial interspecies communication mechanisms (ABC.PE.P), host-microbiome interaction (TSTA3) and genetic information processes (RP-L35). In our population, selection based on abundances of the 30 most informative microbial genes provided a mitigation potential of 17% of mean CH4 emissions per generation, which is higher than for selection based on measured CH4 using respiration chambers (13%), indicating the high potential of microbiome-driven breeding to cumulatively reduce CH4 emissions and mitigate climate change.

Passive breath monitoring of livestock: Using factor analysis to deconvolve the cattle shed

Respiratory and metabolic diseases in livestock cost the agriculture sector billions each year, with delayed diagnosis a key exacerbating factor. Previous studies have shown the potential for breath analysis to successfully identify incidence of disease in a range of livestock. However, these techniques typically involve animal handling, the use of nasal swabs or fixing a mask to individual animals to obtain a sample of breath. Using a cohort of 26 cattle as an example, we show how the breath of individual animals within a herd can be monitored using a passive sampling system, where no such handling is required. These benefits come at the cost of the desired breath samples unavoidably mixed with the complex cocktail of odours that are present within the cattle shed. Data were analysed using positive matrix factorisation (PMF) to identify and remove non-breath related sources of VOC. In total three breath factors were identified (endogenous-, non-endogenous breath and rumen) and seven factors related to other sources within and around the cattle shed (e.g. foodcattle feed, traffic, urine and faeces). Simulation of a respiratory disease within the herd showed that the abnormal change in breath composition were captured in the residuals of the 10 factor PMF solution, highlighting the importance of their inclusion as part of the breath fraction. Increasing the number of PMF factors to 17 saw the identification of a "diseased" factor, which coincided with the visits of the three "diseased" cattle to the breath monitor platform. This work highlights the important role that factor analysis techniques can play in analysing passive breath monitoring data.

Identification of Microbial Genetic Capacities and Potential Mechanisms Within the Rumen Microbiome Explaining Differences in Beef Cattle Feed Efficiency

In this study, Bos Taurus cattle offered one high concentrate diet (92% concentrate-8% straw) during two independent trials allowed us to classify 72 animals comprising of two cattle breeds as "Low" or "High" feed efficiency groups. Digesta samples were taken from individual beef cattle at the abattoir. After metagenomic sequencing, the rumen microbiome composition and genes were determined. Applying a targeted approach based on current biological evidence, 27 genes associated with host-microbiome interaction activities were selected. Partial least square analysis enabled the identification of the most significant genes and genera of feed efficiency (VIP > 0.8) across years of the trial and breeds when comparing all potential genes or genera together. As a result, limited number of genes explained about 40% of the variability in both feed efficiency indicators. Combining information from rumen metagenome-assembled genomes and partial least square analysis results, microbial genera carrying these genes were determined and indicated that a limited number of important genera impacting on feed efficiency. In addition, potential mechanisms explaining significant difference between Low and High feed efficiency animals were analyzed considering, based on the literature, their gastrointestinal location of action. High feed efficiency animals were associated with microbial species including several Eubacterium having the genetic capacity to form biofilm or releasing metabolites like butyrate or propionate known to provide a greater contribution to cattle energy requirements compared to acetate. Populations associated with fucose sensing or hemolysin production, both mechanisms specifically described in the lower gut by activating the immune system to compete with pathogenic colonizers, were also identified to affect feed efficiency using rumen microbiome information. Microbial mechanisms associated with low feed efficiency animals involved potential pathogens within Proteobacteria and Spirochaetales, releasing less energetic substrates (e.g., acetate) or producing sialic acid to avoid the host immune system. Therefore, this study focusing on genes known to be involved in host-microbiome interaction improved the identification of rumen microbial genetic capacities and potential mechanisms significantly impacting on feed efficiency in beef cattle fed high concentrate diet.

Identification of Complex Rumen Microbiome Interaction Within Diverse Functional Niches as Mechanisms Affecting the Variation of Methane Emissions in Bovine

A network analysis including relative abundances of all ruminal microbial genera (archaea, bacteria, fungi, and protists) and their genes was performed to improve our understanding of how the interactions within the ruminal microbiome affects methane emissions (CH4). Metagenomics and CH4 data were available from 63 bovines of a two-breed rotational cross, offered two basal diets. Co-abundance network analysis revealed 10 clusters of functional niches. The most abundant hydrogenotrophic Methanobacteriales with key microbial genes involved in methanogenesis occupied a different functional niche (i.e., "methanogenesis" cluster) than methylotrophic Methanomassiliicoccales (Candidatus Methanomethylophylus) and acetogens (Blautia). Fungi and protists clustered together and other plant fiber degraders like Fibrobacter occupied a seperate cluster. A Partial Least Squares analysis approach to predict CH4 variation in each cluster showed the methanogenesis cluster had the best prediction ability (57.3%). However, the most important explanatory variables in this cluster were genes involved in complex carbohydrate degradation, metabolism of sugars and amino acids and Candidatus Azobacteroides carrying nitrogen fixation genes, but not methanogenic archaea and their genes. The cluster containing Fibrobacter, isolated from other microorganisms, was positively associated with CH4 and explained 49.8% of its variability, showing fermentative advantages compared to other bacteria and fungi in providing substrates (e.g., formate) for methanogenesis. In other clusters, genes with enhancing effect on CH4 were related to lactate and butyrate (Butyrivibrio and Pseudobutyrivibrio) production and simple amino acids metabolism. In comparison, ruminal genes negatively related to CH4 were involved in carbohydrate degradation via lactate and succinate and synthesis of more complex amino acids by γ-Proteobacteria. When analyzing low- and high-methane emitters data in separate networks, competition between methanogens in the methanogenesis cluster was uncovered by a broader diversity of methanogens involved in the three methanogenesis pathways and larger interactions within and between communities in low compared to high emitters. Generally, our results suggest that differences in CH4 are mainly explained by other microbial communities and their activities rather than being only methanogens-driven. Our study provides insight into the interactions of the rumen microbial communities and their genes by uncovering functional niches affecting CH4, which will benefit the development of efficient CH4 mitigation strategies.

Temporal stability of the rumen microbiota in beef cattle, and response to diet and supplements

BACKGROUND

Dietary intake is known to be a driver of microbial community dynamics in ruminants. Beef cattle go through a finishing phase that typically includes very high concentrate ratios in their feed, with consequent effects on rumen metabolism including methane production. This longitudinal study was designed to measure dynamics of the rumen microbial community in response to the introduction of high concentrate diets fed to beef cattle during the finishing period. A cohort of 50 beef steers were fed either of two basal diet formulations consisting of approximately 10:90 or 50:50 forage:concentrate ratios respectively. Nitrate and oil rich supplements were also added either individually or in combination. Digesta samples were taken at time points over ~ 200 days during the finishing period of the cattle to measure the adaptation to the basal diet and long-term stability of the rumen microbiota.

RESULTS

16S rRNA gene amplicon libraries were prepared from 313 rumen digesta samples and analysed at a depth of 20,000 sequences per library. Bray Curtis dissimilarity with analysis of molecular variance (AMOVA) revealed highly significant (p < 0.001) differences in microbiota composition between cattle fed different basal diets, largely driven by reduction of fibre degrading microbial groups and increased relative abundance of an unclassified Gammaproteobacteria OTU in the high concentrate fed animals. Conversely, the forage-based diet was significantly associated with methanogenic archaea. Within basal diet groups, addition of the nitrate and combined supplements had lesser, although still significant, impacts on microbiota dissimilarity compared to pre-treatment time points and controls. Measurements of the response and stability of the microbial community over the time course of the experiment showed continuing adaptation up to 25 days in the high concentrate groups. After this time point, however, no significant variability was detected.

CONCLUSIONS

High concentrate diets that are typically fed to finishing beef cattle can have a significant effect on the microbial community in the rumen. Inferred metabolic activity of the different microbial communities associated with each of the respective basal diets explained differences in methane and short chain fatty acid production between cattle. Longitudinal sampling revealed that once adapted to a change in diet, the rumen microbial community remains in a relatively stable alternate state.

Correction to: The rumen microbiome as a reservoir of antimicrobial resistance and pathogenicity genes is directly affected by diet in beef cattle

Following publication of the original article [1], the authors reported an error in the Additional file 1.

Identification of Rumen Microbial Genes Involved in Pathways Linked to Appetite, Growth, and Feed Conversion Efficiency in Cattle

The rumen microbiome is essential for the biological processes involved in the conversion of feed into nutrients that can be utilized by the host animal. In the present research, the influence of the rumen microbiome on feed conversion efficiency, growth rate, and appetite of beef cattle was investigated using metagenomic data. Our aim was to explore the associations between microbial genes and functional pathways, to shed light on the influence of bacterial enzyme expression on host phenotypes. Two groups of cattle were selected on the basis of their high and low feed conversion ratio. Microbial DNA was extracted from rumen samples, and the relative abundances of microbial genes were determined via shotgun metagenomic sequencing. Using partial least squares analyses, we identified sets of 20, 14, 17, and 18 microbial genes whose relative abundances explained 63, 65, 66, and 73% of the variation of feed conversion efficiency, average daily weight gain, residual feed intake, and daily feed intake, respectively. The microbial genes associated with each of these traits were mostly different, but highly correlated traits such as feed conversion ratio and growth rate showed some overlapping genes. Consistent with this result, distinct clusters of a coabundance network were enriched with microbial genes identified to be related with feed conversion ratio and growth rate or daily feed intake and residual feed intake. Microbial genes encoding for proteins related to cell wall biosynthesis, hemicellulose, and cellulose degradation and host–microbiome crosstalk (e.g., aguA, ptb, K01188, and murD) were associated with feed conversion ratio and/or average daily gain. Genes related to vitamin B12 biosynthesis, environmental information processing, and bacterial mobility (e.g., cobD, tolC, and fliN) were associated with residual feed intake and/or daily feed intake. This research highlights the association of the microbiome with feed conversion processes, influencing growth rate and appetite, and it emphasizes the opportunity to use relative abundances of microbial genes in the prediction of these performance traits, with potential implementation in animal breeding programs and dietary interventions.

Evaluation of reticuloruminal pH measurements from individual cattle: Sampling strategies for the assessment of herd status

The application of pH observations to clinical practice in dairy cattle is based on criteria derived primarily from single time-point observations more than 20 years ago. The aims of this study were to evaluate these criteria using data collected using continuous recording methods; to make recommendations that might improve their interpretation; and to determine the relationship between the number of devices deployed in a herd and the accuracy of the resulting estimate of the herd-mean reticuloruminal pH. The study made use of 815,475 observations of reticuloruminal pH values obtained from 75 cattle in three herds (one beef and two twice-daily milking herds) to assess sampling strategies for the diagnosis of sub-acute rumen acidosis (SARA), and to evaluate the ability of different numbers of bolus devices to accurately estimate the true herd-mean reticuloruminal pH value at any time. The traditional criteria for SARA provide low diagnostic utility, the probability of detection of animals with pH values below specified thresholds being affected by a strong effect of time of day and herd. The analysis suggests that regardless of time of feeding, sampling should be carried out in the late afternoon or evening to obtain a reasonable probability of detection of animals with pH values below the threshold level. The among-cow variation varied strongly between herds, but for a typical herd, if using reticuloruminal pH boluses to detect a predisposition to fermentation disorders while feeding a diet that is high in rapidly fermentable carbohydrates, it is recommended to use a minimum of nine boluses.

Fat accretion measurements strengthen the relationship between feed conversion efficiency and Nitrogen isotopic discrimination while rumen microbial genes contribute little

The use of biomarkers for feed conversion efficiency (FCE), such as Nitrogen isotopic discrimination (Δ15N), facilitates easier measurement and may be useful in breeding strategies. However, we need to better understand the relationship between FCE and Δ15N, particularly the effects of differences in the composition of liveweight gain and rumen N metabolism. Alongside measurements of FCE and Δ15N, we estimated changes in body composition and used dietary treatments with and without nitrates, and rumen metagenomics to explore these effects. Nitrate fed steers had reduced FCE and higher Δ15N in plasma compared to steers offered non-nitrate containing diets. The negative relationship between FCE and Δ15N was strengthened with the inclusion of fat depth change at the 3rd lumbar vertebrae, but not with average daily gain. We identified 1,700 microbial genes with a relative abundance >0.01% of which, 26 were associated with Δ15N. These genes explained 69% of variation in Δ15N and showed clustering in two distinct functional networks. However, there was no clear relationship between their relative abundances and Δ15N, suggesting that rumen microbial genes contribute little to Δ15N. Conversely, we show that changes in the composition of gain (fat accretion) provide additional strength to the relationship between FCE and Δ15N.

Identification, Comparison, and Validation of Robust Rumen Microbial Biomarkers for Methane Emissions Using Diverse Bos Taurus Breeds and Basal Diets

Previous shotgun metagenomic analyses of ruminal digesta identified some microbial information that might be useful as biomarkers to select cattle that emit less methane (CH4), which is a potent greenhouse gas. It is known that methane production (g/kgDMI) and to an extent the microbial community is heritable and therefore biomarkers can offer a method of selecting cattle for low methane emitting phenotypes. In this study a wider range of Bos Taurus cattle, varying in breed and diet, was investigated to determine microbial communities and genetic markers associated with high/low CH4 emissions. Digesta samples were taken from 50 beef cattle, comprising four cattle breeds, receiving two basal diets containing different proportions of concentrate and also including feed additives (nitrate or lipid), that may influence methane emissions. A combination of partial least square analysis and network analysis enabled the identification of the most significant and robust biomarkers of CH4 emissions (VIP > 0.8) across diets and breeds when comparing all potential biomarkers together. Genes associated with the hydrogenotrophic methanogenesis pathway converting carbon dioxide to methane, provided the dominant biomarkers of CH4 emissions and methanogens were the microbial populations most closely correlated with CH4 emissions and identified by metagenomics. Moreover, these genes grouped together as confirmed by network analysis for each independent experiment and when combined. Finally, the genes involved in the methane synthesis pathway explained a higher proportion of variation in CH4 emissions by PLS analysis compared to phylogenetic parameters or functional genes. These results confirmed the reproducibility of the analysis and the advantage to use these genes as robust biomarkers of CH4 emissions. Volatile fatty acid concentrations and ratios were significantly correlated with CH4, but these factors were not identified as robust enough for predictive purposes. Moreover, the methanotrophic Methylomonas genus was found to be negatively correlated with CH4. Finally, this study confirmed the importance of using robust and applicable biomarkers from the microbiome as a proxy of CH4 emissions across diverse production systems and environments.

Changes in feed intake during isolation stress in respiration chambers may impact methane emissions assessment

Respiration chambers are considered the ‘gold standard’ technique for measuring in vivo methane (CH4) emissions in live animals. However, the imposed isolation required may alter feeding behaviour and intake, which ultimately impact CH4 emissions. The aim of this study was to assess the impact of isolation within respiration chambers on feed intake and CH4 emissions with two different diets and breeds of beef cattle. In addition, a routine stressor (transport) was used to examine the relationship between individual stress responsiveness and changes in feed intake during isolation. Eighty-four steers (castrated males) (569 ± 5.7 kg bodyweight, BW) were divided into two groups and each group fed with one of two basal diets consisting of (g/kg dry matter, DM) either 50 : 50 (Mixed) or 8 : 92 (Concentrate) forage to concentrate ratios. Within each basal diet there were three supplementation treatments: (1) control, (2) calcium nitrate, and (3) rapeseed cake. The stress biomarkers plasma cortisol, creatine kinase (CK), and free fatty acids (FFA) were determined before (0 h) and after (30 min, 3 h, 6 h and 9 h) a 30-min journey, when steers were transported to the respiration chamber facilities. Methane emissions were measured over a 3-day period using individual respiration chambers. Dry matter intake (DMI) was assessed within the group-housed pens (4 weeks before entry to training pen), in the training pens and the chambers. Cortisol, FFA and CK increased (P < 0.05) after transport confirming a stress response. DMI (g/kg BW) decreased (P < 0.001) during isolation in the training pens (14.7 ± 0.28) and the chambers (14.3 ± 0.26) compared with that of the same animals in the group pens (16.8 ± 0.23). DMI during isolation decreased more in those animals which had an increased (P < 0.05) stress response during transport as measured by cortisol, FFA and CK. With the Mixed diet, the decline in DMI was estimated to result in an increase in CH4 (g/kg DMI) (r = 0.25, P = 0.001), which did not occur with the Concentrate diet. According to the results of this experiment, the stress associated with isolation reduces the DMI resulting in an increase in g CH4/kg DMI in fibrous diets. Habituation to isolation needs refinement in order to reduce the impact of stress on intake and therefore achieve more accurate estimates of CH4 emissions. Alternatively, modelling CH4 estimations according to behavioural and physiological changes associated with isolation stress would improve accuracy of CH4 estimations.

The rumen microbiome as a reservoir of antimicrobial resistance and pathogenicity genes is directly affected by diet in beef cattle

BACKGROUND

The emergence and spread of antimicrobial resistance is the most urgent current threat to human and animal health. An improved understanding of the abundance of antimicrobial resistance genes and genes associated with microbial colonisation and pathogenicity in the animal gut will have a major role in reducing the contribution of animal production to this problem. Here, the influence of diet on the ruminal resistome and abundance of pathogenicity genes was assessed in ruminal digesta samples taken from 50 antibiotic-free beef cattle, comprising four cattle breeds receiving two diets containing different proportions of concentrate.

RESULTS

Two hundred and four genes associated with antimicrobial resistance (AMR), colonisation, communication or pathogenicity functions were identified from 4966 metagenomic genes using KEGG identification. Both the diversity and abundance of these genes were higher in concentrate-fed animals. Chloramphenicol and microcin resistance genes were dominant in samples from forage-fed animals (P < 0.001), while aminoglycoside and streptomycin resistances were enriched in concentrate-fed animals. The concentrate-based diet also increased the relative abundance of Proteobacteria, which includes many animal and zoonotic pathogens. A high ratio of Proteobacteria to (Firmicutes + Bacteroidetes) was confirmed as a good indicator for rumen dysbiosis, with eight cases all from concentrate-fed animals. Finally, network analysis demonstrated that the resistance/pathogenicity genes are potentially useful as biomarkers for health risk assessment of the ruminal microbiome.

CONCLUSIONS

Diet has important effects on the complement of AMR genes in the rumen microbial community, with potential implications for human and animal health.

Association of Temperament and Acute Stress Responsiveness with Productivity, Feed Efficiency, and Methane Emissions in Beef Cattle: An Observational Study

The aim of this study was to assess individual differences in temperament and stress response and quantify their impact on feed efficiency, performance, and methane (CH₄) emissions in beef cattle. Eighty-four steers (castrated males) (Charolais or Luing) were used. Temperament was assessed using two standardized tests: restlessness when restrained [crush score (CS)] and flight speed (FS) on release from restraint. Over a 56-day period individual animal dry matter intake (DMI) and weekly body weight was measured. Ultrasound fat depth was measured at the end of 56 days. Average daily gain (ADG), feed conversion ratio (FCR), and residual feed intake (RFI) were calculated. After the 56-day test period, animals were transported in groups of six/week to respiration chamber facilities. Blood samples were taken before and 0, 3, 6, and 9 h after transport. Plasma cortisol, creatine kinase (CK), glucose, and free fatty acids (FFA) were determined to assess physiological stress response. Subsequently, CH₄ emissions were measured over a 3-day period in individual respiration chambers. CS (1.7 ± 0.09) and FS (1.6 ± 0.60 m/s) were repeatable (0.63 and 0.51, respectively) and correlated (r = 0.36, P < 0.001). Plasma cortisol, CK, and FFA concentrations increased after transport (P = 0.038, P = 0.006, and P < 0.001, respectively). Temperament (CS) and CK concentration were correlated (r = 0.29; P = 0.015). The extreme group analysis reveals that excitable animals (FS; P = 0.032) and higher stress response (cortisol, P = 0.007; FFA, P = 0.007; and CK, P = 0.003) were associated with lower DMI. ADG was lower in more temperamental animals (CS, P = 0.097, and FS, P = 0.030). Fat depth was greater in steers showing calmer CS (P = 0.026) and lower plasma CK (P = 0.058). Temperament did not show any relationship with RFI or CH₄ emissions. However, steers with higher cortisol showed improved feed efficiency (lower FCR and RFI) (P < 0.05) and greater CH₄ emissions (P = 0.017). In conclusion, agitated temperament and higher stress responsiveness is detrimental to productivity. A greater stress response is associated with a reduction in feed intake that may both increase the efficiency of consumed feed and the ratio of CH₄ emissions/unit of feed. Therefore, temperament and stress response should be considered when designing strategies to improve efficiency and mitigate CH₄ emissions in beef cattle.

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

Methane produced by methanogenic archaea in ruminants contributes significantly to anthropogenic greenhouse gas emissions. The host genetic link controlling microbial methane production is unknown and appropriate genetic selection strategies are not developed. We used sire progeny group differences to estimate the host genetic influence on rumen microbial methane production in a factorial experiment consisting of crossbred breed types and diets. Rumen metagenomic profiling was undertaken to investigate links between microbial genes and methane emissions or feed conversion efficiency. Sire progeny groups differed significantly in their methane emissions measured in respiration chambers. Ranking of the sire progeny groups based on methane emissions or relative archaeal abundance was consistent overall and within diet, suggesting that archaeal abundance in ruminal digesta is under host genetic control and can be used to genetically select animals without measuring methane directly. In the metagenomic analysis of rumen contents, we identified 3970 microbial genes of which 20 and 49 genes were significantly associated with methane emissions and feed conversion efficiency respectively. These explained 81% and 86% of the respective variation and were clustered in distinct functional gene networks. Methanogenesis genes (e.g. mcrA and fmdB) were associated with methane emissions, whilst host-microbiome cross talk genes (e.g. TSTA3 and FucI) were associated with feed conversion efficiency. These results strengthen the idea that the host animal controls its own microbiota to a significant extent and open up the implementation of effective breeding strategies using rumen microbial gene abundance as a predictor for difficult-to-measure traits on a large number of hosts. Generally, the results provide a proof of principle to use the relative abundance of microbial genes in the gastrointestinal tract of different species to predict their influence on traits e.g. human metabolism, health and behaviour, as well as to understand the genetic link between host and microbiome.

Methane is a highly potent greenhouse gas and ruminants are the major source of methane emissions from anthropogenic activities. Here we show in an experiment with cattle that genetic selection of low-emitting animals is a viable option based on a newly developed selection criterion. The experimental data provided a comprehensive insight into the host additive genetic influence on the microbiome, the impact of nutrition on genetics and the microbiome, and the effect of metagenomic microbial genes on the analysed traits. We developed a selection criterion to change those traits by evaluation of hosts based on the relative abundance of microbial genes. This criterion is shown to be highly informative and it is therefore suggested to be used in studies analysing different traits and species. This study provides a proof of principle that there is an additive genetic influence of the host on its microbiome and that selection for the desired host can be based on the abundance of a suite of genes in the ruminal metagenome associated with the trait. The use of this criterion will allow genetic analysis on a large number of hosts, previously a significant barrier to determination of host genetic effects on such traits.

The rumen microbial metagenome associated with high methane production in cattle

Background

Methane represents 16 % of total anthropogenic greenhouse gas emissions. It has been estimated that ruminant livestock produce ca. 29 % of this methane. As individual animals produce consistently different quantities of methane, understanding the basis for these differences may lead to new opportunities for mitigating ruminal methane emissions. Metagenomics is a powerful new tool for understanding the composition and function of complex microbial communities. Here we have applied metagenomics to the rumen microbial community to identify differences in the microbiota and metagenome that lead to high- and low-methane-emitting cattle phenotypes.

Methods

Four pairs of beef cattle were selected for extreme high and low methane emissions from 72 animals, matched for breed (Aberdeen-Angus or Limousin cross) and diet (high or medium concentrate). Community analysis was carried out by qPCR of 16S and 18S rRNA genes and by alignment of Illumina HiSeq reads to the GREENGENES database. Total genomic reads were aligned to the KEGG genes databasefor functional analysis.

Results

Deep sequencing produced on average 11.3 Gb per sample. 16S rRNA gene abundances indicated that archaea, predominantly Methanobrevibacter, were 2.5× more numerous (P = 0.026) in high emitters, whereas among bacteria Proteobacteria, predominantly Succinivibrionaceae, were 4-fold less abundant (2.7 vs. 11.2 %; P = 0.002). KEGG analysis revealed that archaeal genes leading directly or indirectly to methane production were 2.7-fold more abundant in high emitters. Genes less abundant in high emitters included acetate kinase, electron transport complex proteins RnfC and RnfD and glucose-6-phosphate isomerase. Sequence data were assembled de novo and over 1.5 million proteins were annotated on the subsequent metagenome scaffolds. Less than half of the predicted genes matched matched a domain within Pfam. Amongst 2774 identified proteins of the 20 KEGG orthologues that correlated with methane emissions, only 16 showed 100 % identity with a publicly available protein sequence.

Conclusions

The abundance of archaeal genes in ruminal digesta correlated strongly with differing methane emissions from individual animals, a finding useful for genetic screening purposes. Lower emissions were accompanied by higher Succinovibrionaceae abundance and changes in acetate and hydrogen production leading to less methanogenesis, as similarly postulated for Australian macropods. Large numbers of predicted protein sequences differed between high- and low-methane-emitting cattle. Ninety-nine percent were unknown, indicating a fertile area for future exploitation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2032-0) contains supplementary material, which is available to authorized users.

Novel insight into the genomic architecture of feed and nitrogen efficiency measured by residual energy intake and nitrogen excretion in growing pigs

BACKGROUND

Improvement of feed efficiency in pigs is of great economical and environmental interest and contributes to use limited resources efficiently to feed the world population. Genome scans for feed efficiency traits are of importance to reveal the underlying biological causes and increase the rate of genetic gain. The aim of this study was to determine the genomic architecture of feed efficiency measured by residual energy intake (REI), in association with production, feed conversion ratio (FCR) and nitrogen excretion traits through the identification of quantitative trait loci (QTL) at different stages of growth using a three generation full-sib design population which originated from a cross between Pietrain and a commercial dam line.

RESULTS

Six novel QTL for REI were detected explaining 2.7-6.1% of the phenotypic variance in REI. At growth from 60-90 kg body weight (BW), a QTL with a significant dominance effect was identified for REI on SSC14, at a similar location to the QTL for feed intake and nitrogen excretion traits. At growth from 90-120 kg BW, three QTL for REI were detected on SSC2, SSC4 and SSC7 with significant additive, imprinting and additive effects, respectively. These QTL (except for the imprinted QTL) were positionally overlapping with QTL for FCR and nitrogen excretion traits. During final growth (120-140 kg BW), a further QTL for REI was identified on SSC8 with significant additive effect, which overlapped with QTL for nitrogen excretion. During entire analysed growth (60-140 kg BW), a novel additive QTL for REI on SSC4 was observed, with no overlapping with QTL for any other traits considered.

CONCLUSIONS

The occurrence of only one overlapping QTL of REI with feed intake suggests that only a small proportion of the variance in REI was explained by change in feed intake, whereas four overlapping QTL of REI with those of nitrogen excretion traits suggests that mostly underlying factors of feed utilisation such as metabolism and protein turnover were the reason for change in REI. Different QTL for REI were identified at different growth stages, indicating that different genes are responsible for efficiency in feed utilisation at different stages of growth.

Quantitative trait loci for meat quality traits in pigs considering imprinting and epistatic effects

The aim of the research was to gain a better understanding of the genomic regulation of meat quality by investigating individual and epistatic QTL in a three-generation full-sib population (Pietrain x crossbred dam line). In total, 386 animals were genotyped for 96 markers. Analysed traits included pH, reflectance value, conductivity, and meat colour. Thirteen significant individual QTL were identified. The most significant QTL were detected on SSC1 and SSC9 for pH, on SSC4 for meat colour, and on SSC8 for conductivity, accounting for 3.4% to 4.7% of the phenotypic variance. Nine significant epistatic QTL pairs were detected accounting for between 5.7% and 10.9% of the phenotypic variance. Epistatic QTL pairs showing the largest effects were for reflectance value between two locations of SSC4, and for pH between SSC10 and SSC13, explaining 9.5% and 10.9% of the phenotypic variance, respectively. This study indicates that meat quality traits are influenced by numerous QTL as well as a complex network of interactions.

Genomic scan for quantitative trait loci of chemical and physical body composition and deposition on pig chromosome X including the pseudoautosomal region of males

A QTL analysis of pig chromosome X (SSCX) was carried out using an approach that accurately takes into account the specific features of sex chromosomes i.e. their heterogeneity, the presence of a pseudoautosomal region and the dosage compensation phenomenon. A three-generation full-sib population of 386 animals was created by crossing Pietrain sires with a crossbred dam line. Phenotypic data on 72 traits were recorded for at least 292 and up to 315 F2 animals including chemical body composition measured on live animals at five target weights ranging from 30 to 140 kg, daily gain and feed intake measured throughout growth, and carcass characteristics obtained at slaughter weight (140 kg). Several significant and suggestive QTL were detected on pig chromosome X: (1) in the pseudoautosomal region of SSCX, a QTL for entire loin weight, which showed paternal imprinting, (2) closely linked to marker SW2456, a suggestive QTL for feed intake at which Pietrain alleles were found to be associated with higher feed intake, which is unexpected for a breed known for its low feed intake capacity, (3) at the telomeric end of the q arm of SSCX, QTL for jowl weight and lipid accretion and (4) suggestive QTL for chemical body composition at 30 kg. These results indicate that SSCX is important for physical and chemical body composition and accretion as well as feed intake regulation.