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.

Identification of eQTLs using different sets of single nucleotide polymorphisms associated with carcass and body composition traits in pigs.

BACKGROUND: Mapping expression quantitative trait loci (eQTLs) in skeletal muscle tissue in pigs is crucial for understanding the relationship between genetic variation and phenotypic expression of carcass traits in meat animals. Therefore, the primary objective of this study was to evaluate the impact of different sets of single nucleotide polymorphisms (SNP), including scenarios removing SNPs pruned for linkage disequilibrium (LD) and SNPs derived from SNP chip arrays and RNA-seq data from liver, brain, and skeletal muscle tissues, on the identification of eQTLs in the Longissimus lumborum tissue, associated with carcass and body composition traits in Large White pigs. The SNPs identified from muscle mRNA were combined with SNPs identified in the brain and liver tissue transcriptomes, as well as SNPs from the GGP Porcine 50 K SNP chip array. Cis- and trans-eQTLs were identified based on the skeletal muscle gene expression level, followed by functional genomic analyses and statistical associations with carcass and body composition traits in Large White pigs. RESULTS: The number of cis- and trans-eQTLs identified across different sets of SNPs (scenarios) ranged from 261 to 2,539 and from 29 to 13,721, respectively. Furthermore, 6,180 genes were modulated by eQTLs in at least one of the scenarios evaluated. The eQTLs identified were not significantly associated with carcass and body composition traits but were significantly enriched for many traits in the "Meat and Carcass" type QTL. The scenarios with the highest number of cis- (n = 304) and trans- (n = 5,993) modulated genes were the unpruned and LD-pruned SNP set scenarios identified from the muscle transcriptome. These genes include 84 transcription factor coding genes. CONCLUSIONS: After LD pruning, the set of SNPs identified based on the transcriptome of the skeletal muscle tissue of pigs resulted in the highest number of genes modulated by eQTLs. Most eQTLs are of the trans type and are associated with genes influencing complex traits in pigs, such as transcription factors and enhancers. Furthermore, the incorporation of SNPs from other genomic regions to the set of SNPs identified in the porcine skeletal muscle transcriptome contributed to the identification of eQTLs that had not been identified based on the porcine skeletal muscle transcriptome alone.

Transcriptome profile of skeletal muscle using different sources of dietary fatty acids in male pigs.

Pork is of great importance in world trade and represents the largest source of fatty acids in the human diet. Lipid sources such as soybean oil (SOY), canola (CO), and fish oil (FO) are used in pig diets and influence blood parameters and the ratio of deposited fatty acids. In this study, the main objective was to evaluate changes in gene expression in porcine skeletal muscle tissue resulting from the dietary oil sources and to identify metabolic pathways and biological process networks through RNA-Seq. The addition of FO in the diet of pigs led to intramuscular lipid with a higher FA profile composition of C20:5 n-3, C22:6 n-3, and SFA (C16:0 and C18:0). Blood parameters for the FO group showed lower cholesterol and HDL content compared with CO and SOY groups. Skeletal muscle transcriptome analyses revealed 65 differentially expressed genes (DEG, FDR 10%) between CO vs SOY, and 32 DEG for CO vs FO, and 531 DEG for SOY vs FO comparison. Several genes, including AZGP1, PDE3B, APOE, PLIN1, and LIPS, were found to be down-regulated in the diet of the SOY group compared to the FO group. The enrichment analysis revealed DEG involved in lipid metabolism, metabolic diseases, and inflammation between the oil groups, with specific gene functions in each group and altered blood parameters. The results provide mechanisms to help us understand the behavior of genes according to fatty acids.

Brain fatty acid and transcriptome profiles of pig fed diets with different levels of soybean oil.

BACKGROUND: The high similarity in anatomical and neurophysiological processes between pigs and humans make pigs an excellent model for metabolic diseases and neurological disorders. Lipids are essential for brain structure and function, and the polyunsaturated fatty acids (PUFA) have anti-inflammatory and positive effects against cognitive dysfunction in neurodegenerative diseases. Nutrigenomics studies involving pigs and fatty acids (FA) may help us in better understanding important biological processes. In this study, the main goal was to evaluate the effect of different levels of dietary soybean oil on the lipid profile and transcriptome in pigs' brain tissue. RESULTS: Thirty-six male Large White pigs were used in a 98-day study using two experimental diets corn-soybean meal diet containing 1.5% soybean oil (SOY1.5) and corn-soybean meal diet containing 3.0% soybean oil (SOY3.0). No differences were found for the brain total lipid content and FA profile between the different levels of soybean oil. For differential expression analysis, using the DESeq2 statistical package, a total of 34 differentially expressed genes (DEG, FDR-corrected p-value < 0.05) were identified. Of these 34 DEG, 25 are known-genes, of which 11 were up-regulated (log2 fold change ranging from + 0.25 to + 2.93) and 14 were down-regulated (log2 fold change ranging from - 3.43 to -0.36) for the SOY1.5 group compared to SOY3.0. For the functional enrichment analysis performed using MetaCore with the 34 DEG, four pathway maps were identified (p-value < 0.05), related to the ALOX15B (log2 fold change - 1.489), CALB1 (log2 fold change - 3.431) and CAST (log2 fold change + 0.421) genes. A "calcium transport" network (p-value = 2.303e-2), related to the CAST and CALB1 genes, was also identified. CONCLUSION: The results found in this study contribute to understanding the pathways and networks associated with processes involved in intracellular calcium, lipid metabolism, and oxidative processes in the brain tissue. Moreover, these results may help a better comprehension of the modulating effects of soybean oil and its FA composition on processes and diseases affecting the brain tissue.

RNA-seq transcriptome profiling of pigs' liver in response to diet with different sources of fatty acids.

Pigs (Sus scrofa) are an animal model for metabolic diseases in humans. Pork is an important source of fatty acids (FAs) in the human diet, as it is one of the most consumed meats worldwide. The effects of dietary inclusion of oils such as canola, fish, and soybean oils on pig gene expression are mostly unknown. Our objective was to evaluate FA composition, identify changes in gene expression in the liver of male pigs fed diets enriched with different FA profiles, and identify impacted metabolic pathways and gene networks to enlighten the biological mechanisms' variation. Large White male pigs were randomly allocated to one of three diets with 18 pigs in each; all diets comprised a base of corn and soybean meal to which either 3% of soybean oil (SOY), 3% canola oil (CO), or 3% fish oil (FO) was added for a 98-day trial during the growing and finishing phases. RNA sequencing was performed on the liver samples of each animal by Illumina technology for differential gene expression analyses, using the R package DESeq2. The diets modified the FA profile, mainly in relation to polyunsaturated and saturated FAs. Comparing SOY vs. FO, 143 differentially expressed genes (DEGs) were identified as being associated with metabolism, metabolic and neurodegenerative disease pathways, inflammatory processes, and immune response networks. Comparing CO vs. SOY, 148 DEGs were identified, with pathways related to FA oxidation, regulation of lipid metabolism, and metabolic and neurodegenerative diseases. Our results help explain the behavior of genes with differential expression in metabolic pathways resulting from feeding different types of oils in pig diets.

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.

Differential Gene Expression Associated with Soybean Oil Level in the Diet of Pigs.

The aim of this study was to identify the differentially expressed genes (DEG) from the skeletal muscle and liver samples of animal models for metabolic diseases in humans. To perform the study, the fatty acid (FA) profile and RNA sequencing (RNA-Seq) data of 35 samples of liver tissue (SOY1.5, n = 17 and SOY3.0, n = 18) and 36 samples of skeletal muscle (SOY1.5, n = 18 and SOY3.0, n = 18) of Large White pigs were analyzed. The FA profile of the tissues was modified by the diet, mainly those related to monounsaturated (MUFA) and polyunsaturated (PUFA) FA. The skeletal muscle transcriptome analysis revealed 45 DEG (FDR 10%), and the functional enrichment analysis identified network maps related to inflammation, immune processes, and pathways associated with oxidative stress, type 2 diabetes, and metabolic dysfunction. For the liver tissue, the transcriptome profile analysis revealed 281 DEG, which participate in network maps related to neurodegenerative diseases. With this nutrigenomics study, we verified that different levels of soybean oil in the pig diet, an animal model for metabolic diseases in humans, affected the transcriptome profile of skeletal muscle and liver tissue. These findings may help to better understand the biological mechanisms that can be modulated by the diet.

Effect of dietary soybean oil inclusion on liver-related transcription factors in a pig model for metabolic diseases.

Dietary fatty acids (FA) are components of the lipids, which contribute to membrane structure, energy input, and biological functions related to cellular signaling and transcriptome regulation. However, the consumers still associate dietary FA with fat deposition and increased occurrence of metabolic diseases such as obesity and atherosclerosis. Previous studies already demonstrated that some fatty acids are linked with inflammatory response, preventing metabolic diseases. To better understand the role of dietary FA on metabolic diseases, for the first time, a study to identify key transcription factors (TF) involved in lipid metabolism and inflammatory response by transcriptome analysis from liver samples of animal models was performed. The key TF were identified by functional enrichment analysis from the list of differentially expressed genes identified in liver samples between 35 pigs fed with 1.5% or 3.0% soybean oil. The functional enrichment analysis detected TF linked to lipid homeostasis and inflammatory response, such as RXRA, EGFR, and SREBP2 precursor. These findings demonstrated that key TF related to lipid metabolism could be modulated by dietary inclusion of soybean oil. It could contribute to nutrigenomics research field that aims to elucidate dietary interventions in animal and human health, as well as to drive food technology and science.

Effect of different selenium sources and concentrations on glutathione peroxidase activity and cholesterol metabolism of beef cattle.

The objective of this study was to investigate the effects of different Se sources and concentrations on glutathione forms and cholesterol metabolism in beef cattle. Sixty-three Nellore bulls (412 ± 19 kg body weight (BW); 24 mo old) were randomly assigned to a completely randomized design in a 2 × 3 + 1 factorial arrangement (63 pens; one animal/pen) with two Se sources (sodium selenite, ING and Se-yeast, ORG), three concentrations (0.3, 0.9, and 2.7 mg supplemental Se/kg dry matter (DM)), and control treatment (without Se supplementation) fed for 90 d. Blood samples were collected on day 0, 28, 56, and 84. Muscle and liver samples were collected at harvest. Hepatic GSSG (P = 0.004), GSH/GSSG ratio (P = 0.030), and GSH-Px (P = 0.004) were affected by Se source × concentration interaction. Oxidized glutathione was higher in the ORG group vs. ING at concentration 2.7 mg supplemental Se/kg DM, but at 0.3 mg supplemental Se/kg DM the ING group was higher than ORG. The liver GSH-Px activity was higher in the ORG group vs. ING at concentration 0.9 and 2.7 mg supplemental Se/kg DM. The GSH/GSSG ratio was the highest in animals fed 0.3 mg supplemental Se/kg DM of ORG. Selenium liver concentration increased linearly with the supplemental Se concentration in the diet (y = 0.0583 + 0.4254x, R2 = 0.92, P < 0.0001), regardless of source. Total meat cholesterol was greater (P < 0.001) in CON (control) vs. SUP (supplemented, regardless source) group. The muscle GSH-Px activity was higher (P < 0.001) in SUP vs. CON and increased (P < 0.004) with increasing supplemental Se concentrations. There was an increase on very low-density lipoprotein (VLDL), glucose, and triglycerides in ORG vs. ING (P ≤ 0.035). In general, serum Se was higher (P < 0.001) in SUP vs. CON and increased with increasing supplemental Se concentration. Lastly, the 3-hydroxy-3-methyl-glutaryl-coenzyme A reductase (HMGCR) concentration was lower (P = 0.002) in SUP (0.39 ng/mL) vs. CON (0.55 ng/mL). Selenium supplementation with different sources and concentrations has the potential to affect cholesterol metabolism by affecting GSH/GSSG ratio, GSH-Px, and the HMGCR.

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.

Gene co-expression networks associated with carcass traits reveal new pathways for muscle and fat deposition in Nelore cattle.

BACKGROUND: Positively correlated with carcass weight and animal growth, the ribeye area (REA) and the backfat thickness (BFT) are economic important carcass traits, which impact directly on producer's payment. The selection of these traits has not been satisfactory since they are expressed later in the animal's life and multigene regulated. So, next-generation technologies have been applied in this area to improve animal's selection and better understand the molecular mechanisms involved in the development of these traits. Correlation network analysis, performed by tools like WGCNA (Weighted Correlation Network Analysis), has been used to explore gene-gene interactions and gene-phenotype correlations. Thus, this study aimed to identify putative candidate genes and metabolic pathways that regulate REA and BFT by constructing a gene co-expression network using WGCNA and RNA sequencing data, to better understand genetic and molecular variations behind these complex traits in Nelore cattle. RESULTS: The gene co-expression network analysis, using WGCNA, were built using RNA-sequencing data normalized by transcript per million (TPM) from 43 Nelore steers. Forty-six gene clusters were constructed, between them, three were positively correlated (p-value< 0.1) to the BFT (Green Yellow, Ivory, and Light Yellow modules) and, one cluster was negatively correlated (p-value< 0.1) with REA (Salmon module). The enrichment analysis performed by DAVID and WebGestalt (FDR 5%) identified eight Gene Ontology (GO) terms and three KEGG pathways in the Green Yellow module, mostly associated with immune response and inflammatory mechanisms. The enrichment of the Salmon module demonstrated 19 GO terms and 21 KEGG pathways, related to muscle energy metabolism, lipid metabolism, muscle degradation, and oxidative stress diseases. The Ivory and Light yellow modules have not shown significant results in the enrichment analysis. CONCLUSION: With this study, we verified that inflammation and immune response pathways modulate the BFT trait. Energy and lipid metabolism pathways, highlighting fatty acid metabolism, were the central pathways associated with REA. Some genes, as RSAD2, EIF2AK2, ACAT1, and ACSL1 were considered as putative candidate related to these traits. Altogether these results allow us to a better comprehension of the molecular mechanisms that lead to muscle and fat deposition in bovine.

Data from proteomic analysis of bovine Longissimus dorsi muscle associated with intramuscular fat content.

The proteomic data presented in this article are associated with the research article entitled "Longissimus dorsi muscle label-free quantitative proteomic reveals biological mechanisms associated with intramuscular fat deposition" published in Journal of Proteomics [1]. In this article, we characterized the proteomic profile of bovine Longissimus dorsi muscle from Nelore steers and identified differentially abundant proteins associated with the intramuscular fat (IMF) content. An integrated transcriptome-assisted label-free quantitative proteomic approach by High Definition Mass Spectrometry (HDMSE) was employed to identify and quantify the proteins. A functional enrichment analysis using the differentially abundant proteins list was performed to understand the biological processes involved in IMF deposition. Moreover, to explore and clarify the biological mechanisms that influence IMF content, the mRNA data for the same trait from Cesar and collaborators [2] obtained by RNA-sequencing technology was compared with proteomic data. The mRNA data is deposited in the European Nucleotide Archive (ENA) repository (EMBL-EBI), under accession PRJEB13188.

Label-free quantitative proteomic analysis reveals muscle contraction and metabolism proteins linked to ultimate pH in bovine skeletal muscle.

The purpose of this research was to investigate the causes and consequences of pHu variations in beef cattle. A group of 176 Nellore beef cattle was evaluated and classified into two different pHu groups: High (≥6.0, N = 17) and Normal (<5.8, N = 159). Plasma concentrations of cortisol and adrenocorticotropic hormone, lactate and glycogen muscular content, meat color, shear force and Longissimus thoracis muscle proteomic profile were evaluated and compared between pHu groups. Muscle glycogen content, meat color and shear force statistically differed between pHu groups. Label-free quantitative proteomic analysis revealed ten differentially abundant proteins between pHu groups, involved in metabolic processes and muscle contraction, which also were significantly correlated with pHu. Thirty-six and 31 proteins were exclusively present in Normal and High pHu group, respectively, which were related to TCA cycle, cortisol production, calcium regulation, and antioxidant function. The MYH7, UGP2, H2AFJ and VDAC3 were identified as potential indicators of pHu variations. CALM and NNT appeared to be interesting proteins to understand the metabolic pathways behind pHu. Data are available via ProteomeXchange with identifier PXD009320.

Longissimus dorsi muscle label-free quantitative proteomic reveals biological mechanisms associated with intramuscular fat deposition.

The pathways involved in intramuscular fat (IMF) deposition in Longissimus dorsi muscle were investigated using an integrated transcriptome-assisted label-free quantitative proteomic approach by High Definition Mass Spectrometry. We quantified 1582 proteins, of which 164 were differentially abundant proteins (DAPs, p < 0.05) between animals with high (H) and low (L) genomic estimated breeding values (GEBV) for IMF content. Ingenuity pathway analysis (IPA) revealed that these DAPs were mainly involved in glycolysis metabolism, actin cytoskeleton signaling, cell-cell adherens junction and pathways for MAPK and insulin. A comparative study between transcriptomic (mRNA) and proteomic data showed 17 differentially expressed genes corresponding to DAPs, of which three genes/proteins did not agree on the direction of the fold change between groups. Moreover, we investigated microRNAs data to explain these differences in fold change direction, being able to unravel two of the three unexpected mRNA/protein relationships. Results demonstrated that changes in protein/mRNA levels of sarcomere organization, intracellular signal transduction and regulation of actin cytoskeleton, are involved in IMF deposition. These findings provide a deeper understanding of the highly complex regulatory mechanisms involved in IMF deposition in cattle and indicate target pathways for future studies. SIGNIFICANCE: Intramuscular fat is the amount of fat deposited inside muscle and plays an important role in human health and meat quality attributes, influencing energy metabolism of skeletal muscle, as well as, tenderness, flavor, and juiciness of beef. We performed for the first time the utilization of integrated transcriptome-assisted label-free quantitative proteomic approach using High Definition Mass Spectrometry for characterization of the changes in the proteomic profile of the Longissimus dorsi muscle associated with intramuscular fat deposition in cattle. Furthermore, we compared the muscle proteome with the muscle transcriptome (mRNA and microRNAs), obtained by RNA-sequencing, to better understand the relationship between expression of mRNAs and proteins and to unravel essential biological mechanisms involved in bovine skeletal muscle IMF deposition.

Comparative muscle transcriptome associated with carcass traits of Nellore cattle.

BACKGROUND: Commercial cuts yield is an important trait for beef production, which affects the final value of the products, but its direct determination is a challenging procedure to be implemented in practice. The measurement of ribeye area (REA) and backfat thickness (BFT) can be used as indirect measures of meat yield. REA and BFT are important traits studied in beef cattle due to their strong implication in technological (carcass yield) and nutritional characteristics of meat products, like the degree of muscularity and total body fat. Thus, the aim of this work was to study the Longissimus dorsi muscle transcriptome of Nellore cattle, associated with REA and BFT, to find differentially expressed (DE) genes, metabolic pathways, and biological processes that may regulate these traits. RESULTS: By comparing the gene expression level between groups with extreme genomic estimated breeding values (GEBV), 101 DE genes for REA and 18 for BFT (false discovery rate, FDR 10%) were identified. Functional enrichment analysis for REA identified two KEGG pathways, MAPK (Mitogen-Activated Protein Kinase) signaling pathway and endocytosis pathway, and three biological processes, response to endoplasmic reticulum stress, cellular protein modification process, and macromolecule modification. The MAPK pathway is responsible for fundamental cellular processes, such as growth, differentiation, and hypertrophy. For BFT, 18 biological processes were found to be altered and grouped into 8 clusters of semantically similar terms. The DE genes identified in the biological processes for BFT were ACHE, SRD5A1, RSAD2 and RSPO3. RSAD2 has been previously shown to be associated with lipid droplet content and lipid biosynthesis. CONCLUSION: In this study, we identified genes, metabolic pathways, and biological processes, involved in differentiation, proliferation, protein turnover, hypertrophy, as well as adipogenesis and lipid biosynthesis related to REA and BFT. These results enlighten some of the molecular processes involved in muscle and fat deposition, which are economically important carcass traits for beef production.