The inter-cattle growth variations stem from the interaction of many metabolic processes making animal selection difficult. We hypothesized that growth could be predicted using metabolomics. Urinary biomarkers of cattle feed efficiency were explored using mass spectrometry-based untargeted and targeted metabolomics. Feed intake and weight-gain was measured in steers (n=75) on forage-based growing rations (stage-1, 84 d) followed by high-concentrate finishing rations (stage-2, 84 d). Urine from days 0, 21, 42, 63, and 83 in each stage were analyzed from steers with the greater (n = 14) and least (n = 14) average-daily-gain (ADG) and comparable dry-matter-intake (DMI; within 0.32 SD of the mean). Steers were slaughtered after stage-2. Adjusted fat-thickness and carcass-yield-grade increased in greater-ADG-cattle selected in stage-1, but carcass traits did not differ between ADG-selected in stage-2. Overall 85 untargeted metabolites segregated greater- and least-ADG animals, with overlap across diets (both stages) and breed type, despite sampling time effects. Total 18-BAs and 5-steroids were quantified and associated with performance and carcass quality across ADG-classification depending on the stage. Stepwise logistic regression of urinary BA and steroids had >90% accuracy identifying efficient-ADG-steers. Urine metabolomics provides new insight into the physiological mechanisms and potential biomarkers for feed efficiency.
Research Article
Non-Invasive Metabolomics Biomarkers of Production Efficiency and Beef Carcass Quality Traits
https://doi.org/10.21203/rs.3.rs-381369/v1
This work is licensed under a CC BY 4.0 License
Journal Publication
published 07 Jan, 2022
You are reading this latest preprint version
Efficient and sustainable cattle farming requires producing quality beef while reducing operating costs and production wastes to meet demand from both producers and consumers in the complex U.S. beef production system 1. The industry uses a wide range of nutritional inputs and animal types that span a diversity of geographies and climates. Cereal grains and grain by-products are the predominant dietary components during the finishing period while forage traditionally constitutes the main fraction in cattle feed during the growing period2. Breeding and nutrition programs during the last decades notably improved animal performance and carcass composition including cattle growth rates and fat deposition 3,4. These types of economically relevant traits are controlled by many genes and environmental factors that are variably expressed contributing to the moderate prediction of high priority phenotypes selection for the beef industry 5. Next-generation phenotyping using metabolomics is becoming a fundamental approach to refine trait description to improve the prediction of breeding values of the animals aligned with the objectives of selection programs 6-8. This analytical technique provides a distinctive insight into the biochemical activity caused by environmental and genetic factors, and therefore can robustly correlate with complex phenotypes and identify biomarkers of specific physiological states 9. A compiled composition of the bovine metabolome that sets the reference of metabolite values in the bovine tissues is already freely available 10. Recent systematic reviews on livestock metabolomics applied in both research and industry summarized the importance of biomarker discovery for prognostic strategies in animal health/performance as well as increasing prediction accuracies 11.
We previously identified biomarkers of the complex expression of average daily gain (ADG) in ruminal-fluid/plasma and multiple tissues by multivariate analysis during the finishing stage in steers 12-14. These observations raise the possibility of using metabolomics for status monitoring of important production traits in beef cattle. Identifying metabolites of physiological expression for production efficiency and carcass quality may be a useful tool for cattle selection to improve sustainable beef production. Therefore, a non-invasive biofluid easy to collect overtime on a routine basis without the need for special training is crucial. In the current study, we aimed to examine urine comprehensive metabolomic profiling using untargeted/targeted LC-MS technology, with the goal of identifying urinary metabolite markers that can be used to predict feed efficiency and beef carcass quality traits on different breed types of steers on forage- and concentrate-based diets.
Study Design and Samples Details
Research protocols were approved and monitored by the USDA, ARS, U.S. Meat Animal Research Center Institutional and Animal Care Committee in accordance with the Guide for the Care and Use of Agricultural Animals in Agricultural Research and Teaching (2010). All authors complied with the ARRIVE guidelines. Calves (n=75) used in this experiment were the steer progeny of Red Angus, Charolais, and Simmental sires bred by artificial insemination to four dam breeds: MARC II Composite (¼ Hereford, ¼ Angus, ¼ Gelbvieh, ¼ Simmental), MARC III Composite (¼ Hereford, ¼ Angus, ¼ Pinzgauer, ¼ Red Poll), ½ MARC II X ½ MARC III, and ½ MARC II X (¼Red Angus X ¼ Angus).
Stage- 1 (forage based-diet)
After weaning, steers received an implant containing 80 mg of trenbolone acetate and 16mg estradiol (Revalor®-IS; Merck), and were housed in a facility with Calan Broadbent electronic headgates (American Calan, Inc.,Northwood, NH) to measure individual feed intake. Steers were trained to use Calan Headgates during an adaptation period (approximately 21 d).
At the beginning of the study (0 d), mean age and weight of cattle was 271.4 ± 5.0 d and 301.7 ± 31.9 kg, respectively. The diet consisted of 30.0% chopped alfalfa, 69.8% corn silage, and 0.2% salt on a DM basis.
Feed intake was measured for 84 d with cattle weighed on d 0, 1, 21, 42, 63, 83 and 84 of the experiment. Daily feed offered was recorded and orts determined weekly. A urine sample was collected at d 0, 42, and 83 during the scheduled weigh dates.
Stage- 2 (concentrate based-diet)
Following the completion of stage 1, steers were transitioned to a concentrate-based ration over a 21-day period. On a DM basis, the diet consisted of 8.00% chopped alfalfa, 67.75% rolled corn,
20.00% wet distillers grain with solubles, and 4.25% vitamin/mineral supplement with Rumensin® (Elanco, Greenfild, IN). About 49 days after the steers received an implant containing 120 mg of trenbolone acetate and 24 mg of estradiol (Revalor®-S, Merck Animal Health, Madison, NJ) and had acclimated to the finishing ration, feed intake and body weight gain data were collected for 84 days. Body weight measurements were determined at 0, 1, 21, 42, 63, 83 and 84 days on study. Daily feed offered were recorded and orts determined weekly. At the beginning of stage-2 (0 d), mean age and weight of cattle was 397.3 ± 5.0 d and 448.6 ± 41.8 kg, respectively. A urine sample was collected at d 0, 42, and 83 during the scheduled weigh dates.
Urine collection
On days 0, 42, and 83 for each stage-1/2 urine samples (~50 mL) were collected by preputial stimulation and immediately centrifuge 10,000g for 10 min to remove particulates (4 °C) before storage at -80 °C. No urine sample could be collected by stimulation 23% (n =19) of the selected steers in stage-1 and 11% (n =9) of the selected steers in stage-2.
Steer selection by divergent growth for each stage
Steers were selected according to differences in ADG at each stage with those with the greatest and least ADG, whose dry matter intake was within 0.32 SD of the mean intake (Fig. 7). After the end of intake study (stage-2), selected steers in both stages received the same ration ad libitum and remained in the same pen until slaughter (5-8 d) at an abattoir. Carcass data was collected using the VBG2000 beef carcass grading camera (Vision For You, L.L.C., Dakota Dunes, SD).
Sample preparation for untargeted metabolomics analysis.
All urine sample extractions and LC-MS analysis were randomized across time sampling to avoid batch effects. Duplicate samples of each urine were extracted for metabolomics analysis within batches. Both extraction and chromatographic solvents were UPLC-MS Optima Grade (Fisher Chemical Ltd., Waltham, MA). Sample preparation was adapted from previous reports 45. Briefly, 50 µL urine aliquots were diluted with 150 µL of water, vortexed and centrifuged at 16,000 × g for 10 min at 4ºC. The supernatant was filtered at 0.22 µm by centrifugation in Spin-X™ LC (Fisher Chemical Ltd., Waltham, MA) membrane-tube filters (5,000 × g for 5 min at 4ºC) before transfer to LC-vials.
UPLC-qTOF spectral acquisition of urine samples.
The UPLC/MS analysis was performed using a Waters ACQUITY ultra-performance liquid-chromatography (UPLC) system (Waters Corporation, Milford, MA) equipped with an auto sampler and coupled with a hybrid triple quadrupole-time of-flight mass spectrometry (XEVO-G2-S-qTOF; Waters). Instrument calibration was performed before running the samples using 0.5 nM sodium formate. To obtain information regarding system suitability and stability, quality control (QC) samples were injected at regular intervals (every ten samples) throughout the analytical run. Quality control samples were prepared by mixing ~100 urine aliquots (10 µL), producing separate QC samples.
For analyses, samples were held at 4ºC, and analytes in 10 µL of urine extracts were separated on a 2.1 x 50 mm x 1.7 µm Acquity UPLC BEH C18 column (Waters Corp.) held at 40 ºC using a 0.4 mL/min flow of 0.1% formic acid in water (mobile phase A) and 0.1% formic acid in acetonitrile (mobile phase B) with the following gradient: 99% A from 0 to 2 min; 85% A at 4 min; 50% A at 8 min; 5% A at 12 min; 99% A at for 16 min followed by 2 min of re-equilibration. All gradient ramps were linear.
Mass spectrometry was performed in both positive and negative modes. The capillary voltage was 3.2 kV and 2.4 kV for positive and negative mode, respectively. The system parameters were set as follows: source of temperature 120ºC, desolvation temperature 350ºC, and cone gas flow (nitrogen) 25 L/h and desolvation gas flow (nitrogen) 900 L/h. Data were collected in a centroid mode using the lockspray to ensure accuracy and reproducibility. Leucine enkephalin was used as lock-mass in a 2 ng/mL concentration solution. The lockspray frequency was set at 15 s, and the lock mass data were average over 15 scans for correction. The scan mass range was from 50-1200 m/z using an extended dynamic range. The MS/MS analysis was carried out by ramping the collision energy from 10 to 50 V using argon as a collision gas.
Untargeted-metabolomics data processing. Raw data obtained from UPLC-qTOF analysis was analyzed using Progenesis QI v1.0 software (Waters Corp.). The data was aligned, deconvoluted, and normalized using total ion intensity. Solvent blanks were run between samples and each mass was checked against the blank run to exclude possible sources of contaminations. The coefficient of variation (CV) of the mass abundance for all biological replicates, was calculated across pooled repeated samples for each feature and those with a CV above a 20% cutoff were removed. Metabolites were identified by comparison with online Bovine Metabolome Database (http://www.cowmetdb.ca/) using exact m/z values and retention times. The identities of selected metabolites were confirmed by MS/MS fragment ion analysis using Mass-Fragment application manager software (Water MassLynk v4.1, Waters Corp.). The MS/MS fragmentation of the molecules was compared, with ChemSpider database (www.chemspider. com), by way of chemically intelligent peak-matching algorithms.
Urine bile acid and steroid quantification
Endogenous urinary steroids and bile acids were quantified using isotopic enrichment and authentic calibration standards detected by LC-MS/MS with electrospray ionization and multiple reaction monitoring (MRM) on an API 6500QTRAP (Sciex, Framingham, MA, USA) as previously reported for plasma 47. Urine extraction was performed by solid-phase extraction (SPE) using Oasis-HLB SPE 96-wellplates 30 mg (Milford, MA, USA). Briefly, 100 µL urine aliquots were spiked with 10 µL of methanol containing a suite of deuterated bile acids and steroids at 12.5 nM, vortexed, and loaded onto SPE cartridges pre-conditioned with 1 mL MeOH, followed by 1 mL H2O. Loaded cartridges were then washed with 1 mL H2O and eluted with 2 mL MeOH. The eluate was evaporated under vacuum at room temperature and reconstituted in 100 µL of 50:50 MeOH:H2O (v/v). Extracts were separated on a 2.1 × 100 mm, 1.7μm Acquity C18 BEH column using a Shimadzu Nexera X2 UPLC (Shimadzu, Kyoto, Japan). Analytes were quantified using internal standard methods with surrogate/analyte associations for response ratio calculations. Data was processed with AB Sciex MultiQuant v 3.0.1.
Statistical analysis
Untargeted urine metabolomics analysis
The resulting data matrix was imported into MetaboAnalyst version 5.0 48 for subsequent univariate/multivariate analysis. To account for repeated measures structure of the data, the focus of our methodology for analysis was to explain the metabolic changes across steers over time and stage. ANOVA-simultaneous component analysis (ASCA) was conducted to identify the effect of time or stage (diet) on the metabolome urine profile of an ADG-classification (least-ADG vs. greater-ADG steers) using a significance level of P <0.05. Permutation testing (1000 times) was performed to minimize the possibility that observed effects were by chance. Leverage was used to evaluate the importance of the metabolite to the model, and square prediction error (SPE) was used as a test of the fitness of the model for the particular metabolite. Variables with low SPE and higher leverage had a significant contribution to the model and were selected as influentially affected compounds. Two-way ANOVA univariate analysis was conducted to evaluate the effect of breed type on metabolomic profiles of selected ADG-groups. Post hoc analysis with false discovery rate (FDR) correction was performed on both main effects level and interaction of the metabolites identified. Hierarchical cluster analyses were used to group metabolite by ADG-classification time-dependent changes (FDR <0.05). The interactive principal components analysis (iPCA) was conducted in both two-way –ANOVO and ASCA to identify and visualize stage and breed type effects of least-ADG vs. greater-ADG steers. Pathway analysis was performed using a Bos taurus pathway library, which integrate global pathway enrichment analysis and relative between centrality pathway topology analysis from MetaboAnalyst 3.0 software. The identification and visualization of the top altered pathways were based on KEGG (http://www.genome.jp/kegg/) database sources 49. Receiver operating characteristic curve (ROC) analysis was performed to evaluate the minimum number of untargeted metabolites that have the sensitivity and specificity to distinguish between ADG- groups using a linear support vector machine (SVM) algorithm. These models were assessed by area under the curve (AUC) and prediction accuracies by cross-validation. The significance of the model was further validated using permutation test.
Targeted urine bile acids and steroids
BA and steroid concentrations were log-transformed to achieve normality. The generalized linear mixed model (GLIMMIX) procedure of SAS (SAS Institute, Cary, North Carolina) were used to examine the bile acid and steroid responses to experimental variables. The model included time, stage, breed type, ADG-classification and interactions as fixed effect, residual of steers nested by stage as a random effect, with time as the repeated measure, and spatial power covariance as structure of sampling interval.
Pearson’s correlations were performed to identify the associations between urine bile acids and steroid concentrations with performance and carcass traits. Correlation strength was considered strong (|r| ≥ 0.7), moderate (0.7 > |r| ≥ 0.4) or weak (0.4 > |r| ≥ 0.1) (15). A predictive equation for ADG-classification using bile acids and steroids log transform concentrations was performed by stepwise nominal logistic regression analysis with a JMP Pro version 15.0 (SAS Institute, Cary, North Carolina). Basically, for each ADG-category, multiple regression is run to predict the odds ratio into that category. The odds ratios are then converted into probabilities, and the most likely class is chosen. For ADG-prediction model building, data were divided into a 75/25 training:validation set splits to avoid over-fitting and minimize model reliance on sample selection. Validation set was stratified by ADG-classification to balance training set representation and used to test model fit probabilities.
Availability of Data and Materials
All data generated or analyzed during this study are included in this published article and its Supplementary files.
Characteristics of steers in the selected ADG classification. In Table 1, the performance and carcass traits of selected ADG-classification is provided. The dry matter intake (DMI) did not differ between ADG groups in both stages; ADG was greater in the greater-ADG group (0.99 ± 0.08 kg/d; 2.25 ± 0.08 kg/d) than the least-ADG group (0.69 ± 0.05 kg/d; 1.90 ± 0.08 kg/d) for stage-1 and stage-2, respectively (P < 0.01). The feed conversion ratio was lower in the greater- ADG than least-ADG in both stages (P <0.01). The residual feed intake (RFI) differed between ADG groups in stage-1 (P =0.01) but not in stage-2. Final body weight of the steers did not differ (P =0.42) within ADG-classification on 1/2-stages. Adjusted-fat-thickness and image analysis yield-grade improved in the greater-ADG in stage-1 (P <0.05) while in stage-2 carcass traits did not differ (Table 1). There was no breed type effect on performance and carcass characteristics for ADG-groups for either stage (Supplementary Table S1). Steers classified as greater- or least-ADG in stage-1 were not the same steers classified in stage-2.
Urine metabolomic characterization of ADG-classification
Overall, 85 metabolites were identified between the greater/least-ADG groups in both stages (1-forage diet; and 2- concentrate diet) and the six time samplings during growth (Supplementary table S2). In order to analyze the overall metabolic shift over the entire timeframe of the study, ANOVA and ASCA models were constructed with the resulting scores graphically displayed and grouped according to the time component of the study design. A scree plot was used to graphically display the relationship between eigenvalues and factors. The major patterns associated with factor A (time or stage variation), factor B (greater/least-ADG) and their interaction were shown in the score scatter plots (Supplementary Figure S1) based on PC1 of the corresponding sub models. In a follow-up process, a significance test was performed by a permutation approach to validate the model for the ADG-classification, time or stage, and their interaction. Based on Supplementary Fig. S1, it is evident that although ADG-classification and time both have a dramatic impact (P<0.01) on the urine metabolome, these changes were independent as indicated by no interaction (P>0.8). However, stage and the interaction of stage by ADG-classification was significant (Fig. 1a). The significant variables associated with a specific factor were identified based on the leverage/squared prediction error (SPE) plots (Fig. 1b). A total of 12 metabolites were greatly altered by ADG-classification, 10 metabolites by stage and 16 metabolites by interactions between ADG-classification and stage including bile acids, steroids and fatty acids (Fig. 1b, Supplementary Table S3).
Considering time as an independent variable, a univariate two-way ANOVA analysis was performed to define the impact of ADG-classification for each stage, time, and their interactive effects. As shown in Figure 2a, 59 metabolites showed a significant ADG-classification/stage effect, 54 metabolites showed a significant time effect, and 4 metabolites showed a significant interaction effect for each time point differing by ADG-classification (Supplementary Table S4). In total 59 metabolites differentially expressed between least/greater ADG were consistently expressed over time such that they could be considered as suitable model metabolites to reflect the ADG-specific phenotype. A heatmap (Fig. 3) was constructed to present the metabolic alterations of these identified metabolites of ADG-classification by stage (FDR <0.05). Pathway enrichment analysis showed significant enrichment (P-values < 0.05, Supplementary Table S5) of several pathways in the ADG-classification, including steroid biosynthesis, primary bile acid metabolism, biosynthesis of unsaturated fatty acids alpha‑linolenic acid metabolism, vitamin B6 metabolism, cysteine and methionine metabolism.
The iPCA represents the untargeted metabolites separation by ADG-classification and the result interaction with stage and breed type. (Fig. 4). There were 16 metabolites with identified breed type effects and 39 metabolites showing interaction of breed type that differentiated ADG-classification by each stage in the two-way ANOVA (Fig 2b; Supplementary table S5). There were no identified untargeted metabolites in the MARC 3-Charolais cross in low-ADG phenotype in either stage (Fig. 4c).
As a result of the ROC analysis of the untargeted metabolites, potential biomarkers were identified as androsterone sulfate, 11-hydroxyprogesterone-glucuronide deoxycholic acid glycine conjugated and taurochenodeoxycholate in stage-1 (AUC=0.93, CI:0.83-0.99) while dexocholic acid glycine conjugated, progesterone and taurochenodeoxycholate were identified in stage-2 (AUC=0.99, CI:0.96-1). Permutation test in both models were P <0.01 (Supplementary Fig. S2).
Urine bile acid and steroidconcentrations of ADG-groups
Most bile acid (BA) and all detected steroid concentrations were affected by time of sampling (P <0.05). Sampling time was again identified to be significantly and independently associated with ADG-groups (Supplementary Table S7). However, there was an interaction of ADG-groups with stage effect for the urine concentration of corticosterone (P =0.05) unconjugated BAs (P =0.04), and, the sum of all bile acids (P =0.05) (Supplementary Table S7). Of the total BAs in urine, unconjugated BA (%) were higher in least-ADG groups as compared to greater-ADG in stage-1, and higher than both ADG-phenotypes at stage-2 (Fig. 5). As a percent of total detected steroids, glucocorticoids tended to be lower (P =0.10) and sex-steroids to be higher in least-ADG compare to greater- ADG in stage-1 but were not different in stage-2.
There were distinctive associations of BA and steroids with performance and carcass quality across ADG-classification depending on stage (Fig. 6). In both stages, the feed conversion ratio (FCR; i.e. the DMI/ADG) FCR was positively associated with marbling (r> 0.34; P =0.05) in least-ADG-groups while FCR was negatively associated with adjusted-fat-thickness (r> -0.31; P =0.05) in greater-ADG-groups. In stage-2, the concentration of unconjugated-BA was associated with marbling (r =0.48; P <0.01), ribeye-area (r =0.36; P <0.01), and carcass-weight (r =0.53; P <0.001) in greater-ADG but were not associated in least-ADG-group. The sex-steroids (sum of the concentration of testosterone and progesterone) were negatively associated (P <0.01) with carcass traits in least-ADG-groups in both stages. Glucocorticoids (sum of the concentration of cortisol, cortisone and corticosterone) were associated with marbling differently in stage-1 in least-ADG (r =-0.45; P =0.01), and stage-2 in greater-ADG-group (r =0.34; P =0.05).
A predictive equation for least/greater-ADG was developed using a stepwise nominal logistic regression using log transformed urine bile acid and steroid concentrations (R2=0.61; P <0.01; AUC=0.97). The standard errors of the estimates are in supplementary table S8. Cutoff values of bile acids and steroids to predict least/greater-ADG are in supplementary table S9. The accuracy of ADG prediction by urine bile acids and steroids was 91.2% for greater-ADG and 96.7% for least-ADG.
Improving production efficiency while obtaining quality meat is an increasing agricultural need from farm-to-table. Livestock are selected for performance by measuring the amount of feed consumed for individual animals using high-cost special facilities and intensive labor 15. Multiple individual records are needed to identify an efficient phenotype that is expected to be moderately heritable 15,16. With the rapid development of molecular biology and related disciplines, animal breeding has evolved from conventional to molecular breeding 17. While marker-assisted selection using genomics has become a reliable mainstream practice, the SNPs and GWAS do not fully explain the molecular basis of phenotypic variations in complex traits 18. In this study, we identify novel biomarkers and biochemical pathways in urine that have the potential to improve cattle selection for feed efficiency and carcass quality promising both, a deeper metabolic understanding of ADG-expression and, ultimately, a high sensitivity screening tool to refine cattle selection. The advantage of sampling urine, compared to other biofluids, it is noninvasive, low-cost, requires minimal time- and minimal training. The dependability is relatively good with only ~15% of steers not producing a sample.
The ADG-phenotypic expression represents one of the metrics to evaluate feed efficiency. The inter-animal ADG -variation stems from the interaction of many metabolic processes influenced, in turn, by physiological status and herd management 19. In our study, extreme ADG phenotypes were selected independently from DMI in two different diets in growing animals. The selection criteria in both stages (diets) were the same. Steers selected in stage-1 for greater ADG had higher yield grades at slaughter. In stage-2 carcass traits were not different between selected ADG-groups likely due to an increase in the growth rate of the smaller animals promoted by high-energy diets 20. The association of feed efficiency with carcass traits was previously evaluated based on RFI selection on growing beef cattle on a high-concentrate diet 21-23. While selection for RFI may result in slight negative impacts on growth and reduce carcass quality 21 others observed a reduction of the proportion of visceral organs which impact energy expenditure and hence protein turnover with a greater muscle content in the carcass 22. The evidence in the scientific literature is certainly variable in the breed, sex and stage of physiological maturity of the cattle employed. This emphasizes the relevance of understanding the biological basis of animal traits to better define sustained progress on selection for feed efficiency and its consequences on carcass traits of beef cattle.
In our study, a distinctive urine metabolic fingerprint of the ADG-phenotypic groups was identified using non-targeted metabolomics and multivariate analysis, suggesting metabolic differences between ADG-phenotypical selection similarly to previous results in metabolomics analysis in rumen 13 and multiple-tissues in beef-cattle 14. These previous differences were involved in the metabolism of alpha-linolenic acid, unsaturated fatty acids, glycerophospholipids, taurine, and hypotaurine, creatinine, cysteine and methionine as well as bile acids and steroids, which were all in common with urinary factors characterized here. We identified for the first time metabolomic differences in sampling-time, diet, and breed type but only diet and breed type were both related to ADG-metabolome differences. Breed type had an effect in almost half of the untargeted metabolites identified for selected ADG-groups and only MARC-3xCharolais metabolome were not in least-ADG in both stages. Prior metabolomics studies reported differential molecular mechanisms of breeds to adapt to the environment 9 and diet, 24 and association between feeding system and breed on fatty acids composition in beef muscle 25. Selection for feed efficiency (ADG-groups) and carcass quality is a complex trait of metabolic diverse expression in common across different tissues and biofluids, and some of these metabolic mechanisms may excel by breed and environment.
Using biomarkers analysis, we identified a small number of urine untargeted BA and steroid metabolites with excellent sensitivity and specificity for ADG-classification distinctive for each diet. Then, a predictive equation for ADG-(least/greater) -independently of diet and time sampling- was develop with a >90 % of accuracy using urine bile acid and steroid concentrations selected by stepwise analysis. A potential screening tool of feed efficiency that could be implement after weaning. Nevertheless, this algorithm requires further validations, including a larger cohort of animals, to warrants the accuracy of suggested cutoff values.
Urine steroids and BA are cholesterol-derived metabolites involved in the regulation of different biological pathways mainly in metabolic homeostasis and maintenance 26. Total BA concentration is traditionally used as a clinical index of liver pathology in cattle27, although the advance of current quantification techniques of BA and steroid profiling makes it possible to elucidate their roles in diet intervention 28, breeding 29 and metabolic functions 30. Here a total 18 BAs and 5 steroids were detected in cattle urine, with conjugated-BAs positively associated with sex-steroids and glucocorticoids in both diets. Previously, similar associations of urine bile acid and glucocorticoid concentration were shown to be related to the level of dietary cholesterol in monogastric (humans) 31. In ruminants, cholesterol is exclusively obtained by hepatic de novo synthesis 32, which is associated with cattle efficiency by reductions in both in gene expression of hepatic lipid synthesis 33 and ultimately cholesterol concentrations 14,34,35. In our results, the urine concentration of conjugated-BAs was higher in a concentrate-diet while differences by ADG-groups in conjugated-BA were only observed in a forage-based diet. In high-grain diets, conjugated-BA were increased by signature microbial population changes in the small intestine of Angus-cattle 28. Differences in urine conjugated-BA concentrations in relation to cattle efficiency may reflect the hepatic capacity of cholesterol synthesis, transport, and oxidation, which in turn depends on the diet. An increased conjugated-BA, furthermore, triggers metabolic activities of microbes (e.g. glycerophospholipids) in the lower GI tract of beef cattle, and is suggested to be one of the links of increased unhealthy fatty acids in beef resulting from high-corn diets 28. Thus, the association of urine BA concentrations with carcass quality may reflect differences in nutrient absorption at the small intestine and/or changes in the intestinal microbiome.
The endogenous concentration of sex-steroids and glucocorticoids were previously suggested as a useful tool for the prediction of beef carcass quality 36,37. Endogenous sex-hormones exhibit anabolic effects. For instance, testosterone regulates body composition by increasing skeletal muscle mass while decreasing fat mass 38. Similar to our results, carcass leanness was associated with testosterone in growing beef bulls 36,39. The relationship between testosterone concentration and cattle performance -within the same animal category- has not been extensively reported. In twenty “Red Danish” bulls, plasma testosterone measured 3-times from 4 to 10 months of age was not associated with the rate of gain or linear growth 40. Nevertheless, glucocorticoid concentrations were related to both productive performance and carcass quality in our results, and the direction of their correlation depended on ADG-group and diet. Plasma cortisol and fecal corticosterone previously measured in beef cattle on a high-corn diet during the finishing-period were indicators of production traits in a contradictory manner between studies 41-43. Glucocorticoids are steroids produced by the adrenal cortex and play a key role in energy homeostasis by peripheral tissues including skeletal muscle 44. Physiologically concentration of endogenous glucocorticoids has major stimulatory effects on food intake and energy balance, and is well understood in monogastric models 45, which are a major adaptive regulatory mechanism for phenotype and environment.
Taken together, based on our data, urine analysis offers metabolic phenotypes that are both quantitative and highly targeted for feed-efficient animals without compromising growth and meat quality. Noninvasive, low-cost, and minimal-training time of urine sampling compared to other biofluids makes it an especially attractive screening tool for performance and carcass quality that could potentially be included in breeding programs to refine phenotypes and improve the accuracy of predictions of economically important traits. Although high-throughput metabolomics approaches for biomarker discovery at the molecular level are available, challenges remain to translate these technologies as affordable analytical techniques coupled to transparent and translatable data mining to farm programs. Target analysis of urinary bile acid and steroid concentrations may represent a potential screening tool to implement in breeding programs.
Author contributions
A.P.F., H.C.F. and V.M.A. designed and planned the experiment. A.P.F. conducted the feeding and growth study. V.M.A. collected urine samples. J.W.N. developed LC-MS urine bile acid and steroid targeted protocol and quantification analysis. S.D.S., D.A.K, and T.L.W. collected carcass traits. R.M.L consulted on data analysis and revised the manuscript. V.M.A performed the LC-MS analysis, data analysis, and wrote the manuscript. All authors have read and agreed to the published version of the manuscript.
Competing interests
The authors declare no competing interests.
Additional Information
Mention of trade names or commercial products in this article is solely for the purpose of providing specific information and does not imply recommendations or endorsement by the USDA. The USDA is an equal opportunity provider and employer.
- Rotz, C. A., Asem-Hiablie, S., Place, S. & Thoma, G. Environmental footprints of beef cattle production in the United States. Agric. Syst. 169, 1–13 (2019).
- Drouillard, J. S. Current situation and future trends for beef production in the United States of America - A review. Asian-australas. J. Anim. Sci. 31, 1007–1016 (2018).
- Asem-Hiablie, S., Rotz, C. A., Stout, R. & Stackhouse-Lawson, K. Management characteristics of beef cattle production in the Northern Plains and Midwest regions of the United States. Prof. Anim. Sci. 32, 736–749 (2016).
- Terry, S. A., Basarab, J. A., Guan, L. & McAllister, T. A. Invited Review: Strategies to improve the efficiency of beef cattle production. Can. J. Anim. Sci. 101 (1), 1–19 (2020).
- Retallick, K. J. et al. Genetic variance and covariance and breed differences for feed intake and average daily gain to improve feed efficiency in growing cattle. J. Anim. Sci. 95, 1444–1450 (2017).
- Novais, F. J. et al. Identification of a metabolomic signature associated with feed efficiency in beef cattle. BMC Genomics. 20, 8 (2019).
- Wang, X. & Kadarmideen, H. N. Metabolomics analyses in high-low feed efficient dairy cows reveal novel biochemical mechanisms and predictive biomarkers. Metabolites. 9, 151 (2019).
- Goldansaz, S. A. et al. Candidate serum metabolite biomarkers of residual feed intake and carcass merit in sheep. J. Anim. Sci. 98, 298 (2020).
- Liao, Y. et al. Metabolomics profiling of serum and urine in three beef cattle breeds revealed different levels of tolerance to heat stress. J. Agric. Food Chem. 66, 6926–6935 (2018).
- Foroutan, A. et al. The Bovine Metabolome. Metabolites. 10, 233 (2020).
- Goldansaz, S. A. et al. Livestock metabolomics and the livestock metabolome: A systematic review. PLoS One. 12, e0177675 (2017).
- Artegoitia, V. et al. Endocannabinoids concentrations in plasma associated with feed efficiency and carcass composition of beef steers. J. Anim. Sci. 94 (12), 5177–5181 (2016).
- Artegoitia, V. M., Foote, A. P., Lewis, R. M. & Freetly, H. C. Rumen Fluid Metabolomics Analysis Associated with Feed Efficiency on Crossbred Steers. Sci. Rep. 7, 2864 (2017).
- Artegoitia, V. M., Foote, A. P., Lewis, R. M. & Freetly, H. C. Metabolomics Profile and Targeted Lipidomics in Multiple Tissues Associated with Feed Efficiency in Beef Steers. ACS omega. 4, 3973–3982 (2019).
- Koch, R. M., Swiger, L. A., Chambers, D. & Gregory, K. E. Efficiency of feed use in beef cattle. J. Anim. Sci. 22, 486–494 (1963).
- Rolfe, K. M. et al. Genetic and phenotypic parameter estimates for feed intake and other traits in growing beef cattle, and opportunities for selection. J. Anim. Sci. 89, 3452–3459 (2011).
- Rexroad, C. et al. Genome to phenome: improving animal health, production, and well-being–a new USDA blueprint for animal genome research 2018–2027. Front. Genet. 10, 327 (2019).
- Schweer, K. R. et al. Genome-wide association study for feed efficiency traits using SNP and haplotype models. J. Anim. Sci. 96, 2086–2098 (2018).
- Berry, D. P., Crowley, J. J. & CELL BIOLOGY, S. Y. M. P. O. S. I. U. M. Genetics of feed efficiency in dairy and beef cattle. J. Anim. Sci. 91, 1594–1613(2013).
- Peng, Q. et al. Effects of different pomace and pulp dietary energy density on growth performance and intramuscular fat deposition relating mRNA expression in beef cattle. J. Food. Agric. Environ. 10, 404–407 (2012).
- Torres-Vázquez, J. A., van der Werf, J. H. J. & Clark, S. A. Genetic and phenotypic associations of feed efficiency with growth and carcass traits in Australian Angus cattle1. J. Anim. Sci. 96, 4521–4531 (2018).
- Taussat, S. et al. Relationship between feed efficiency and slaughter traits of French Charolais bulls. J. Anim. Sci. 97, 2308–2319 (2019).
- Elolimy, A. A., Abdelmegeid, M. K., McCann, J. C., Shike, D. W. & Loor, J. J. Residual feed intake in beef cattle and its association with carcass traits, ruminal solid-fraction bacteria, and epithelium gene expression. J. Anim. Sci. Biotech. 9, 67 (2018).
- Yang, Y. et al. Rumen and plasma metabolomics profiling by UHPLC-QTOF/MS revealed metabolic alterations associated with a high-corn diet in beef steers. PLoS One. 13, e0208031 (2018).
- Nuernberg, K. et al. Effect of a grass-based and a concentrate feeding system on meat quality characteristics and fatty acid composition of longissimus muscle in different cattle breeds. Livest. Prod. Sci. 94, 137–147 (2005).
- Sjövall, J. Fifty years with bile acids and steroids in health and disease. Lipids. 39, 703–722 (2004).
- Kalaitzakis, E. et al. Clinicopathologic evaluation of hepatic lipidosis in periparturient dairy cattle. J. Vet. Intern. Med. 21, 835–845 (2007).
- Liu, J. et al. Diet-induced changes in bacterial communities in the jejunum and their associations with bile acids in Angus beef cattle. Anim. Microbiome. 2, 1–14 (2020).
- Blaschka, C., Sánchez-Guijo, A., Wudy, S. A. & Wrenzycki, C. Profile of bile acid subspecies is similar in blood and follicular fluid of cattle. Vet Med Sci. 6, 167–176 (2020).
- Sánchez-Guijo, A., Blaschka, C., Hartmann, M. F., Wrenzycki, C. & Wudy, S. A. Profiling of bile acids in bovine follicular fluid by fused-core-LC–MS/MS. J. Steroid. Biochem. Mol. Biol. 162, 117–125 (2016).
- Anderson, G. W., Kenyon, C. J. & Al-Dujaili, E. A. S. Post-Prandial Changes in Salivary Glucocorticoids: Effects of Dietary Cholesterol and Associations with Bile Acid Excretion. Nutrients. 11, 360 (2019).
- Nestel, P. et al. The effect of dietary fat supplements on cholesterol metabolism in ruminants. J. Lipid Res. 19, 899–909 (1978).
- Mukiibi, R. et al. Transcriptome analyses reveal reduced hepatic lipid synthesis and accumulation in more feed efficient beef cattle. Sci. Rep. 8, 7303 (2018).
- Montanholi, Y. R. et al. Liver morphometrics and metabolic blood profile across divergent phenotypes for feed efficiency in the bovine. Acta Vet. Scand. 59, 39 (2017).
- Bourgon, S., de Amorim, M. D., Miller, S. & Montanholi, Y. Associations of blood parameters with age, feed efficiency and sampling routine in young beef bulls. Livest. Sci. 195, 27–37 (2017).
- Anderson, P. T. et al. The relationship between composition of gain and circulating hormones in growing beef bulls fed three dietary crude protein levels. J. Anim. Sci. 66, 3059–3067 (1988).
- Brandt, M. M., Keisler, D. H., Meyer, D. L., Schmidt, T. B. & Berg, E. P. Serum hormone concentrations relative to carcass composition of a random allotment of commercial-fed beef cattle12. J. Anim. Sci. 85, 267–275 (2007).
- Singh, R., Artaza, J. N., Taylor, W. E., Gonzalez-Cadavid, N. F. & Bhasin, S. Androgens stimulate myogenic differentiation and inhibit adipogenesis in C3H 10T1/2 pluripotent cells through an androgen receptor-mediated pathway. Endocrinology. 144, 5081–5088 (2003).
- Gortsema, S. R., Jacobs, J. A., Sasser, R. G., Gregory, T. L. & Bull, R. C. Effects of Endogenous Testosterone on Production and Carcass Traits in Beef Cattle. J. Anim. Sci. 39, 680–686 (1974).
- Lund-Larsen, T. R., Sundby, A., Kruse, V. & Velle, W. Relation between growth rate, serum somatomedin and plasma testosterone in young bulls. J. Anim. Sci. 44, 189–194 (1977).
- Montanholi, Y. R. et al. Assessing feed efficiency in beef steers through feeding behavior, infrared thermography and glucocorticoids. Anim. 4, 692–701 (2010).
- Montanholi, Y. R. et al. On the relationships between glucocorticoids and feed efficiency in beef cattle. LivesT. Sci. 155, 130–136 (2013).
- Foote, A. P. et al. Relationship of glucocorticoids and hematological measures with feed intake, growth, and efficiency of finishing beef cattle. J. Anim. Sci. 94, 275–283 (2016).
- Braun, T. P. & Marks, D. L. The regulation of muscle mass by endogenous glucocorticoids.Front. Physiol.6 (2015).
- Dallman, M. F. et al. Minireview: Glucocorticoids—Food Intake, Abdominal Obesity, and Wealthy Nations in 2004. Endocrinology. 145, 2633–2638 (2004).
- Chan, E. C. Y., Pasikanti, K. K. & Nicholson, J. K. Global urinary metabolic profiling procedures using gas chromatography–mass spectrometry. Nat. Protoc. 6, 1483–1499 (2011).
- Pedersen, T. L., Gray, I. J. & Newman, J. W. Plasma and serum oxylipin, endocannabinoid, bile acid, steroid, fatty acid and nonsteroidal anti-inflammatory drug quantification in a 96-well plate format. Anal. Chim. Acta. 1143, 189–200 (2021).
- Chong, J., Wishart, D. S., Xia, J. & Using MetaboAnalyst 4.0 for Comprehensive and Integrative Metabolomics Data Analysis. Curr. Protoc. Bioinformatics. 68, e86 (2019).
- Xia, J. & Wishart, D. S. Web-based inference of biological patterns, functions and pathways from metabolomic data using MetaboAnalyst. Nat. Protoc. 6, 743–760 (2011).
Table 1. Performance characteristics (means ± SE) of the steers selected at the end of stage-1(forage) and stage-2(concentrate) and carcass traits at slaughter according to differences in ADG with those with the greatest ADG and least ADG, whose dry matter intake was within 0.32 SD of the mean intake.
|
Variables |
Stage 1 |
Stage 2 |
||||||
|
Least-ADG (n=5) |
Greater-ADG (n=7) |
SE |
P-value |
Least-ADG (n=8) |
Greater-ADG (n=8) |
SE |
P-value |
|
|
Performance |
||||||||
|
DMI (kg/d) |
8.46 |
8.62 |
0.29 |
0.39 |
10.8 |
11.1 |
0.11 |
0.16 |
|
ADG (kg/d) |
0.69 |
0.99 |
0.07 |
<0.01 |
1.87 |
2.24 |
0.04 |
<0.01 |
|
FCR |
12.2 |
8.69 |
0.25 |
<0.01 |
5.82 |
4.94 |
0.16 |
<0.01 |
|
RFI |
0.30 |
-0.16 |
0.10 |
0.01 |
-0.10 |
-0.16 |
0.14 |
0.90 |
|
Initial body weight (kg) |
313 |
303 |
8.87 |
0.42 |
482 |
438 |
8.63 |
<0.01 |
|
Final body weight (kg) |
372 |
386 |
11.2 |
0.44 |
627 |
635 |
13.2 |
0.55 |
|
Age (d) |
274 |
270 |
1.45 |
0.57 |
398 |
395 |
1.55 |
0.78 |
|
Carcass traits |
||||||||
|
Marbling score |
535 |
532 |
24.5 |
0.95 |
488 |
468 |
17.8 |
0.43 |
|
Hot-carcass weight (kg) |
395 |
392 |
6.94 |
0.77 |
405 |
384 |
8.6 |
0.10 |
|
Ribeye-area (cm) |
13.3 |
12.5 |
0.42 |
0.22 |
13.9 |
13.8 |
0.52 |
0.86 |
|
Adjusted-fat-thickness (cm) |
0.46 |
0.62 |
0.55 |
0.03 |
0.43 |
0.47 |
0.04 |
0.60 |
|
Yield grade |
3.08 |
3.79 |
0.16 |
0.02 |
3.05 |
2.71 |
0.26 |
0.38 |
|
Fat thickness (cm) |
0.47 |
0.60 |
0.05 |
0.12 |
4.43 |
0.44 |
0.04 |
0.87 |
*Daily average dry matter intake (DMI); Average daily gain (ADG), Feed conversion ratio (FCR) Residual feed intake (RFI).
No competing interests reported.