Transcriptome analysis reveals the synergistic involvement of MGLL and LPIN1 in fatty acid synthesis in broiler pectoral muscles

Fatty acids (FAs) are one of the most important bioactive compounds affecting the quality of meat. In this study, we compared the expression profiles of genes involved in FA production in the breast muscle of Jingxing Yellow chickens at different days of age determined by transcriptomic analysis to identify key genes and pathways regulating the FA composition of the breast muscle. Through clustering analysis of gene expression data, the growth process of broiler chickens can be divided into two stages, namely the growth and development stage at the 35th and 63rd days of age (D35, D63), and the mature stage at the 119th day of age (D119). The content of some important unsaturated fatty acids (UFAs), such as C18:2n6c, C20:4n6, and C22:6n3, in the pectoral muscles, differed significantly between these two stages (p < 0.05). Therefore, we compared the gene expression profiles at D35 and D63 with those at D119, and identified differentially expressed genes (DEGs). The gene modules related to the five UFAs with significant changes were identified by weighted gene co-expression network analysis (WGCNA), and then 150 crossover genes were identified by crossover analysis of the detected DEGs and WGCNA. The results of the pathway enrichment analysis revealed the glycerolipid metabolism pathway related to lipid metabolism, in which the MGLL and LPIN1 genes were particularly enriched. In this study, the expression levels of MGLL and LPIN1 showed an increasing trend during the growth process of broilers, with a negative regulatory effect on the significantly reduced content of C18:2n6c in the pectoral muscle, and a positive regulatory effect on the significantly increased content of C20:4n6. These findings indicated that MGLL and LPIN1 synergistically promote the deposition of FAs, which may further promote the conversion of linoleic acid (C18:2n6c) to arachidonic acid (C20:4n6). Therefore, screening and identifying FA production-related functional genes are key to elucidate the regulatory molecular mechanism of production of FAs in chicken muscle, aiming to provide a theoretical basis for improving chicken meat quality.

Evaluation of genomic prediction considering non-additive genetic effects on fatty acid traits of Japanese Black cattle

Genomic prediction was conducted using 2494 Japanese Black cattle from Hiroshima Prefecture and both single-nucleotide polymorphism information and phenotype data on monounsaturated fatty acid (MUFA) and oleic acid (C18:1) analyzed with gas chromatography. We compared the prediction accuracy for four models (A, additive genetic effects; AD, as for A with dominance genetic effects; ADR, as for AD with the runs of homozygosity (ROH) effects calculated by ROH-based relationship matrix; and ADF, as for AD with the ROH-based inbreeding coefficient of the linear regression). Bayesian methods were used to estimate variance components. The narrow-sense heritability estimates for MUFA and C18:1 were 0.52-0.53 and 0.57, respectively; the corresponding proportions of dominance genetic variance were 0.04-0.07 and 0.04-0.05, and the proportion of ROH variance was 0.02. The deviance information criterion values showed slight differences among the models, and the models provided similar prediction accuracy.

Genetic architecture and key regulatory genes of fatty acid composition in Gushi chicken breast muscle determined by GWAS and WGCNA

BACKGROUND

Fatty acids composition in poultry muscle is directly related to its tenderness, flavour, and juiciness, whereas its genetic mechanisms have not been elucidated. In this study, the genetic structure and key regulatory genes of the breast muscle fatty acid composition of local Chinese chicken, Gushi-Anka F2 resource population by integrating genome-wide association study (GWAS) and weighted gene co-expression network analysis (WGCNA) strategies. GWAS was performed based on 323,306 single nucleotide polymorphisms (SNPs) obtained by genotyping by sequencing (GBS) method and 721 chickens from the Gushi-Anka F2 resource population with highly variable fatty acid composition traits in the breast muscle. And then, according to the transcriptome data of the candidate genes that were obtained and phenotypic data of fatty acid composition traits in breast muscle of Gushi chickens at 14, 22, and 30 weeks of age, we conducted a WGCNA.

RESULTS

A total of 128 suggestive significantly associated SNPs for 11 fatty acid composition traits were identified and mapped on chromosomes (Chr) 2, 3, 4, 5, 13, 17, 21, and 27. Of these, the two most significant SNPs were Chr13:5,100,140 (P = 4.56423e-10) and Chr13:5,100,173 (P = 4.56423e-10), which explained 5.6% of the phenotypic variation in polyunsaturated fatty acids (PUFA). In addition, six fatty acid composition traits, including C20:1, C22:6, saturated fatty acid (SFA), unsaturated fatty acids (UFA), PUFA, and average chain length (ACL), were located in the same QTL intervals on Chr13. We obtained 505 genes by scanning the linkage disequilibrium (LD) regions of all significant SNPs and performed a WGCNA based on the transcriptome data of the above 505 genes. Combining two strategies, 9 hub genes (ENO1, ADH1, ASAH1, ADH1C, PIK3CD, WISP1, AKT1, PANK3, and C1QTNF2) were finally identified, which could be the potential candidate genes regulating fatty acid composition traits in chicken breast muscle.

CONCLUSION

The results of this study deepen our understanding of the genetic mechanisms underlying the regulation of fatty acid composition traits, which is helpful in the design of breeding strategies for the subsequent improvement of fatty acid composition in poultry muscle.

Effects of the Expressions and Variants of the CAST Gene on the Fatty Acid Composition of the Longissimus Thoracis Muscle of Grazing Sonid Sheep

Fatty acid (FA) composition has an important impact on the nutrition and flavor of meat, and on consumer health, and is receiving more attention in the sheep industry. This study aimed to evaluate the relationship between the expression levels of the CAST gene and the FA composition in the longissimus thoracis (LL) muscle, to identify novel variants of CAST, and to perform association analysis with the FA composition in grazing Sonid lambs. The correlation results showed that high expression levels of CAST are correlated with better FA compositions and classes in LL. For association studies, the results showed that c.1210C>T and c.1437G>A in LD-M, and c.2097C>T mutations are associated with some compositions and classes of FA in the LL of grazing Sonid sheep. Two missense c.646G>C (G216R) and c.1210C>T (R404C) mutations were predicted to influence the Calpain_inhib domains of CAST. Thus, the correlation results and associated mutations are expected to be genetic selection markers for the FA composition and meat quality of grazing Sonid lamb muscle and provide new insights into sheep meat quality traits influenced by the ovine CAST gene.

Application of linear and machine learning models to genomic prediction of fatty acid composition in Japanese Black cattle

We collected 3180 records of oleic acid (C18:1) and monounsaturated fatty acid (MUFA) measured using gas chromatography (GC) and 6960 records of C18:1 and MUFA measured using near-infrared spectroscopy (NIRS) in intermuscular fat samples of Japanese Black cattle. We compared genomic prediction performance for four linear models (genomic best linear unbiased prediction [GBLUP], kinship-adjusted multiple loci [KAML], BayesC, and BayesLASSO) and five machine learning models (Gaussian kernel [GK], deep kernel [DK], random forest [RF], extreme gradient boost [XGB], and convolutional neural network [CNN]). For GC-based C18:1 and MUFA, KAML showed the highest accuracies, followed by BayesC, XGB, DK, GK, and BayesLASSO, with more than 6% gain of accuracy by KAML over GBLUP. Meanwhile, DK had the highest prediction accuracy for NIRS-based C18:1 and MUFA, but the difference in accuracies between DK and KAML was slight. For all traits, accuracies of RF and CNN were lower than those of GBLUP. The KAML extends GBLUP methods, of which marker effects are weighted, and involves only additive genetic effects; whereas machine learning methods capture non-additive genetic effects. Thus, KAML is the most suitable method for breeding of fatty acid composition in Japanese Black cattle.

A selected population study reveals the biochemical mechanism of intramuscular fat deposition in chicken meat

BACKGROUND

Increasing intramuscular fat (IMF) is an important strategy to improve meat quality, but the regulation mechanism of IMF deposition needs to be systematically clarified.

RESULTS

A total of 520 chickens from a selected line with improved IMF content and a control line were used to investigate the biochemical mechanism of IMF deposition in chickens. The results showed that the increased IMF would improve the flavor and tenderness quality of chicken meat. IMF content was mainly determined both by measuring triglyceride (TG) and phospholipid (PLIP) in muscle tissue, but only TG content was found to be decisive for IMF deposition. Furthermore, the increase in major fatty acid (FA) components in IMF is mainly derived from TGs (including C16:0, C16:1, C18:1n9c, and C18:2n6c, etc.), and the inhibition of certain very-long-chain FAs would help to IMF/TG deposition.

CONCLUSIONS

Our study elucidated the underlying biochemical mechanism of IMF deposition in chicken: Prevalent accumulation of long-chain FAs and inhibitions of medium-chain FAs and very long chain FA would jointly result in the increase of TGs with the FA biosynthesis and cellular uptake ways. Our findings will guide the production of high-quality chicken meat.

The accuracy of genomic prediction for meat quality traits in Hanwoo cattle when using genotypes from different SNP densities and preselected variants from imputed whole genome sequence

Context Genomic prediction is the use of genomic data in the estimation of genomic breeding values (GEBV) in animal breeding. In beef cattle breeding programs, genomic prediction increases the rates of genetic gain by increasing the accuracy of selection at earlier ages. Aims The objectives of the study were to examine the effect of single-nucleotide polymorphism (SNP) density and to evaluate the effect of using SNPs preselected from imputed whole-genome sequence for genomic prediction. Methods Genomic and phenotypic data from 2110 Hanwoo steers were used to predict GEBV for marbling score (MS), meat texture (MT), and meat colour (MC) traits. Three types of SNP densities including 50k, high-density (HD), and whole-genome sequence data and preselected SNPs from genome-wide association study (GWAS) were used for genomic prediction analyses. Two scenarios (independent and dependent discovery populations) were used to select top significant SNPs. The accuracy of GEBV was assessed using random cross-validation. Genomic best linear unbiased prediction (GBLUP) was used to predict the breeding values for each trait. Key results Our result showed that very similar prediction accuracies were observed across all SNP densities used in the study. The prediction accuracy among traits ranged from 0.29±0.05 for MC to 0.46±0.04 for MS. Depending on the studied traits, up to 5% of prediction accuracy improvement was obtained when the preselected SNPs from GWAS analysis were included in the prediction analysis. Conclusions High SNP density such as HD and the whole-genome sequence data yielded a similar prediction accuracy in Hanwoo beef cattle. Therefore, the 50K SNP chip panel is sufficient to capture the relationships in a breed with a small effective population size such as the Hanwoo cattle population. Preselected variants improved prediction accuracy when they were included in the genomic prediction model. Implications The estimated genomic prediction accuracies are moderately accurate in Hanwoo cattle and for searching for SNPs that are more productive could increase the accuracy of estimated breeding values for the studied traits.

Hepatic transcriptome analysis identifies genes, polymorphisms and pathways involved in the fatty acids metabolism in sheep

Fatty acids (FA) in ruminants, especially unsaturated FA (USFA) have important impact in meat quality, nutritional value, and flavour quality of meat, and on consumer’s health. Identification of the genetic factors controlling the FA composition and metabolism is pivotal to select sheep that produce higher USFA and lower saturated (SFA) for the benefit of sheep industry and consumers. Therefore, this study was aimed to investigate the transcriptome profiling in the liver tissues collected from sheep with divergent USFA content in longissimus muscle using RNA deep-sequencing. From sheep (n = 100) population, liver tissues with higher (n = 3) and lower (n = 3) USFA content were analysed using Illumina HiSeq 2500. The total number of reads produced for each liver sample were ranged from 21.28 to 28.51 million with a median of 23.90 million. Approximately, 198 genes were differentially regulated with significance level of p-adjusted value <0.05. Among them, 100 genes were up-regulated, and 98 were down-regulated (p<0.01, FC>1.5) in the higher USFA group. A large proportion of key genes involved in FA biosynthesis, adipogenesis, fat deposition, and lipid metabolism were identified, such as APOA5, SLC25A30, GFPT1, LEPR, TGFBR2, FABP7, GSTCD, and CYP17A. Pathway analysis revealed that glycosaminoglycan biosynthesis- keratan sulfate, adipokine signaling, galactose metabolism, endocrine and other factors-regulating calcium metabolism, mineral metabolism, and PPAR signaling pathway were playing important regulatory roles in FA metabolism. Importantly, polymorphism and association analyses showed that mutation in APOA5, CFHR5, TGFBR2 and LEPR genes could be potential markers for the FA composition in sheep. These polymorphisms and transcriptome networks controlling the FA variation could be used as genetic markers for FA composition-related traits improvement. However, functional validation is required to confirm the effect of these SNPs in other sheep population in order to incorporate them in the sheep breeding program.

Estimation of genetic parameters for carcass grading traits, image analysis traits, and monounsaturated fatty acids in Japanese Black cattle from Hyogo Prefecture

Genetic parameters for carcass grading traits, image analysis traits, and monounsaturated fatty acid (MUFA) percentages were estimated in 29,942 Japanese Black cattle from Hyogo Prefecture. The analyzed traits included five carcass grading traits, two image analysis traits, fat area ratio and fineness index, and two MUFA traits, one measured in intermuscular fat using near-infrared spectroscopy (NIRS) and the other in intramuscular fat using gas chromatography (GC). The heritability estimates of image analysis traits and MUFA were moderate to high, ranging from 0.395 to 0.740, and it was considered that they could be improved simultaneously with carcass grading traits because no severe genetic antagonism was observed. Although the heritability of the NIRS-based intermuscular MUFA was slightly lower than that of the GC-based intramuscular MUFA, the genetic correlation between the two methods was as high as 0.804. These results indicate that the NIRS method can be used as an alternative evaluation procedure to predict MUFA in intramuscular fat in the longissimus muscle.

Immortalizing Cells for Human Consumption

The need to produce immortal, food-relevant cell lines is one of the most pressing challenges of cellular agriculture, the field which seeks to produce meat and other animal products via tissue engineering and synthetic biology. Immortal cell lines have a long and complicated story, from the first recognized immortal human cell lines taken from Henrietta Lacks, to today, where they are used to assay toxicity and produce therapeutics, to the future, where they could be used to create meat without harming an animal. Although work in immortal cell lines began more than 50 years ago, there are few existing cell lines made of species and cell types appropriate for cultured meat. Cells in cultured meat will be eaten by consumers; therefore, cultured meat cell lines will also require unique attributes not selected for in other cell line applications. Specifically, cultured meat cell lines will need to be approved as safe for consumption as food, proliferate and differentiate efficiently at industrial scales, and have desirable taste, texture, and nutrition characteristics for consumers. This paper defines what cell lines are needed, the existing methods to produce new cell lines and their limitations, and the unique considerations of cell lines used in cultured meat.

Are Marbling and the Prediction of Beef Eating Quality Affected by Different Grading Sites?

For the European abattoirs, the preferred carcass grading site is at the fifth rib, and cutting at the tenth rib as in Australia could lead to a lower economic value of the carcass. Therefore, the objective of this study was to compare the grading scores of marbling and the meat and fat color on Musculus longissimus thoracis et lumborum (LTL) at the fifth and the tenth thoracic vertebrae. The consequences on the prediction of beef eating quality using the Meat Standards Australia (MSA) grading scheme were also evaluated for cull cows, which produce the majority of beef consumed in France. Carcasses from 208 French cattle, mainly Limousine cows, were graded according to the Australian Beef Carcase Chiller Assessment System (ABCAS) used for the implementation of the MSA system. The results indicate that there was no significant difference in the marbling score, between the fifth and the tenth ribs and hence in the MSA index and in the Global Quality [meat quality (MQ4)] scores calculated from marbling values from either the fifth rib or the tenth rib. However, the meat color at the tenth rib was significantly darker than that at the fifth rib (p < 0.01), and the fat color at the tenth rib was significantly yellower than that at the fifth rib (p < 0.001). The results of this study suggest that the grading of marbling can be conducted on M. LTL at the fifth thoracic vertebrae for routine use of the MSA system in France and, more generally, in Europe. However, further investigation and adjustment would be needed for other critical MSA scores (such as rib fat thickness) while respecting the European carcass quartering practices.

Genomic imprinting variances of beef carcass traits and physiochemical characteristics in Japanese Black cattle

The objective of this study was to estimate variance components related to imprinting for carcass traits and physiochemical characteristics in Japanese Black cattle. The carcass records obtained from 4,220 Japanese Black feedlot cattle included carcass weight (CW), rib eye area (REA), rib thickness, subcutaneous fat thickness, and beef marbling score (BMS), and the physiochemical characteristics were fat, moisture, glycogen per proportion of moisture content, oleic acid, and monounsaturated fatty acids (MUFA). To detect gametic effects, an imprinting model was fitted. High additive heritabilities were estimated for all traits (from 0.516 for glycogen to 0.853 for fat) and were reduced in Mendelian heritability. The range of the differences was from 0.002 (CW) to 0.331 (fat and moisture), and the reductions were due to their imprinting variances. The ratio of the imprinting variance to the total additive genetic variance for REA (0.374), BMS (0.291), fat (0.387), moisture (0.388), and MUFA (0.337) were large (p < 0.05). These imprinting variances were due to the maternal contribution and suggested the existence of maternally expressed genomic imprinting effects on the traits in Japanese Black cattle. Therefore, maternal gametic effects should be considered in breeding programs for Japanese Black cattle.

Genetic associations between reproductive traits for first calving and growth curve characteristics of Japanese Black cattle

The objective of this study was to estimate genetic parameters for first calving reproductive traits and growth curve characteristics in Japanese Black cattle. The Gompertz growth function was fitted to body weight-age data to obtain the mature weight (MWT) and rate of maturing (ROM) of cows. Data of reproductive traits including the first service conception rate (CR) for heifers, age at the first calving (AFC), and gestation length for the first calving were collected. Records of 3,204 animals were used for analysis. Genetic parameters were estimated using a linear uni- and bivariate animal model. The heritability estimates were moderate (0.29 for ROM) and high (0.57 for MWT) for growth curve parameters and low (0.03-0.11) for reproductive traits. There was a negative genetic correlation between MWT and ROM (-0.26), suggesting that an animal with a faster ROM would show a lower MWT. CR was negatively correlated with MWT (-0.42) but significantly and positively correlated with ROM (0.91). There was a negative genetic correlation between AFC and MWT (-0.49). These results suggest that a heifer with a faster ROM and lower MWT would show a higher CR. Meanwhile, a heifer with a lower MWT would show a higher AFC.

Relationships Between Fat and Cholesterol Contents and Fatty Acid Composition in Different Meat-Producing Animal Species

The objective of this study was to determine the relationships between intramuscular fat, cholesterol contents and fatty acid composition in the muscles of different animal species. Intramuscular fat, cholesterol and fatty acid composition in 207 muscle samples from 129 animals of different species (pigs, beef cattle, farmed red deer, horses and geese) were determined and analysed. The obtained results indicated unequal relations between intramuscular fat and cholesterol contents and fatty acid proportions in the muscles of different animal species. The increase of intramuscular fat content resulted in higher monounsaturated and lower polyunsaturated fatty acid contents in most muscles of meat producing animals. In all the species higher fatness did not show any increase in cholesterol content and also cholesterol contents were lower as fat increased in m. semimembranosus of pigs and m. pectoralis profundus of horses. The cholesterol content positively correlated with saturated and monounsaturated fatty acids in the longissimus muscle with the lowest fat content found in red deer and beef cattle, whereas the correlations between these measures were negative in m. pectoralis profundus of horses and the breast of goose containing high fat levels. Negative correlations between polyunsaturated fatty acids and cholesterol content were found in the longissimus muscle of red deer and cattle, whereas these correlations in goose breast and horse meat were positive.

Effects of parent-of-origin models with different pedigree information on beef carcass traits and fatty acid composition in Japanese Black cattle

Genomic imprinting should be considered in animal breeding systems to avoid lead in bias in genetic parameter estimation. The objective of this study was to clarify the effects of pedigree information on imprinting variances for carcass traits and fatty acid composition in Japanese Black cattle. Carcass records [carcass weight, rib eye area, rib thickness (RT), subcutaneous fat thickness and beef marbling score (BMS)] and fatty acid composition were obtained for 11,855 Japanese Black feedlot cattle. To estimate and compare the imprinting variances for the traits, two imprinting models with different pedigree information [the sire-dam gametic relationship matrix (Model 1) and the sire-maternal grandsire (MGS) numerator relationship matrix (Model 2)] were fitted. The ratio of the imprinting variance to the total additive genetic variance for RT (6.33%) and BMS (19.00%) was significant in Model 1, but only that for BMS (21.09%) was significant in Model 2. This study revealed that fitting the sire-MGS model could be useful in estimating imprinting variance under certain conditions, such as when restricted pedigree information is available. Furthermore, the present result suggested that the maternal gametic effects on BMS should be included in breeding programmes for Japanese Black cattle to avoid selection bias caused by imprinting effects.

Genome-Wide Association Study of Meat Quality Traits in Hanwoo Beef Cattle Using Imputed Whole-Genome Sequence Data

The discovery of single nucleotide polymorphisms (SNP) and the subsequent genotyping of large numbers of animals have enabled large-scale analyses to begin to understand the biological processes that underpin variation in animal populations. In beef cattle, genome-wide association studies using genotype arrays have revealed many quantitative trait loci (QTL) for various production traits such as growth, efficiency and meat quality. Most studies regarding meat quality have focused on marbling, which is a key trait associated with meat eating quality. However, other important traits like meat color, texture and fat color have not commonly been studied. Developments in genome sequencing technologies provide new opportunities to identify regions associated with these traits more precisely. The objective of this study was to estimate variance components and identify significant variants underpinning variation in meat quality traits using imputed whole genome sequence data. Phenotypic and genomic data from 2,110 Hanwoo cattle were used. The estimated heritabilities for the studied traits were 0.01, 0.16, 0.31, and 0.49 for fat color, meat color, meat texture and marbling score, respectively. Marbling score and meat texture were highly correlated. The genome-wide association study revealed 107 significant SNPs located on 14 selected chromosomes (one QTL region per selected chromosome). Four QTL regions were identified on BTA2, 12, 16, and 24 for marbling score and two QTL regions were found for meat texture trait on BTA12 and 29. Similarly, three QTL regions were identified for meat color on BTA2, 14 and 24 and five QTL regions for fat color on BTA7, 10, 12, 16, and 21. Candidate genes were identified for all traits, and their potential influence on the given trait was discussed. The significant SNP will be an important inclusion into commercial genotyping arrays to select new breeding animals more accurately.

Identifying labelling and marketing advantages of nutrients in minced beef meat: A case study

In 2006, the European Commission approved nutrition and health claim regulations of foods to stimulate healthier choices. To document how commercial, minced beef meat complied with regulations, meat samples from 72 carcasses were analysed. These samples were a source of niacin (B3), pyridoxine (B6), phosphorus and iron (cow meat only), and a rich source of protein, monounsaturated fat, vitamin B12 and zinc. A potential exists for establishing beef meat as a source of vitamin K, iron and selenium. The meat's nutrient relevance for young women when ingesting 150 g of raw beef mince/day was estimated. Increased levels of riboflavin (B2), pantothenic acid (B5), iron and selenium beyond presently observed, would better support this group's recommended nutrient intake. If the bioactivity of 25-OH-vitamin D3 could be used in calculations, vitamin D3 in the minced meat would add positively to the intake of vitamin D3 that was 49% of the recommended intake.

Genome wide association study identifies SNPs associated with fatty acid composition in Chinese Wagyu cattle

Background

Fatty acids are important traits that affect meat quality and nutritive values in beef cattle. Detection of genetic variants for fatty acid composition can help to elucidate the genetic mechanism underpinning these traits and promote the improvement of fatty acid profiles. In this study, we performed a genome-wide association study (GWAS) on fatty acid composition using high-density single nucleotide polymorphism (SNP) arrays in Chinese Wagyu cattle.

Results

In total, we detected 15 and 8 significant genome-wide SNPs for individual fatty acids and fatty acid groups in Chinese Wagyu cattle, respectively. Also, we identified nine candidate genes based on 100 kb regions around associated SNPs. Four SNPs significantly associated with C14:1 cis-9 were embedded with stearoyl-CoA desaturase (SCD), while three SNPs in total were identified for C22:6 n-3 within Phospholipid scramblase family member 5 (PLSCR5), Cytoplasmic linker associated protein 1 (CLASP1), and Chymosin (CYM). Notably, we found the top candidate SNP within SCD can explain ~ 7.37% of phenotypic variance for C14:1 cis-9. Moreover, we detected several blocks with high LD in the 100 kb region around SCD. In addition, we found three significant SNPs within a 100 kb region showing pleiotropic effects related to multiple FA groups (PUFA, n-6, and PUFA/SFA), which contains BAI1 associated protein 2 like 2 (BAIAP2L2), MAF bZIP transcription factor F (MAFF), and transmembrane protein 184B (TMEM184B).

Conclusions

Our study identified several significant SNPs and candidate genes for individual fatty acids and fatty acid groups in Chinese Wagyu cattle, and these findings will further assist the design of breeding programs for meat quality in cattle.