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

Sarah Jade Meale, Marc D. Auffret, Mick Watson, Diego Morgavi, Gonzalo Cantalapiedra-Hijar, Carol-Anne Duthie, Rainer Roehe, Richard J. Dewhurst

► To cite this version:

Sarah Jade Meale, Marc D. Auffret, Mick Watson, Diego Morgavi, Gonzalo Cantalapiedra-Hijar, et al.. Fat accretion measurements strengthen the relationship between feed conversion efficiency and Nitrogen isotopic discrimination while rumen microbial genes contribute little. Scientific Reports, 2018, 8, 10.1038/s41598-018-22103-4. hal-01778640

HAL Id: hal-01778640

https://hal.science/hal-01778640

Submitted on 25 Apr 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

<figure><img src="image_2.png" /></figure>

OPEN

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

Sarah J. Meale<fnref n="1" />,⁵, Marc D. Auffret<fnref n="2" />, Mick Watson<fnref n="1" />,³,⁴, Diego P. Morgavi<fnref n="1" />, Gonzalo Cantalapiedra-Hijar<fnref n="1" />, Carol-Anne Duthie<fnref n="2" />, Rainer Roehe<fnref n="1" />,<fnref n="2" /> & Richard J. Dewhurst<fnref n="1" />,<fnref n="2" />

Received: 4 October 2017

Accepted: 16 February 2018

Published online: 01 March 2018

The use of biomarkers for feed conversion efficiency (FCE), such as Nitrogen isotopic discrimination (Δ¹⁵N), facilitates easier measurement and may be useful in breeding strategies. However, we need to better understand the relationship between FCE and Δ¹⁵N, particularly the effects of differences in the composition of liveweight gain and rumen N metabolism. Alongside measurements of FCE and Δ¹⁵N, 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 Δ¹⁵N in plasma compared to steers offered non-nitrate containing diets. The negative relationship between FCE and Δ¹⁵N 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 Δ¹⁵N. These genes explained 69% of variation in Δ¹⁵N and showed clustering in two distinct functional networks. However, there was no clear relationship between their relative abundances and Δ¹⁵N, suggesting that rumen microbial genes contribute little to Δ¹⁵N. Conversely, we show that changes in the composition of gain (fat accretion) provide additional strength to the relationship between FCE and Δ¹⁵N.

Improving the efficiency of ruminants is an important goal due to the rising human population, growing competition for feed and food resources and increased concerns about environmental consequences of ruminant production. Feed efficiency is central to this issue, but can be difficult to measure, so there has been recent interest in potential biomarkers using samples that are easy to collect. One such biomarker is the phenomenon of Nitrogen isotopic discrimination (Δ¹⁵N), or the difference in natural ¹⁵N abundance (δ¹⁵N) between animal proteins and the diet. Variability in Δ¹⁵N has been linked to the efficiency of N metabolism<fnref n="1" />,<fnref n="2" /> and rumen microbial N metabolism³. At the mechanistic level, differences in the rates of biochemical reactions causing N isotopic discrimination contribute to the overall observed discrimination (Δ¹⁵N). Previous studies have reported negative relationships between FCE and Δ¹⁵N⁴,⁵, suggesting that the relationship with FCE may stem from its strong association with N intake and growth, thus reflecting the incorporation of N into tissue versus its excretion in urine⁶ – hence current interest in estimating the composition of gain.

<footnote n="1">¹Université Clermont Auvergne, INRA, VetAgro Sup, UMR Herbivores, F-63122, Saint-Genès-Champanelle, France.</footnote>

<footnote n="2">²SRUC, West Mains Road, Edinburgh, EH9 3JG, United Kingdom. ³Division of Genetics and Genomics, The Roslin Institute and R(D)SVS, University of Edinburgh, Edinburgh, EH25 9RG, United Kingdom. ⁴Edinburgh Genomics, The Roslin Institute and R(D)SVS, University of Edinburgh, Edinburgh, EH25 9RG, United Kingdom. ⁵Present address: School of Agriculture and Food Sciences, The University of Queensland, Gatton, QLD 4343, Australia.</footnote>

<footnote>Correspondence and requests for materials should be addressed to R.J.D. (email: richard.dewhurst@sruc.ac.uk)</footnote>
Early work by Wattiaux and Reed³ indicated that rumen microbial metabolism can cause N isotopic discrimination, with effects particularly related to the incorporation and emission of ammonia into and out of microbial amino acids, respectively. Some previous studies suggested that rumen effects may be significant when there is relatively low variation in hepatic deamination and transamination⁷. Studies in other ecosystems such as soils⁸, oceans⁹ and in pure bacterial cultures¹⁰ have identified N isotopic discrimination associated with denitrification and it is possible that supplementation with dietary nitrate could enhance this pathway in ruminants¹¹. Additionally, supplementation with nitrate linearly increases rumen ammonia concentrations, which is more energetically favourable than the reduction of carbon dioxide to methane¹²,¹³, theoretically increasing the availability of dietary energy. Yet, the effect of nitrate supplementation on N isotopic discrimination has not been studied and should be considered when determining its suitability as a biomarker for the feed efficiency phenotype. Furthermore, dietary lipids can have large effects on the number of rumen microbes¹⁴, which may influence the processes of denitrification and conversion of nitrate to ammonia, warranting investigation of its inclusion alone, and in combination with nitrate, on N isotopic discrimination. Overall, this study aimed to improve understanding of the relationship between FCE and Δ¹⁵N by (i) using ultrasound measurements to describe the composition of liveweight gain during a FCE test; and (ii) using dietary treatments and rumen metagenomics to explore the effects of rumen N metabolism across different diets.

## Materials and Methods

This experiment was conducted at Scotland's Rural College (SRUC) Beef and Sheep Research Centre in Edinburgh in 2014. The experimental protocol was approved by SRUC's Animal Welfare and Ethical Review Body, the Animal Experiments Committee and was conducted in accordance with the requirements of the UK Animals (Scientific Procedures) Act, 1986.

**Experimental design, animals and diets.** Details of the experimental design, and growth and performance results have been reported previously¹⁵. Briefly, eighty 13–15 month old Aberdeen Angus (AAx) and Limousin (LIMx) sired steers, from an established rotational crossing study, were tested for feed efficiency. In a 2 × 4 factorial design (breed × diet), animals were offered a basal diet of 55:45 forage (grass and whole crop barley silages) to concentrate (DM basis), in which the main protein source (rapeseed meal), was replaced with 1) no additive (CONT; n = 20); 2) 128 g/kg DM maize dark grains containing ~12% lipids (LIPID; n = 20); 3) calcium nitrate at a maximum of 25 g/kg diet DM (NIT; n = 20; Calcinit, Yara, Oslo, Norway); or, 4) a combination of LIPID and NIT at the listed inclusion rates (COMB; n = 20). A full description of dietary ingredients and chemical composition is presented in Duthie *et al.*¹⁵. Steers in the NIT and COMB groups were progressively adapted to the calcium nitrate supplement over a 4 week period receiving 25, 50, 75 and 100% of the dietary nitrate each week. Treatments were balanced for sire within each breed and live-weight (BW) at the start of the experiment. Fresh water was provided for *ad libitum* intake, animals were bedded on sawdust and diets were offered using 32 electronic feeders (8 per pen, Insentec BV, Marknesse, The Netherlands) to record individual daily feed intakes.

**Feed efficiency test.** Feed efficiency was estimated over a 56 d test period. Individual dry matter intake (DMI, kg/d) was recorded daily and BW was measured weekly, before fresh feed was offered, using a calibrated weigh scale. Hyslop *et al.*¹⁶ evaluated the precision of average daily gain (ADG) measurements over alternative test lengths and demonstrated that a 56 d measurement period, with weekly weighing is sufficient for characterising ADG of finishing beef cattle. For all steers, ultrasonic muscle and fat depth was obtained at the 10th and 12th/13th ribs, and at the 3rd lumbar vertebrae at the start and end of the 56 d test using industry-standard equipment (Aloka 500, BCF Technology Ltd, Bellshill, UK). Images were analysed using Matrox Inspector 8 software (Matrox Video and Imaging Technology Europe Ltd., Middlesex, UK). Changes in body composition were calculated as the difference between the start and final depth (in mm). Following the feed efficiency test, 18 steers with extreme high or low FCE, balanced for breed type and diet were selected for further metagenomic analysis. Following the feed efficiency test, steers remained on the same experimental diets for a period of up to 14 weeks before slaughter, to allow measurement of methane emissions in all animals (once per animal in batches of six, as reported by Duthie *et al.*¹⁵). Animals were slaughtered in larger batches, comprising multiple measurement batches. Transportation to the abattoir took 1 hour and there was less than 3 hours between the animals having *ad libitum* access to feed and slaughter.

**Rumen sampling and analysis.** Post-mortem digesta samples collected immediately after the rumen was opened to be drained (<45 min post-slaughter) in a commercial abattoir, were taken from the selected 18 steers, following our previous discovery that rumen samples from slaughter or live cattle showed high microbial community similarities¹⁷. Digesta was filtered through two layers of muslin cloth and samples of rumen liquid were transferred on dry ice then stored at -20 °C until analysis.

**DNA extractions.** DNA was extracted from the rumen samples following the protocol described in Rooke *et al.*¹⁸. Briefly, the method involved a combination of repeated bead beating followed by column filtration based on the protocol of Yu and Morrison¹⁹.

**Metagenomic analysis.** Illumina TruSeq libraries were prepared from genomic DNA and sequenced on Illumina HiSeq systems 4000 by Edinburgh Genomics. Further analyses of the data set followed the same procedure previously described in Wallace *et al.*²⁰ and Roehe *et al.*¹⁷. Briefly, paired-end reads (2 × 100 bp) were generated, resulting in between 8.6 and 14.5 GB per sample (between 43.4 and 72.7 million paired reads). The microbial functional genes encoding for proteins and including the genes detailed in this study were identified using the KEGG genes database²¹. Parameters were adjusted such that all hits were reported that were equal in quality to the best hit for each genomic read. The read and best hits have to be more than 90% identical and have
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Effect of dietary additive on (a) feed conversion efficiency and (b) isotopic nitrogen discrimination in Aberdeen Angus x and Limousin x steers. <sup>a-b</sup>Boxes with different superscripts have differing means (P < 0.05).</figcaption></figure>

to be belonging to a single KEGG orthologue group to be kept in the data. If the best hits are spread over more than one KEGG orthologue group, the reads were disregarded. Read counts for KEGG orthologues were summed and normalised to the total number of hits. Then, annotated genes with a relative abundance of more than 0.01% were selected for further analysis. These data can be downloaded from the European Nucleotide Archive under accession PRJEB21624.

**Nitrogen isotopic discrimination of plasma proteins.** Samples of feeds were taken weekly during the FCE test period and blood was sampled from the caudal vein (9 ml lithium heparin evacuated tubes; Greiner Bio One Ltd., Gloucestershire, UK) on the final day of the test period (d 56). Plasma was separated by centrifugation at 2 500 × g for approximately 15 min (Thermo Scientific, Heraeus Multifuge 3SR + centrifuge, Waltham, MA, USA) and stored at -20 °C until analysis. Weekly samples of each diet were dried at 60 °C and ground. Four milligram sub-samples of each diet, and plasma samples (20 µl) were prepared for ¹⁵N analysis by isotope ratio mass spectrometry (Iso-Analytical Ltd., Crewe, UK). Nitrogen isotopic discrimination (Δ¹⁵N) was calculated as the difference between the natural abundance of ¹⁵N in plasma proteins minus the natural abundance of ¹⁵N in the diet, such that Δ¹⁵N = δ¹⁵N<sub>animal</sub> - δ¹⁵N<sub>diet</sub>. Results were expressed using the delta notation in parts per 1000 (‰) related to the international standard (air).

**Calculations and statistical analyses.** Data from one steer from the 56 d test period was discarded as the steer was removed from the trial for health reasons unconnected to the diets and treatments imposed, leaving data from 79 steers available for analysis. Across all 79 steers, growth was modelled by linear regression of BW against test date, to obtain average daily gain (ADG). Feed conversion efficiency (FCE) was calculated as the ratio of ADG/DMI.

For effects of diet, breed and their interaction on Δ¹⁵N and FCE, least squares means (LSM) were estimated using a general linear model analysis (GLM, Version 9.1 for Windows, SAS Institute Inc., Cary, NC, USA). Regression analysis to determine the contribution of body composition traits (muscle and fat depth at the 10th or 12/13th rib, or 3rd lumbar vertebrae) and ADG on the relationship between Δ¹⁵N and FCE were also examined using the GLM procedure of SAS, with diet and breed as fixed effects. However, breed was never significant and was therefore, removed from the model. The contribution of each variable to the model was cross validated using stepwise regression and Variance Inflation Factors (VIF) to check for multicollinearity of selected variables. Significance was declared at P < 0.05.

Preliminary GLM analysis was carried out to estimate the influence of the KEGG annotated genes, from the selected 18 steers, on Δ¹⁵N by fitting the significant effects (diet) as well as the relative abundance of one gene each time¹⁷. The 18 steers comprised of 4 CONT, 4 OIL, 4 NIT, and 6 COMB. The residuals of each model were normally distributed. We used partial least squares analysis (PLS, Version 9.1 for Windows, SAS Institute Inc., Cary, NC, USA) to identify the most significant genes associated with Δ¹⁵N. The PLS analysis accounts for multiple testing and the correlation between microbial genes. The model selection was based on the variable importance for projection (VIP) criterion²², whereby microbial genes with a VIP < 0.8 contributed little to the prediction.

# Results and Discussion

**Dietary nitrate increases isotopic N discrimination.** Growth and production traits of these animals have been described previously¹⁵. Briefly, no breed effect was observed on FCE. However, FCE was reduced (P = 0.02) in steers supplemented with NIT, compared to CONT and LIPID (Fig. 1). This reduction in feed efficiency with the NIT diet was mitigated with the inclusion of lipids to the nitrate diet (COMB).

The natural abundance of isotopic N (δ¹⁵N, ‰) was higher in the COMB diet, compared to NIT, LIPID or CONT diets (4.40, 4.18, 4.10 and 4.08, respectively). The δ¹⁵N of the NIT diet was however, numerically higher than non-nitrate containing diets. The elevated δ¹⁵N could be expected to reduce N isotopic discrimination in plasma proteins (Δ¹⁵N, ‰) of steers fed diets containing calcium nitrate, yet this was not observed here. Nitrogen
isotopic discrimination was not affected by breed. This was not unexpected because animals came from rota-
tional crosses of two breeds so that AAx animals had an expectation of 2/3 AA and 1/3 LIM genes and LIMx
vice versa. In contrast, changing the dietary protein source altered N isotopic discrimination such that, steers
fed diets containing calcium nitrate (NIT and COMB) had higher isotopic N discrimination, compared to those
fed non-nitrate diets (P < 0.001; Fig. 1). Wattiaux and Reed³ showed that discrimination of N isotopes can occur
during the incorporation of ammonia into amino acids and microbial proteins, and as dietary nitrate supplemen-
tation linearly increases rumen ammonia concentrations this may explain the increase in N isotopic discrimina-
tion. Similarly, the inclusion of dietary nitrate may increase nitrous oxide emissions²³,²⁴ and this may be a source
of additional N isotopic discrimination.

To explore the contribution of rumen N metabolism on N isotopic discrimination we conducted metagen-
omics analysis on rumen content samples. We revealed 4 827 microbial genes, 1 700 of which had a relative
abundance >0.01%, and we have used their relative abundances to explore possible effects of N metabolism on
overall Δ¹⁵N. The effect of dietary nitrate is of interest in relation to possible effects of denitrification on N iso-
topic discrimination. Complete denitrification of NO₃⁻ to N₂ involves four reduction steps which are catalyzed
by four main types of functional enzymes, those being nitrate reductases Nar or Nap; NO-generating nitrite
reductases NirK or NirS; N₂O generating nitric oxide reductases cNor, qNor, or qCuNor; and the two types of
nitrous oxide reductases N₂OR²⁵,²⁶. Despite deep metagenomic sequencing, we only observed two genes listed
in the KEGG database, K00370 and K00371 that are part of the denitrification pathway, the first of which corre-
sponds to the nitrate reductase, *narG*²⁷. Yet, despite these two genes also being present in the dissimilatory nitrate
reduction pathway, they were only observed in one animal at a relative abundance of 0.01%, suggesting that their
contribution to N isotopic discrimination is likely negligible in the rumen. Investigations of denitrification genes
in the rumen are scarce. The absence of denitrification genes in the rumen does not provide direct support for
the suggestion that the effect of nitrate on Δ¹⁵N can be explained by denitrification. However, it remains possible
that the KEGG representatives we mapped to, do not reflect well the rumen enzymes for denitrification and/or
denitrification may be occurring in the mouth. Peterson *et al*.²⁴ measured nitrous oxide emissions from cattle in
respiration chambers and, based on diurnal patterns of emissions relative to methane, suggested that it may be
produced by denitrifying bacteria in the mouth. Studies in dentistry literature²⁷ show nitrous oxide emissions
generated by oral bacteria utilising nitrate which is concentrated and recycled via saliva. The large volumes of
saliva produced in ruminants, as well as the regular regurgitation and chewing of partially digested feed during
rumination suggests that this process could be more important in ruminants; this would be consistent with our
Δ¹⁵N observations, but requires verification.

Furthermore, the process of dissimilatory nitrate reduction to ammonia (DNRA) is considered to be favoured over denitrification in anaerobic, reduced environments²⁸, such as the rumen. Despite the lack of evidence for discrimination during DNRA, discrimination of N isotopes may occur during the incorporation of ammonia into amino acids and microbial proteins³. Here we observed that unlike denitrification genes, enzymes involved in the biosynthesis of amino acids, including the small and large chain of glutamate synthase (NADPH/NADH; K00266 and K00266, respectively) were present in all animals at an average relative abundance of 0.354%. This is in agreement with previous work, which reported a low presence of denitrifying genes in the bovine rumen compared to those contributing to nitrate and nitrite ammonification (85 vs. 636 genes, respectively, out of a total 1233 involved in nitrogen metabolism)²⁹. Similarly, a limited number of rumen bacterial species have been proposed as denitrifiers including *Pseudomonas aeruginosa* and species of *Propionibacterium* and *Nitrosomonas*, which have been considered transient or minor rumen colonizers in the absence of dietary nitrate³⁰,³¹ further suggesting that denitrification is not a primary pathway of N metabolism in the bovine rumen.

As part of a larger study³², a network analysis was carried out based on the relative abundance of KEGG genes and found 15 distinct functional clusters of gene networks. To further investigate the influence of rumen microbial genes on N isotopic discrimination we performed a partial least squares analysis, firstly, on all genes whose relative abundance significantly differed with Δ¹⁵N. From this, 50 genes were identified belonging to either cluster 2 or 5. Subsequently, all genes from these 2 clusters (n = 195 and 92, respectively) where examined via partial least squares analysis, revealing that 26 rumen microbial genes accounted for 68.6% of variation in Δ¹⁵N (model variation was 67.6%; Fig. 2). Surprisingly, no genes from the Nitrogen metabolism pathway were identified as influential in the prediction of N isotopic discrimination. The identified genes have primary roles in carbohydrate metabolism, energy metabolism and amino acid metabolism. Similarly, no clear pattern of N isotopic discrimination in association with the relative abundance of genes was observed (Fig. 2) indicating that ruminal genes may have a very limited contribution to the process of N isotopic discrimination in ruminants under these conditions. It is possible, however, that gene expression is not intimately linked to gene presence, and as such, a metatranscriptomic analysis may reveal further insights into the role of rumen N metabolism on N isotopic discrimination. An important consideration is the time of sampling, as blood was sampled at d 56 of the feed efficiency test, and rumen contents were sampled at slaughter at d 154. Despite this temporal separation, the animals remained on the same experimental diets during the interim period and as reported by Cantalapiedra-Hijar *et al*.³³, isotopic equilibrium is achieved after 45 d on an experimental diet, indicating that sampling blood at 56 d would be expected to yield similar N isotopic discrimination as if the blood samples were collected at the time of rumen sample collection.

**Relationship between Nitrogen isotopic discrimination and FCE.** To better understand the relationship between FCE and N isotopic discrimination we conducted a regression analysis which revealed a strong negative correlation (R² = 0.534; P < 0.001; Fig. 3a) between the two variables, with the inclusion of diet. This correlation was reduced to 28.7% (P < 0.001) when excluding the effect of dietary additive. The general form of this relationship is commonly observed across different species and breeds of ruminants⁴,⁵, indicating the potential of Δ¹⁵N<sub>animal-diet</sub> as a biomarker for FCE. Previous work has shown that a significant proportion of N
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Heat map of the relative abundance of microbial genes associated with N isotopic discrimination as identified in the partial least squares analysis. The relative abundance of microbial genes (aqua = low to purple = high) changed depending on N isotopic discrimination. The coloured bars on the horizontal axis represent diet (CONT = yellow, LIPID = green, NIT = blue, COMB = red) and FCE (low = light grey, high = black). Coloured bars on the vertical axis represent the cluster in which genes were located (red = cluster 2, and orange = cluster 5). The values on the x-axis represent the corresponding Δ¹⁵N values for each animal arranged in ascending order.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Effect of Feed Conversion Efficiency on N isotopic discrimination in plasma proteins (Δ¹⁵N) determined by regression analysis. The predicted relationship between Δ¹⁵N vs FCE including diet effect ((a) R² = 0.534; P < 0.001) was strengthened with the inclusion of changes in fat depth at the 3rd lumbar vertebrae during the feed efficiency test period ((b) R² = 0.552; P < 0.001).</figcaption></figure>

isotopic discrimination results from processes involved in N metabolism – particularly in the liver¹,², and rumen microbial metabolism³. Whilst there is a relationship between these processes and overall FCE, it seems likely that there would be a stronger relationship with the efficiency of N utilization (NUE). A few studies have linked N isotopic discrimination to NUE in sheep³⁴ and dairy cattle³⁵,³⁶, but this is not so easy to accomplish in growing animals where daily N gain is difficult to measure. To explore this relationship further we included body composition parameters in our model, specifically the change in back fat depth at the 3rd lumbar vertebrae over the
FCE test period, which accounted for an additional (P < 0.001) 1.80% of variation in N isotopic discrimination (R² = 0.552; Fig. 3b). Furthermore, we considered the overall change in body composition, i.e., the average change in fat depth across the 3 measurement sites. However, the predictive ability of the model did not increase beyond that of the single measurement at the 3rd lumbar vertebrae (R² = 0.556; P = 0.801). Interestingly, the inclusion of ADG over the test period, in the basal model with diet and FCE, accounted for less variation in Δ¹⁵N (R² = 0.539; P < 0.001), compared to the inclusion of change in fat depth. The lack of contribution of ADG in the model to estimate Δ¹⁵N, may result from its high correlation with FCE. As when FCE was excluded as a predictive variable, ADG accounted for 43.9% of variability in Δ¹⁵N, including the diet effect (P < 0.001). Conversely, Δ¹⁵N and diet explained 30.2% (29.7% without diet) of the variation in FCE. The inclusion of both ADG and Δ¹⁵N explained 58.9% of variability in FCE. However, the inclusion of body composition parameters only showed significance for the change in fat depth at the 10th rib, which increased the amount of variability accounted for to 60.4% (P < 0.001; Supplementary Fig. S1). Conversely, fat depth at the 12-13th rib did not increase the predictive power of the model. This may be due to the short duration of the test period (i.e. the relatively small changes in measurements of muscle and back fat depth over a 56 d period). Nevertheless, our results suggest that the composition of gain, i.e., partitioning between fat and muscle may account for variation in Δ¹⁵N in growing steers to a greater extent than the total amount of gain (ADG). Whereas, FCE is more susceptible to changes in total gain (ADG).

## Conclusion

The small percentage of variation in Δ¹⁵N explained by rumen microbial genes and the lack of association with increasing Δ¹⁵N levels suggests that rumen microbes have very little influence on the N isotopic discrimination observed in plasma proteins of beef cattle. However, the effect of FCE on Δ¹⁵N was strengthened with the inclusion of body compositional changes, specifically fat accretion at the 3rd lumbar vertebrae during the feed efficiency test.

## References

1. Poupin, N. et al. The nature of the dietary protein impacts the tissue-to-diet 15N discrimination factors in laboratory rats. *PLoS ONE* **6**, e28046, https://doi.org/10.1371/journal.pone.0028046 (2011).

2. Martínez del Rio, C. & Carleton, S. A. How fast and how faithful: the dynamics of isotopic incorporation into animal tissues. *J Mammal* **93**, 353–359, https://doi.org/10.1644/11-MAMM-S-165.1 (2012).

3. Wattiaux, M. A. & Reed, J. D. Fractionation of nitrogen isotopes by mixed ruminal bacteria. *J Anim Sci* **73**, 257–266, https://doi.org/10.2527/1995.731257x (1995).

4. Wheadon, N. M., McGee, M., Edwards, G. R. & Dewhurst, R. J. Plasma nitrogen isotopic fractionation and feed efficiency in growing beef heifers. *Br J Nutr* **111**, 1705–1711, https://doi.org/10.1017/S0007114513004078 (2014).

5. Cantalapiedra-Hijar, G. et al. Relationship between efficiency of nitrogen utilization and isotopic nitrogen fractionation in dairy cows: contribution of digestion v. metabolism? *Animal* **10**, 221–229, https://doi.org/10.1017/S1751731115002025 (2015).

6. Cheng, L. et al. Effects of feed intake and genetics on tissue nitrogen-15 enrichment and feed conversion efficiency in sheep. *J Anim Sci* **93**, 5849–5855, https://doi.org/10.2527/jas.2015-9638 (2015).

7. Cabrita, A. R. J., Fonseca, A. J. M. & Dewhurst, R. J. Short communication: Relationship between the efficiency of utilization of feed nitrogen and ¹⁵N enrichment in casein from lactating dairy cows. *J Dairy Sci* **97**, 7225–7229, https://doi.org/10.3168/jds.2014-8297 (2014).

8. Böttcher, J., Strebel, O., Voerkelius, S. & Schmidt, H. L. Using isotope fractionation of nitrate-nitrogen and nitrate-oxygen for evaluation of microbial denitrification in a sandy aquifer. *J Hydrol* **114**, 413–424, https://doi.org/10.1016/0022-1694(90)90068-9 (1990).

9. Cline, J. D. & Kaplan, I. R. Isotopic fractionation of dissolved nitrate during denitrification in the eastern tropical north pacific ocean. *Mar Chem* **3**, 271–299, https://doi.org/10.1016/0304-4203(75)90009-2 (1975).

10. Mariotti, A. et al. Experimental determination of nitrogen kinetic isotope fractionation: Some principles; illustration for the denitrification and nitrification processes. *Plant Soil* **62**, 413–430, https://doi.org/10.1007/bf02374138 (1981).

11. Tiso, M. & Schechter, A. N. Nitrate reduction to nitrite, nitric oxide and ammonia by gut bacteria under physiological conditions. *PLoS ONE* **10**, e0119712, https://doi.org/10.1371/journal.pone.0119712 (2015).

12. Olijhoek, D. W. et al. Effect of dietary nitrate level on enteric methane production, hydrogen emission, rumen fermentation, and nutrient digestibility in dairy cows. *J Dairy Sci* **99**, 6191–6205, https://doi.org/10.3168/jds.2015-10691 (2016).

13. Ungerfeld, E. M. & Kohn, R. A. The role of thermodynamics in the control of ruminal fermentation in: *Ruminant physiology: digestion, metabolism and impact of nutrition on gene expression, immunology and stress* (eds K. Sejrsen, T. Hvelplund, & M.O. Nielsen) 55-85 (Wageningen Academic Publishers, Wageningen, the Netherlands, 2006).

14. Toral, P. G. et al. Comparison of ruminal lipid metabolism in dairy cows and goats fed diets supplemented with starch, plant oil, or fish oil. *J Dairy Sci* **99**, 301–316, https://doi.org/10.3168/jds.2015-10292 (2016).

15. Duthie, C. A. et al. The effect of dietary addition of nitrate or increase in lipid concentrations, alone or in combination, on performance and methane emissions of beef cattle. *Animal* **12**, 280–287, https://doi.org/10.1017/S175173111700146X (2018).

16. Hyslop, J. J., Duthie, C.-A., Ross, D. W., Rooke, J. A. & Roehe, R. An assessment of alternative test length periods when measuring liveweight change in finishing cattle during feed efficiency studies. *Adv Anim Biosci* **3**, 46 (2012).

17. Roehe, R. 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* **12**, e1005846, https://doi.org/10.1371/journal.pgen.1005846 (2016).

18. Rooke, J. A. et al. Hydrogen and methane emissions from beef cattle and their rumen microbial community vary with diet, time after feeding and genotype. *Br J Nutr* **112**, https://doi.org/10.1017/s0007114514000932 (2014).

19. Yu, Z. & Morrison, M. Improved extraction of PCR-quality community DNA from digesta and fecal samples. *Biotechniques* **36** (2004).

20. Wallace, R. J. et al. The rumen microbial metagenome associated with high methane production in cattle. *BMC Genomics* **16**, 839, https://doi.org/10.1186/s12864-015-2032-0 (2015).

21. KEGG. The Kyoto Encyclopedia of Genes and Genomes - KEGG. *Yeast* **17**, 48–55 (2000).

22. Wold, S. PLS for multivariate linear modelling in: *Chemometric Methods in Molecular Design* (ed H. Waterbeemd) (VCH, 1994).

23. Stackhouse, K. R., Pan, Y., Zhao, Y. & Mitloehner, F. M. Greenhouse gas and alcohol emissions from feedlot steers and calves. *J Environ Qu* **40**, 899–906, https://doi.org/10.2134/jeq.2010.0354 (2011).

24. Petersen, S. O. et al. Dietary nitrate for methane mitigation leads to nitrous oxide emissions from dairy cows. *J Environ Qu* **44**, 1063–1070, https://doi.org/10.2134/jeq.2015.02.0107 (2015).

25. Lycus, P. et al. Phenotypic and genotypic richness of denitrifiers revealed by a novel isolation strategy. *ISME J.* https://doi.org/10.1038/ismej.2017.82 (2017).
26. Shapleigh, J. Denitrifying Prokaryotes in *The prokaryotes*, Fourth edn (ed E. Rosenbers) 405–425 (Springer-Verlag, 2013).

27. Schreiber, F. et al. Denitrification in human dental plaque. *BMC Biol* **8**, 24, https://doi.org/10.1186/1741-7007-8-24 (2010).

28. Mohan, S. B. & Cole, J. A. The dissimilatory reduction of nitrate to ammonia by anaerobic bacteria in *The biology of the nitrogen cycle* (eds H. Bothe, S. J. Ferguson, & W. E. Newton) (Elsevier, 2007).

29. Brulc, J. M. et al. Gene-centric metagenomics of the fiber-adherent bovine rumen microbiome reveals forage specific glycoside hydrolases. *Proc Natl Acad Sci USA* **106**, https://doi.org/10.1073/pnas.0806191110 (2009).

30. Mitsumori, M., Ajisaka, N., Tajima, K., Kajikawa, H. & Kurihara, M. Detection of Proteobacteria from the rumen by PCR using methanotroph-specific primers. *Lett. Appl. Microbiol.* **35**, 251–255, https://doi.org/10.1046/j.1472-765X.2002.01172.x (2002).

31. Arai, H., Mizutani, M. & Igarashi, Y. Transcriptional regulation of the nos genes for nitrous oxide reductase in Pseudomonas aeruginosa. *Microbiology* **149**, 29–36, https://doi.org/10.1099/mic.0.25936-0 (2003).

32. Auffret, M. D. et al. Identification, comparison and validation of robust rumen microbial biomarkers for methane emissions using diverse Bos Taurus breeds and basal diets. *Front. Microbiol.* **8**, e02548.

33. Cantalapiedra-Hijar, G. et al. Diet–animal fractionation of nitrogen stable isotopes reflects the efficiency of nitrogen assimilation in ruminants. *Br J Nutr* **113**, 1158–1169, https://doi.org/10.1017/S0007114514004449 (2015).

34. Cheng, L., Nicol, A. M., Dewhurst, R. J. & Edwards, G. R. The effects of dietary nitrogen to water-soluble carbohydrate ratio on isotopic fractionation and partitioning of nitrogen in non-lactating sheep. *Animal* **7**, 1274–1279, https://doi.org/10.1017/S1751731113000311 (2013).

35. Cheng, L. et al. Technical note: Nitrogen isotopic fractionation can be used to predict nitrogen-use efficiency in dairy cows fed temperate pasture. *J Anim Sci* **91**, 5785–5788, https://doi.org/10.2527/jas.2012-5378 (2013).

36. Cantalapiedra-Hijar, G. et al. Natural N abundance in key amino acids from lamb muscle: exploring a new horizon in diet authentication and assessment of feed efficiency in ruminants. *J Agric Food Chem* **64**, 4058–4067, https://doi.org/10.1021/acs.jafc.6b00967 (2016).

## Acknowledgements

SJM received funding from the European Union's Horizon 2020 research and innovation program under the Marie Skłodowska-Curie grant agreement #658126. The research was funded by the Scottish Government (RESAS) as part of the 2011–2016 and 2016–2021 Programmes, as well as Department for Environment Food & Rural Affairs (DEFRA) and the devolved administrations through the UK Agricultural Greenhouse Gas Inventory Research Platform (http://www.ghgplatform.org.uk).

## Author Contributions

M.D.A., C.-A.D., R.J.D. and R.R. made substantial contributions to the conception or design of the work. M.D.A. and M.W. performed sequencing and bioinformatics. S.J.M. and M.D.A. conducted data analysis and all authors were involved in the interpretation of data. S.J.M. drafted the manuscript and prepared all figures and tables, and M.D.A., M.W., C.-A.D., G.C.-H., D.P.M., R.R. and R.J.D. revised it critically for important intellectual content. All authors approved the submitted versions and agree to be accountable for all aspects of the work.

## Additional Information

Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-018-22103-4.

**Competing Interests:** The authors declare no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_6.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2018