Meat quality is a complex trait that is influenced by genetic and environmental factors, which includes mineral concentration. However, the association between mineral concentration and meat quality, and the specific molecular pathways underlying this association, are not well explored. We therefore analyzed gene expression as measured with RNA-seq in Longissimus thoracis muscle of 194 Nelore steers for association with three meat quality traits (intramuscular fat, meat pH, and tenderness) and the concentration of 13 minerals (Ca, Cr, Co, Cu, Fe, K, Mg, Mn, Na, P, S, Se, and Zn). We identified seven sets of co-expressed genes (modules) associated with at least two traits, which indicates that common pathways influence these traits. From pathway analysis of module hub genes, we further found an over-representation for energy and protein metabolism (AMPK and mTOR signaling pathways) in addition to muscle growth, and protein turnover pathways. Among the identified hub genes FASN, ELOV5, and PDE3B are involved with lipid metabolism and were affected by previously identified eQTLs associated to fat deposition. The reported hub genes and over-represented pathways provide evidence of interplay among gene expression, mineral concentration, and meat quality traits. Future studies investigating the effect of different levels of mineral supplementation in the gene expression and meat quality traits could help us to elucidate the regulatory mechanism by which the genes/pathways are affected.

Feed efficiency helps to reduce environmental impacts from livestock production, improving beef cattle profitability. We identified potential biomarkers (hub genes) for feed efficiency, by applying co-expression analysis in Longissimus thoracis RNA-Seq data from 180 Nelore steers. Six co-expression modules were associated with six feed efficiency-related traits (p-value ≤ 0.05). Within these modules, 391 hub genes were enriched for pathways as protein synthesis, muscle growth, and immune response. Trait-associated transcription factors (TFs) ELF1, ELK3, ETS1, FLI1, and TCF4, were identified with binding sites in at least one hub gene. Gene expression of CCDC80, FBLN5, SERPINF1, and OGN was associated with multiple feed efficiency-related traits (FDR ≤ 0.05) and were previously related to glucose homeostasis, oxidative stress, fat mass, and osteoblastogenesis, respectively. Potential regulatory elements were identified, integrating the hub genes with previous studies from our research group, such as the putative cis-regulatory elements (eQTLs) inferred as affecting the PCDH18 and SPARCL1 hub genes related to immune system and adipogenesis, respectively. Therefore, our analyses contribute to a better understanding of the biological mechanisms underlying feed efficiency in bovine and the hub genes disclosed can be used as biomarkers for feed efficiency-related traits in Nelore cattle.

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.

Mineral contents in bovine muscle can affect meat quality, growth, health, and reproductive traits. To better understand the genetic basis of this phenotype in Nelore (Bosindicus) cattle, we analysed genome-wide mRNA and miRNA expression data from 114 muscle samples. The analysis implemented a new application for two complementary algorithms: the partial correlation and information theory (PCIT) and the regulatory impact factor (RIF), in which we included the estimated genomic breeding values (GEBVs) for the phenotypes additionally to the expression levels, originally proposed for these methods. We used PCIT to determine putative regulatory relationships based on significant associations between gene expression and GEBVs for each mineral amount. Then, RIF was adopted to determine the regulatory impact of genes and miRNAs expression over the GEBVs for the mineral amounts. We also investigated over-represented pathways, as well as pieces of evidences from previous studies carried in the same population and in the literature, to determine regulatory genes for the mineral amounts. For example,NOX1expression level was positively correlated to Zinc and has been described as Zinc-regulated in humans. Based on our approach, we were able to identify genes, miRNAs and pathways not yet described as underlying mineral amount. The results support the hypothesis that extracellular matrix interactions are the core regulator of mineral amount in muscle cells. Putative regulators described here add information to this hypothesis, expanding the knowledge on molecular relationships between gene expression and minerals.

Mineral content affects the biological processes underlying beef quality. Muscle mineral concentration depends not only on intake-outtake balance and muscle type, but also on age, environment, breed, and genetic factors. To unveil the genetic factors involved in muscle mineral concentration, we applied a pairwise differential gene expression analysis in groups of Nelore steers genetically divergent for nine different mineral concentrations. Here, based on significant expression differences between contrasting groups, we presented candidate genes for the genetic regulation of mineral concentration in muscle. Functional enrichment and protein-protein interaction network analyses were carried out to search for gene regulatory processes concerning each mineral. The core genetic regulation for all minerals studied, except Zn, seems to rest on interactions between components of the extracellular matrix. Regulation of adipogenesis-related pathways was also significant in our results. Antagonistic patterns of gene expression for fatty acid metabolism-related genes may explain the Cu and Zn antagonistic effect on fatty acid accumulation. Our results shed light on the role of these minerals on cell function.

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.