Skip to main content
Research Article

Untargeted Metabolomics for Beef Flavor Beyond Fat in Ground Beef With Different Lean Sources and Different Fat Content From a Common Fat Source

Authors
  • Zena M. Hicks (University of Nebraska)
  • Chris R. Kerth (Texas A&M University)
  • Kayley R. Wall (Cargill Innovation Center)
  • Stephen B. Smith (Texas A&M University)
  • Rhonda K. Miller (Texas A&M University)

Abstract

To explore the effects of lean source on the flavor of ground beef patties, beef inside rounds (n=9) were procured from each of the following treatment groups: commodity upper two-thirds choice (HC; USDA modest and higher marbling), heart-healthy-branded (HEART), natural grass-fed (NATURAL), and commodity USDA Select (SELECT) beef. Rounds from each source were ground and supplemented with commercially sourced, pre-ground commodity fat trim to form treatment batches containing10% or 20% total fat. Batches were then fine-ground and formed into 113.5-g patties. Patties were vacuum packaged and frozen until analyzed. Trained sensory panel, fatty acid profile, volatile compound composition, and metabolomic features subsequently were analyzed. Lean source had no impact (P>0.05) on any major trained sensory traits, but patties with 20% fat had higher (P<0.05) fat-like, buttery, and juiciness scores. Patties made with HEART lean had the highest (P<0.05) monounsaturated fatty acids, and those from SELECT lean had the lowest. Those patties made with NATURAL lean and 20% fat tended (0.05<P<0.10) to have higher volatile concentrations of alcohol, aldehyde, ketones, and especially terpenoid compounds. Using discriminate analysis, metabolites (n=64 metabolites) from raw samples were accurately segregated by lean source only, while the cooked patties showed that the 138 metabolites were able to discriminate lean source for HC, HEART, and NATURAL within both 10% and 20% fat treatments. Patties made with SELECT lean clustered by themselves and generally had the opposite reaction to metabolite concentration as the other lean sources. Overall, while lean sources did not impact flavor, patties with different lean sources impacted the fatty acid content, volatile aroma compounds, and metabolite distribution in ground beef patties. Metabolomics may be another valuable tool to help describe meat quality, and it could be used to determine these traits in lean prior to sensory testing.

Keywords: fat content, flavor, ground beef, HPLC, lean source, metabolomics

How to Cite:

Hicks, Z. M., Kerth, C. R., Wall, K. R., Smith, S. B. & Miller, R. K., (2023) “Untargeted Metabolomics for Beef Flavor Beyond Fat in Ground Beef With Different Lean Sources and Different Fat Content From a Common Fat Source”, Meat and Muscle Biology 7(1): 16080, 1-17. doi: https://doi.org/10.22175/mmb.16080

Rights: © 2023 Hicks, et al. This is an open access article distributed under the CC BY license.

Funding

Name
Beef Checkoff

739 Views

230 Downloads

Published on
2023-11-15

Peer Reviewed

Important findings:

  • Lean sources did not impact any major sensory traits of ground beef patties.

  • Monounsaturated fatty acids were highest in patties made with heart-healthy lean and lowest in pattie made with USDA Select lean.

  • Ground beef patties made with natural, grass-fed beef and 20% fat had higher terpenoid volatile aroma compounds.

  • Both raw and cooked metabolites metabolomics can discriminate lean source and fat content of ground beef patties.

Introduction

Flavor is the most important attribute for consumer acceptance of beef as long as tenderness is acceptable (Killinger et al., 2004; O’Quinn et al., 2012; Kerth and Miller, 2015). “Fat is flavor” is a well-known adage among the beef industry (Wasserman and Talley, 1968; Huang et al., 2010). Due to this line of thinking, the USDA quality grade system, which evaluates carcass quality based on maturity and marbling, has been a long-used tool for consumers to predict quality and for processors and retails to set prices on products. Thus, numerous studies (Corbin et al., 2014; Stewart et al., 2021) have been conducted examining marbling and fat deposition as they relate to consumer satisfaction and determining logical slaughter endpoints to increase profitability by producers. In turn, the degree of marbling has been shown to be a factor that impacts tenderness and consumer acceptability of beef products (Platter et al., 2003; Killinger et al., 2004).

In contrast, marbling has been only moderately correlated with flavor differences among carcasses from different breeds of beef cattle (Arshad et al., 2018). Research by Kerth and Miller (2015) reported that “beef identity” and “brown/roasted” flavor attributes as defined by the beef lexicon (Adhikari et al., 2011) are highly correlated to Maillard reaction products and consumer liking. The Maillard reaction, in which reducing sugars and amino acids react during high dry heat conditions, produces the browning and caramelization present in beef steaks (Dashdorj et al., 2015). However, there remains a lack of information about how particular Maillard products are produced. Recent research by Dinh et al. (2018) has shown the role of sulfur-containing amino acids, certain sugars, and nucleic compounds such as adenosine monophosphate to play a larger role in the generation of beef and brown flavor compounds than was initially thought.

Untargeted metabolomics focuses on the detection of as many groups of metabolites as possible to obtain patterns or fingerprints of biological phenomena (Cevallos-Cevallos et al., 2009). Hundreds or even thousands of metabolites can be measured, utilizing advanced chemistry detection techniques such as high-performance liquid chromatography–quadrupole time of flight (HPLC-qTOF). Statistical analyses are utilized to narrow the number of metabolites from thousands to just those that are statistically impactful on the applied treatments. This smaller number of compounds (100 to 200) may then be analyzed statistically using various multivariate approaches to then identify those influenced by applied treatments. Many pre-harvest factors have been shown to have an influence on metabolites that act as substrates in flavor-producing reactions. Specifically, Arshad et al. (2018) noted that breed variations due to genetic differences in metabolites resulted in over 40 different Maillard reaction products. The role of small sugar molecules, peptide chains, and free amino acids in the development of beef flavor is largely unknown, outside of impacts on basic tastes (Yoo et al., 2020). This leads to the possibility that the regulation of different metabolite concentrations may result in flavor differences across lean source.

In the beef industry, numerous branded programs exist based on differences in animal biological type and meat quality. These claims have led to premium products that are based on high degrees of marbling and a guaranteed tender product. Consequently, consumers would pay premium prices for these products. These claims have not been confirmed by sensory and basic meat science research. However, novel instrumentation and research methods allow us to explore the water-soluble metabolites in the lean portion of meat separate from the lipid portion (Mottram and Edwards, 1983; Huang et al., 2010). Furthermore, Legako et al. (2015) reported that the polar lipid fraction of fatty acids is largely responsible for differences in flavor as the neutral fraction undergoes little change and degradation during the cooking process, thus contributing less to flavor differences. Since the polar lipid fraction is primarily made up of phospholipids found in the muscle cell membranes and only to a minor degree in adipose tissue (Mottram and Edwards, 1983), we propose that formulating ground beef with a common fat source, but creating differential metabolomic profiles from ground beef made from differently marketed lean beef sources, would discriminate differences in flavor contribution of lean aside from adipose tissue. Defining differences in metabolites due to differences in biological type could provide evidence for differences in flavor from fatty acid profiles and up- or down-regulation of muscle metabolism. The objective of this study was to investigate differences in flavor formation in ground beef due only to the lean source by formulating ground beef with a common fat source. We hypothesize that differences in flavor and aroma compounds in the lean portion of ground beef will be the driving force for flavor differences across ground beef patty types.

Materials and Methods

Product procurement and patty formation

Beef inside rounds (n = 9/lean soure) were purchased from a commercial meat distributor that advertised or qualified for 4 different lean source groups: a USDA premium branded program (upper 2/3 USDA Choice as advertised for the Certified Angus Beef Program; HC); a heart-healthy-branded program (from Akaushi cattle containing a high degree of marbling; HEART); a commodity, grass-fed all-natural program (meeting specifications for all-natural in addition to consuming a grass-based diet; NATURAL); and a USDA commodity USDA Select control group (SELECT). All rounds were processed within 7 d of packing plant fabrication and boxing. Because of the leanness of the rounds, fat content was tested for each round and found to be similar across all batches; therefore, each round was then used to make 2 batches (10% and 20% fat) of ground beef, the experimental unit was a batch of ground beef within a round, and each experimental unit was replicated 9 times. Rounds were trimmed of any visible exterior fat and ground individually (with a grinder cleanout in between rounds) using a 12.7 mm plate. Commercial, pre-ground beef fat from carcasses of grain-fed steers less than 30 mo of age was added to the ground rounds and was mixed to get an assigned targeted fat percentage (10% fat or 20% fat). Actual fat percentage was validated using the Foss FoodScan2 Meat (FOSS Global, Hilleroed, Denmark). After each batch of ground beef was formulated to the appropriate fat percentage, the blend was reground using a 4.76 mm plate. Patties were weighed and hand-formed using a handheld hamburger patty press (Oneida Hospitality Group, Lincolnshire, IL) to achieve a 113.5-g patty with a 12.7 mm thickness and 11.43 cm diameter. Six patties were formed per treatment per replication (9 replications x 2 fat contents) and were randomly designated to sensory analysis (3 patties per round per treatment), gas chromatography–mass spectrometry (GC/MS) analysis (1 patty per round per treatment), fatty acid analysis (1 patty per round per treatment), or HPLC analysis (1 patty per round per treatment). Patty paper was placed on either side of each patty, and patties were individually crust frozen at −10°C, vacuum packaged, and stored at −20°C until analysis.

Cookery and trained sensory analysis

The patties were thawed for 24 h in a cooler at 4°C for each day of the evaluation. Patties were cooked on a commercial flat top grill (Star Max 536TGF 91.4 cm Countertop Electric Griddle, Star International Holdings Inc. Company, St. Louis, MO) set at 177°C to an internal temperature of 35°C, flipped, and removed when the internal temperature reached 71°C (AMSA, 2016). Internal temperatures were monitored by iron-constantan thermocouples (Omega Engineering, Stanford, CT) inserted into the geometric center of each patty and monitored using a digital thermometer (Omega Engineering, model HH501BT type T, Stamford, CT) with a type T thermocouple (Omega Engineering, model TMQSS).

After cooking, patties were cut into 6 pie-shaped wedges. Two wedges per sample were served in 59 mL clear, plastic soufflé cups (translucent plastic 2 oz. portion cups, Georgia-Pacific, Asheboro, NC) tested to ensure that they did not impart flavors. Samples were identified with random three-digit codes and served in random order. Samples were cut and served immediately to ensure samples were approximately 37°C upon time of serving. During evaluation, panelists were seated in individual breadbox-style booths separated from the preparation area, and samples were evaluated under red lights (44.2 lux). To prevent taste fatigue, each evaluation day was divided into 2 sessions, with a 10-min break between sessions, and samples were served 4 min apart.

Human subject use was approved by the Texas A&M Human Research Protection Program prior to the study (IRB2019-0596M). Three patties (5 wedges/patty—each panelist was served 2 random wedges) were evaluated by a 5-member, expert trained beef flavor and texture descriptive attribute panel that helped develop and validate the beef lexicon. This panel was retrained using the beef lexicon for 14 d using 16-point intensity scales where 0 = none and 15 = extremely intense flavor (Adhikari et al., 2011). After training was complete, panelists were presented 12 samples per day, divided into 2 sessions. Prior to the start of each trained panel evaluation day, panelists were calibrated using one orientation or “warm up” sample that was evaluated and discussed orally. Double-distilled, deionized water, unsalted saltine crackers, and fat-free ricotta cheese were available for cleansing the palette between samples.

Fatty acid analysis

Total lipids of raw patties were extracted by a modification of the method of Folch et al. (1957). Five grams of homogenized beef was extracted in chloroform:methanol (2:1, v/v), and fatty acid methyl esters (FAME) were prepared as described by Morrison and Smith (1964), modified to include an additional saponification step described by Archibeque et al. (2005). The FAME were analyzed using an FID detector Varian gas chromatograph (model CP-3800 fixed with a CP-8200 auto sampler, Varian Inc., Walnut Creek, CA). Separation of FAME was accomplished on a fused silica capillary column CP-Sil88 (100 m × 0.25 mm [i.d.]) (Chrompack Inc., Middleburg, The Netherlands), with hydrogen as the carrier gas (flow rate = 35 mL/min; injection split ratio 20:1). Initial oven temperature was 150°C; oven temperature increased at 5°C/min to 220°C and were held for 22 min. Injector and detector temperatures were set at 270°C and 300°C, respectively. Individual fatty acids were identified using genuine external standard GLC-68D (Nu-Chek Prep, Inc., Elysian, MN).

Volatile compound analysis—GC/MS

Immediately after cooking for sensory analysis described earlier, leftover patty wedges from the sensory patties were frozen in liquid nitrogen and stored in −80°C until collection of volatile compounds. Frozen patties were powdered in liquid nitrogen, and 5 g of powdered sample was placed in a 20-mL glass vial with a Teflon lid and placed on a heating block (Block analog 2 120V with block modular 28M, VWR) held at 65°C. The volatile compounds present in the headspace were collected using a solid-phase micro-extraction (SPME) Portable Field Sampler (Supelco 504831, 75 μm carboxen/polydimethylsiloxane, Sigma-Aldrich, St. Louis, MO) for 20 min. Volatile aroma compounds were eluted from the SPME and separated using GC (Agilent Technologies 7920 series GC, Santa Clara, CA). The sample was desorbed at 280°C for 3 min. The sample was then loaded onto the gas chromatograph column (Agilent VF 5MS 30 m × 0.25 mm ID/1 μ film thickness, SGE Analytical Sciences, Austin, TX). Through the column, the temperature started at 40°C (held for 1 min) and increased at a rate of 20°C/min until reaching 250°C. Compounds were identified and quantified (using a 1,3-dichlorobenzene internal standard) with a mass spectrometer (Agilent Technologies 5975 series MSD, Santa Clara, CA) using the NIST/Wiley Chemical Library (Palisade, Ithaca, NY). Only compounds with total ion counts of 500 or more (about 10 ng/g) were kept for analyses.

High-performance liquid chromatography–quadrupole time-of-flight

Frozen (powdered in liquid nitrogen in a blender), 2-g samples (both raw and cooked taken from sensory patties) were placed in a 15 mL centrifuge tube containing 8 mL acidified acetonitrile (2.0% formic acid) and then centrifuged for 5 min at 4,000 × g (4°C). Next, 5 mL of supernate was transferred to a new tube with dSPE Enhanced Matrix Removal (Agilent Technologies; hydrated with 5 mL of water). The samples were centrifuged at 4,000 × g (4°C) for 3 min. The supernate was transferred to a new centrifuge tube with 3.5 g MgSO4 and centrifuged at 4,000 × g (4°C) for 3 min. Then 200 μL of the supernate was added to a 2.5 mL sample vial with 800 μL of water. Each sample was run in duplicate using an Agilent 6545 LC/MS-QTOF using a 3.0 x 100 mm, 2.7 μm LC column (35°C; Agilent Poroshell 120 EC-C18, Santa Clara, CA) with a 1.0 μL injection volume, 0.4 mL/min flow rate in positive ion MS mode. The mobile phase consisted of acidified (0.1% formic acid) HPLC-grade water for solvent A and acidified (0.1% formic acid) methanol for solvent B with the following gradient: 0 min 97% A and 3% B; 2 min 97% A and 3% B; 8 min 0% A and 100% B; 14 min 0% A and 100% B; 15 min 97% A and 3% B. The liquid was injected into an Agilent qTOF (model 6545, Agilent Technologies, Santa Clara, CA) using dual electro-spray injection with the drying gas set at 12 L/min at 320°C, the nebulizer at 241 kpa, sheath gas flow at 11 L/min at 350°C. The mass range was 50 to 1400 m/z with a spectral acquisition rate of 5 spectra/s and 1553 transients/spectrum in centroid mode. Quality control samples were created by pooling equivalent volumes of all samples within a treatment, and a blank containing HPLC-grade water was used to subtract background noise. Compounds with an abundance greater than 1.0 were kept for analyses.

Statistical analyses

Trained sensory panel, volatile aroma compound, and fatty acid profile data were analyzed using ANOVA as a 4 (HC, NATURAL, HEART, or SELECT lean source) by 2 (10% or 20% fat inclusion) factorial arrangement of a completely randomized design with 9 replications, with 1 round within each lean source serving as the experimental unit and each experimental unit being subjected to both 10% and 20% fat. Lean source and fat inclusion served as fixed effects, and replication across treatments served as a random effect. Least-squares means of significant (P < 0.05; trends with 0.05 < P < 0.10 were also discussed) F-tests for main and interaction effects were separated using Fisher’s protected LSD (a two-sample t-test) using JMP Pro version 16 (SAS Institute Inc., Cary, NC).

Metabolomic data from the HPLC-qTOF were integrated using Agilent MassHunter Workstation Workflows (version B.08.00, Agilent Technologies, Inc., Santa Clara, CA). The extracted compounds were then imported into Mass Profiler Professional (version 14.9.1, Agilent Technologies, Inc., Santa Clara, CA). That subset of compounds was then filtered for those present in 75% of samples in at least one treatment, with a coefficient of variation within raw data of 25% or less, normality with Shapiro-Wilk cutoff of 0.1, a fold-change of at least 2.0, and significance in an ANOVA (P < 0.05) conducted on log10-normalized data for the main or interaction effects. Compounds that were kept in the model for further analyses were identified by comparing the mass spectrum and retention time to the METLIN compound database (Guijas et al., 2018) using MassHunter ID Browser (version B.08.00, Agilent Technologies, Inc., Santa Clara, CA). Since only MS spectra were collected, we are only reporting relative differences among unidentified, untargeted compounds. Abundance values for compounds were exported to JMP where they were analyzed for hierarchical clustering and discriminate analysis.

Results and Discussion

Trained sensory panel

No interaction effects were found between lean source and fat percentage (P > 0.05; data not shown). Trained sensory green descriptor was higher (P < 0.05) in HC patties than in NATURAL patties, whereas HEART and SELECT patties were similar (P > 0.05) to all lean source treatments (Table 1). Textural hardness scores were higher in HEART and SELECT patties than in NATURAL patties (P < 0.05). Additionally, trained panel fat-like and rancid scores tended (P = 0.075 and 0.081, respectively) to be higher and springiness scores tended (P = 0.069) to be lower for NATURAL patties than in all other lean sources. Trained sensory panel scores for fat-like, burnt, buttery, heated oil, and juiciness were higher (P < 0.05) for patties with 20% fat than in those with 10% fat. Brown/roasted sensory scores tended (P = 0.094) to be higher in patties with 20% fat than in patties containing 10% fat.

Table 1.

Least-squares means and SEM for the main effects of lean source and fat content of trained sensory panel scores of ground beef patties

Lean source main effect Fat percentage main effect
Sensorya HCb HEART NAT SELECT SEM P > F 10% Fat 20% Fat SEM P > F
Beef identity 6.8 6.7 6.7 7.0 0.18 0.34 6.7 6.9 0.13 0.59
Brown/roast 8.8 8.5 8.6 8.7 0.22 0.74 8.5 8.8 0.16 0.093
Bloody 0.8 0.7 0.8 1.1 0.19 0.60 0.8 0.9 0.13 0.49
Fat-like 3.5 3.5 4.0 3.5 0.16 0.075 3.4d 3.9c 0.11 <0.001
Bitter 2.1 2.2 2.0 2.2 0.12 0.69 2.2 2.1 0.09 0.78
Salty 1.3 1.2 1.2 1.4 0.14 0.87 1.2 1.3 0.10 0.49
Sweet 0.3 0.2 0.3 0.1 0.06 0.12 0.2 0.2 0.04 0.65
Sour 2.4 2.0 2.0 2.4 0.15 0.080 2.2 2.2 0.11 0.90
Umami 4.0 3.9 4.0 4.1 0.20 0.94 3.9 4.2 0.14 0.16
Metallic 2.5 2.4 2.4 2.6 0.10 0.37 2.5 2.5 0.07 0.80
Musty earthy 0.8 0.8 0.7 0.6 0.10 0.28 0.8 0.6 0.07 0.13
Cardboardy 1.2 1.0 0.9 1.0 0.17 0.66 1.1 1.0 0.12 0.64
Burnt 0.5 0.4 0.5 0.4 0.15 0.93 0.3d 0.6c 0.11 0.047
Buttery 0.4 0.5 0.5 0.4 0.10 0.69 0.2d 0.6c 0.07 <0.001
Green 0.4c 0.1cd 0.0d 0.2cd 0.08 0.029 0.2 0.2 0.06 0.89
Heated oil 0.8 1.0 1.1 0.8 0.16 0.45 0.7d 1.2c 0.11 0.005
Liver-like 0.5 0.4 0.3 0.4 0.16 0.89 0.4 0.4 0.11 0.84
Painty 0.3c 0.0d 0.1d 0.0d 0.07 0.050 0.1 0.1 0.05 0.34
Rancid 0.0 0.0 0.1 0.0 0.04 0.081 0.1 0.0 0.03 0.31
Refrigerator 0.2 0.1 0.2 0.3 0.07 0.56 0.2 0.2 0.05 0.72
Juiciness 6.9 7.0 7.7 7.3 0.25 0.10 6.9d 7.6c 0.17 0.016
Springiness 5.2 5.5 4.9 5.4 0.17 0.069 5.3 5.2 0.12 0.41
Hardness 5.3cd 5.4c 4.9d 5.6c 0.15 0.012 5.4 5.2 0.11 0.16
Cohesiveness 5.5 5.8 5.7 6.0 0.18 0.36 5.8 5.8 0.12 0.81
Cook loss, % 29.6 28.3 27.5 27.9 1.03 0.52 26.9d 29.8c 0.73 0.007
  • Sensory descriptor scores are on a scale of 0 to 15 with 0 = absent and 15 = extremely intense.

  • HC = high choice lean; HEART = heart-healthy-branded lean; NAT = naturally branded lean; SELECT = USDA Select lean.

  • Means in a row within a main effect category lacking a common superscript differ (P < 0.05).

We hypothesized that we would find differences in flavor and aroma compounds in ground beef patties with a common fat source but differing lean sources. However, the major trained sensory panel attributes (beef identity, brown/roasted, bloody, fat-like, and umami) were not impacted by lean source in ground beef patties regardless of fat content. Legako et al. (2015) indicated that differences in fatty acids, and in particular fatty acids from phospholipids, were responsible for much of the differences in whole muscle strip steaks from differing quality grades. Kerth et al. (2015) reported that when ground beef patties were formulated with a common lean source, but fat sources with different fatty acid profiles, very few differences were found in consumer acceptance scores or volatile aroma compounds. It was concluded that, because a large portion of any differences in sensory attributes would be from the impact of the polar lipids found in the lean, and the lean source was common to all treatments, no differences would be found (Kerth et al., 2015). However, when both lean and fat from sources with divergent fatty acid profiles were used to make ground beef patties (Blackmon et al., 2015), differences in fatty acid profiles as well as some trained sensory panel scores (e.g., fat-like) were improved in patties that had a higher percentage of oleic acid and total monounsaturated fatty acids (MUFA). In this present study, only the green flavor descriptor significantly differed between lean sources. Surprisingly, this descriptor was perceived as higher intensity in the HC patties compared to NATURAL patties. Frank et al. (2016) similarly reported that steaks from Angus cattle finished on grain had higher intensity of grassy flavors compared to Angus cattle finished on grass-fed diets.

Fatty acid profile

Myristic acid (C14:0), palmitic (C16:0) acid, aracidic acid (20:0), eicosenoic acid (C20:1), eicosadienoic acid (C20:2), docosenoic acid (C22:1), and docosahexanoic acid (C22:6n-3) were not significantly impacted by lean source or fat percentage (P > 0.05; Table 2). Furthermore, myristoleic acid (C14:1), vaccenic acid (C18:1n-11), linoleic acid (C18:2n-6), behenic acid (C22:0), other, and total polyunsaturated fatty acids (PUFA) were not impacted by lean source (P > 0.05). Palmitic oleic acid (16:1) was higher (P < 0.05) in HC compared to either NATURAL or SE but was similar (P > 0.05) to HEART. Stearic acid (C18:0) was lowest (P < 0.05) in HEART and highest (P < 0.05) in SELECT compared to both HC and NATURAL, which were similar (P > 0.05). Oleic acid (C18:1n-9) was lowest (P < 0.05) in patties with SELECT lean compared to all other lean source treatments. Alpha-linolenic acid (C18:3) was higher (P < 0.05) in HC and SELECT compared to HEART and NATURAL, while arachidonic acid (C20:4) was higher (P < 0.05) in patties that used HEART lean compared to those that used either HC or SELECT. Erucic acid (C22:1) was highest in patties that had lean sourced from HC or SELECT compared to those from either HEART or NATURAL (P < 0.05). Lignoceric acid (24:0) was highest (P < 0.05) in HEART compared to all other lean source treatments, and nervonic acid (24:1) was highest (P < 0.05) in both HC and SELECT compared to the other lean sources. Patties with SELECT lean had the highest level of total saturated fatty acids (SFA), whereas HEART had the lowest (P < 0.05). Finally, patties with HEART lean had 2.3% to 4.9% more (P < 0.05) total MUFA compared to NATURAL and SELECT lean sources; however, HEART patties had similar (P > 0.05) total MUFA compared to HC.

Table 2.

Least-squares means and SEM for the main effects of lean source and fat content of fatty acid profiles of ground beef patties

Lean source main effect Fat percentage main effect
Fatty acida HCb HEART NAT SELECT SEM P > F 10% Fat 20% Fat SEM P > F
14:0 3.19 2.73 2.89 3.03 0.140 0.13 2.91 3.01 0.099 0.50
14:1n-5 0.53 0.55 0.49 0.51 0.041 0.39 0.58c 0.49d 0.029 0.031
16:0 24.84 24.66 24.51 24.51 0.094 0.054 24.55 24.71 0.068 0.093
16:1 2.37c 2.27cd 1.96de 1.71e 0.140 0.006 2.27c 1.88d 0.099 0.008
18:0 21.96d 19.16e 21.97d 24.52c 0.603 <0.001 20.15d 23.66c 0.408 <0.001
18:1n-9 37.32c 38.88c 36.71c 34.34d 0.784 0.001 37.50 36.13 0.546 0.084
18:1n-11 0.70 0.44 0.71 0.66 0.108 0.26 0.76c 0.49d 0.077 0.016
18:2n-6 3.66 3.41 3.66 3.49 0.212 0.78 3.82c 3.29d 0.150 0.014
18:3 0.29c 0.22d 0.23d 0.28c 0.285 <0.001 0.24 0.26 0.009 0.16
20:0 0.24 0.19 0.23 0.24 0.024 0.42 0.22 0.23 0.017 0.55
20:1 0.22 0.18 0.20 0.19 0.019 0.34 0.21 0.19 0.013 0.36
20:2 0.04 0.05 0.03 0.04 0.006 0.53 0.05 0.04 0.004 0.096
20:4 0.35d 0.47c 0.43cd 0.34d 0.031 0.013 0.52c 0.27d 0.022 <0.001
22:0 0.03 0.03 0.02 0.03 0.004 0.089 0.02d 0.03c 0.003 0.001
22:1 0.04c 0.01d 0.01d 0.04c 0.008 0.001 0.03 0.02 0.005 0.21
24:0 0.04d 0.07c 0.05d 0.04d 0.006 0.003 0.07c 0.03d 0.004 <0.001
24:1 0.09c 0.05d 0.04d 0.08c 0.009 0.001 0.08c 0.05d 0.006 0.001
22:6n-3 0.03 0.02 0.02 0.03 0.004 0.20 0.02d 0.03c 0.003 <0.001
Other 6.60 6.30 6.26 6.25 0.363 0.89 6.66 6.05 0.26 0.10
Total SFA 49.93d 46.83e 49.59d 52.30c 0.577 <0.001 47.77d 51.55c 0.403 <0.001
Total MUFA 41.04cd 42.36c 40.04d 37.49e 0.814 0.001 41.27c 39.20d 0.576 0.013
Total PUFA2 4.37 4.16 4.35 4.18 0.260 0.90 4.64c 3.89d 0.181 0.004
  • Fatty acids are expressed as a percentage of the total fatty acids; fatty acids are defined as carbon chain length:number of double bonds with n-5, n-9, n-11, n-6, and n-3 identifying which carbon contains the double bond; SFA = saturated fatty acids; MUFA = monounsaturated fatty acids; PUFA = polyunsaturated fatty acids.

  • HC = high choice lean; HEART = heart-healthy-branded lean; NAT = naturally branded lean; SELECT = USDA Select lean.

  • Means in a row within a main effect lacking a common superscript differ (P < 0.05).

Patties with 10% fat were 0.09% higher in myristoleic, 0.39% lower in palmitoleic acid, 3.51% lower in stearic acid, 0.49% higher in oleic acid, 0.25% higher in arachidonic acid, 0.01% lower in behenic acid, 0.04% higher in lignoceric acid, 0.03% higher in nervonic acid, and 0.01% lower in docosahexanoic acid compared to their counterparts with 20% fat (P < 0.05). Finally, patties with 10% fat had 3.85% lower total SFA and 2.07% and 0.75% higher levels of total MUFA and PUFA, respectively (P < 0.05).

Turk and Smith (2009) reported that stearic acid decreased and oleic acid was greater in ground beef from either Angus- or Wagyu-branded ground beef than in pasture-fed ground beef. Furthermore, the MUFA:SFA ratio was significantly greater in Angus and Wagyu ground beef compared to not only pasture-fed ground beef but all other commodity and fast-food ground beef. Overall, the addition of the grain-fed fat source resulted in an increase in the total SFA for all lean sources and decrease of the MUFA for all lean sources (Westerling and Hedrick, 1979). This is unfavorable from a human health perspective as unsaturated fatty acids have been shown to have more beneficial aspects compared to saturated fatty acids in food products (Lunn and Theobald, 2006; Margină et al., 2020; Djuricic and Calder, 2021). High-density lipoprotein-C (HDL-C), which is considered the “good” cholesterol, has been shown to be elevated in women (Gilmore et al., 2011, 2013) and in men (Adams et al., 2010) when they consume beef that is high in mono-unsaturated fatty acids. This is particularly true when consuming beef with high levels of oleic acid and as well as beef lower in total, saturated, and trans-fatty acids. Furthermore, in an extensive review, Smith et al. (2020) discussed producing high-oleic acid beef and beef consumption on the risk factors for cardiovascular disease.

Legako et al. (2015) reported that the percentage of MUFA increased as the USDA quality grade increased from standard to low choice to prime but decreased in cooked samples compared to raw samples regardless of quality grade. Furthermore, the percentage of PUFA from phospholipids decreased with an increase in USDA quality grade but significantly increased in cooked samples compared to raw samples, regardless of quality grade. This may indicate that the disappearance of MUFA is from their participation in lipid thermal degradation and their contribution to volatile aroma and other flavor compounds. Alternatively, MUFA could be lost disproportunately as drip loss during cooking.

Since the ground beef batches were formulated with such high percentages of fat compared to the fat inherently present from the lean source in the round, it was expected that the percentages of fatty acids would not be impacted by the lean source. There are two different theories for why differences were found. The first is that the batches formulated from HC and HEART lean contained a high enough concentration of intramuscular fat that they were not impacted as heavily by the additional fat, especially in the lower fat formulation (Lunt et al., 1993). Conversely, the additional fat source was purchased in a large amount and thus contained fat from several cattle. Therefore, the fat source itself could have been so diverse that it still led to differences between the treatment batches. Analyzing the lean sources prior to the addition of the fat and analyzing the fat separately could also help prove or disprove that theory. In future research, creating a baseline fatty acid composition from the lean source prior to formulating could help explain the differences in final composition and could help form a better theory for how fat is impacting flavor.

Volatile aroma compounds

The ng/g of volatile aroma compounds are reported in Table 3 for the main effects of lean source, fat percentage, and their interaction. In evaluating the volatiles affected by the fat% main effect, acetic acid was higher (P = 0.002) in patties with 20% fat than in those with 10% fat, while 1-octen-3-ol (mushroom-like aroma; Kerth and Miller, 2015) was higher in patties with HEART than in those with 20% fat (P = 0.031). In fact, 1-octen-3-ol was only present in 10% fat patties and absent in 20% fat patties. Patties with 20% fat had higher levels of 2-methyl butanal (a malty, fruity arom), heptane, octane (gas-like aroma), and 3-hydroxy-2-butanone (buttery, creamy aroma; P < 0.032) compared to patties with 10% fat. Blackmon et al. (2015) found that as the percentage fat in ground beef patties increases, the concentration of volatiles derived from the Maillard reaction decreases. This is likely because phospholipids and their degradation products inhibit important reactions involved in the formation of heterocyclic aroma compounds in the Maillard reaction.

Table 3.

Least-squares means and SEM for the main effects of lean source and fat content of volatile aroma compounds (ng/g) of ground beef patties

P > F
10% Fat 20% Fat Fat% x
Compound, ng/g HCa HEART NAT SELECT HC HEART NAT SELECT SEM Fat% Lean Lean LRI
Acids
 Acetic acid 0.21 0.00 0.03 0.25 0.85 0.31 0.47 0.65 0.218 0.002 0.23 0.88 579
Alcohols
 1-octen-3-ol 0.08 0.27 0.04 0.07 0.00 0.00 0.00 0.00 0.082 0.031 0.41 0.41 983
 1-pentanol 1.40c 0.97c 2.58c 0.78c 1.44c 1.37c 9.66b 0.79c 1.470 0.012 0.001 0.001 769
 1-penten-3-ol 0.26b 0.00c 0.00c 0.00c 0.01c 0.00c 0.04c 0.00c 0.048 0.11 0.017 0.008 685
Aldehydes
 2-heptenal 0.00k 0.01jk 0.02j 0.00k 1.04k 0.05jk 0.18j 0.00k 0.044 0.086 0.017 0.087 964
 2-methyl-butanal 4.95 5.32 0.67 3.33 7.21 10.53 11.14 5.03 2.192 0.001 0.37 0.060 665
 2-methyl propanal 1.58cd 2.30cd 0.25d 1.35cd 3.20bc 2.15cd 5.17b 2.29cd 1.013 0.007 0.79 0.021 556
 3-methyl-butanal 6.79cd 7.26cd 1.27d 3.80cd 10.17c 10.94bc 19.28b 7.52cd 3.294 0.001 0.44 0.012 658
 Acetaldehyde 0.81c 1.59c 0.35c 0.97c 2.09c 1.99c 4.87b 1.24c 0.994 0.012 0.29 0.028 450
 Benzaldehyde 1.17j 0.70jk 0.09k 0.78j 0.60j 0.97jk 0.47k 0.46j 0.291 0.74 0.048 0.18 978
 Heptanal 6.46c 4.60cd 1.36d 3.66cd 6.77c 5.43cd 15.57b 6.58c 2.228 0.002 0.20 0.001 963
 Hexanal 104.62c 77.07c 34.75c 53.36c 107.70c 93.83c 319.13b 48.09c 56.354 0.040 0.052 0.003 704
 Octanal 2.78c 2.27cd 0.64d 1.95cd 3.54c 3.00c 6.95b 3.15c 0.886 0.001 0.30 0.001 1008
 Pentanal 7.99bc 4.00c 2.49c 3.88c 1.66c 3.71c 15.78b 0.00c 3.621 0.76 0.10 0.005 702
Hydrocarbons
 4-methyl octane 0.89bc 0.35cd 0.02d 1.21b 0.37cd 0.25d 0.11d 0.23d 0.22 0.01 0.002 0.035 864
 Decane 1.35bc 0.74c 0.13d 1.55b 1.00bc 0.77c 1.02bc 0.81c 0.241 0.78 0.006 0.001 999
 Heptane 1.13 1.06 0.20 1.16 5.51 4.54 5.58 5.82 1.225 0.001 0.91 0.83 701
 Nonane 0.30c 0.29c 0.06d 0.42bc 0.69b 0.24cd 0.44bc 0.45bc 0.107 0.006 0.021 0.035 901
 Octane 2.71e 3.01de 0.70f 2.90de 6.22c 4.92cde 9.57b 5.54cd 0.978 0.001 0.49 0.001 864
 Pentane 1.00 6.15 0.53 0.00 8.53 6.16 20.86 4.24 4.736 0.009 0.17 0.053 503
Ketones
 2-butanone 0.05k 0.00k 0.14k 3.46j 0.00k 2.61k 0.00k 6.61j 1.750 0.26 0.009 0.55 596
 2-heptanone 0.37jk 0.57j 0.13k 0.18k 0.40jk 0.44j 0.25k 0.18k 0.143 0.96 0.031 0.79 593
 2-propanone 2.23c 0.77c 1.40c 2.08c 3.79c 1.22c 17.72b 3.03c 2.621 0.005 0.001 0.001 500
 2,3-butanedione 2.67cd 8.76b 0.30d 6.36bc 4.24bcd 6.35bc 7.68bc 5.03bc 2.212 0.36 0.17 0.034 587
 2,3-pentanedione 0.14j 0.12j 0.01k 0.13j 0.06j 0.08j 0.00k 0.23j 0.045 0.80 0.001 0.16 699
 2,5-octanedione 1.84k 2.95j 0.24k 0.62k 0.10k 3.56j 1.75k 1.44k 0.769 0.54 0.003 0.097 985
 3-hydroxy-2-butanone 16.57k 40.26j 1.94k 21.67jk 24.01k 56.19j 37.13k 39.56jk 10.665 0.006 0.011 0.43 716
 5-methyl-2-hexanone 0.08bc 0.00c 0.03c 0.05c 0.00c 0.00c 0.27b 0.05c 0.076 0.41 0.11 0.049 893
Pyrazines
 2,5-dimethyl-pyrazine 1.06k 1.09j 0.09k 0.95jk 0.94k 3.37j 0.94k 1.10jk 0.638 0.056 0.020 0.20 924
 3-ethyl-2,5-dimethyl-pyrazine 2.20c 2.30c 0.59d 2.41c 3.14c 2.84c 5.28b 3.46bc 0.741 0.001 0.92 0.002 1087
 Methyl-pyrazine 0.24jk 0.36j 0.04k 0.25j 0.33jk 0.35j 0.02k 0.58j 0.153 0.32 0.014 0.59 834
Sulfur
 Carbon disulfide 0.70b 0.17cde 0.06e 0.40bcd 0.20cde 0.09de 0.36bcd 0.49bc 0.135 0.58 0.019 0.008 541
 Thiobis-methane 0.11k 0.08k 0.02k 0.33j 0.05k 0.00k 0.00k 0.0j 0.077 0.13 0.002 0.86 521
Terpenoids
 1-octene 0.17cd 0.25cd 0.05d 0.03d 0.13cd 0.44bc 0.74b 0.19cd 0.140 0.006 0.051 0.010 791
 3-carene 0.21cd 0.09cd 0.00d 0.00d 0.05c 0.43bc 0.70b 0.05d 0.131 0.007 0.023 0.001 1020
 Alpha-pinene 0.73de 0.80d 0.16e 0.63de 1.04cd 1.60c 2.68b 1.03cd 0.285 0.001 0.050 0.001 944
 DL-limonene 0.17cd 0.34c 0.08d 0.16cd 0.28c 0.74b 0.78b 0.29c 0.088 0.001 0.001 0.001 1040
 Methyl benzene 0.00c 0.27c 0.05c 0.98b 0.00c 0.40c 0.27c 0.28c 0.195 0.47 0.007 0.045 773
 Toluene 1.09j 1.33j 0.24k 0.30jk 1.49j 1.00j 0.29k 0.90jk 0.439 0.52 0.020 0.69 772
  • 10% fat = calculated 10% fat and 90% lean in the ground beef; 20% fat = calculated 20% fat and 80% lean in the ground beef; HC = high choice lean; HEART = heart-healthy-branded lean; NAT = naturally branded lean; SELECT = USDA Select lean; LRI = linear retention index (Raza et al., 2019).

  • Interaction means in a row lacking a common superscript differ (P < 0.05).

  • Lean source main effect means in a row lacking a common superscript differ (P < 0.05).

Patties made from NATURAL lean had slightly higher (P = 0.017) 2-heptenal compared to patties made with HC or SELECT but lower (P = 0.048) benzaldehyde (nutty, woody aroma) compared to patties made with either HC or SELECT lean. The fruity, green aroma (Kerth and Miller, 2015) of 2-butanone was by far higher in patties made in the present study with SELECT lean compared to all the other lean sources, but 2-heptanone (banana aroma) was highest for patties made with HC lean compared to those made with either NATURAL or SELECT (P < 0.05). The ketone 2,3-pentanedione was lowest (P < 0.05) for patties made with NATURAL lean compared to the other lean sources, which did not differ (P > 0.05). Volatiles 2,5-octanedione, 3-hydroxy-2-butanone (buttery, creamy aroma), 2,5-dimethyl-pyrazine (musty potato aroma), and methyl-pyrazine were highest in patties that were made with HEART lean compared to the other lean sources (P < 0.05), although volatile levels for the latter three were similar (P < 0.05) to patties made with SELECT. The fishy aroma thiobis-methane was highest in patties made with SELECT lean compared to the other three sources, and toluene was lower in patties made with NATURAL lean compared to either HC or HEART (P < 0.05).

1-pentanol was significantly higher in NATURAL patties with 20% fat, while 1-penten-3-ol (grassy/green aroma; Burdock, 2009) was higher in HC 10% patties compared to all other fat% and lean source combinations (P < 0.05). All 7 of the aldehyde volatiles significant for a fat% by lean source interaction were higher in patties made with 20% fat and the NATURAL lean source (P < 0.03). However, patties made with 20% fat and HEART lean were similar (P > 0.05) to the 20% NATURAL patties in 3-methyl-butanal (malty), and both heptanal (fatty, harsh aroma) and octanal (citrus-like) volatiles were lower (P < 0.05) in 10% fat NATURAL patties compared to 10% fat HC patties. The volatile 4-methyl-octane was higher (P < 0.05) in patties with 10% fat and SELECT lean except 10% HC patties, which were similar (P > 0.05).

Dimethyl disulfide (onion; Burdock, 2009) was higher in HC and HEART than in NATURAL patties (P < 0.05). Patties with 10% fat and SELECT lean along with HC and NATURAL patties with 20% fat are higher in the volatile decane compared to 10% HEART, 20% HC, and 20% SELECT patties, and 10% NATURAL was the lowest compared to all other treatment combinations (P < 0.05). Lean source did not impact (P > 0.05) nonane within the patties made with 20% fat, but patties made from NATURAL lean was the lowest (P < 0.05) within those patties with 10% fat. Interestingly, octane is lowest for patties made with NATURAL lean within those with 10% fat and is the highest for patties made with NATURAL lean within those patties with 20% fat. The volatile 2-propanone (acetone) was unaffected by lean source within the 10% fat patties, but patties made with NATURAL lean with 20% fat was 13.93% to 17.00% higher than all other treatment combinations (P < 0.05). The buttery-smelling volatile 2,3-butanedione was unaffected (P > 0.05) by lean source within patties with 20% fat, but within patties with 10% fat those with HEART lean were higher and NATURAL lean were lower (P < 0.05). Patties with 20% fat and NATURAL lean had the highest (P < 0.05) 5-methyl-2-hexanone compared to others, except for patties with 10% fat and HC lean (P > 0.05). The caramel/coffee-smelling volatile of 3-ethyl-2,5-dimethyl-pyrazine was highest for patties made with NATURAL lean and 20% fat and lowest for those patties with NATURAL lean made with 10% fat (P < 0.05). Patties with 10% fat had higher (P < 0.05) carbon disulfide volatile in patties made with either HC or SELECT, while patties with 20% fat had higher (P < 0.05) carbon disulfide in patties with SELECT lean compared to those with HEART lean. The terpenoids 1-octene and 3-carene had similar (P > 0.05) concentrations for all lean sources within the patties made with 10% fat, while patties made with 20% fat had higher (P < 0.05) levels in patties made with NATURAL lean compared to those made with HC or SELECT. Alpha-pinene was similar (P > 0.05) among lean sources in 10% fat patties, but in patties with 20% fat and NATURAL lean were higher (P < 0.05) than the other lean sources. The lemon-like DL-limonene volatile was lowest in the NATURAL lean patties with 10% fat, while patties with HEART and NATURAL lean within 20% fat had higher levels compared to patties made with either HC or SELECT (P < 0.05). Finally, patties made with HEART lean and 10% fat had the highest levels of methyl-benzene compared to all other treatment combinations (P < 0.05).

It is not surprising that the NATURAL patties exhibited more differences in volatile aroma compounds, as this was the only treatment that claimed on the label to be finished on a grass-based diet. It is assumed that cattle for all other treatments were finished on a concentrate diet. Additionally, 2-heptenal could also be a contributing factor to the grassy flavor as it has been described as herbaceous. Dimethyl disulfide was found in greater concentrations in patties made with lean from higher USDA quality grades. Gardner and Legako (2018) reported the highest concentration of dimethyl disulfide in USDA Prime steaks and decreasing concentration with decreasing quality grades. Therefore, dimethyl disulfide could be partially responsible for increase in beefy, brown/roasted flavors that drive consumer liking of premium branded projects. Mottram (1998) indicated that sulfur-containing compounds formed in the Maillard reaction from sulfur-containing amino acids like cysteine and reducing sugars like ribose (from inosine monophosphate and other ribonucleotides) seem to be particularly important for the characteristic aroma of meat. Elmore and Mottram (2006) reported greater concentrations of 2-octene, an isomer of 4-octene in cattle finished on a silage diet. Additionally, they reported that cattle fed on concentrate diets had greater concentrations of 1-octen-3-ol and pentanal compared to silage finished cattle.

Mottram (1998) did a thorough job of detailing the formation of flavor in meat, noting the water-soluble products from the Maillard reaction as well as the lipid-derived volatiles that are a result of lipid thermal degradation. He also noted that these two categories of volatile compounds do not act independently of each other and that, in particular, lipid oxidation products like aldehydes and carbonyls react readily with Maillard-derived intermediates. This gives rise to additional aroma compounds and may limit those compounds normally seen in the Maillard reaction. The current research does not necessarily follow this theory as many of the lipid-derived compounds like aldehydes, alcohols, and acids were found in higher concentrations in those patties that had 20% fat, and many of the Maillard products like the pyrazines increased in the patties with the higher fat percentage, but then some products like the sulfur-containing compounds and terpenoids were lower.

Numerous investigators have reported on the resulting differences in flavors and volatiles resulting from beef fed different diets (Melton, 1990, 1999; Maruri and Larick, 1992; Cox et al., 2006; Kerth et al., 2007). In particular, Maruri and Larick (1992) found that terpenoids and long-chain hydrocarbons were moderately and positively correlated with a gamey/stale off-flavor found in ground beef made from cattle that have been finished on grass. Brown et al. (1979) reported that ground beef from steers fed on forage had higher amounts of free fatty acids and that steers fed grain had higher free sugar content compared to their forage-fed counterparts. In general, cattle that are finished on grain diets induce higher levels of branched-chain fatty acids, some aldehydes, and lactones, whereas cattle finished on grass tend to contain higher levels of phenols, terpenes, indoles, and sulphur compounds (Vasta and Priolo, 2006). As possible precursors, Koutsidis et al. (2008) found that animals fed grass had higher levels of free amino acids, whereas those fed concentrates had a higher total reducing sugar content.

HPLC—Metabolomics

Zhang et al. (2021) described the workflow of metabolomics analyses that included biological questions, sample preparation, metabolomics analyses, and data interpretation. In interpreting data, the two major steps are data preprocessing or pretreatment and biological interpretation. Due to the complexity of the data as well as the relatively large data sets, multivariate analyses have been the most chosen statistical pathway. This includes partial least squares discriminate analysis, principal component analysis, and machine learning. We will report hierarchical clustering and discriminate analyses.

Ground beef patties were analyzed using HPLC-qTOF for untargeted metabolomics in both the raw (Figures 1 to 3) as well as the fully cooked (Figures 4 to 6) state. A total of 64 molecules were significantly affected (P < 0.05) by treatment combinations in raw ground beef samples. A two-way hierarchical cluster analysis for these metabolites is illustrated in Figure 1. Different colored bars represent differences in the log10 relative abundance of each metabolite across each treatment combination. The raw cluster is the presence of 4 different compounds that are likely carnitines which have been correlated to consumer acceptance, tenderness, and juiciness (Antonelo et al., 2020a, 2020b). Hierarchacal clusters using the Ward method were identified both across treatment combinations but also across metabolite. Whereas the patties with either 10% or 20% fat clustered together across HEART, HC, and NATURAL lean sources, both the heat map as well as the constellation plot (Figure 2) shows that SELECT patties with 20% fat clustered with the remaining 10% fat lean sources and vice versa for the SELECT patties with 10% fat. Using partial least squares discriminate analyses (Figure 3), treatment combination means with 95% confidence interval ellipses and 50% contours showed that only the lean source main effect accurately clustered metabolites into the 4 lean source categories.

Figure 1.
Figure 1.

Two-way hierarchical cluster analysis of 64 raw small molecule metabolites. Each color bar represents the log10 relative abundance of each annotated metabolite. Those molecules not identified are expressed as mass:charge (m/z) @ the column retention time. 10 = 10% fat ground beef patties, 20 = 20% fat ground beef patties, HC = high choice lean; HEART = heart-healthy-branded lean; NAT = naturally branded lean; SELECT = USDA Select lean.

Figure 2.
Figure 2.

Constellation plot for 64 raw small-molecule metabolites. 10 = 10% fat ground beef patties, 20 = 20% fat ground beef patties, HC = high choice lean; HEART = heart-healthy-branded lean; NAT = naturally branded lean; SELECT = USDA Select lean.

Figure 3.
Figure 3.

Partial least squares discriminative analysis of 64 raw small-molecule metabolites sampled from ground beef with four different sources of lean with 95% confidence ellipses and 50% countours around each treatment mean. HC = high choice lean; HEART = heart-healthy-branded lean; NAT = naturally branded lean; SELECT = USDA Select lean.

Figure 4.
Figure 4.

Two-way hierarchical cluster analysis of 138 cooked small-molecule metabolites. Each color bar represents the log10 relative abundance of each annotated or unidentified metabolite. Those molecules not identified are expressed as mass:charge (m/z) @ the column retention time. 10 = 10% fat ground beef patties, 20 = 20% fat ground beef patties, HC = high choice lean; HEART = heart-healthy-branded lean; NAT = naturally branded lean; SELECT = USDA Select lean.

Figure 5.
Figure 5.

Constellation plot for 138 cooked small-molecule metabolites. 10 = 10% fat ground beef patties, 20 = 20% fat ground beef patties, HC = high choice lean; HEART = heart-healthy-branded lean; NAT = naturally branded lean; SELECT = USDA Select lean.

Figure 6.
Figure 6.

Partial least squares discriminative analysis of 64 raw small molecule metabolites sampled from ground beef with four different sources of lean with 95% confidence ellipses and 50% countours around each treatment mean. 10 = 10% fat ground beef patties, 20 = 20% fat ground beef patties, HC = high choice lean; HEART = heart-healthy-branded lean; NAT = naturally branded lean; SELECT = USDA Select lean.

The cooked patty two-way hierarchical cluster analyses (Figure 4) showed 3 main treatment clusters. Again, patties with both 10% and 20% fat in the SELECT lean source clustered together and apart from the other treatment combinations. The remaining lean sources with 10% fat and those with 20% fat each clustered together. Furthermore, the metabolites (as listed in Table 4) appeared to cluster into 3 primary clusters across the treatment combinations. These 3 treatment combination clusters are also illustrated in the constellation plot found in Figure 5. Interestingly, when looking at the fatty acid profile (Table 2), it appears that the percentage of oleic acid parallels the metabolites from the cooked patties. Perhaps many of the metabolites may be regulating the proportions of fatty acids in the ground beef. Several investigators (Kim et al., 2018; Setyabrata et al., 2021, 2022) have used metabolomics to describe the flavor precursors and flavor changes as a result of aging, and Jeong et al. (2020) even described the impact of specific carnitines present in highly marbled beef that had high sensory scores. The discriminate analysis showed very tight clustering of treatments for cooked ground beef patty small-molecule metabolites across all fat content by lean source interaction means.

Table 4.

List of likely annotated small metabolite compounds or unknown (m/z @ retention time) from raw and cooked ground beef samples. Each compound that was found matches a colored bar across all fat percentage and lean source treatments representing the log10 relative abundance of each compound in the heat map found in Figures 1 and 4.

Raw metabolomic compound n = 64
Unknown 25
Metabolites 18
Amino acids/peptides 5
Carnitines 5
Organic acids 7
Fatty acids/derivatives 4
Cooked metabolomic compound n
Unknown 89
Metabolites 22
Amino acids/peptides 10
Carnitines 4
Organic acids 4
Fatty acids/derivatives 3
Heteroaromatic compound 3
Hormones 3

Researchers have shown that beef high in MUFA, and especially oleic acid, is correlated to positive trained sensory and consumer acceptance scores (Blackmon et al., 2015; Kerth et al., 2015; Listrat et al., 2020). There is some discussion that the source of these positive fatty acids may be sourced in the polar phospholipids of the lean cell membrane and to a lesser extent in the lipids from adipose tissue (Legako et al., 2015; Hunt et al., 2016). This would support the interaction of lean source and fat content found in the present study as fatty acid differences were found among lean sources in the 10% fat patties, but it seems that the added fat in the 20% fat patties diluted the MUFA and palmitoleic acid and oleic acid in particular. It should be noted that the MUFA content of the SELECT 10% patties was similarly low as that found in all of the 20% fat patties. When evaluating the small-molecule metabolite clustering in both raw (Figure 1) and cooked (Figure 4), HEART, NATURAL, and HC all three cluster together in the 10% fat patties. Furthermore, in the cooked heat map, 10% and 20% fat SELECT ground beef patties cluster together and generally have opposite reactions in the metabolites compared to all other treatment combinations.

Jiang and Bratcher (2016) reported 22 metabolic compounds that were significant in a principal components analysis of metabolites from ground beef sourced from different lean sources. In their analysis, they demonstrated that grass-fed beef differentiated from natural beef and that oleic acid was identified as one of the important compounds influencing beef flavor. Ma et al. (2017) found that acyl carnitines were dominant in the psoas major muscle compared to the longissimus lumborum or semimembranosus and in muscles that were aged 9 d compared to those aged 16 or 23 d. They also reported that L-carnitine was negatively correlated to meat color traits during display. They hypothesized that higher levels of carnitine may indicate a higher level of energy metabolism and mitochondrial enzyme activity in antemortem muscles. This condition may cause an increase in the oxygen consumption rate of aged muscle and eliminate free oxygen penetrating the muscle and reacting with myoglobin, resulting in a dark purplish color from the presence of deoxymyoglobin.

Cho et al. (2017) demonstrated a positive relationship between short-chain fatty acid-carnitine (e.g., proprionyl carnitine) levels and insulin resistance. Other research (Koves et al., 2005) indicated that rats that were fed high-fat diets throughout the study had consistently higher levels of acylcarnitine levels, particularly proprionyl, butanoyl, and hexanoyl carnitine, whether fasted or fed, indicating the inability of the muscle to metabolize short-chain fatty acids and amino acids with beta-oxidation. In a study comparing Korean cattle with low marbling (11.6% lipid) and high marbling (16.9% lipid), NMR analysis showed that carnitine along with other metabolites like creatine, lactate, and carnosine increased in longissimus muscle from high marbled beef (Jeong et al., 2020). Conversely, Mateescu et al. (2012) reported that when they tested over 2,000 head of Angus cattle, not only did carnitine not correlate with any sensory panel descriptive ratings, but carnitine was lowly heritable (h = 0.015). It is worth considering in the latter case whether the low heritability coefficient is due to a relatively homogeneous group of cattle within one breed. It would be interesting to see whether heritability would change with an increase in the variability of genetics.

Conclusions

While the lean source and fat content of the ground beef patties used in this study had little impact on trained sensory panel scores, several differences in fatty acid profiles, volatile aroma compounds, and non-volatile flavor precursors and flavor compounds were found. It was particularly interesting to see that both raw and cooked samples analyzed by metabolic profiling had compounds that contributed to differences in ground beef lean and fat content treatments. Certainly, metabolic factors and(or) small molecules in the raw and cooked patties may be those that impact traditional meat quality parameters, and additional work to positively annotate those compounds is necessary. It appears that with these metabolomic findings, future research should focus on correlating these compounds with sensory and flavor traits. The ultimate goal, in considering the entire production system from growth to consumption, would be able to track these metabolites further back in the production of beef that would allow us to be able to determine future beef quality. This could include sampling serum from cattle destined for harvest to determine quality and sensory traits.

Literature Cited

AMSA. 2016. Research guidelines for cookery, sensory evaluation, and instrumental tenderness measurements of meat. 2nd ed. American Meat Science Association, Champaign, IL.

Adams, T. H., R. L. Walzem, D. R. Smith, S. Tseng, and S. B. Smith. 2010. Hamburger high in total, saturated and trans-fatty acids decreases HDL cholesterol and LDL particle diameter, and increases TAG, in mildly hypercholesterolaemic men. Brit. J. Nutr. 103:91–98. doi: https://doi.org/10.1017/S0007114509991516.

Adhikari, K., E. Chambers IV, R. Miller, L. Vazquez-Araujo, N. Bhumiratana, and C. Philip. 2011. Development of a lexicon for beef flavor in intact muscle. J. Sens. Stud. 26:413–420. doi: https://doi.org/10.1111/j.1745-459X.2011.00356.x.

Antonelo, D. S., N. R. Cônsolo, J. F. Gómez, M. Beline, R. S. Goulart, R. R. P. d. S. Corte, L. A. Colnago, M. W. Schilling, D. E. Gerrard, and S. L. Silva. 2020a. Metabolite profile and consumer sensory acceptability of meat from lean Nellore and Angus× Nellore crossbreed cattle fed soybean oil. Food Res. Int. 132:109056. doi: https://doi.org/10.1016/j.foodres.2020.109056.

Antonelo, D. S., J. F. Gómez, N. R. Cônsolo, M. Beline, L. A. Colnago, W. Schilling, X. Zhang, S. P. Suman, D. E. Gerrard, and J. C. Balieiro. 2020b. Metabolites and metabolic pathways correlated with beef tenderness. Meat Muscle Biol. 4. doi: https://doi.org/10.22175/mmb.10854.

Archibeque, S., D. Lunt, C. Gilbert, R. Tume, and S. Smith. 2005. Fatty acid indices of stearoyl-CoA desaturase do not reflect actual stearoyl-CoA desaturase enzyme activities in adipose tissues of beef steers finished with corn-, flaxseed-, or sorghum-based diets. J. Anim. Sci. 83:1153–1166. doi: https://doi.org/10.2527/2005.8351153x.

Arshad, M. S., M. Sohaib, R. S. Ahmad, M. T. Nadeem, A. Imran, M. U. Arshad, J.-H. Kwon, and Z. Amjad. 2018. Ruminant meat flavor influenced by different factors with special reference to fatty acids. Lipids Health Dis. 17:1–13. doi: https://doi.org/10.1186/s12944-018-0860-z.

Blackmon, T., R. K. Miller, C. Kerth, and S. B. Smith. 2015. Ground beef patties prepared from brisket, flank and plate have unique fatty acid and sensory characteristics. Meat Sci. 103:46–53. doi: https://doi.org/10.1016/j.meatsci.2015.01.004.

Brown, H., S. Melton, M. Riemann, and W. Backus. 1979. Effects of energy intake and feed source on chemical changes and flavor of ground beef during frozen storage. J. Anim. Sci. 48:338–347. doi: https://doi.org/10.2527/jas1979.482338x.

Burdock, G. A. 2009. Fernoli’s handbook of flavor ingredients. 6th ed. CRC Press/Taylor and Francis, Boca Raton, FL. doi: https://doi.org/10.1201/9781439847503.

Cevallos-Cevallos, J. M., J. I. Reyes-De-Corcuera, E. Etxeberria, M. D. Danyluk, and G. E. Rodrick. 2009. Metabolomic analysis in food science: A review. Trends Food Sci. Tech. 20:557–566. doi: https://doi.org/10.1016/j.tifs.2009.07.002.

Cho, J., Y. Zhang, S.-Y. Park, A.-M. Joseph, C. Han, H.-J. Park, S. Kalavalapalli, S.-K. Chun, D. Morgan, and J.-S. Kim. 2017. Mitochondrial ATP transporter depletion protects mice against liver steatosis and insulin resistance. Nat. Commun. 8:1–12. doi: https://doi.org/10.1038/ncomms14477.

Corbin, C., M. Miller, T. O’Quinn, T. Dinh, J. Legako, A. Garmyn, M. Hunt, and C. Brooks. 2014. Consumer assessment and fatty acid analysis of beef strip steaks of similar tenderness with varying marbling levels. Meat Sci. 1:469. doi: https://doi.org/10.1016%2Fj.meatsci.2013.07.100.

Cox, R. B., C. R. Kerth, J. G. Gentry, J. W. Prevatt, K. W. Braden, and W. R. Jones. 2006. Determining acceptance of domestic forage-or grain-finished beef by consumers from three southeastern US states. J. Food Sci. 71:S542–S546. doi: https://doi.org/10.1111/j.1750-3841.2006.00124.x.

Dashdorj, D., T. Amna, and I. Hwang. 2015. Influence of specific taste-active components on meat flavor as affected by intrinsic and extrinsic factors: An overview. Eur. Food Res. Technol. 241:157–171. doi: https://doi.org/10.1007/s00217-015-2449-3.

Dinh, T., J. Legako, M. Miller, and J. Brooks. 2018. Effects of USDA quality grade and cooking on water-soluble precursors of beef flavor. Meat Sci. 146:122–130. doi: https://doi.org/10.1016/j.meatsci.2018.08.008.

Djuricic, I., and P. C. Calder. 2021. Beneficial outcomes of omega-6 and omega-3 polyunsaturated fatty acids on human health: An update for 2021. Nutrients 13:2421. doi: https://doi.org/10.3390/nu13072421.

Elmore, J. S., and D. S. Mottram. 2006. The role of lipid in the flavour of cooked beef. Dev. Food Sci. 43:375–378.

Folch, J., M. Lees, and G. H. Sloane Stanley. 1957. A simple method for the isolation and purification of total lipids from animal tissues. J. Biol. Chem. 226:497–509.

Frank, D., A. Ball, J. Hughes, R. Krishnamurthy, U. Piyasiri, J. Stark, P. Watkins, and R. Warner. 2016. Sensory and flavor chemistry characteristics of Australian beef: influence of intramuscular fat, feed, and breed. J. Agr. Food Chem. 64:4299–4311. doi: https://doi.org/10.1021/acs.jafc.6b00160.

Gardner, K., and J. F. Legako. 2018. Volatile flavor compounds vary by beef product type and degree of doneness. J. Anim. Sci. 96:4238–4250. doi: https://doi.org/10.1093/jas/sky287.

Gilmore, L. A., S. F. Crouse, A. Carbuhn, J. Klooster, J. A. E. Calles, T. Meade, and S. B. Smith. 2013. Exercise attenuates the increase in plasma monounsaturated fatty acids and high-density lipoprotein cholesterol but not high-density lipoprotein 2b cholesterol caused by high-oleic ground beef in women. Nutr. Res. 33:1003–1011. doi: https://doi.org/10.1016/j.nutres.2013.09.003.

Gilmore, L. A., R. L. Walzem, S. F. Crouse, D. R. Smith, T. H. Adams, V. Vaidyanathan, X. Cao, and S. B. Smith. 2011. Consumption of high-oleic acid ground beef increases HDL-cholesterol concentration but both high-and low-oleic acid ground beef decrease HDL particle diameter in normocholesterolemic men. J. Nutr. 141:1188–1194. doi: https://doi.org/10.3945/jn.110.136085.

Guijas, C., J. R. Montenegro-Burke, X. Domingo-Almenara, A. Palermo, B. Warth, G. Hermann, G. Koellensperger, T. Huan, W. Uritboonthai, and A. E. Aisporna. 2018. METLIN: A technology platform for identifying knowns and unknowns. Anal. Chem. 90:3156–3164. doi: https://doi.org/10.1021/acs.analchem.7b04424.

Huang, Y.-C., H.-J. Li, Z.-F. He, T. Wang, and G. Qin. 2010. Study on the flavor contribution of phospholipids and triglycerides to pork. Food Sci. Biotechnol. 19:1267–1276. doi: https://doi.org/10.1007/s10068-010-0181-0.

Hunt, M., J. Legako, T. Dinh, A. Garmyn, T. O’quinn, C. Corbin, R. Rathmann, J. Brooks, and M. Miller. 2016. Assessment of volatile compounds, neutral and polar lipid fatty acids of four beef muscles from USDA Choice and Select graded carcasses and their relationships with consumer palatability scores and intramuscular fat content. Meat Sci. 116:91–101. doi: https://doi.org/10.1016/j.meatsci.2016.02.010.

Jeong, J. Y., M. Kim, S.-Y. Ji, Y.-C. Baek, S. Lee, Y. K. Oh, K. E. Reddy, H.-W. Seo, S. Cho, and H.-J. Lee. 2020. Metabolomics analysis of the beef samples with different meat qualities and tastes. Food Science of Animal Resources 40:924–937. doi: https://doi.org/10.5851/kosfa.2020.e59.

Jiang, T., and C. L. Bratcher. 2016. Differentiation of commercial ground beef products and correlation between metabolites and sensory attributes: A metabolomic approach. Food Res. Int. 90:298–306. doi: https://doi.org/10.1016/j.foodres.2016.11.002.

Kerth, C., K. Braden, R. Cox, L. Kerth, and D. Rankins Jr. 2007. Carcass, sensory, fat color, and consumer acceptance characteristics of Angus-cross steers finished on ryegrass (Lolium multiflorum) forage or on a high-concentrate diet. Meat Sci. 75:324–331. doi: https://doi.org/10.1016/j.meatsci.2006.07.019.

Kerth, C. R., A. L. Harbison, S. B. Smith, and R. K. Miller. 2015. Consumer sensory evaluation, fatty acid composition, and shelf-life of ground beef with subcutaneous fat trimmings from different carcass locations. Meat Sci. 104:30–36. doi: https://doi.org/10.1016/j.meatsci.2015.01.014.

Kerth, C. R., and R. K. Miller. 2015. Beef flavor: A review from chemistry to consumer. J. Sci. Food Agr. 95:2783–2798. doi: https://doi.org/10.1002/jsfa.7204.

Killinger, K. M., C. R. Calkins, W. J. Umberger, D. M. Feuz, and K. M. Eskridge. 2004. Consumer sensory acceptance and value for beef steaks of similar tenderness, but differing in marbling level. J. Anim. Sci. 82:3294–3301. doi: https://doi.org/10.2527/2004.82113294x.

Kim, Y. H. B., D. Ma, D. Setyabrata, M. M. Farouk, S. M. Lonergan, E. Huff-Lonergan, and M. C. Hunt. 2018. Understanding postmortem biochemical processes and post-harvest aging factors to develop novel smart-aging strategies. Meat Sci. 144:74–90. doi: https://doi.org/10.1016/j.meatsci.2018.04.031.

Koutsidis, G., J. Elmore, M. J. Oruna-Concha, M. M. Campo, J. D. Wood, and D. Mottram. 2008. Water-soluble precursors of beef flavour: I. Effect of diet and breed. Meat Sci. 79:124–130.  http://doi.org/doi:/10.1016/j.meatsci.2007.08.008.

Koves, T. R., R. C. Noland, A. L. Bates, S. T. Henes, D. M. Muoio, and R. N. Cortright. 2005. Subsarcolemmal and intermyofibrillar mitochondria play distinct roles in regulating skeletal muscle fatty acid metabolism. Am. J. Physiol.-Cell Ph. 288:C1074–C1082. doi: https://doi.org/10.1152/ajpcell.00391.2004.

Legako, J., T. Dinh, M. Miller, and J. Brooks. 2015. Effects of USDA beef quality grade and cooking on fatty acid composition of neutral and polar lipid fractions. Meat Sci. 100:246–255. doi: https://doi.org/10.1016/j.meatsci.2014.10.013.

Listrat, A., M. Gagaoua, D. Andueza, D. Gruffat, J. Normand, G. Mairesse, B. Picard, and J.-F. Hocquette. 2020. What are the drivers of beef sensory quality using metadata of intramuscular connective tissue, fatty acids and muscle fiber characteristics? Livest. Sci. 240:104209. doi: https://doi.org/10.1016/j.livsci.2020.104209.

Lunn, J., and H. Theobald. 2006. The health effects of dietary unsaturated fatty acids. Nutrition Bulletin 31:178–224. doi: https://doi.org/10.1111/j.1467-3010.2006.00571.x.

Lunt, D., R. Riley, and S. Smith. 1993. Growth and carcass characteristics of Angus and American Wagyu steers. Meat Sci. 34:327–334. doi: https://doi.org/10.1016/0309-1740(93)90081-R.

Ma, D., Y. H. B. Kim, B. Cooper, J.-H. Oh, H. Chun, J.-H. Choe, J. P. Schoonmaker, K. Ajuwon, and B. Min. 2017. Metabolomics profiling to determine the effect of postmortem aging on color and lipid oxidative stabilities of different bovine muscles. J. Agr. Food Chem. 65:6708–6716. doi: https://doi.org/10.1021/acs.jafc.7b02175.

Margină, D., A. Ungurianu, C. Purdel, G. M. Niţulescu, D. Tsoukalas, E. Sarandi, M. Thanasoula, T. I. Burykina, F. Tekos, and A. Buha. 2020. Analysis of the intricate effects of polyunsaturated fatty acids and polyphenols on inflammatory pathways in health and disease. Food Chem. Toxicol. 143:111558. doi: https://doi.org/10.1016/j.fct.2020.111558.

Maruri, J., and D. Larick. 1992. Volatile concentration and flavor of beef as influenced by diet. J. Food Sci. 57:1275–1281. doi: https://doi.org/10.1111/j.1365-2621.1992.tb06836.x.

Mateescu, R., A. Garmyn, M. O’neil, R. G. Tait Jr., A. Abuzaid, M. Mayes, D. Garrick, A. Van Eenennaam, D. VanOverbeke, and G. Hilton. 2012. Genetic parameters for carnitine, creatine, creatinine, carnosine, and anserine concentration in longissimus muscle and their association with palatability traits in Angus cattle. J. Anim. Sci. 90:4248–4255. doi: https://doi.org/10.2527/jas.2011-5077.

Melton, S. L. 1990. Effects of feeds on flavor of red meat: A review. J. Anim. Sci. 68:4421–4435. doi: https://doi.org/10.2527/1990.68124421x.

Melton, S. L. 1999. Current status of meat flavor. Quality attributes of muscle foods. Springer. p 115–133.

Morrison, W. R., and L. M. Smith. 1964. Preparation of fatty acid methyl esters and dimethylacetals from lipids with boron fluoride–methanol. J. Lipid Res. 5:600–608. doi: https://doi.org/10.1016/S0022-2275(20)40190-7.

Mottram, D. S. 1998. Flavour formation in meat and meat products: A review. Food Chem. 62:415–424. doi: https://doi.org/10.1016/S0308-8146(98)00076-4.

Mottram, D. S., and R. A. Edwards. 1983. The role of triglycerides and phospholipids in the aroma of cooked beef. J. Sci. Food Agr. 34:517–522. doi: https://doi.org/10.1002/jsfa.2740340513.

O’Quinn, T., J. Brooks, R. Polkinghorne, A. Garmyn, B. Johnson, J. Starkey, R. Rathmann, and M. Miller. 2012. Consumer assessment of beef strip loin steaks of varying fat levels. J. Anim. Sci. 90:626–634. doi: https://doi.org/10.2527/jas.2011-4282.

Platter, W. J., J. D. Tatum, K. E. Belk, P. L. Chapman, J. A. Scanga, and G. C. Smith. 2003. Relationships of consumer sensory ratings, marbling score, and shear force value to consumer acceptance of beef strip loin steaks. J. Anim. Sci. 81:2741–2750.

Raza, A., N. Begum, H. Song, K. Li, and P. Li. 2019. Optimization of headspace solid-phase microextraction (HS-SPME) parameters for the analysis of pyrazines in yeast extract via gas chromatography mass spectrometry (GC-MS). J. Food Sci. 84:2031–2041. doi: https://doi.org/10.1111/1750-3841.14694.

Setyabrata, D., B. R. Cooper, T. J. Sobreira, J. F. Legako, S. Martini, and Y. H. B. Kim. 2021. Elucidating mechanisms involved in flavor generation of dry-aged beef loins using metabolomics approach. Food Res. Int. 139:109969. doi: https://doi.org/10.1016/j.foodres.2020.109969.

Setyabrata, D., K. Vierck, T. R. Sheets, J. F. Legako, B. R. Cooper, T. A. Johnson, and Y. H. B. Kim. 2022. Characterizing the flavor precursors and liberation mechanisms of various dry-aging methods in cull beef loins using metabolomics and microbiome approaches. Metabolites 12:472. doi: https://doi.org/10.3390/metabo12060472.

Smith, S. B., D. K. Lunt, D. R. Smith, and R. L. Walzem. 2020. Producing high-oleic acid beef and the impact of ground beef consumption on risk factors for cardiovascular disease: A review. Meat Sci. 163:108076. doi: https://doi.org/10.1016/j.meatsci.2020.108076.

Stewart, S., G. Gardner, P. McGilchrist, D. Pethick, R. Polkinghorne, J. Thompson, and G. Tarr. 2021. Prediction of consumer palatability in beef using visual marbling scores and chemical intramuscular fat percentage. Meat Sci. 181:108322. doi: https://doi.org/10.1016/j.meatsci.2020.108322.

Turk, S., and S. Smith. 2009. Carcass fatty acid mapping. Meat Sci. 81:658–663. doi: https://doi.org/10.1016/j.meatsci.2008.11.005.

Vasta, V., and A. Priolo. 2006. Ruminant fat volatiles as affected by diet. A review. Meat Sci. 73:218–228. doi: https://doi.org/10.1016/j.meatsci.2005.11.017.

Wasserman, A. E., and F. Talley. 1968. Organoleptic identification of roasted beef, veal, lamb and pork as affected by fat. J. Food Sci. 33:219–223. doi: https://doi.org/10.1111/j.1365-2621.1968.tb01353.x.

Westerling, D. B., and H. Hedrick. 1979. Fatty acid composition of bovine lipids as influenced by diet, sex and anatomical location and relationship to sensory characteristics. J. Anim. Sci. 48:1343–1348. doi: https://doi.org/10.2527/jas1979.4861343x.

Yoo, J. H., J. W. Kim, H. I. Yong, K. H. Baek, H. J. Lee, and C. Jo. 2020. Effects of searing cooking on sensory and physicochemical properties of beef steak. Food Science of Animal Resources 40:44–54. doi: https://doi.org/10.5851/kosfa.2019.e80.

Zhang, T., C. Chen, K. Xie, J. Wang, and Z. Pan. 2021. Current state of metabolomics research in meat quality analysis and authentication. Foods 10:2388. doi: https://doi.org/10.3390/foods10102388.