Genome scan for meat quality traits in Nelore beef cattle

Novel putative causal mutations associated with fat traits in Nellore cattle uncovered by eQTLs located in open chromatin regions

Intramuscular fat (IMF) and backfat thickness (BFT) are critical economic traits impacting meat quality. However, the genetic variants controlling these traits need to be better understood. To advance knowledge in this area, we integrated RNA-seq and single nucleotide polymorphisms (SNPs) identified in genomic and transcriptomic data to generate a linkage disequilibrium filtered panel of 553,581 variants. Expression quantitative trait loci (eQTL) analysis revealed 36,916 cis-eQTLs and 14,408 trans-eQTLs. Association analysis resulted in three eQTLs associated with BFT and 24 with IMF. Functional enrichment analysis of genes regulated by these 27 eQTLs revealed noteworthy pathways that can play a fundamental role in lipid metabolism and fat deposition, such as immune response, cytoskeleton remodeling, iron transport, and phospholipid metabolism. We next used ATAC-Seq assay to identify and overlap eQTL and open chromatin regions. Six eQTLs were in regulatory regions, four in predicted insulators and possible CCCTC-binding factor DNA binding sites, one in an active enhancer region, and the last in a low signal region. Our results provided novel insights into the transcriptional regulation of IMF and BFT, unraveling putative regulatory variants.

Genome-wide scans identify biological and metabolic pathways regulating carcass and meat quality traits in beef cattle

Genome association studies (GWAS) provides knowledge about the genetic architecture of beef-related traits that allow linking the target phenotype to genomic information aiding breeding decision. Thus, the present study aims to uncover the genetic mechanism involved in carcass (REA: rib eye area, BF: backfat thickness, and HCW: hot carcass weight) and meat quality traits (SF: shear-force, MARB: marbling score, and IMF: intramuscular fat content) in Nellore cattle. For this, 6910 young bulls with phenotypic information and 23,859 animals genotyped with 435 k markers were used to perform the weighted single-step GBLUP (WssGBLUP) approach, considering two iterations. The top 10 genomic regions explained 8.13, 11.81, and 9.58% of the additive genetic variance, harboring a total of 119, 143, and 95 positional candidate genes for REA, BF, and HCW, respectively. For meat quality traits, the top 10 windows explained a large proportion of the total genetic variance for SF (14.95%), MARB (17.56%), and IMF (21.41%) surrounding 92, 155, and 111 candidate genes, respectively. Relevant candidate genes (CAST, PLAG1, XKR4, PLAGL2, AQP3/AQP7, MYLK2, WWOX, CARTPT, and PLA2G16) are related to physiological aspects affecting growth, carcass, meat quality, feed intake, and reproductive traits by signaling pathways controlling muscle control, key signal metabolic molecules INS / IGF-1 pathway, lipid metabolism, and adipose tissue development. The GWAS results provided insights into the genetic control of the traits studied and the genes found are potential candidates to be used in the improvement of carcass and meat quality traits.

Identification of genomic regions and pathways associated with traits related to rumen acidosis in feedlot Nellore cattle

There may be an increased risk of metabolic disorders, such as rumen acidosis, in cattle fed high-concentrate diets, particularly those from Bos taurus indicus genotypes, which have shown to be more sensitive to ruminal acidification. Therefore, this study aimed to estimate (co)variance components and identify genomic regions and pathways associated with ruminal acidosis in feedlot Nellore cattle fed high-concentrate diets. It was utilized a dataset containing a total of 642 Nellore bulls that were genotyped from seven feedlot nutrition studies. The GGP Indicus 35k panel was used with the single step genome-wide association study methodology in which the effects of the markers were obtained from the genomic values estimated by the GBLUP model. A bivariate model to estimate genetic correlations between the economically important traits and indicator traits for acidosis was used. The traits evaluated in this study that were nutritionally related to rumen acidosis included average daily gain (ADG), final body weight, time spent eating (TSE), time spent ruminating, rumenitis score (RUM), rumen absorptive surface area (ASA), rumen keratinized layer thickness (KER) and hot carcass weight (HCW). The identified candidate genes were mainly involved in the negative or non-regulation of the apoptotic process, salivary secretion, and transmembrane transport. The genetic correlation between HCW and ASA was low positive (0.27 ± 0.23), and between ADG and ASA was high moderate (0.58 ± 0.59). A positive genetic correlation between RUM and all performance traits was observed, and TSE correlated negatively with HCW (-0.33 ± 0.21), ASA (-0.75 ± 0.48), and KER (-0.40 ± 0.27). The genetic association between economically important traits and indicator traits for acidosis suggested that Nellore cattle may be more sensitive to acidosis in feedlot systems.

Whole-genome SNP markers reveal runs of homozygosity in indigenous cattle breeds of Pakistan

The runs of homozygosity (ROH) were identified in 14 Pakistani cattle breeds (n = 105) by genotyping with the Illumina 50 K SNP BeadChip. These breeds were categorized into Dairy, Dual, and Draft breeds based on their utility and production performance. We identified a total of 10,936 ROHs which mainly consisted of a high number of shorter segments (1-4 Mb). Dairy group exhibited the highest level of inbreeding (FROH: 0.078 ± 0.028) while the lowest (FROH: 0.002 ± 0.008) was observed in Dual group. In 48 genomic regions identified with a high frequency of ROH, 207 genes were detected in the three breed groups. A substantially higher number of ROH islands detected in dairy breeds indicated the impact of the positive selection pressure over the years. Important candidate genes and QTL were detected in the ROH islands associated with economic traits like milk production, reproduction, meat, carcass, and health traits in dairy cattle.

Estimating genetic parameters of reproductive, carcass, and meat quality traits in Polled Nellore cattle

Considering the economic and commercial efficiency of the beef production chain, the yield and quality of the meat produced must also be included in breeding programs. For the Nellore breed, including the polled herd, these aspects have not been much studied. The aim of this study was to estimate genetic parameters for scrotal circumference adjusted to 365 (SC365) and 450 (SC450) days of age, age at first calving (AFC), accumulated productivity (AP), stayability (STAY), longissimus muscle area (LMA), thickness of subcutaneous fat over the 12th-13th ribs (BF), thickness of subcutaneous fat over the rump (RF), and shear force measured by Warner-Bratzler shear force (WBSF) of polled Nellore cattle. Bayesian analyses were performed by adopting a linear animal model, whereas STAY analyses used the linear threshold model. Heritability estimates were 0.31 (SC365), 0.37 (SC450), 0.16 (AFC), 0.25 (AP), 0.16 (STAY), 0.30 (LMA), 0.13 (BF), 0.24 (RF), and 0.15 (WBSF), indicating moderate response to selection. Genetic and residual correlations between SC365 and SC450 were high (0.91 and 0.74, respectively), as well as the genetic correlations of AP with SC365, SC450, AFC, and STAY (0.61, 0.62, - 0.69, and 0.83, respectively). Genetic and residual correlations of WBSF with reproductive and carcass characteristics exhibited high standard deviations, however favorable. Based on the results, it is expected that in the medium term, animals with greater sexual precocity will also have greater accumulated productivity and longer permanence of females in the herd, along with superior carcass traits. However, due to the low heritabilities and small genetic associations with reproductive traits, fat thickness characteristics (BF and RF) will still require direct selection.

Bivariate GWAS reveals pleiotropic regions among feed efficiency and beef quality-related traits in Nelore cattle

Feed-efficient cattle selection is among the most leading solutions to reduce cost for beef cattle production. However, technical difficulties in measuring feed efficiency traits had limited the application in livestock. Here, we performed a Bivariate Genome-Wide Association Study (Bi-GWAS) and presented candidate biological mechanisms underlying the association between feed efficiency and meat quality traits in a half-sibling design with 353 Nelore steers derived from 34 unrelated sires. A total of 13 Quantitative Trait Loci (QTL) were found explaining part of the phenotypic variations. An important transcription factor of adipogenesis in cattle, the TAL1 (rs133408775) gene located on BTA3 was associated with intramuscular fat and average daily gain (IMF-ADG), and a region located on BTA20, close to CD180 and MAST4 genes, both related to fat accumulation. We observed a low positive genetic correlation between IMF-ADG (r = 0.30 ± 0.0686), indicating that it may respond to selection in the same direction. Our findings contributed to clarifying the pleiotropic modulation of the complex traits, indicating new QTLs for bovine genetic improvement.

METTL21C mediates lysine trimethylation of IGF2BP1 to regulate chicken myoblast proliferation

1. Methyltransferase-like 21C (METTL21C) and insulin-like growth factor 2 mRNA-binding protein 1 (IGF2BP1) play important roles in the proliferation of chicken myoblasts. However, it remains unclear whether there is protein-protein interaction between METTL21C and IGF2BP1 to regulate proliferation of chicken myoblasts.2. In this study, the Igf2bp1 gene was amplified from cDNA of liver tissue of Lueyang black-bone chicken to construct the overexpression vector HA-Igf2bp1. The HA-Igf2bp1 and Flag-Mettl21c vectors were individually transfected and co-transfected into HEK293T, respectively. Co-immunoprecipitation (Co-IP) assay indicated a protein-protein interaction between METTL21C and IGF2BP1.3. Using the Western blotting and LC-MS/MS, it was found that METTL21C could mediate the lysine methylation modification of IGF2BP1. Furthermore, the His-tagged overexpression vector HA-Igf2bp1-His was constructed, transfected and co-transfected with Flag-Mettl21c into HEK293T. His-tagged IGF2BP1 was purified by nickel ion affinity chromatography. Western blotting revealed that IGF2BP1 was successfully purified, and the trimethylation modification level of co-transfection group was significantly elevated compared with the single-transfection Igf2bp1 group.4. Mettl21c and Igf2bp1 overexpression vectors were transfected and co-transfected into primary chicken myoblasts, respectively. The results of 5-ethynyl-2'-deoxyuridine assay and the expression level of Pax7 and MyoD indicated that overexpression of Igf2bp1 alone inhibited the chicken myoblast proliferation, whereas co-expression of Mettl21c and Igf2bp1 eliminated the inhibitory effects of Igf2bp1, thereby favouring cell proliferation and differentiation.5. The results, for the first time, revealed that METTL21C mediated the lysine trimethylation modification of IGF2BP1 to regulate the proliferation of chicken myoblasts, which provided a new insight into in-depth analysis of the molecular mechanism of METTL21C methylation involved in regulating the growth and development of skeletal muscle in Lueyang black-bone chicken.

Genomic prediction for meat and carcass traits in Nellore cattle using a Markov blanket algorithm

This study was carried out to evaluate the advantage of preselecting SNP markers using Markov blanket algorithm regarding the accuracy of genomic prediction for carcass and meat quality traits in Nellore cattle. This study considered 3675, 3680, 3660 and 524 records of rib eye area (REA), back fat thickness (BF), rump fat (RF), and Warner-Bratzler shear force (WBSF), respectively, from the Nellore Brazil Breeding Program. The animals have been genotyped using low-density SNP panel (30 k), and subsequently imputed for arrays with 777 k SNPs. Four Bayesian specifications of genomic regression models, namely Bayes A, Bayes B, Bayes Cπ and Bayesian Ridge Regression methods were compared in terms of prediction accuracy using a five folds cross-validation. Prediction accuracy for REA, BF and RF was all similar using the Bayesian Alphabet models, ranging from 0.75 to 0.95. For WBSF, the predictive ability was higher using Bayes B (0.47) than other methods (0.39 to 0.42). Although the prediction accuracies using Markov blanket of SNP markers were lower than those using all SNPs, for WBSF the relative gain was lower than 13%. With a subset of informative SNPs markers, identified using Markov blanket, probably, is possible to capture a large proportion of the genetic variance for WBSF. The development of low-density and customized arrays using Markov blanket might be cost-effective to perform a genomic selection for this trait, increasing the number of evaluated animals, improving the management decisions based on genomic information and applying genomic selection on a large scale.

Differential Allele-Specific Expression Revealed Functional Variants and Candidate Genes Related to Meat Quality Traits in B. indicus Muscle

Traditional transcriptomics approaches have been used to identify candidate genes affecting economically important livestock traits. Regulatory variants affecting these traits, however, remain under covered. Genomic regions showing allele-specific expression (ASE) are under the effect of cis-regulatory variants, being useful for improving the accuracy of genomic selection models. Taking advantage of the better of these two methods, we investigated single nucleotide polymorphisms (SNPs) in regions showing differential ASE (DASE SNPs) between contrasting groups for beef quality traits. For these analyses, we used RNA sequencing data, imputed genotypes and genomic estimated breeding values of muscle-related traits from 190 Nelore (Bos indicus) steers. We selected 40 contrasting unrelated samples for the analysis (N = 20 animals per contrasting group) and used a beta-binomial model to identify ASE SNPs in only one group (i.e., DASE SNPs). We found 1479 DASE SNPs (FDR ≤ 0.05) associated with 55 beef-quality traits. Most DASE genes were involved with tenderness and muscle homeostasis, presenting a co-expression module enriched for the protein ubiquitination process. The results overlapped with epigenetics and phenotype-associated data, suggesting that DASE SNPs are potentially linked to cis-regulatory variants affecting simultaneously the transcription and phenotype through chromatin state modulation.

Bibliometric Analysis of Research on the Main Genes Involved in Meat Tenderness

Simple Summary

A bibliometric analysis was carried out to know the evolution of research on genes associated with meat tenderness, of interest for the development of selection programs. Since 1993, studies have been limited to a few researchers in high-income countries due to the costs associated with the techniques. The main findings showed that the scientific production had a discontinuous growth because science experienced a significant change since approximately 2010. Marker-assisted selection was widely used, evaluating mainly CAPN (calpain) and CAST (calpastatin) genes for their contribution to meat tenderness, especially in cattle. However, the effects are small; therefore, genomic selection was implemented by genotyping thousands of single nucleotide polymorphisms (SNPs) for further explanation of genetic variation. The results shown are important for scholars to identify emerging methodologies and gaps in the literature and to know who the prolific authors and institutions in the field for possible collaborations, etc., are.

Abstract

Tenderness is one of the main characteristics of meat because it determines its price and acceptability. This is the first bibliometric study on the trend of research on the role of genes in meat tenderness. A total of 175 original and English-language articles published up to 2021 were retrieved from Scopus. The bibliometric analysis was carried out with VOSviewer (version 1.6.18, Eck and Waltman, Leiden, Netherlands) and complemented with the Analyze search results service from Scopus. Erroneous and duplicate data were eliminated, and incomplete information was added to standardize the results. Scientific production was evaluated by means of quantity, quality and structure indicators. As a first glance, 8.816% of authors have published more than 50% of papers mainly related to genes encoding the calpain (CAPN)-calpastatin (CAST) system and single nucleotide polymorphisms (SNPs). Among other findings, a strong link was found between the contribution of the main countries (led by the United States with) and their institutions (led by the USDA Agricultural Research Service with) to their gross domestic product. Most studies on the topic are published in the Journal of Animal Science, and other journals with high impact according to the number of citations and different metrics. Finally, when evaluating the most cited articles, the occurrence and association of the main keywords, it was confirmed that research is focused on the role of CAPN and CAST genes and of SNPs in beef tenderness. The change in science was emphasized; although marker-assisted selection is still used, genes have an infinitesimal effect on complex traits. Therefore, since about 2010, new research groups adopted genomic selection to evaluate dense panels of SNPs and better explain genetic variation in meat tenderness.

Multi-omic data integration for the study of production, carcass, and meat quality traits in Nellore cattle

Data integration using hierarchical analysis based on the central dogma or common pathway enrichment analysis may not reveal non-obvious relationships among omic data. Here, we applied factor analysis (FA) and Bayesian network (BN) modeling to integrate different omic data and complex traits by latent variables (production, carcass, and meat quality traits). A total of 14 latent variables were identified: five for phenotype, three for miRNA, four for protein, and two for mRNA data. Pearson correlation coefficients showed negative correlations between latent variables miRNA 1 (mirna1) and miRNA 2 (mirna2) (−0.47), ribeye area (REA) and protein 4 (prot4) (−0.33), REA and protein 2 (prot2) (−0.3), carcass and prot4 (−0.31), carcass and prot2 (−0.28), and backfat thickness (BFT) and miRNA 3 (mirna3) (−0.25). Positive correlations were observed among the four protein factors (0.45–0.83): between meat quality and fat content (0.71), fat content and carcass (0.74), fat content and REA (0.76), and REA and carcass (0.99). BN presented arcs from the carcass, meat quality, prot2, and prot4 latent variables to REA; from meat quality, REA, mirna2, and gene expression mRNA1 to fat content; from protein 1 (prot1) and mirna2 to protein 5 (prot5); and from prot5 and carcass to prot2. The relations of protein latent variables suggest new hypotheses about the impact of these proteins on REA. The network also showed relationships among miRNAs and nebulin proteins. REA seems to be the central node in the network, influencing carcass, prot2, prot4, mRNA1, and meat quality, suggesting that REA is a good indicator of meat quality. The connection among miRNA latent variables, BFT, and fat content relates to the influence of miRNAs on lipid metabolism. The relationship between mirna1 and prot5 composed of isoforms of nebulin needs further investigation. The FA identified latent variables, decreasing the dimensionality and complexity of the data. The BN was capable of generating interrelationships among latent variables from different types of data, allowing the integration of omics and complex traits and identifying conditional independencies. Our framework based on FA and BN is capable of generating new hypotheses for molecular research, by integrating different types of data and exploring non-obvious relationships.

Integrative Analysis Between Genome-Wide Association Study and Expression Quantitative Trait Loci Reveals Bovine Muscle Gene Expression Regulatory Polymorphisms Associated With Intramuscular Fat and Backfat Thickness

Understanding the architecture of gene expression is fundamental to unravel the molecular mechanisms regulating complex traits in bovine, such as intramuscular fat content (IMF) and backfat thickness (BFT). These traits are economically important for the beef industry since they affect carcass and meat quality. Our main goal was to identify gene expression regulatory polymorphisms within genomic regions (QTL) associated with IMF and BFT in Nellore cattle. For that, we used RNA-Seq data from 193 Nellore steers to perform SNP calling analysis. Then, we combined the RNA-Seq SNP and a high-density SNP panel to obtain a new dataset for further genome-wide association analysis (GWAS), totaling 534,928 SNPs. GWAS was performed using the Bayes B model. Twenty-one relevant QTL were associated with our target traits. The expression quantitative trait loci (eQTL) analysis was performed using Matrix eQTL with the complete SNP dataset and 12,991 genes, revealing a total of 71,033 cis and 36,497 trans-eQTL (FDR < 0.05). Intersecting with QTL for IMF, we found 231 eQTL regulating the expression levels of 117 genes. Within those eQTL, three predicted deleterious SNPs were identified. We also identified 109 eQTL associated with BFT and affecting the expression of 54 genes. This study revealed genomic regions and regulatory SNPs associated with fat deposition in Nellore cattle. We highlight the transcription factors FOXP4, FOXO3, ZSCAN2, and EBF4, involved in lipid metabolism-related pathways. These results helped us to improve our knowledge about the genetic architecture behind important traits in cattle.

EEF1A1 transcription cofactor gene polymorphism is associated with muscle gene expression and residual feed intake in Nelore cattle

Cis-acting effects of noncoding variants on gene expression and regulatory molecules constitute a significant factor for phenotypic variation in complex traits. To provide new insights into the impacts of single-nucleotide polymorphisms (SNPs) on transcription factors (TFs) and transcription cofactors (TcoF) coding genes, we carried out a multi-omic analysis to identify cis-regulatory effects of SNPs on these genes' expression in muscle and describe their association with feed efficiency-related traits in Nelore cattle. As a result, we identified one SNP, the rs137256008C > T, predicted to impact the EEF1A1 gene expression (β = 3.02; P-value = 3.51E-03) and the residual feed intake trait (β = - 3.47; P-value = 0.02). This SNP was predicted to modify transcription factor sites and overlaps with several QTL for feed efficiency traits. In addition, co-expression network analyses showed that animals containing the T allele of the rs137256008 SNP may be triggering changes in the gene network. Therefore, our analyses reinforce and contribute to a better understanding of the biological mechanisms underlying gene expression control of feed efficiency traits in bovines. The cis-regulatory SNP can be used as biomarker for feed efficiency in Nelore cattle.

Genome-Wide Association Study of Potential Meat Quality Trait Loci in Ducks

With continuously increasing living standards and health requirements of consumers, meat quality is becoming an important consideration while buying meat products. To date, no genome-wide association study (GWAS) for copy number variants (CNVs) and single nucleotide polymorphisms (SNPs) has been conducted to reveal the genetic effects on meat quality in ducks. This study analyzed the phenotypic correlation and heritability of fat, water, collagen, and protein content of duck breast muscle. To identify the candidate variants for meat quality, we performed a GWAS using 273 ducks from an F2 population. The results of the SNP GWAS showed that the BARHL2, COPS7B, and CCDC50 genes were associated with fat content; BLM, WDR76, and EOMES with water content; CAMTA1, FGD5, GRM7, and RAPGEF5 with collagen production; and RIMS2, HNRNPU, and SPTBN1 with protein content. Additionally, 3, 7, 1, and 3 CNVs were associated with fat, water, collagen, and protein content, respectively, in duck breast muscle. The genes identified in this study can serve as markers for meat quality. Furthermore, our findings may help devise effective breeding plans and selection strategies to improve meat quality.

DNA methylation may affect beef tenderness through signal transduction in Bos indicus

BACKGROUND

Beef tenderness is a complex trait of economic importance for the beef industry. Understanding the epigenetic mechanisms underlying this trait may help improve the accuracy of breeding programs. However, little is known about epigenetic effects on Bos taurus muscle and their implications in tenderness, and no studies have been conducted in Bos indicus.

RESULTS

Comparing methylation profile of Bos indicus skeletal muscle with contrasting beef tenderness at 14 days after slaughter, we identified differentially methylated cytosines and regions associated with this trait. Interestingly, muscle that became tender beef had higher levels of hypermethylation compared to the tough group. Enrichment analysis of predicted target genes suggested that differences in methylation between tender and tough beef may affect signal transduction pathways, among which G protein signaling was a key pathway. In addition, different methylation levels were found associated with expression levels of GNAS, PDE4B, EPCAM and EBF3 genes. The differentially methylated elements correlated with EBF3 and GNAS genes overlapped CpG islands and regulatory elements. GNAS, a complex imprinted gene, has a key role on G protein signaling pathways. Moreover, both G protein signaling pathway and the EBF3 gene regulate muscle homeostasis, relaxation, and muscle cell-specificity.

CONCLUSIONS

We present differentially methylated loci that may be of interest to decipher the epigenetic mechanisms affecting tenderness. Supported by the previous knowledge about regulatory elements and gene function, the methylation data suggests EBF3 and GNAS as potential candidate genes and G protein signaling as potential candidate pathway associated with beef tenderness via methylation.

Genetic diversity and population structure of six autochthonous pig breeds from Croatia, Serbia, and Slovenia

BACKGROUND

The importance of local breeds as genetic reservoirs of valuable genetic variation is well established. Pig breeding in Central and South-Eastern Europe has a long tradition that led to the formation of several local pig breeds. In the present study, genetic diversity parameters were analysed in six autochthonous pig breeds from Slovenia, Croatia and Serbia (Banija spotted, Black Slavonian, Turopolje pig, Swallow-bellied Mangalitsa, Moravka and Krskopolje pig). Animals from each of these breeds were genotyped using microsatellites and single nucleotide polymorphisms (SNPs). The results obtained with these two marker systems and those based on pedigree data were compared. In addition, we estimated inbreeding levels based on the distribution of runs of homozygosity (ROH) and identified genomic regions under selection pressure using ROH islands and the integrated haplotype score (iHS).

RESULTS

The lowest heterozygosity values calculated from microsatellite and SNP data were observed in the Turopolje pig. The observed heterozygosity was higher than the expected heterozygosity in the Black Slavonian, Moravka and Turopolje pig. Both types of markers allowed us to distinguish clusters of individuals belonging to each breed. The analysis of admixture between breeds revealed potential gene flow between the Mangalitsa and Moravka, and between the Mangalitsa and Black Slavonian, but no introgression events were detected in the Banija spotted and Turopolje pig. The distribution of ROH across the genome was not uniform. Analysis of the ROH islands identified genomic regions with an extremely high frequency of shared ROH within the Swallow-bellied Mangalitsa, which harboured genes associated with cholesterol biosynthesis, fatty acid metabolism and daily weight gain. The iHS approach to detect signatures of selection revealed candidate regions containing genes with potential roles in reproduction traits and disease resistance.

CONCLUSIONS

Based on the estimation of population parameters obtained from three data sets, we showed the existence of relationships among the six pig breeds analysed here. Analysis of the distribution of ROH allowed us to estimate the level of inbreeding and the extent of homozygous regions in these breeds. The iHS analysis revealed genomic regions potentially associated with phenotypic traits and allowed the detection of genomic regions under selection pressure.

Sustained Effects of Muscle Calpain System Genotypes on Tenderness Phenotypes of South African Beef Bulls during Ageing up to 20 Days

Simple Summary

When searching for genetic markers for the selection of more tender beef, it is important to maintain minimal environmental variation from pre-slaughter, right through to the ageing process, to ensure the accuracy of the obtained phenotypes. This is because beef quality traits have a large environmental component that can greatly alter the characteristics of the meat, which would not reflect a true genetic effect. We propose that variable ageing times are especially important in determining whether markers are associated with tenderization or not. Our analyses included candidate genes for the protein degrading enzyme system for calpains, because they contribute the most to tenderization. We were able to validate these markers in South African beef cattle, where they could be useful for selection. The timing of the collection of tenderness data was critical, as only a few (6/134) genetic markers sustained their association with tenderization over ageing to 20 days. A larger tenderization effect earlier in ageing, as shown here for the capn1_187 and capn1_4751 markers, would decrease the length of ageing. This would not only increase profits, but also decrease the energy needed during the storage and refrigeration of aged beef, decreasing the carbon footprint of beef production.

Abstract

The most important factor that determines beef tenderness is its proteolytic activity, and the balance between calpain-1 protease activity and calpastatin inhibition is especially important, while contributions can also arise from calpain-2 and, possibly, calpain-3. The meat ageing process itself affects these processes. To determine whether genotypes in the calpain–calpastatin system can enhance tenderness through a 20-day ageing period, South African purebred beef bulls (n = 166) were genotyped using the Illumina BovineHD SNP BeadChip through a gene-based association analysis targeting the cast, capn3, capn2 and capn1 genes. The Warner–Bratzler shear force (WBSF) and myofibril fragment length (MFL) of Longissimus thoracis et lumborum (LTL) steaks were evaluated between d 3 and d 20 of ageing, with protease enzyme activity in the first 20 h post-mortem. Although several of the 134 SNPs are associated with tenderness, only seven SNP in the cast, capn2 and capn1 genes sustained genetic associations, additive to the ageing-associated increases in tenderness for at least three of the four ageing periods. While most genomic associations were relatively stable over time, some genotypes within the SNP responded differently to ageing, resulting in altered genomic effects over time. The level of ageing at which genomic associations are performed is an important factor that determines whether SNPs affect tenderness phenotypes.

Genome-wide association study for meat tenderness in beef cattle identifies patterns of the genetic contribution in different post-mortem stages

The beef tenderization process during the post-mortem period is one of the most important sensorial attributes and it is well-established. The aim of this study was to identify the genetic contribution pattern to meat tenderness at 7-(LMD7), 14-(LMD14), and 21-(LMD21) days post-mortem. The heritabilities for LMD7 (0.194), LMD14 (0.142) and LMD21 (0.048) are well established in the population evaluated here. However, its genetic contribution in terms of genomic candidate regions is still poorly understood. Tenderness was measured in the Longissiums thoracis using Warner-Bratzler shear force in the three post-mortem periods. A total of 4323 crossbred beef cattle were phenotyped and genotyped using the Illumina BovineSNP50K. The percentage of the total genetic variance was estimated using the weighted single-step genomic best linear unbiased prediction method. The main candidate windows for LMD7 were associated with proteolysis of myofibrillar structures and the weakening endomysium and perimysium. Candidate windows for LMD14 and LMD21 were mapped in bovine QTLs for body composition, height and growth. Results presented herein highlight, the largest contribution of proteolysis related processes before 14-days post-mortem and body composition characteristics in later stages for meat tenderness.

Expression of calpain system transcripts responds inversely to beef tenderization after vitamin D3 supplementation in Nellore cattle

ABSTRACT The calpain system is the central player for meat tenderization and the calpastatin activity plays an important role in beef tenderness of Bos taurus indicus cattle. This study investigated whether dietary vitamin D3 induced gene expression of μ-calpain (CAPN1), total calpastatin (CAST T), and their variants (CAST I and II). Animals received none or 2 × 106 IU of vitamin D3 for either 2 or 8 days before slaughter and were submitted to different conditions during feedlot: sun exposure or artificial shade. Steaks from Longissimus lumborum were fabricated, aged for 1, 7, and 21 days post-mortem, and later used for the analyses of shear force and the myofibrillar fragmentation index. Vitamin D3 did not influence mRNA abundance; however, it induced a greater CAST II transcript in animals supplemented 8 days before slaughter. There was a negative association between CAST II abundance and the shear force, which revealed an important modulation of the calpain system expression due to vitamin D supplementation. This result is an important tool for strategies to improve beef tenderness.

Genes and proteins associated with ribeye area and meat tenderness in a commercial Nellore cattle population

Despite several studies on genetic markers and differently expressed genes related to ribeye area (REA) and tenderness traits in beef cattle, there is divergence in the results regarding the genes associated with these traits. Thirteen genes that had been associated or have biological functions that may influence such phenotypes were included in this study. A total of five genes for REA (IGF-1, IGF-2, MSTN, NEDD4, and UBE4A) and eight genes for meat tenderness (CAPN1, CAPN2, CAST, HSPB1, DNAJA1, FABP4, SCD, and PRKAG3) were selected from previously studies in beef cattle. Genes and its respective proteins expression were validated in a commercial population of Nellore cattle using quantitative real-time PCR (RT-qPCR) and advanced mass spectrometry (LC / MS-MS) techniques, respectively. MSTN gene was upregulated in animals with low REA. CAPN1, CAPN2, CAST, HSPB1, and DNAJA1 genes were upregulated in animals with tougher meat. The proteins translated by these genes were not differentially expressed. Our results could confirm the potential of some studied genes as biomarkers for carcass and meat quality in Nellore cattle.

Genome-wide association study of resistance/susceptibility to infectious bovine keratoconjunctivitis in Brazilian Hereford cattle

Genome-wide association studies were conducted to identify the more informative genomic regions and SNPs, as well as to identify candidate genes associated with infectious bovine keratoconjunctivitis (IBK) resistance/susceptibility in Hereford cattle. A Bayes B statistical approach was initially applied in genome-wide association studies by using deregressed estimated breeding values for IBK resistance/susceptibility. To estimate the combined effect of a genomic region that is potentially associated with QTL, 2504 non-overlapping 1-Mb windows that varied in SNP number were defined, with the most informative 24 windows including 427 SNPs and explaining more than 20% of the estimated genetic variance for IBK resistance/susceptibility. These regions were explored with respect to their biological functions through functional analysis to map potential candidate genes. The significant SNPs were mapped on chromosomes 1, 3, 5, 6, 7, 8, 10, 11, 12, 13, 14, 15, 18, 20, 23, and 28, and candidate genes were detected as related to the IBK. Most informative SNPs in term of genetic variance were located in proximity of genes related to phenotypic expression of lesions and biological processes associated to the IBK. Knowledge about phenotypic and genomic variation generated in the present study can be used to on design selection strategies to improve the resistance to IBK of Hereford cattle herds.

Genome-wide association and genomic prediction for scrotal circumference in Hereford and Braford bulls

Scrotal circumference (SC) is widely used as a selection criterion for bulls in breeding programs, since it is easily assessed and correlated with several desirable reproductive traits. The objectives of this study were: to perform a genome-wide association study (GWAS) to identify genomic regions associated with SC adjusted for age (SCa) and for both age and weight (SCaw); to select Tag SNPs from GWAS to construct low-density panel for genomic prediction; and to compare the prediction accuracy of the SC through different methods for Braford and Hereford bulls from the same genetic breeding program. Data of SC from 18,172 bulls (30.4 ± 3.7 cm) and of genotypes from 131 sires and 3,545 animals were used. From GWAS, the top 1% of 1-Mb windows were observed on chromosome (BTA) 2, 20, 7, 8, 15, 3, 16, 27, 6 and 8 for SCa and on BTA 8, 15, 16, 21, 19, 2, 6, 5 and 10 for SCaw, representing 17.4% and 18.8% of the additive genetic variance of SCa and SCaw, respectively. The MeSH analysis was able to translate genomic information providing biological meanings of more specific gene functions related to the SCa and SCaw. The genomic enhancement methods, especially single step GBLUP, that combined phenotype and pedigree data with direct genomic values generated gains in accuracy in relation to pedigree BLUP, suggesting that genomic predictions should be applied to improve genetic gain and to narrow the generation interval compared to traditional methods. The proposed Tag-SNP panels may be useful for lower-cost commercial genomic prediction applications in the future, when the number of bulls in the reference population increases for SC in Hereford and Braford breeds.

Introgression contributes to distribution of structural variations in cattle

Structural variations (SVs) are an important source of phenotypic diversity in cattle. Here, 72 whole genome sequences representing taurine and zebu cattle were used to identify SVs. Applying multiple approaches, 16,738 SVs were identified. A comparison against the Database of Genomic Variants archives revealed that 1575 SVs were novel in our data. A novel duplication covering the entire GALNT15 gene, was observed only in N'Dama. A duplication, which was previously reported only in zebu and associated with navel length, was also observed in N'Dama. Investigation of a novel deletion located upstream of CAST13 gene and identified only in Italian cattle and zebu, revealed its introgressed origin in the former. Overall, our data highlights how the SVs distribution in cattle is also shaped by forces such as demographical differences and gene flow. The cattle SVs of this study and its meta-data can be visualized on an interactive genome browser at https://tinyurl.com/svCowArs.

Transcriptome profiling analysis of muscle tissue reveals potential candidate genes affecting water holding capacity in Chinese Simmental beef cattle

Water holding capacity (WHC) is an important sensory attribute that greatly influences meat quality. However, the molecular mechanism that regulates the beef WHC remains to be elucidated. In this study, the longissimus dorsi (LD) muscles of 49 Chinese Simmental beef cattle were measured for meat quality traits and subjected to RNA sequencing. WHC had significant correlation with 35 kg water loss (r = − 0.99, p < 0.01) and IMF content (r = 0.31, p < 0.05), but not with SF (r = − 0.20, p = 0.18) and pH (r = 0.11, p = 0.44). Eight individuals with the highest WHC (H-WHC) and the lowest WHC (L-WHC) were selected for transcriptome analysis. A total of 865 genes were identified as differentially expressed genes (DEGs) between two groups, of which 633 genes were up-regulated and 232 genes were down-regulated. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment revealed that DEGs were significantly enriched in 15 GO terms and 96 pathways. Additionally, based on protein–protein interaction (PPI) network, animal QTL database (QTLdb), and relevant literature, the study not only confirmed seven genes (HSPA12A, HSPA13, PPARγ, MYL2, MYPN, TPI, and ATP2A1) influenced WHC in accordance with previous studies, but also identified ATP2B4, ACTN1, ITGAV, TGFBR1, THBS1, and TEK as the most promising novel candidate genes affecting the WHC. These findings could offer important insight for exploring the molecular mechanism underlying the WHC trait and facilitate the improvement of beef quality.

Genome-Wide Association Study Provides Insights into Important Genes for Reproductive Traits in Nelore Cattle

Simple Summary

In this study, we investigated the association between single nucleotide polymorphisms (SNPs) and reproductive traits in order to identify candidate genes and biological pathways associated with these traits in Nelore beef cattle. The genome-wide association analysis revealed genomic regions that could explain part of the genetic variance of the studied traits. The results revealed genes with important functions for reproductive traits, such as fertility and precocity. Some genes were associated with more than one trait, being important for reproductive efficiency. The identification of candidate genes that were associated with the studied traits as well as genes enriched in the functional terms and pathways may be useful for exploring the genetic architecture underlying reproductive traits and may be used in Nelore breeding programs.

Abstract

The identification of genomic regions associated with reproductive traits as well as their biological processes allows a better understanding of the phenotypic variability of these traits. This information could be applied to animal breeding programs to accelerate genetic gain. The aim of this study was to evaluate the association between single nucleotide polymorphisms (SNP) with a scrotal circumference at 365 days of age (SC365) and at 450 days of age (SC450), gestation length (GL) as a calf trait, age at first calving (AFC), accumulated productivity (ACP), heifer early calving until 30 months (HC30), and stayability (STAY) traits, in order to identify candidate genes and biological pathways associated with reproductive traits in Nelore cattle. The data set consisted of pedigree, phenotypes, and genotypes of Nelore cattle from the “Associação Nacional de Criadores e Pesquisadores” (ANCP). The association analyses were performed using the Weighted Single-Step Genome-Wide Association method; the regions, consisting of 10 consecutive SNP, which explained more than 0.5% of additive genetic variance, were considered as a significant association. A total of 3, 6, 7, 5, 10, 25, and 12 windows were associated with SC355, SC450, GL, AFC, ACP, HC30, and STAY, respectively. The results revealed genes with important functions for reproductive traits, such as fertility and precocity. Some genes were associated with more than one trait, among them CAMK1D, TASP1, ACOXL, RAB11FIP5, and SFXN5. Moreover, the genes were enriched in functional terms, like negative regulation of fat cell differentiation, fatty acid alpha-oxidation, and sphingolipids signaling pathway. The identification of the genes associated with the traits, as well as genes enriched in the terms and pathway mentioned above, should contribute to future biological validation studies and may be used as candidate genes in Nelore breeding programs.

Allelic and genotypic frequencies for loci associated with meat quality in Mexican Braunvieh cattle

The objective was to estimate allelic and genotypic frequencies for loci associated with meat quality in a Mexican population of Braunvieh cattle. Information was obtained from 300 animals genotyped with the Genomic Profile Bovine LD chip of 30K and 50K SNPs. After the final edition, including quality control, the data contained information for 12 loci of the CAPN1, CAPN3, CAPN5, CAPN14, DGAT1, DGAT2, TG, ANK1, and MADH3 genes. Allelic and genotypic frequencies and Hardy-Weinberg equilibrium were estimated with the Cervus 3.0.7 software. The studied population markers were in Hardy-Weinberg equilibrium, except for those associated with CAPN1, DGAT1, and MADH3. Frequencies higher than those reported for other breeds were found for genotypes associated with meat softness, higher marbling score, lower quantity of saturated fatty acids, and lower shear force (CAPN1 and DGAT2). There were similarities with frequencies reported for Bos taurus breeds for the CAPN3 and TG genes. For the DGAT1 and ANK1 genes, the frequencies of the desired genotypes were low. A marker for DGAT1 and another for MADH3 were monomorphic. The results of this study are encouraging in terms of the potential of the Braunvieh population studied for breeding programs aiming to increase meat quality. The breed has strengths that could be used either by crossbreeding to generate heterozygous animals or by selection to increase frequencies of valuable alleles.

Muscle allele-specific expression QTLs may affect meat quality traits in Bos indicus

Single nucleotide polymorphisms (SNPs) located in transcript sequences showing allele-specific expression (ASE SNPs) were previously identified in the Longissimus thoracis muscle of a Nelore (Bos indicus) population consisting of 190 steers. Given that the allele-specific expression pattern may result from cis-regulatory SNPs, called allele-specific expression quantitative trait loci (aseQTLs), in this study, we searched for aseQTLs in a window of 1 Mb upstream and downstream from each ASE SNP. After this initial analysis, aiming to investigate variants with a potential regulatory role, we further screened our aseQTL data for sequence similarity with transcription factor binding sites and microRNA (miRNA) binding sites. These aseQTLs were overlapped with methylation data from reduced representation bisulfite sequencing (RRBS) obtained from 12 animals of the same population. We identified 1134 aseQTLs associated with 126 different ASE SNPs. For 215 aseQTLs, one allele potentially affected the affinity of a muscle-expressed transcription factor to its binding site. 162 aseQTLs were predicted to affect 149 miRNA binding sites, from which 114 miRNAs were expressed in muscle. Also, 16 aseQTLs were methylated in our population. Integration of aseQTL with GWAS data revealed enrichment for traits such as meat tenderness, ribeye area, and intramuscular fat . To our knowledge, this is the first report of aseQTLs identification in bovine muscle. Our findings indicate that various cis-regulatory and epigenetic mechanisms can affect multiple variants to modulate the allelic expression. Some of the potential regulatory variants described here were associated with the expression pattern of genes related to interesting phenotypes for livestock. Thus, these variants might be useful for the comprehension of the genetic control of these phenotypes.

Genomic Approaches Reveal Pleiotropic Effects in Crossbred Beef Cattle

Carcass and meat quality are two important attributes for the beef industry because they drive profitability and consumer demand. These traits are of even greater importance in crossbred cattle used in subtropical and tropical regions for their superior adaptability because they tend to underperform compared to their purebred counterparts. Many of these traits are challenging and expensive to measure and unavailable until late in life or after the animal is harvested, hence unrealistic to improve through traditional phenotypic selection, but perfect candidates for genomic selection. Before genomic selection can be implemented in crossbred populations, it is important to explore if pleiotropic effects exist between carcass and meat quality traits. Therefore, the objective of this study was to identify genomic regions with pleiotropic effects on carcass and meat quality traits in a multibreed Angus-Brahman population that included purebred and crossbred animals. Data included phenotypes for 10 carcass and meat quality traits from 2,384 steers, of which 1,038 were genotyped with the GGP Bovine F-250. Single-trait genome-wide association studies were first used to investigate the relevance of direct additive genetic effects on each carcass, sensory and visual meat quality traits. A second analysis for each trait included all other phenotypes as covariates to correct for direct causal effects from identified genomic regions with pure direct effects on the trait under analysis. Five genomic windows on chromosomes BTA5, BTA7, BTA18, and BTA29 explained more than 1% of additive genetic variance of two or more traits. Moreover, three suggestive pleiotropic regions were identified on BTA10 and BTA19. The 317 genes uncovered in pleiotropic regions included anchoring and cytoskeletal proteins, key players in cell growth, muscle development, lipid metabolism and fat deposition, and important factors in muscle proteolysis. A functional analysis of these genes revealed GO terms directly related to carcass quality, meat quality, and tenderness in beef cattle, including calcium-related processes, cell signaling, and modulation of cell-cell adhesion. These results contribute with novel information about the complex genetic architecture and pleiotropic effects of carcass and meat quality traits in crossbred beef cattle.

Breed Ancestry, Divergence, Admixture, and Selection Patterns of the Simbra Crossbreed

In this study, we evaluated an admixed South African Simbra crossbred population, as well as the Brahman (Indicine) and Simmental (Taurine) ancestor populations to understand their genetic architecture and detect genomic regions showing signatures of selection. Animals were genotyped using the Illumina BovineLD v2 BeadChip (7K). Genomic structure analysis confirmed that the South African Simbra cattle have an admixed genome, composed of 5/8 Taurine and 3/8 Indicine, ensuring that the Simbra genome maintains favorable traits from both breeds. Genomic regions that have been targeted by selection were detected using the linkage disequilibrium-based methods iHS and Rsb. These analyses identified 10 candidate regions that are potentially under strong positive selection, containing genes implicated in cattle health and production (e.g., TRIM63, KCNA10, NCAM1, SMIM5, MIER3, and SLC24A4). These adaptive alleles likely contribute to the biological and cellular functions determining phenotype in the Simbra hybrid cattle breed. Our data suggested that these alleles were introgressed from the breed's original indicine and taurine ancestors. The Simbra breed thus possesses derived parental alleles that combine the superior traits of the founder Brahman and Simmental breeds. These regions and genes might represent good targets for ad-hoc physiological studies, selection of breeding material and eventually even gene editing, for improved traits in modern cattle breeds. This study represents an important step toward developing and improving strategies for selection and population breeding to ultimately contribute meaningfully to the beef production industry.

Multi-Omics Approach Reveals miR-SNPs Affecting Muscle Fatty Acids Profile in Nelore Cattle

MicroRNAs (miRNAs) are key regulators of gene expression, potentially affecting several biological processes, whose function can be altered by sequence variation. Hence, the integration of single nucleotide polymorphisms (SNP) and miRNAs can explain individual differences in economic traits. To provide new insights into the effects of SNPs on miRNAs and their related target genes, we carried out a multi-omic analysis to identify SNPs in miRNA mature sequences (miR-SNPs) associated with fatty acid (FA) composition in the Nelore cattle. As a result, we identified 3 miR-SNPs in different miRNAs (bta-miR-2419-3p, bta-miR-193a-2, and bta-miR-1291) significantly associated with FA traits (p-value < 0.02, Bonferroni corrected). Among these, the rs110817643C>T, located in the seed sequence of the bta-miR-1291, was associated with different ω6 FAs, polyunsaturated FA, and polyunsaturated:saturated FA ratios. Concerning the other two miR-SNPs, the rs43400521T>C (located in the bta-miR-2419-3p) was associated with C12:0 and C18:1 cis-11 FA, whereas the rs516857374A>G (located in the bta-miR-193a-2) was associated with C18:3 ω6 and ratio of ω6/ω3 traits. Additionally, to identify potential biomarkers for FA composition, we described target genes affected by these miR-SNPs at the mRNA or protein level. Our multi-omics analysis outlines the effects of genetic polymorphism on miRNA, and it highlights miR-SNPs and target candidate genes that control beef fatty acid composition.

Genome-wide association study between copy number variation regions and carcass- and meat-quality traits in Nellore cattle

Context Indicine breeds are the main source of beef products in tropical and subtropical regions. However, genetic improvement for carcass- and meat-quality traits in zebu cattle have been limited and genomics studies concerning structural variations that influence these traits are essential. Aim The aim of this study was to perform a genome-wide association study between copy number variation regions (CNVRs) and carcass- and meat quality-traits in Nellore cattle. Methods In total, 3794 animals, males and females included, were genotyped using a 777962 single-nucleotide polymorphism platform of BovineHD BeadChip (777k; Illumina Inc.). Of these, 1751 Nellore bulls were slaughtered at 24 months of age for further carcass beef analysis. The following traits were studied: beef tenderness, marbling, rib-eye area, backfat thickness and meat colour (lightness, redness and yellowness). The CNV detection was conducted through PennCNV software. The association analyses were performed using CNVRuler software. Key results Several identified genomic regions were linked to quantitative trait loci associated with fat deposition (FABP7) and lipid metabolism (PPARA; PLA2 family; BCHE), extracellular matrix (INS; COL10A1), contraction (SLC34A3; TRDN) and muscle development (CAPZP). The gene-enrichment analyses highlighted biological mechanisms directly related to the metabolism and synthesis of lipids and fatty acids. Conclusions The large number of potential candidate genes identified within the CNVRs, as well as the functions and pathways identified, should help better elucidate the genetic mechanisms involved in the expression of beef and carcass traits in Nellore cattle. Several CNVRs harboured genes that might have a functional impact to improve the beef and carcass traits. Implications The results obtained contribute to upgrade the sensorial and organoleptic attributes of Nellore cattle and make feasible the genetic improvement of carcass- and meat-quality traits.

Genome-wide association study and pathway analysis for fat deposition traits in nellore cattle raised in pasture-based systems

Genome-wide association study (GWAS) is a powerful tool to identify candidate genes and genomic regions underlying key biological mechanisms associated with economically important traits. In this context, the aim of this study was to identify genomic regions and metabolic pathways associated with backfat thickness (BFT) and rump fat thickness (RFT) in Nellore cattle, raised in pasture-based systems. Ultrasound-based measurements of BFT and RFT (adjusted to 18 months of age) were collected in 11,750 animals, with 39,903 animals in the pedigree file. Additionally, 1,440 animals were genotyped using the GGP-indicus 35K SNP chip, containing 33,623 SNPs after the quality control. The single-step GWAS analyses were performed using the BLUPF90 family programs. Candidate genes were identified through the Ensembl database incorporated in the BioMart tool, while PANTHER and REVIGO were used to identify the key metabolic pathways and gene networks. A total of 18 genomic regions located on 10 different chromosomes and harbouring 23 candidate genes were identified for BFT. For RFT, 22 genomic regions were found on 14 chromosomes, with a total of 29 candidate genes identified. The results of the pathway analyses showed important genes for BFT, including TBL1XR1, AHCYL2, SLC4A7, AADAT, VPS53, IDH2 and ETS1, which are involved in lipid metabolism, synthesis of cellular amino acids, transport of solutes, transport between Golgi Complex membranes, cell differentiation and cellular development. The main genes identified for RFT were GSK3β, LRP1B, EXT1, GRB2, SORCS1 and SLMAP, which are involved in metabolic pathways such as glycogen synthesis, lipid transport and homeostasis, polysaccharide and carbohydrate metabolism. Polymorphisms located in these candidate genes can be incorporated in commercial genotyping platforms to improve the accuracy of imputation and genomic evaluations for carcass fatness. In addition to uncovering biological mechanisms associated with carcass quality, the key gene pathways identified can also be incorporated in biology-driven genomic prediction methods.

Integrating Single-Step GWAS and Bipartite Networks Reconstruction Provides Novel Insights into Yearling Weight and Carcass Traits in Hanwoo Beef Cattle

Simple Summary

Hanwoo is an indigenous cattle breed in Korea and popular for meat production owing to its rapid growth and high-quality meat. Its yearling weight and carcass traits (backfat thickness, carcass weight, eye muscle area, and marbling score) are economically important for the selection of young and proven bulls. In recent decades, the advent of high throughput genotyping technologies has made it possible to perform genome-wide association studies (GWAS) for the detection of genomic regions associated with traits of economic interest in different species. In this study, we conducted a weighted single-step genome-wide association study which combines all genotypes, phenotypes and pedigree data in one step (ssGBLUP). It allows for the use of all SNPs simultaneously along with all phenotypes from genotyped and ungenotyped animals. Our results revealed 33 relevant genomic regions related to the traits of interest. Gene set enrichment analysis indicated that the identified candidate genes were related to biological processes and functional terms that were involved in growth and lipid metabolism. In conclusion, these results suggest that the incorporation of GWAS results and network analysis can help us to better understand the genetic bases underlying growth and carcass traits.

Abstract

In recent years, studies on the biological mechanisms underlying complex traits have been facilitated by innovations in high-throughput genotyping technology. We conducted a weighted single-step genome-wide association study (WssGWAS) to evaluate backfat thickness, carcass weight, eye muscle area, marbling score, and yearling weight in a cohort of 1540 Hanwoo beef cattle using BovineSNP50 BeadChip. The WssGWAS uncovered thirty-three genomic regions that explained more than 1% of the additive genetic variance, mostly located on chromosomes 6 and 14. Among the identified window regions, seven quantitative trait loci (QTL) had pleiotropic effects and twenty-six QTL were trait-specific. Significant pathways implicated in the measured traits through Gene Ontology (GO) term enrichment analysis included the following: lipid biosynthetic process, regulation of lipid metabolic process, transport or localization of lipid, regulation of growth, developmental growth, and multicellular organism growth. Integration of GWAS results of the studied traits with pathway and network analyses facilitated the exploration of the respective candidate genes involved in several biological functions, particularly lipid and growth metabolism. This study provides novel insight into the genetic bases underlying complex traits and could be useful in developing breeding schemes aimed at improving growth and carcass traits in Hanwoo beef cattle.

Genome-wide association study identified genomic regions and putative candidate genes affecting meat color traits in Nellore cattle

Single and multiple-trait GWAS were conducted to detect genomic regions and candidate genes associated with meat color traits (L*, lightness; a*, redness; b*, yellowness) in Nellore cattle. Phenotypic records of 5000 animals, and 3794 genotypes for 614,274 SNPs were used. The BLUPF90 family programs were used through single step GWAS approach. The top 10 genomic regions from single-trait GWAS explained 13.64%, 15.12% and 13% of genetic variance of L*, a* and b*, which harbored 129, 70, and 84 candidate genes, respectively. Regarding multiple-trait GWAS, the top 10 SNP windows explained 17.46%, 18.98% and 13.74% of genetic variance of L*, a* and b*, and harbored 124, 86, and 82 candidate genes, respectively. Pleiotropic effects were evidenced by the overlapping regions detected on BTA 15 and 26 associated with L* and a* (genetic correlation of -0.53), and on BTA 18 associated with a* and b* (genetic correlation of 0.60). Similar genomic regions located on BTA 2, 5, 6, and 18 were detected through single and multi-trait GWAS. Overlapped regions harbored a total of 30 functional candidate genes involved in mitochondrial activity, structural integrity of muscles, lipid oxidation, anaerobic metabolism, and muscular pH.

Genetic parameters associated with meat quality of Nellore cattle at different anatomical points of longissimus: Brazilian standards

The present study estimated genetic parameters and evaluated the genetic and phenotypic correlations between meat quality characteristics of Nellore cattle evaluated at different anatomical points of the longissimus. Data from 1329 Nellore young bulls were used to evaluate, in the 5th and 12th ribs, marbling score (MAR), shear force (SF), cooking weight losses (CWL) and intramuscular fat (IMF). In addition, the subcutaneous fat thickness was measured at the 12th rib (SFT₁₂) and between the last lumbar and the first sacral vertebrae (SFT_LR), in the separation of loin and round. Results yielded moderate heritability coefficients for evaluated characteristics, except CWL. High genetic correlations (0.61) were found between measurements of SFT₁₂ and SFT_LR. MAR, IMF and SF were evaluated at the 5th and 12th rib. Meat quality and subcutaneous fat thickness measured at different anatomical points of the longissimus are genetically correlated and can be used in genetic selection programs to improve meat quality characteristics in Nellore cattle.

The application of omics in ruminant production: a review in the tropical and sub-tropical animal production context

The demand for animal products (e.g. dairy and beef) in tropical regions is expected to increase in parallel with the public demand for sustainable practices, due to factors such as population growth and climate change. The necessity to increase animal production output must be achieved with better management and production technologies. For this to happen, novel research methodologies, animal selection and postgenomic tools play a pivotal role. Indeed, improving breeder selection programs, the quality of meat and dairy products as well as animal health will contribute to higher sustainability and productivity. This would surely benefit regions where resource quality and quantity are increasingly unstable, and research is still very incipient, which is the case of many regions in the tropics. The purpose of this review is to demonstrate how omics-based approaches play a major role in animal science, particularly concerning ruminant production systems and research associated to the tropics and developing countries. SIGNIFICANCE: Environmental conditions in the tropics make livestock production harder, compared to temperate regions. Due to global warming, the sustainability of livestock production will become increasingly problematic. The use of novel omics technologies could generate useful information to understand adaptation mechanisms of resilient breeds and/or species. The application of omics to tropical animal production is still residual in the currently available literature. With this review, we aim to summarize the most notable results in the field whilst encouraging further research to deal with the future challenges that animal production in the tropics will need to face.

Allele-specific expression is widespread in Bos indicus muscle and affects meat quality candidate genes

Differences between the expression of the two alleles of a gene are known as allele-specific expression (ASE), a common event in the transcriptome of mammals. Despite ASE being a source of phenotypic variation, its occurrence and effects on genetic prediction of economically relevant traits are still unexplored in bovines. Furthermore, as ASE events are likely driven by cis-regulatory mutations, scanning them throughout the bovine genome represents a significant step to elucidate the mechanisms underlying gene expression regulation. To address this question in a Bos indicus population, we built the ASE profile of the skeletal muscle tissue of 190 Nelore steers, using RNA sequencing data and SNPs genotypes from the Illumina BovineHD BeadChip (770 K bp). After quality control, 820 SNPs showed at least one sample with ASE. These SNPs were widespread among all autosomal chromosomes, being 32.01% found in 3'UTR and 31.41% in coding regions. We observed a considerable variation of ASE profile among individuals, which highlighted the need for biological replicates in ASE studies. Functional analysis revealed that ASE genes play critical biological functions in the development and maintenance of muscle tissue. Additionally, some of these genes were previously reported as associated with beef production and quality traits in livestock, thus indicating a possible source of bias on genomic predictions for these traits.

Signatures of positive selection underlying beef production traits in Korean cattle breeds

The difference in the breeding programs and population history may have diversely shaped the genomes of Korean native cattle breeds. In the absence of phenotypic data, comparisons of breeds that have been subjected to different selective pressures can aid to identify genomic regions and genes controlling qualitative and complex traits. In this study to decipher genetic variation and identify evidence of divergent selection, 3 Korean cattle breeds were genotyped using the recently developed high-density GeneSeek Genomic Profiler F250 (GGP-F250) array. The three Korean cattle breeds clustered according to their coat color phenotypes and breeding programs. The Heugu breed reliably showed smaller effective population size at all generations considered. Across the autosomal chromosomes, 113 and 83 annotated genes were identified from Hanwoo-Chikso and Hanwoo-Heugu comparisons, respectively of which 16 genes were shared between the two pairwise comparisons. The most important signals of selection were detected on bovine chromosomes 14 (24.39–25.13 Mb) and 18 (13.34–15.07 Mb), containing genes related to body size, and coat color (XKR4, LYN, PLAG1, SDR16C5, TMEM68, CDH15, MC1R, and GALNS). Some of the candidate genes are also associated with meat quality traits (ACSF3, EIF2B1, BANP, APCDD1, and GALM) and harbor quantitative trait locus (QTL) for beef production traits. Further functional analysis revealed that the candidate genes (DBI, ACSF3, HINT2, GBA2, AGPAT5, SCAP, ELP6, APOB, and RBL1) were involved in gene ontology (GO) terms relevant to meat quality including fatty acid oxidation, biosynthesis, and lipid storage. Candidate genes previously known to affect beef production and quality traits could be used in the beef cattle selection strategies.

Proteome alterations associated with the oleic acid and cis-9, trans-11 conjugated linoleic acid content in bovine skeletal muscle

Oleic acid (OA) and cis-9, trans-11 conjugated linoleic acid (c9t11-CLA) are fatty acids found in beef with beneficial effects in human health. This study investigated differentially abundant proteins (DAPs) in skeletal muscle of bovines with extreme values of OA, and c9t11-CLA. For each one of the fatty acids, twenty muscle samples were divided into two groups (N = 10_High; N = 10_Low) and analyzed by high definition mass spectrometry. We identified 103 and 133 DAPs between the groups for each fatty acid. We found 64 and 45 up-regulated and 39 and 68 down-regulated proteins for OA and c9t11-CLA, respectively. Comparative analysis between proteomic and transcriptomic data revealed eight and ten genes with a consistent between mRNA expression levels and protein abundance for OA and c9t11-CLA, respectively. Unconventional myosin-Id (MYO1D), mineralocorticoid receptor (NR3C2), geranylgeranyl transferase type-2 subunit-alpha (RABGGTA), and uveal autoantigen with coiled-coil domains and ankyrin repeats (UACA) were found as putative candidate proteins for OA content. Fatty acid synthase (FASN), tubulin alpha-4A chain (TUBA4A), vinculin (VCL), NADH dehydrogenase 1 alpha subcomplex 5 (NDUFA5), and prefoldin subunit 6 (PFDN6) for c9t11-CLA. Our findings contribute to a deeper understanding of the molecular mechanisms behind the regulation of the OA and c9t11-CLA content in cattle skeletal muscle. SIGNIFICANCE: Questions about the association between meat intake and disease incidence in humans has driven animal scientist to pursue a better understanding of the biological processes associated with differences in the intramuscular fat composition. The beneficial effects of oleic acid and conjugated linoleic acid in human health have been demonstrated by improving the immune system and preventing atherosclerosis, different types of cancers, hypertension, and diabetes. Previous genome-wide association and gene expression studies identified genomic regions and differentially expressed genes associated with the fatty acid profile in skeletal muscle. In this work, differences were evaluated at the protein level. The use of a label-free quantitative proteomic approach, compared with muscle transcriptome results obtained by RNA-sequencing, allowed us to earn new insights into the variability in fatty acid deposition in skeletal muscle of farm animals. This study opens new avenues to explore the effect of the fatty acids in the skeletal muscle of livestock animals, which is associated with nutritional values of the meat, and perhaps to understand the mechanisms correlated with metabolic diseases in other species.

Use of gene expression profile to identify potentially relevant transcripts to myofibrillar fragmentation index trait

The myofibrillar fragmentation index (MFI) is an indicative trait for meat tenderness. Longissimus thoracis muscle samples from the 20 most extreme bulls (out of 80 bulls set) for MFI (high (n = 10) and low (n = 10) groups) trait were used to perform transcriptomic analysis, using RNA Sequencing (RNA-Seq). An average of 24.616 genes was expressed in the Nellore muscle transcriptome analysis. A total of 96 genes were differentially expressed (p value ≤ 0.001) between the two groups of divergent bulls for MFI. The HEBP2 and BDH1 genes were overexpressed in animals with high MFI. The MYBPH and MYL6, myosin encoders, were identified. The differentially expressed genes were related to increase mitochondria efficiency, especially in cells under oxidative stress conditions, and these also were related to zinc and calcium binding, membrane transport, and muscle constituent proteins, such as actin and myosin. Most of those genes were involved in metabolic pathways of oxidation-reduction, transport of lactate in the plasma membrane, and muscle contraction. This is the first study applying MFI phenotypes in transcriptomic studies to identify and understand differentially expressed genes for beef tenderness. These results suggest that differences detected in gene expression between high and low MFI animals are related to reactive mechanisms and structural components of oxidative fibers under the condition of cellular stress. Some genes may be selected as positional candidate genes to beef tenderness, MYL6, MYBPH, TRIM63, TRIM55, TRIOBP, and CHRNG genes. The use of MFI phenotypes could enhance results of meat tenderness studies.

Genetic analysis of carcass and meat quality traits in Nelore cattle1

The objective of this study was to estimate genetic parameters for carcass and meat quality traits, as well as their genetic correlations using pedigree and genomic information. A total of 3,716; 3,702; 3,439; 3,705; and 3,714 records of 12th-13th rib LM area (LMA), backfat thickness (BF), HCW, marbling score (MARB), and Warner-Bratzler peak shear force (WBSF), respectively, were used. Animals were genotyped with BovineHD BeadChip and GeneSeek Genomic Profiler Indicus HD - GGP75Ki panel. The (co)variance components were estimated by Bayesian inference using a multitrait ssGBLUP analysis. The animal model included fixed effects of contemporary group (defined by the combination of farm and year of birth, and management group at yearling) and age of animal at slaughtering as a covariate (linear). Direct additive genetic and residual effects were fitted as random. The posterior means and SD of heritabilities for LMA, BF, HCW, MARB, and WBSF were 0.28 (0.03), 0.21 (0.04), 0.21 (0.04), 0.12 (0.04), and 0.11 (0.03), respectively. The posterior means for genetic correlations between LMA and meat quality were positive and moderate with MARB (0.38 ± 0.12) and negative with WBSF (-0.47 ± 0.12). Low genetic correlations were estimated between BF and WBSF (-0.03 ± 0.16) and between HCW and MARB (-0.04 ± 0.14), indicating that these traits are not controlled by the same set or linked genes. Carcass traits (LMA, BF, and HCW) presented moderate heritability providing quick response to the selection purpose. The estimates of heritability for meat quality traits (MARB and WBSF) were low and indicate that the rate of genetic improvement for these traits would be slow. Genetic correlations indicated that selection for carcass traits would not be strongly antagonistic for improving meat quality.

A Gene-Set Enrichment and Protein-Protein Interaction Network-Based GWAS with Regulatory SNPs Identifies Candidate Genes and Pathways Associated with Carcass Traits in Hanwoo Cattle

Non-synonymous SNPs and protein coding SNPs within the promoter region of genes (regulatory SNPs) might have a significant effect on carcass traits. Imputed sequence level data of 10,215 Hanwoo bulls, annotated and filtered to include only regulatory SNPs (450,062 SNPs), were used in a genome-wide association study (GWAS) to identify loci associated with backfat thickness (BFT), carcass weight (CWT), eye muscle area (EMA), and marbling score (MS). A total of 15, 176, and 1 SNPs were found to be significantly associated (p < 1.11 × 10-7) with BFT, CWT, and EMA, respectively. The significant loci were BTA4 (CWT), BTA6 (CWT), BTA14 (CWT and EMA), and BTA19 (BFT). BayesR estimated that 1.1%~1.9% of the SNPs contributed to more than 0.01% of the phenotypic variance. So, the GWAS was complemented by a gene-set enrichment (GSEA) and protein-protein interaction network (PPIN) analysis in identifying the pathways affecting carcass traits. At p < 0.005 (~2,261 SNPs), 25 GO and 18 KEGG categories, including calcium signaling, cell proliferation, and folate biosynthesis, were found to be enriched through GSEA. The PPIN analysis showed enrichment for 81 candidate genes involved in various pathways, including the PI3K-AKT, calcium, and FoxO signaling pathways. Our finding provides insight into the effects of regulatory SNPs on carcass traits.

Association of growth traits with a structural variation downstream of the KCNJ11 gene: a large population-based study in chickens

1. The potassium voltage-gated channel subfamily J member 11 gene (KCNJ11) is involved in the insulin secretion pathway. Studies have shown that mutation in this gene is associated with muscle weakness. The objective of the present study was to establish the association between KCNJ11 gene polymorphism and chicken growth performance and to analyse its expression pattern. 2. A novel 163-bp insertion/deletion (indel) polymorphism was identified in the region downstream of the KCNJ11 gene in 2330 individuals from ten populations by polymerase chain reaction (PCR). An F₂ resource population was used to investigate the genetic effects of the chicken KCNJ11 gene. Association analysis showed that the indel was significantly associated with chicken growth traits and that the phenotypic value of the ins-ins (II) genotype is higher than that of the ins-del (ID) and del-del (DD) genotypes. 3. Gene expression for different genotypes showed that birds carrying the II allele had a higher expression level than the DD genotypes. Analysis of tissue and spatiotemporal expression patterns indicated that the KCNJ11 gene was highly expressed in muscle tissues, with the highest levels in muscle tissue at one week of age, and that a 10% crude protein diet reduced the expression of this gene, average daily gain and muscle fibre diameter. 4. The results suggested that this novel 163-bp indel has the potential to become a new target for marker-assisted selection.

Aggregation of Omic Data and Secretome Prediction Enable the Discovery of Candidate Plasma Biomarkers for Beef Tenderness

Beef quality is a complex phenotype that can be evaluated only after animal slaughtering. Previous research has investigated the potential of genetic markers or muscle-derived proteins to assess beef tenderness. Thus, the use of low-invasive biomarkers in living animals is an issue for the beef sector. We hypothesized that publicly available data may help us discovering candidate plasma biomarkers. Thanks to a review of the literature, we built a corpus of articles on beef tenderness. Following data collection, aggregation, and computational reconstruction of the muscle secretome, the putative plasma proteins were searched by comparison with a bovine plasma proteome atlas and submitted to mining of biological information. Of the 44 publications included in the study, 469 unique gene names were extracted for aggregation. Seventy-one proteins putatively released in the plasma were revealed. Among them 13 proteins were predicted to be secreted in plasma, 44 proteins as hypothetically secreted in plasma, and 14 additional candidate proteins were detected thanks to network analysis. Among these 71 proteins, 24 were included in tenderness quantitative trait loci. The in-silico workflow enabled the discovery of candidate plasma biomarkers for beef tenderness from reconstruction of the secretome, to be examined in the cattle plasma proteome.