Accuracy of genomic breeding values for meat tenderness in Polled Nellore cattle

Leveraging Functional Genomics for Understanding Beef Quality Complexities and Breeding Beef Cattle for Improved Meat Quality

Consumer perception of beef is heavily influenced by overall meat quality, a critical factor in the cattle industry. Genomics has the potential to improve important beef quality traits and identify genetic markers and causal variants associated with these traits through genomic selection (GS) and genome-wide association studies (GWAS) approaches. Transcriptomics, proteomics, and metabolomics provide insights into underlying genetic mechanisms by identifying differentially expressed genes, proteins, and metabolic pathways linked to quality traits, complementing GWAS data. Leveraging these functional genomics techniques can optimize beef cattle breeding for enhanced quality traits to meet high-quality beef demand. This paper provides a comprehensive overview of the current state of applications of omics technologies in uncovering functional variants underlying beef quality complexities. By highlighting the latest findings from GWAS, GS, transcriptomics, proteomics, and metabolomics studies, this work seeks to serve as a valuable resource for fostering a deeper understanding of the complex relationships between genetics, gene expression, protein dynamics, and metabolic pathways in shaping beef quality.

A network-based approach to understanding gene-biological processes affecting economically important traits of Nelore cattle

This study aimed to build gene-biological process networks with differentially expressed genes associated with economically important traits of Nelore cattle from 17 previous studies. The genes were clustered into three groups by evaluated traits: group 1, production traits; group 2, carcass traits; and group 3, meat quality traits. For each group, a gene-biological process network analysis was performed with the differentially expressed genes in common. For production traits, 37 genes were found in common, of which 13 genes were enriched for six Gene Ontology (GO) terms; these terms were not functionally grouped. However, the enriched GO terms were related to homeostasis, the development of muscles and the immune system. For carcass traits, four genes were found in common. Thus, it was not possible to functionally group these genes into a network. For meat quality traits, the analysis revealed 222 genes in common. CSRP3 was the only gene differentially expressed in all three groups. Non-redundant biological terms for clusters of genes were functionally grouped networks, reflecting the cross-talk between all biological processes and genes involved. Many biological processes and pathways related to muscles, the immune system and lipid metabolism were enriched, such as striated muscle cell development and triglyceride metabolic processes. This study provides insights into the genetic mechanisms of production, carcass and meat quality traits of Nelore cattle. This information is fundamental for a better understanding of the complex traits and could help in planning strategies for the production and selection systems of Nelore cattle.

Assessment of genomic breeding values and their accuracies for carcass traits in Jeju Black cattle using whole-genome SNP chip panels

The objective of the present study was to evaluate the breeding value and accuracy of genomic estimated breeding values (GEBVs) of carcass traits in Jeju Black cattle (JBC) using Hanwoo steers and JBC as a reference population using the single-trait animal model. Our research included genotype and phenotype information on 19,154 Hanwoo steers with 1097 JBC acting as the reference population. Likewise, the test population consisted of 418 genotyped JBC individuals with no phenotypic records for those carcass traits. For estimating the accuracy of GEBV, we divided the entire population into three groups. Hanwoo and JBC make up the first group; Hanwoo and JBC, who has both the genotype and phenotypic records, are referred to as the reference (training) population, and JBC, who lacks phenotypic information is referred to as the test (validation) population. The second group consists of the JBC (without phenotype) as the test population and Hanwoo as a reference population with phenotype and genotypic data. The only JBCs in the third group are those who have genotypic and phenotypic data on them as a reference population but no phenotypic data on them as a test population. The single-trait animal model was used in all three groups for statistical purposes. The reference populations estimated heritabilities for carcass weight (CWT), eye muscle area (EMA), backfat thickness (BF), and marbling score (MS) as 0.30, 0.26, 0.26, and 0.34 for the Hanwoo steer and 0.42, 0.27, 0.26, and 0.48 for JBC. The average accuracy for carcass traits in Group 1 was 0.80 for the Hanwoo and JBC reference population compared with 0.73 for the JBC test population. Although the average accuracy for carcass traits in Group 2 was 0.80, it was 0.80 for the Hanwoo reference population and only 0.56 for the JBC test population. The average accuracy for the JBC reference and test populations was 0.68 and 0.50, respectively, when they were included in the accuracy comparison without the Hanwoo reference population. Groups 1 and 2 used Hanwoo as reference population, which led to a better average accuracy; however, Group 3 only used the JBC reference and test population, which led to a lower average accuracy. This might be due to the fact that Group 3 used a smaller reference size than the group that came before it and that the genetic makeup of the Hanwoo and JBC breeds differed. The GEBV accuracy for MS was higher than that of other traits across all three analysis groups, followed by CWT, EMA, and BF, which may be partially explained by the MS traits' higher heritability. This study suggests that in order to achieve more accuracy, a large reference population particular to a breed should be established. Therefore, to increase the accuracy of GEBV prediction and the genetic benefit from genomic selection in JBC, individual reference breeds, and large populations are required.

Genome-enabled prediction of meat and carcass traits using Bayesian regression, single-step genomic best linear unbiased prediction and blending methods in Nelore cattle

Several methods have been used for genome-enabled prediction (or genomic selection) of complex traits, for example, multiple regression models describing a target trait with a linear function of a set of genetic markers. Genomic selection studies have been focused mostly on single-trait analyses. However, most profitability traits are genetically correlated, and an increase in prediction accuracy of genomic breeding values for genetically correlated traits is expected when using multiple-trait models. Thus, this study was carried out to assess the accuracy of genomic prediction for carcass and meat quality traits in Nelore cattle, using single- and multiple-trait approaches. The study considered 15 780, 15 784, 15 742 and 526 records of rib eye area (REA, cm²), back fat thickness (BF, mm), rump fat (RF, mm) and Warner-Bratzler shear force (WBSF, kg), respectively, in Nelore cattle, from the Nelore Brazil Breeding Program. Animals were genotyped with a low-density single nucleotide polymorphism (SNP) panel and subsequently imputed to arrays with 54 and 777 k SNPs. Four Bayesian specifications of genomic regression models, namely, Bayes A, Bayes B, Bayes Cπ and Bayesian Ridge Regression; blending methods, BLUP; and single-step genomic best linear unbiased prediction (ssGBLUP) methods were compared in terms of prediction accuracy using a fivefold cross-validation. Estimates of heritability ranged from 0.20 to 0.35 and from 0.21 to 0.46 for RF and WBSF on single- and multiple-trait analyses, respectively. Prediction accuracies for REA, BF, RF and WBSF were all similar using the different specifications of regression models. In addition, this study has shown the impact of genomic information upon genetic evaluations in beef cattle using the multiple-trait model, which was also advantageous compared to the single-trait model because it accounted for the selection process using multiple traits at the same time. The advantage of multi-trait analyses is attributed to the consideration of correlations and genetic influences between the traits, in addition to the non-random association of alleles.

Improving genomic prediction accuracy for meat tenderness in Nellore cattle using artificial neural networks

The goal of this study was to compare the predictive performance of artificial neural networks (ANNs) with Bayesian ridge regression, Bayesian Lasso, Bayes A, Bayes B and Bayes Cπ in estimating genomic breeding values for meat tenderness in Nellore cattle. The animals were genotyped with the Illumina Bovine HD Bead Chip (HD, 777K from 90 samples) and the GeneSeek Genomic Profiler (GGP Indicus HD, 77K from 485 samples). The quality control for the genotypes was applied on each Chip and comprised removal of SNPs located on non-autosomal chromosomes, with minor allele frequency <5%, deviation from HWE (p < 10^-6 ), and with linkage disequilibrium >0.8. The FImpute program was used for genotype imputation. Pedigree-based analyses indicated that meat tenderness is moderately heritable (0.35), indicating that it can be improved by direct selection. Prediction accuracies were very similar across the Bayesian regression models, ranging from 0.20 (Bayes A) to 0.22 (Bayes B) and 0.14 (Bayes Cπ) to 0.19 (Bayes A) for the additive and dominance effects, respectively. ANN achieved the highest accuracy (0.33) of genomic prediction of genetic merit. Even though deep neural networks are recognized to deliver more accurate predictions, in our study ANN with one single hidden layer, 105 neurons and rectified linear unit (ReLU) activation function was sufficient to increase the prediction of genetic merit for meat tenderness. These results indicate that an ANN with relatively simple architecture can provide superior genomic predictions for meat tenderness in Nellore cattle.

Genome-enable prediction for health traits using high-density SNP panel in US Holstein cattle

The objective of this study was to compare accuracies of different Bayesian regression models in predicting molecular breeding values for health traits in Holstein cattle. The dataset was composed of 2505 records reporting the occurrence of retained fetal membranes (RFM), metritis (MET), mastitis (MAST), displaced abomasum (DA), lameness (LS), clinical endometritis (CE), respiratory disease (RD), dystocia (DYST) and subclinical ketosis (SCK) in Holstein cows, collected between 2012 and 2014 in 16 dairies located across the US. Cows were genotyped with the Illumina BovineHD (HD, 777K). The quality controls for SNP genotypes were HWE P‐value of at least 1 × 10⁻¹⁰; MAF greater than 0.01 and call rate greater than 0.95. The fimpute program was used for imputation of missing SNP markers. The effect of each SNP was estimated using the Bayesian Ridge Regression (BRR), Bayes A, Bayes B and Bayes Cπ methods. The prediction quality was assessed by the area under the curve, the prediction mean square error and the correlation between genomic breeding value and the observed phenotype, using a leave‐one‐out cross‐validation technique that avoids iterative cross‐validation. The highest accuracies of predictions achieved were: RFM [Bayes B (0.34)], MET [BRR (0.36)], MAST [Bayes B (0.55), DA [Bayes Cπ (0.26)], LS [Bayes A (0.12)], CE [Bayes A (0.32)], RD [Bayes Cπ (0.23)], DYST [Bayes A (0.35)] and SCK [Bayes Cπ (0.38)] models. Except for DA, LS and RD, the predictive abilities were similar between the methods. A strong relationship between the predictive ability and the heritability of the trait was observed, where traits with higher heritability achieved higher accuracy and lower bias when compared with those with low heritability. Overall, it has been shown that a high‐density SNP panel can be used successfully to predict genomic breeding values of health traits in Holstein cattle and that the model of choice will depend mostly on the genetic architecture of the trait.

Development and evaluation of a low-density single-nucleotide polymorphism chip specific to Bos indicus cattle

Context Genomic selection has been of increasing interest in the genetic improvement of Zebu cattle, particularly for quantitative traits that are difficult or expensive to measure, such as carcass traits and meat tenderness. The success of genomic selection depends on several factors, and at its core is the availability of single-nucleotide polymorphism (SNP) chips that are appropriately designed for Bos indicus cattle. However, the currently available commercial bovine SNP chips are mostly designed for Bos taurus cattle. There are two commercial Bos indicus SNP chips; namely, GeneSeek genomic profiler high-density Bos indicus (GGP-HDi) SNP chip and a low-density (LD) Bos indicus SNP chip (Z chip), but these two Bos indicus SNP chips were built with mixed contents of SNPs for Bos indicus and Bos taurus cattle, due to limited availability of genotype data from Bos indicus cattle. Aims To develop a new GGP indicus 35000 SNP chip specifically for Bos indicus cattle, which has a low cost, but high accuracy of imputation to Illumina BovineHD chips. Methods The design of the chip consisted of 34000 optimally selected SNPs, plus 1000 SNPs pre-reserved for those on the Y chromosome, ‘causative’ mutations for a variety of economically relevant traits, genetic health conditions and International Society for Animal Genetics globally recognised parentage markers for those breeds of cattle. Key results The present results showed that this new indicus LD SNP chip had considerably increased minor allele frequencies in indicus breeds than the previous Z-chip. It demonstrated with high imputation accuracy to HD SNP genotypes in five indicus breeds, and with considerable predictability on 14 growth and reproduction traits in Nellore cattle. Conclusions This new indicus LD chip represented a successful effort to leverage existing knowledge and genotype resources towards the public release of a cost-effective LD SNP chip specifically for Bos indicus cattle, which is expected to replace the previous GGP indicus LD chip and to supplement the existing GGP-HDi 80000 SNP chip. Implications A new SNP chip specifically designed for Bos indicus, with high power of imputation to Illumina BovineHD technology and with excellent coverage of the whole genome, is now available on the market for Bos indicus cattle, and Bos indicus and Bos taurus crosses.

Likelihood of obtaining tender meat from confined calf

Abstract The objective was to evaluate the main factors that influence the shear strength of meat from confined steers and the probability of obtaining soft meat. For this purpose, we evaluated the literature on carcass and/or meat of beef steers in Brazil published between January 1999 and April 2019 and extracted the content from three sections for analysis: materials and methods, results, and discussion. Pearson’s correlation was used to analyze the data, and the stepwise statistic was used to determine the proportion of the synchronized effect of variables on shear force. For determining the probability of tenderness, meat with a shear force lower than 4.6 kgf/cm3 was classified as soft; meat with a higher sheer force was classified as hard. Following the classification, logistic regression analysis and odds ratio test were performed. The factors of study location, the proportion of zebu background in the genome, finishing weight, the percentage of concentrate in the diet, and finishing period and meat marbling explained 62.45% of the variability in the shear strength of beef. The following strategies were found to increase the chances of effectively obtaining soft meat from confined steers: starting the termination phase early even in animals with lower weights, prolonging the confinement time, increasing concentrate percentage in the diet, and a higher marbling degree. It is possible to estimate a large proportion of shear force variability using the production variables (ante-mortem), and the process can be adjusted accordingly to considerable increase the possibility of obtaining soft meat.

Beef Tenderness Prediction by a Combination of Statistical Methods: Chemometrics and Supervised Learning to Manage Integrative Farm-To-Meat Continuum Data

This trial aimed to integrate metadata that spread over farm-to-fork continuum of 110 Protected Designation of Origin (PDO)Maine-Anjou cows and combine two statistical approaches that are chemometrics and supervised learning; to identify the potential predictors of beef tenderness analyzed using the instrumental Warner-Bratzler Shear force (WBSF). Accordingly, 60 variables including WBSF and belonging to 4 levels of the continuum that are farm-slaughterhouse-muscle-meat were analyzed by Partial Least Squares (PLS) and three decision tree methods (C&RT: classification and regression tree; QUEST: quick, unbiased, efficient regression tree and CHAID: Chi-squared Automatic Interaction Detection) to select the driving factors of beef tenderness and propose predictive decision tools. The former method retained 24 variables from 59 to explain 75% of WBSF. Among the 24 variables, six were from farm level, four from slaughterhouse level, 11 were from muscle level which are mostly protein biomarkers, and three were from meat level. The decision trees applied on the variables retained by the PLS model, allowed identifying three WBSF classes (Tender (WBSF ≤ 40 N/cm²), Medium (40 N/cm² < WBSF < 45 N/cm²), and Tough (WBSF ≥ 45 N/cm²)) using CHAID as the best decision tree method. The resultant model yielded an overall predictive accuracy of 69.4% by five splitting variables (total collagen, µ-calpain, fiber area, age of weaning and ultimate pH). Therefore, two decision model rules allow achieving tender meat on PDO Maine-Anjou cows: (i) IF (total collagen < 3.6 μg OH-proline/mg) AND (µ-calpain ≥ 169 arbitrary units (AU)) AND (ultimate pH < 5.55) THEN meat was very tender (mean WBSF values = 36.2 N/cm², n = 12); or (ii) IF (total collagen < 3.6 μg OH-proline/mg) AND (µ-calpain < 169 AU) AND (age of weaning < 7.75 months) AND (fiber area < 3100 µm²) THEN meat was tender (mean WBSF values = 39.4 N/cm², n = 30).

Genomewide association mapping and pathway analysis of meat tenderness in Polled Nellore cattle

Brazil is one of the world's largest beef exporters, although the product has a low price due to quality issues. The meat exported by Brazil is considered medium and low quality by international buyers, mainly due to lack of tenderness. The predominant Zebu breeds (80% Nellore) are known for producing tougher beef than taurine breeds. Nonetheless, some studies have shown that there is substantial genetic variability for tenderness within the Nellore breed, although it is a difficult trait to improve by conventional selection methods. Therefore, the aim of this study was to perform a genomewide association study (GWAS) and a gene set enrichment analysis to identify genomic regions and biologically relevant pathways associated with meat tenderness in Polled Nellore cattle. Data consisted of Warner-Bratzler shear force values of LM from 427 Polled Nellore animals divided into 3 experimental slaughters (years 2005, 2008, and 2010). The animals were genotyped with either the Illumina BovineHD BeadChip (777k, on 61 samples) or the GGP Indicus HD chip (77k, on 366 samples). Single nucleotide polymorphisms were excluded when the call rate was <90%, the Hardy-Weinberg proportions -value was <1% (Fisher exact test, Bonferroni adjusted), and the minor allele frequency was <1%. Imputation from the GGP Indicus HD chip to the Illumina BovineHD BeadChip was performed using the FImput program. Genomewide association analysis was performed using the Efficient Mixed Model Association eXpedited (EMMAx) and the population parameters previously determined (P3D) methods. The GWAS was complemented with a gene set enrichment analysis performed using the FatiGO procedure. Significant markers ( < 0.0001) explaining a larger proportion of variation than other significant SNPs were located on chromosomes 3, 13, 17, 20, 21, and 25, indicating QTL associated with meat tenderness throughout the genome. Additionally, gene set analysis identified 22 Gene Ontology functional terms and 2 InterPro entries that showed significant enrichment of genes associated with tenderness. The functional categories included protein tyrosine and serine/threonine kinase activity, calcium ion binding, lipid metabolic process, and growth factors, among others. These results help to elucidate the genetic architecture and metabolic pathways underlying this trait, which is of extreme economic and social importance to Brazil, because Nellore is the dominant beef cattle breed in the country.