- Research article
- Open Access
Differences in the skeletal muscle transcriptome profile associated with extreme values of fatty acids content
BMC Genomics volume 17, Article number: 961 (2016)
Lipids are a class of molecules that play an important role in cellular structure and metabolism in all cell types. In the last few decades, it has been reported that long-chain fatty acids (FAs) are involved in several biological functions from transcriptional regulation to physiological processes. Several fatty acids have been both positively and negatively implicated in different biological processes in skeletal muscle and other tissues. To gain insight into biological processes associated with fatty acid content in skeletal muscle, the aim of the present study was to identify differentially expressed genes (DEGs) and functional pathways related to gene expression regulation associated with FA content in cattle.
Skeletal muscle transcriptome analysis of 164 Nellore steers revealed no differentially expressed genes (DEGs, FDR 10%) for samples with extreme values for linoleic acid (LA) or stearic acid (SA), and only a few DEGs for eicosapentaenoic acid (EPA, 5 DEGs), docosahexaenoic acid (DHA, 4 DEGs) and palmitic acid (PA, 123 DEGs), while large numbers of DEGs were associated with oleic acid (OA, 1134 DEGs) and conjugated linoleic acid cis9 trans11 (CLA-c9t11, 872 DEGs). Functional annotation and functional enrichment from OA DEGs identified important genes, canonical pathways and upstream regulators such as SCD, PLIN5, UCP3, CPT1, CPT1B, oxidative phosphorylation mitochondrial dysfunction, PPARGC1A, and FOXO1. Two important genes associated with lipid metabolism, gene expression and cancer were identified as DEGs between animals with high and low CLA-c9t11, specifically, epidermal growth factor receptor (EGFR) and RNPS.
Only two out of seven classes of molecules of FA studied were associated with large changes in the expression profile of skeletal muscle. OA and CLA-c9t11 content had significant effects on the expression level of genes related to important biological processes associated with oxidative phosphorylation, and cell growth, survival, and migration. These results contribute to our understanding of how some FAs modulate metabolism and may have protective health function.
Lipids are a class of molecules present in all cell types. They contribute to cellular structure, energy storage and several biological functions from transcriptional regulation to physiological processes . Over the last couple decades, it has been reported that consumption of fatty acids is associated with fat deposition and may have metabolic effects, such as altered blood lipid and lipoprotein content . However, there are still contrasting opinions about the role of dietary fatty acids on human health, which can be observed in previous reviews and meta-analysis studies conducted in the last years [3, 4].
Lipid composition in beef has become an important discussion point. Although meat is a significant source of fat in the human diet, it contains high saturated fat content that is associated with some diseases such as obesity, cancer and coronary heart disease [5, 6]. On the other hand, meat has high nutritional value and is an important source of unsaturated fatty acids such as oleic acid (OA) and conjugated linoleic acid (CLA), which have beneficial effects on human health . Thus, there is a need to improve the nutritional value of meat and better understand the biological and molecular processes associated with fatty acid composition in skeletal muscle.
Previous studies have shown that OA as well as other fatty acids such as linoleic acid (LA), docosahexaenoic acid (DHA), eicosapentaenoic acid (EPA) and CLA can regulate gene transcription in tissues such as muscle, liver, adipose, monocytes and blood mononuclear cells [8, 9]. Activation of nuclear receptors such as peroxisome proliferator-activated receptors (PPARs) and coactivators of the sterol response element binding protein (SREBP) can mediate this regulation. Both PPARs and SREBPs play a role in preadipocyte differentiation, lipid homeostasis and peroxisomal beta-oxidation regulation by controlling the transcription of acyl-CoA oxidase and other enzymes involved in adipocyte development [10–12]. While elevated intracellular levels of stearic (SA) and palmitic (PA) acids were associated with apoptotic death  and inflammatory process . However, the effective involvement of these fatty acids in molecular and physiological processes in skeletal muscle remains unclear.
The present investigation was undertaken to identify differentially expressed genes (DEGs) and functional pathways associated with seven FA content in beef cattle. Our hypothesis was that variation in FA content could be associated with difference in gene expression in skeletal muscle. Our study revealed that up regulation of stearoyl-CoA desaturase (SCD) is associated with deposition of unsaturated FA and that increased levels of OA and CLA are associated with expression of genes that have protective function associated with important human diseases.
Phenotypes and sequencing data
Seven different FAs, important to beef quality and human health, were used in this transcriptome study, including: oleic acid (OA, C18:1 cis9), palmitic acid (PA, 16:0), stearic acid (SA, C18:0), linoleic acid (LA, C18:2 cis9,12), conjugated linoleic acid cis9 trans11 (CLA-c9t11, C18:2 cis9 trans11), eicosapentaenoic acid (EPA, C22:5) and docosahexaenoic acid (DHA, C22:6). The skeletal muscle samples analyzed were chosen based on their extreme values for each FA studied (Additional file 1). When ranked by fatty acids content in skeletal muscle, the statistical test of means performed between the high (H) and low (L) groups according to FA had significant differences (p-value < 0.05) for OA, PA, SA, CLA-c9t11, LA, EPA and DHA, but no difference in IMF content or backfat thickness (Additional file 2). However, animals with extreme phenotypes for OA, PA, SA, CLA-c9t11, LA, EPA and DHA also had significant differences for other FAs. For example, animals with high levels of OA had significantly (p < 0.05) lower levels for SA and PA, but high levels for CLA-c9t11. While extreme animals for CLA-c9t11 had significant low levels for SA and PA, but high levels for OA and DHA. The correlation between any two FAs tended to be low (<0.2). However, there were some FAs that high negative correlations, e.g. OA and PA (-0.63) and OA and SA (-0.62) or high positive correlations, e.g. EPA and LA (0.66) and EPA and DHA (0.74) (Table 1).
The average number of mapped read pairs per sample before and after filtering were 14.47 million and 11.5 million respectively (see Additional file 3), which means that on average 75.85% of total read pairs were mapped to the Bos taurus UMD3.1 reference genome assembly. After filtering for transcript expression level, a total of 17,137; 17,070; 17,053; 17,052; 17,038; 16,989 and 17,100 annotated genes were used for differential expression analysis for OA, PA, SA, CLA-c9t11, EPA, DHA and LA content, respectively.
Differential expression analysis
Differential gene expression analysis was performed for each FA by camparing gene expression level between the H and L FA content groups. No DEGs were identified for LA and SA (see Additional files 4 and 5), while 5, 4 and 123 DEGs were identified for EPA, DHA and PA, respectively (see Additional files 6, 7 and 8). In contrast, 1134 and 872 DEGs for OA and CLA-c9t11, respectively, were identified between H and L groups (for more details see Additional file 9 and 10). A Volcano plot of log2 fold-change (x-axis) versus − log10 FDR-corrected p-value (y-axis) for OA is shown in Fig. 1. The volcano plot of other FAs are shown in Additional file 11. Of the 1134 DEGs for OA, 870 genes were up-regulated and 264 were down-regulated and of the 872 DEGs for CLA-c9t11, 536 genes were up-regulated and 336 were down-regulated in the L group relative to the H group. The expression level, fold change, adjusted p-value and annotation of all DEGs for OA and CLA-c9t11 are shown in more detail in Additional files 9 and 10.
Functional enrichment analysis for oleic acid content
In this study, two different tools of pathway analysis were performed: DAVID and Ingenuity Pathway Analysis® (IPA®). According to DAVID enrichment analysis, several ontology terms were significantly enriched (p-value adjusted for multiple testing by Benjamin-Hochberg , BH-adj-pvalue) for the differentially expressed gene lists (see Additional file 12). Cellular component GO terms related to mitochondrial envelope (BH-adj-pvalue = 2.6e-15), mitochondrial membrane (BH-adj-pvalue = 1.0e-14), respiratory chain (BH-adj-pvalue = 6.6e-12), mitochondrial matrix (BH-adj-pvalue = 6.0e-08) and intracellular organelle lumen (BH-adj-pvalue = 5.8e-06) were associated with DEGs for OA content. In addition, six KEGG pathways (BH-adj-p-value < 0.10) that were associated with oxidative phosphorylation (BH-adj-pvalue = 5.8e-10), ribosome (BH-adj-pvalue = 1.6e-02) and proteasome (BH-adj-pvalue = 2.9e-02) were also identified. These pathways are involved with oxidative phosphorylation and oxidative stress diseases such as Huntington (BH-adj-pvalue = 3.0E-10), Parkinson (BH-adj-pvalue = 1.5e-08) and Alzheimer (BH-adj-pvalue = 1.2E-07) (Table 2).
The IPA® software (IPA®, QIAGEN Redwood City, www.qiagen.com/ingenuity) was also performed for functional enrichment using the Global Functional Analysis (GFA) and Global Canonical Pathways (GCP) analysis from the DEGs list for OA content (1134 genes). Global Canonical Pathways analysis estimates the likelihood that the association between a set of DEGs in a pathway is due to random chance using a right-tailed Fishers Exact Test. Significance was set at BH-adj-pvalue ≤ 0.05, which indicated a non-random association. Several canonical pathways (BH-adj-pvalue ≤ 0.05, which indicated a non-random association) were identified by GCP analysis, including oxidative phosphorylation (BH-adj-pvalue = 7.94e-17), mitochondrial dysfunction (BH-adj-pvalue = 3.98e-16), insulin receptor signaling (BH-adj-pvalue = 0.02), LPS-stimulated MAPK signaling (BH-adj-pvalue = 0.035), role of NFAT in cardiac hypertrophy (BH-adj-pvalue = 0.035) and Huntington’s disease signaling (BH-adj-pvalue = 0.04) (see Additional file 13).
Twenty-five networks were enriched as a function of OA content as identified by IPA®GFA analyses (see Additional file 14). One in particular, shown in Fig. 2, was associated with lipid metabolism and energy production. This network contained molecules and genes such as cholesterol (LDL-cholesterol), very-low-density lipoprotein (VLDL-cholesterol), perilipin 5 (PLIN5), SCD, low-density lipoprotein carnitine palmitoyltransferase I (CPT1) and carnitine palmitoyltransferase I (CPT1B), respectively (Fig. 2). Among these only SCD was up-regulated in the H group as shown in red in Fig. 2.
To further explore the observed changes in gene expression, IPA® upstream regulator analyses were performed to identify the cascade of upstream transcriptional regulators that could be involved in the gene expression changes in skeletal muscle due to OA content. The IPA® program estimated the effects between transcriptional regulators and their target genes based on prior knowledge stored in the Ingenuity® Knowledge Base. Upstream regulators were identified from the list of DEGs for OA content (see Additional file 15), which were connected via mechanistic networks detected by IPA®. Peroxisome proliferator-activated receptor gamma, coactivator 1 alpha (PPARGC1A) was among the top predicted upstream regulators. PPARGC1A and ten other transcription factors (TFs) such as forkhead box O3 (FOXO3), peroxisome proliferator-activated receptor gamma, coactivator 1 beta (PPARGC1B), estrogen-related receptor gamma (ESRRG), sterol regulatory element binding transcription factor 1 (SREBF1), tumor protein p53 (TP53), forkhead box O1 (FOXO1), V-Myc avian myelocytomatosis viral oncogene homolog (MYC) and myogenic differentiation 1 (MYOD1) were predicted to be inhibited (blue shapes in Fig. 3), whereas the transcription factors forkhead Box O4 (FOXO4) and sterol regulatory element binding transcription factor 2 (SREBF2) were predicted to be activated (orange shapes in Fig. 3).
The functional enrichment analysis using the IPA® program identified several statically significantly enriched biological functions (BH-adj-p-value ≤ 0.05) from the OA DEGs list (see Additional file 16), e.g. processing of RNA (BH-adj-pvalue = 2.28e-07, mitochondrial disorder (BH-adj-pvalue = 3.46e-07), transcription of RNA (BH-adj-pvalue = 1.70e-06), expression of RNA (BH-adj-pvalue = 3.30e-06), transcription of DNA (BH-adj-pvalue = 2.51e-05), expression of mRNA (BH-adj-pvalue = 4.96e-05), catabolism of fatty acid (BH-adj-pvalue = 1.86e-03), beta-oxidation of palmitic acid (BH-adj-pvalue = 3.20e-03), depletion of triacylglycerol (BH-adj-pvalue = 9.34e-04 and oxidation of long chain fatty acid (BH-adj-pvalue = 2.86e-02). The indirect inhibitory interaction between SCD and β-oxidation and the depletion of triacylglycerol, which was associated with lipid metabolism and energy production are shown in Fig. 4. In the network shown in Fig. 2 it is also possible to observe the indirect interaction between SCD, LDL-cholesterol, VLDL-cholesterol, CPT1, CPT1B, PLIN5 and UCP3.
Functional enrichment analysis for conjugated linoleic acid cis9 trans11 content
The enrichment analysis performed by DAVID software using the list of DEGs for CLA-c9t11 content identified one KEGG pathway, which was related to ribosome function (BH-adj-pvalue = 1.2e-02). Among the molecular functions identified by DAVID associated (BH-adj-pvalue < 0.05) with the CLA-c9t11 DE genes identified between the groups with extreme values of CLA-c9t11 were nucleotide binding (BH-adj-pvalue = 1.7e-03), ATP binding (BH-adj-pvalue = 2.2e-02) and structural constituent of ribosome (BH-adj-pvalue = 2.4e-02) (see Additional file 17).
The top five canonical pathways identified for CLA-c9t11 content by GFA and GCP analysis using IPA® program were protein ubiquitination (nominal pvalue = 4.11e-07), interferon signaling (nominal pvalue = 2.9e-05), hypoxia signaling in the cardiovascular system (nominal pvalue = 8.8e-05), angiopoietin signaling (nominal pvalue = 1.01e-04) and mTOR signaling (nominal pvalue = 1.19e-04) (see Additional file 18).
The top five upstream regulators identified from the list of DEGs for CLA-c9t11, associated network functions and toxicology list (list of molecules/genes that are relevant to causality of the phenotype of interest, which can help to identify potential therapeutic or toxicity targets) for CLA-c9t11 content are shown in Additional file 18. Among the top upstream regulators were: interferon, lambda 1 (IFLN1), tumor necrosis factor (TNF), retinoic acid, signal transducer and activator of transcription 2 (STAT2), and interferon regulatory factor 7 (IRF7), which were predicted as activated in this study (see Additional file 19). However other important cytokines and transcription regulators were identified as activated/inhibited such as IFNK (interferon kappa), CCL5 (chemokine (C-C motif) ligand 5) and, PRL (prolactin), TRIM24 (tripartite motif containing 24), NFKB1A (Nuclear factor NF-kappa-B 1A), SREBF1, FOXO1, TP53 and MYC, respectively (see complete list in Additional file 19). Figure 5 shows the FOXO1 network, which has been associated with the metabolic disorder dyslipidemia.
The top five associated network functions detected from the list of CLA-c9t11 DEGs between the groups with extreme values of CLA-c9t11 were related to RNA damage, repair and post-transcriptional modification; developmental disorder; metabolic disorder; lipid metabolism, molecular transport, RNA trafficking and; cell death and survival (see Additional file 18). Figure 6 illustrates the network associated with lipid and nucleic acid metabolism and small molecule biochemistry, which presents epidermal growth factor receptor (EGFR) as a central gene and which was down-regulated in animals with high CLA-c9t11 content (colored in green). Figure 7 shows a network associated with molecular transport, RNA trafficking, RNA post-transcriptional modification that contained ribonucleoproteins (RNPS), which were up-regulated in animals with high CLA-c9t11 content. The top five Toxicology terms were associated with cardiac necrosis/cell death (nominal p-value = 1.18e-06), cardiac hypertrophy (nominal p-value = 2.4e-04), hypoxia-inducible factor signaling (nominal p-value = 3.2e-03), increases liver steatosis (nominal p-value = 3.3e-03) and, PPARA/RXRA activation (nominal p-value = 5.3e-03) (see Additional file 18).
Validation of differentially expressed genes
The quantitative real-time polymerase chain reaction (qRT-PCR) results and the relative quantification of DEGs selected to validate the present RNA-Seq study are shown in Additional file 20. Nine genes differently expressed between high and low OA groups were tested by qPCR and seven genes (78%) confirmed the RNAseq data (p < 0.1), one gene confirmed but not at a significant level and one gene disagreed with RNAseq data (Additional file 21).
Skeletal muscle is important for energy homeostasis, as it is a site of energy storage and insulin-stimulated glucose. Metabolic diseases such as obesity and coronary disorders could be a consequence of insulin resistance, i.e. the inability of insulin to drive glucose into the skeletal muscle and other tissues as liver, which can be caused by excessive body fat deposition . Fatty acids in skeletal muscle play an important role in the structural and functional processes that influence the fluidity and stability of membrane structures, which impacts membrane functions such as transport, cell signaling, responses to oxidative damage and apoptosis. Studies have reported that the fatty acid profile of skeletal muscle phospholipids and triacylglycerides can be influenced by the fatty acid composition of human and animal diets [17, 18]. Several saturated and unsaturated fatty acids have been implicated, both positively and negatively, in different biological processes in skeletal muscle and other tissues such as liver and cardiac tissue [1, 19, 20]. In contrast, meta-analysis studies reported no evidence that a human diet, which was rich in monounsaturated fatty acids (MUFA) and polyunsaturated fatty acids (PUFA) and depleted for saturated fatty acid (SFA), had beneficial effects on coronary heart disease [3, 4].
In this study, all animals were of the same breed, sex, similar age and fed the same diet. However, different levels of FA content were observed in skeletal muscle. This maybe due to genetic variation among the animals  given that the environmental variables were similar across all animals. Transcriptomic analyses (validated by q RT-PCR [see Additional file 20]) indicated that many genes were differentially expressed (FDR 10%) between animals with high and low OA and CLA-c9t11 content. While, only a few genes were differentially expressed among the extreme groups for PA, SA, EPA and DHA and none for LA and SA. Additionally, previous studies reported only a few genes were affected by different levels of intramuscular fat content  and tenderness  in this same population.
Because the samples used for OA and CLA-c9t11 analysis also differed significantly for other FAs (see Additional file 2), we cannot be certain that the observed DE genes for OA and CLA-c9t11 are associated solely with these FAs. Despite the use of some samples in the analysis of several FA extremes, the lists of differentially expressed genes, and identified pathways were very different. This may not be surprising given the generally low correlation between FAs (Table 1). Thus it is very likely that utilization of only a couple samples in several analyses does not have an negative impact on other analyses. Finally, these findings indicate the biological importance of OA and CLA-c9t11 for skeletal muscle physiology. However, a word of caution is necessary because, a direct relation between gene expression level and protein abundance was not establish . These fatty acids may play an important role in transcription regulation and consequently in energy homeostasis and metabolic diseases such as obesity and diabetes as also presented by Ntambi .
The most abundant MUFA found in skeletal muscle fat is oleic acid (OA, C18:1 cis9) . It is present in membrane phospholipids, triglycerides and cholesterol esters. Human OA consumption has been associated with low levels of low-density lipoprotein (LDL) and also with the potential increase of high-density lipoprotein (HDL) levels in the blood. Furthermore, OA is converted to oleoylethanolamide (OEA), a small molecule in the intestine, which is involved in appetite control . Previous studies have reported an important contribution of OA intake to general human health, which could lead to a decrease in cholesterol levels, atherosclerosis risk and diabetes occurrence. Further, OA has protective effects against viral infection and cancer development [28–30]. Beef has a high unsaturated fatty acid profile, primarily OA . It is considered an important source of OA in human nutrition, and health [31–33]. Another important aspect of the OA content in beef is its relationship to meat quality traits such as tenderness, flavor and shelf-life . Thus, it is appealing to improve our knowledge about biological processes associated with OA content in skeletal muscle, which could both help to improve the nutritional value and quality of beef.
In mammals, the biosynthesis of the OA is dependent on the activity of stearoyl-CoA desaturase (SCD). This soluble enzyme catalyzes the introduction of a cis-double bond at the 9th carbon position of stearic acid (18:0), to form the MUFA, oleic acid (18:1) . SCD is also required for the synthesis of highly unsaturated essential fatty acids or PUFAs such as eicosopentaenoic (EPA, C22:5 n-3) and docosahexaenoic acids (DHA, C22:6 n-6); synthesized from a-linolenic acid (LA, C18:2 cis9,cis12), and arachidonic acid (synthesized from linoleic acid), which are dependent on the presence of oleic acid . In ruminants, such as beef cattle, the concentration of OA in adipose tissue depends more on the concentration of stearic acid hydrogenated by ruminal microorganisms than from stearic acid desaturation by SCD in the tissues .
EPA and DHA are omega-3 (n-3) long-chain PUFA, which are dietary fats associated with health benefits. They are present in cell membranes and play a role in anti-inflammatory processes and in the viscosity of cell membranes . EPA and DHA are also precursors of several metabolites that are potent lipid mediators, considered by many investigators to be beneficial in the prevention or treatment of several diseases such as cardiovascular, obesity and Alzheimer [35–37]. However, in this study both EPA and DHA did not appear to have much of an effect on skeletal muscle gene expression, which may mean that they do not have a biological effect on skeletal muscle. It is interesting to speculate that DEG analyses of other tissues such as liver and adipose may be able to identify tissues responsive to these fatty acids. In contrast, Conjugated linoleic acid cis9, trans11 has a large effect on gene expression. This FA is associated with human health, specifically related with diseases such as cancer, immune and inflammatory responses and obesity [38, 39] and is naturally present in food sources such as fats in milk and meat of ruminant animals . This may indicate that skeletal muscle is a primary target of Conjugated linoleic acid cis9, trans11.
The enriched transcript clusters identified by DAVID based on biological process GO terms (see Additional files 16 and 18) from the list of DEGs were associated with electron transport chain, mRNA translation and generation of precursor metabolites and energy. FAs like OA and CLA-c9t11 are involved in both extracellular and intracellular activities, which are in a partially dissociated state at physiological pH with part of the molecules diffused into the cytosol and other part integrating the membrane structure. These FAs impact several molecular functions in the cell such as: (1) ion transport across membranes; (2) activate intracellular enzymes; (3) modulation of Na+/K + -ATPase activation, (4) modulation of hormone/receptor interactions; (5) gene transcription and; (6) mRNA translation [40–44].
OA can be obtained by either dietary consumption or de novo lipogenesis. De novo lipogenesis of OA can occur from both dietary carbohydrates (will be transformed to triacylglycerol in case of positive energy balance) and fat. The enzyme responsible for MUFA OA biosynthesis is SCD , which is considered as a key enzyme in de novo lipogenesis, a process that is associated with obesity, synthesis of cholesterol and apoptosis. In the present study, SCD was up-regulated in the group of animals with high OA content. This was expected as this gene function is linked to OA synthesis. Hulver et al.  suggested that elevated expression of SCD1 in human skeletal muscle contributes to abnormal lipid metabolism and progression of obesity. However, in this study, the IMF content between the L and H group was not statistically different (p-value = 0.06), as the H group had an average IMF content of 3.1% and the L group had an IMF of 2.7%. In this study, the fatty acid synthase (FAS) and acetyl-CoA carboxylase, responsible for fatty acid metabolism, were not differentially expressed between the groups H an L for any FA studied, considering that the diet was the same for all animals. This may indicate: 1) the increased SCD1 in this study was associated just with higher synthesis of OA and not with higher fat deposition; 2) a difference in how bovine and human skeletal muscles respond to OA content.
Schrauwen , Clapham et al.  reported that mice which overexpressed UCP3 in skeletal muscle had lower body weights and increased AMPK activity, elevated metabolic rate, lower fasting plasma glucose and insulin levels compared to wild type mice. In the first instance UCP3, which was down-regulated in the group of animals with high OA content was associated with regulation of energy homeostasis, but according to Schrauwen et al.  UCP3 played an important role in the protection of mitochondria against lipid-induced oxidative damage by ROS production. This protection was the result of better fatty acid exportation by mitochondria and less production of hydrogen peroxide (H2O2), which caused oxidative damage to proteins and phospholipids. PLIN5, which was down-regulated in the group of animals with high OA content, in turn is a protein associated with lipid droplet that maintains the balance between lipogenesis and lipolysis and regulates fatty acid oxidation in oxidative tissues. These proteins recruit mitochondria to the surface of lipid droplets, which regulates their homeostasis by storing fatty acids in the form of triglycerides and releasing fatty acids for fatty acid oxidation by mitochondria . CPT1, which was up-regulated in the group of animals with high OA content and CPT1B, which was down-regulated in the group of animals with high OA content, are important for the beta-oxidation of long chain fatty acids, which are involved in the long chain fatty acids transport to the mitochondrial matrix and associated with type 2 diabetes and insulin resistance. The inhibition of CPT1 decreases the ability of muscles to oxidize fatty acids due to decreased long chain fatty acids transport into muscle mitochondria that could in turn be responsible for fat accumulation in skeletal muscle . In this study, animals with lower values of OA content presented higher expression levels (Figs. 1 and 4) of UCP3, PLIN5, CPT1 and CPT1B and lower expression level of SCD, corroborating with these previous studies [48–51] where the expression level of these genes were associated with fatty acid metabolism, energy homeostasis and animal health.
Herein, networks associated with reactive oxygen species (ROS) generation such as oxidative phosphorylation and mitochondrial dysfunction were detected, results that corroborate previous studies that reported the role of OA  and fatty acids  in the control of lipid oxidation. Mitochondria, where β-oxidation occurs, are the major source of ROS generation. Fatty acid molecules are broken down to generate acetyl-CoA, which enters the citric acid cycle, and NADH and FADH2 co-enzymes are used in the electron transport chain . Tizioto et al.  performed DEG gene analysis between animals that showed high (inefficient) and low (efficient) residual feed intake from the same population used in this study, and indicated the metabolic processes that underlie oxidative stress are a primarily network related to lipid metabolism and energy required for maintenance. These authors also observed the important role of fatty acids in energy metabolism and mitochondrial action on lipid oxidation process similar to that reported in this study.
One of the major transcription factors identified in our study was PPARGC1A, which was differentially expressed between the H and L groups for CLA-c9t11 and down-regulated in H group. PPARGC1A is a transcriptional coactivator that regulates genes involved in energy metabolism and provides a link between external physiological stimuli and the regulation of mitochondrial biogenesis. PPARGC1A regulates global oxidative metabolism by controlling both the induction of mitochondrial metabolism and the removal of its ROS by-products, which could elevate oxidative metabolism and minimize the impact of ROS on cell physiology in muscle . A review by Puigserver and Spiegelman  reported that the PPARGC1A transcriptional coactivator could be activated or not in response to environmental stimuli or by cellular signals such as cAMP and cytokine pathways, which control energy and nutrient homeostasis. Activated PPARGC1A stimulates mitochondrial oxidative metabolism, fiber-type switching in skeletal muscle, and the fasted response in liver . These effects of PPARGC1A are the result of gene expression regulation by interacting specifically with other transcription factors such as nuclear hormone receptors, nuclear respiratory factors, and muscle-specific transcription factors.
Another important TF associated with lipid metabolism, FOXO1, which regulates the transcriptional cascades of glucose metabolism was differentially expressed in the current study for CLA-c9t11. FOXO1, which was down-regulated in the group of animals with high CLA-c9t11 content, is highly expressed in tissues such as pancreas, liver, skeletal muscle and adipose tissue in response to insulin and is associated with metabolic disorders as dyslipidemia (Fig. 5). Dyslipidemia is a metabolic disorder that contributes to the development of atherosclerosis due to the elevation of plasma cholesterol, triglycerides (TGs) or low high-density lipoprotein level . These results are in agreement with a previous study, which indicated that the level of CLA in human diet can influence atherosclerosis in animal models .
The functional annotation analysis by IPA® for transcripts associated with FAs with many DEG (i.e. OA and CLA-c9t11) included important TFs, such as PPARGC1A, SREBF1, TP53, FOXO1 and MYC as upstream regulators. In addition, biological processes, which can cause damage to human and animals health, were identified as enriched. Several disorders, diseases, molecular and cellular functions were also identified as associated with OA or CLA-c9t11 FA content in skeletal muscle by transcriptome profiling (Additional files 13 and 18).
Two important genes associated with lipid metabolism, gene expression and cancer were identified as DEGs between animals with high and low CLA-c9t11, specifically, epidermal growth factor receptor (EGFR, Fig. 5) and RNPS (Fig. 6). EGFR is a tyrosine protein kinase, which has been implicated in several biological processes such as cell growth, survival, and migration. The overexpression of EGFR was associated with breast, lung and colon cancers [59, 60]. In this study, the EGFR was down-regulated in animals with high CLA-c9t11 content, result that agree with the association between CLA-c9t11 content and cancer prevention. RNPS, which was up-regulated in the group of animals with high CLA-c9t11 content, can affect several basic functions such as protein synthesis, gene expression, and chromosome stability. It is responsible for pseudouridylation of ribosomal and spliceosomal small nuclear RNAs, ribosomal RNAs processing, and RNA telomerase stabilization . MacDonald  and Park et al.  reported that association of CLA-c9t11 transcriptome level can regulate cell proliferation and differentiation.
The effective involvement of OA and CLA-c9t11 in biological processes like transcriptional regulation is still unclear. However, this transcriptome profiling study, which used Longissimus dorsi muscle from cattle with extreme values of OA and CLA-c9t11 content, demonstrated that both OA and CLA-c9t11 content had large effects on the expression level of genes related to important biological processes associated with oxidative phosphorylation, cell growth, survival, and migration. It is tempting to speculate that these fatty acids have direct and maybe indirect physiological effects on skeletal muscle, while other fatty acids do not. Additional nutrigenomic and metabolomic studies are necessary to elucidate the metabolic mechanisms by which these fatty acids influence cell growth, survival, and migration, which are interestingly also associated with several human diseases directly related to human dietary intake of fatty acids.
Animals, samples and phenotypes
Three hundred and ninety Nellore steers from the Brazilian Agricultural Research Corporation (EMBRAPA/Brazil) experimental breeding herd between 2009 and 2011 were included in this study; see Cesar et al.  for additional details. These steers were sired by 34 unrelated sires that represented most of the main genealogies used in Brazil according to the National Summary of Nellore produced by the Brazilian Association of Zebu Breeders (ABCZ) and National Research Center for Beef Cattle. Animals were raised in feedlots under identical nutrition and handling conditions until slaughter at an average age of 25 months. Samples from LD muscle located between the 12th and 13th ribs were collected in two time-points: at slaughter for RNA sequencing analysis to guarantee the RNA integrity, and 24 h after slaughter for the fatty acid content measurement according to research guidelines for cookery, sensory evaluation and instrumental tenderness measurements of fresh meat .
Approximately 100 g samples of LD muscle were lyophilized and ground to a fine powder. Five grams of this ground, lyophilized tissue was used to measure the OA, PA, SA, CLA-c9t11, EPA, DHA and LA content, which were chosen based on their importance in many biological processes and human health [35–40]. The fatty acid content measurement was conducted as described by Hara and Radin , except the hexane to propanol ratio was increased to 3:2. Approximately 4 g of LD muscle were lyophilized, ground in liquid nitrogen, mixed with 28 mL of hexane/propanol (3:2 vol/vol) and homogenized for 1 min. Samples were vacuum filtered and 12 mL sodium sulfate (67 mg mL− 1) solution was added and agitated for 30 s. The supernatant was transferred to a tube with 2 g of sodium sulfate and insufflated with nitrogen gas (N2), after which the tube was sealed and incubated at room temperature for 30 min. Subsequently, the liquid was transferred to 10 mL test tube, insufflated with N2, sealed and kept at -20 °C until dry with N2 for methylation. The extracted lipids were hydrolyzed and methylated as described by Christie , except that hexane and methyl acetate were used instead of hexane:diethyl ether:formic acid (90:10:1). Around 40 mg of lipids were transferred to a tube containing 2 mL of hexane. Subsequently, 40 μL of methyl acetate were added, the sample agitated, and 40 μL of methylation solution (1.75 mL of methanol/0.4 mL of 5.4 mol/L of sodium metoxyde) were added. This mixture was agitated for 2 min and incubated for 10 min at room temperature. Then 60 μL of finishing solution (1 g of oxalic acid/30 mL of diethyl ether) was added and the mixture was agitated for 30 s, after which 200 mg of calcium chloride was added. The sample was then mixed and incubated at room temperature for 1 h. Samples were centrifuged at 3200 rpm, for 5 min at 5 °C. The supernatant was collected for determination of fatty acid. Fatty acid methyl esters were quantified with a gas chromatograph (ThermoFinnigan, Termo Electron Corp., MA, USA) equipped with a flame ionization detector and a 100 m Supelco SP-2560 (Supelco Inc., PA, USA) fused silica capillary column (100 m, 0.25 mm and 0.2 μ m film thickness). The column oven temperature was held at 70 °C for 4 min, then increased to 170 °C at a rate of 13 °C min − 1, and subsequently increased to 250 °C at a rate of 35 °C min − 1, and held at 250 °C for 5 min. The gas fluxes were 1.8 mL min − 1 for carrier gas (He), 45 mL min − 1 for make-up gas (N2), 40 mL min − 1 for hydrogen, and 450 mL min − 1 for synthetic flame gas. One μL sample was analyzed. Injector and detector temperatures were 250 and 300 °C, respectively. The fatty acids were identified by comparison of retention time of methyl esters of the samples with standards of fatty acids butter reference BCR-CRM 164, Anhydrous Milk Fat-Producer (BCR Institute for Materials and Reference Measurements) and also with commercial standard for 37 fatty acids Supelco TM Component FAME Mix (cat 18919, Supelco, Bellefonte, PA). The fatty acid content was quantified by normalizing the area under the curve of methyl esters using Chromquest 4.1 software (Thermo Electron, Italy). The fatty acid content was expressed as a weight percentage (mg/mg). These analyses were performed at the Animal Nutrition and Growth Laboratory at ESALQ, Piracicaba, São Paulo, Brazil.
From a total of 390 animals, we selected 30 animals with either the highest (H) or lowest (L) FA content for each one of the FA chosen (OA, PA, SA, CLA-c9t11, EPA, DHA and LA) to form the groups that were tested for differential expression analysis. Data from some animals were used in more than one analysis, which resulted in a total of 164 animals that were used on this study. The extreme values of FAs were identified based in their adjusted phenotype values, which was estimated via regression as recommended by Naylor et al.  and selected for differential gene expression analysis. SAS PROC MIXED (SAS version 9.1) was used to adjust the phenotypic values and run the regression model, in which the contemporary group (animals with the same origin, birth year and slaughter date) and hot carcass weight were fitted as fixed and covariate, respectively. To verify the difference of FA content between the high and low groups, IMF (%) and BFT (mm) a Student’s t-test was performed using R package. Because of some samples were used in more than one FA extreme group, the statistical test of means was performed between the high (H) and low (L) groups according to FA by Student’s t-test performed in R package. Pairwise correlation analysis was also performed for all FAs. The IMF and BFT measurement were obtained as described by Cesar et al.  and Tizioto et al. .
RNA extraction, quality analysis, library preparation and sequencing
Total RNA was extracted from 100 mg of frozen LD muscle that was collected at slaughter using the TRIzol reagent (Life Technologies, Carlsbad, CA). RNA integrity was verified by Bioanalyzer 2100 (Agilent, Santa Clara, CA, USA). Only samples with RIN score > 8 were used. A total of 2 μg of total RNA from each sample was used for library preparation according to the protocol described in the TruSeq RNA Sample Preparation kit v2 guide (Illumina, San Diego, CA). The average insert size of the libraries was estimated using the Agilent Bioanalyzer 2100 (Agilent, Santa Clara, CA, USA) and quantified using quantitative PCR with the KAPA Library Quantification kit (KAPA Biosystems, Foster City, CA, USA). Quantified, samples were diluted and pooled (three pools of six samples each). Three lanes of a sequencing flowcell, using the TruSeq PE Cluster kit v3-cBot-HS kit (Illumina, San Diego, CA, USA), were clustered and sequenced using HiSeq2500 ultra-high-throughput sequencing system (Illumina, San Diego, CA, USA) with the TruSeq SBS Kit v3-HS (200 cycles), according to manufacturer instructions. The sequencing analyses were performed at the Genomics Center at ESALQ, Piracicaba, São Paulo, Brazil.
Quality control and read alignment
Sequencing adaptors and low-complexity reads were removed in an initial data filtering step. Quality control and reads statistics were estimated with FASTQC version 0.10.1 software [http://www.bioinformatics.bbsrc.ac.uk/projects/fastqc/]. Tophat v. 2.0.11 software  with Bowtie2 (version 2.2.3)  was used to map reads to the UMD3.1 Bos taurus reference assembly available at Ensembl [http://www.ensembl.org/Bos_taurus/Info/Index/]. The abundance (read counts) of mRNAs for all annotated genes was calculated using HTSeq version 0.6.1 software [http://www-huber.embl.de/users/anders/HTSeq/] . Only sequence reads that uniquely mapped to known chromosomes (excluding reads mapped to unassigned contigs) were used in further analyses.
Identification of differential expressed genes and pathway analysis
Differentially expressed genes were identified using the DESeq2 statistical package available at Bioconductor open source software for bioinformatics , using multi-factor design. Prior to statistical analysis, read count data was filtered as follows: i) transcripts with zero counts were removed (unexpressed); ii) genes with less than 1 read per sample on average were removed (very lowly expressed); iii) genes that were not present in at least three samples were removed (rarely expressed). Transcript expression level was fit as a negative binomial distribution. Age (covariate) and contemporary group (animals with the same origin, birth year and slaughter date; class effect) were fit as factors in the multi-factor model. Exploratory diagnostic plots were generated to check the dispersion estimates. The cut-off approach performed to identify the DEGs was Benjamini-Hochberg  methodology was used to control false discovery rate (FDR) at 10%.
Enrichment analysis of gene ontology terms was completed with the Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.7 tool  using the list of genes that presented FDR < 10%. The functional analysis of statistically significant gene expression changes (FDR 10%) between the groups was performed with QIAGEN’s Ingenuity® Pathways Analysis (IPA®, QIAGEN Redwood City, www.qiagen.com/ingenuity), which uses the Fisher’s exact test to calculate a p-value for each biological function assigned to the experimental data. These analyses were performed just for OA and CLA-c9t11 because they presented higher impact on gene expression level.
Quantitative real-time PCR analysis
Total RNAs of three samples of each group (H and L) and nine DE genes were randomly selected for quantitative PCR analysis. The High Capacity cDNA Reverse Transcription Kit (Life Technologies, Grand Island, NY) was performed in reverse transcription (RT) as manufacturer’s procedure. The RT reaction contained 150 ng of total RNA, 1 μL of 10× RT buffer, 0.5 μl of 25× dNTP mixture, 1 μL of 10× random reverse primers, 1 μL of 10× gene-specific reverse primers (1 μM) and 0.5 μLof MultiScribe RT (50 U/μl). The reactions of 10 μL were incubated in thermocycler for 10 min at 25 °C, 2 h at 37 °C and then held at 4 °C. All samples were treated with Ambion TURBO DNA-free® DNase (Life Technologies, Grand Island, NY) to remove contaminating DNA from RNA preparations. Quantitative real-time PCR was performed on an Applied Biosystems QuantStudio™ 12 K Flex System (Thermo Fisher Scientific, Waltham, MA) using the TaqMan OpenArray gene expression plates. The TaqMan assays of nine DE genes selected (NDUFV1, NDUFS6, NDUFB7, COX7A1, CYP4B1, JAM2, PLIN5, COA5 and NDUFS8) for this study was designed and performed by manufacturer (http://www.thermofisher.com/order/genome-database). The PCR reaction (1.2 μLof sample plus 3.8 μL of PCR mix) was prepared for all samples and placed into 384-well sample plate. The PCR reactions were dispensed using the Accufill System onto the corresponding OpenArray plate. The PCR condition was performed as manufacturer’s procedure. The comparative Ct method for quantification was performed to quantify the relative expression of specific genes by QuantStudio 12 K Flex Software (Thermo Fisher Scientific, Waltham, MA) according to Pfaffl . Eukaryotic Translation Elongation Factor 1 Alpha 2 (EEF1A2 - Bt03229236_mH) and Beta-2-Microglobulin (B2M - Bt03251628_m1) were selected as reference genes to normalize gene expression. Analysis of covariance (ANCOVA)  using the same fixed effect and covariate used in RNA-Seq analysis was performed to test the significance (p < 0.10) of FA content (group effect) on gene expression level of the DE genes selected randomly for this study.
The dataset supporting the conclusions of this article is available in the in the European Nucleotide Archive (ENA) repository (EMBL-EBI), under accession PRJEB13188 [http://www.ebi.ac.uk/ena/data/view/PRJEB13188].
Brazilian Association of Zebu Breeders
Conjugated linoleic acid
Conjugated linoleic acid cis 9, trans 11
Database for annotation, visualization and integrated discovery
Differentially expressed gene
False discovery rate
Global canonical pathways
Global functional analysis
Ingenuity pathways analysis
Kyoto encyclopedia of genes and genomes
Monounsaturated fatty acidOA: Oleic acid
Polymerase chain reaction
Polyunsaturated fatty acid
Quantitative real-time polymerase chain reaction
Reactive oxygen species
Saturated fatty acid
Very low-density lipoprotein
Nakamura MT, Nara TY. Essential fatty acid synthesis and its regulation in mammals. Prostaglandins Leukot Essent Fatty Acids. 2003;68:145–50.
Wood JD, Enser M, Fisher AV, Nute GR, Sheard PR, Richardson RI, et al. Fat deposition, fatty acid composition and meat quality: a review. Meat Sci. 2008;78:343–58.
Schwingshackl L, Hoffmann G. Monounsaturated fatty acids and risk of cardiovascular disease: synopsis of the evidence available from systematic reviews and meta-analyses. Nutrients. 2012;4(12):1989–2007.
Schwingshackl L, Hoffmann G. Association of dietary, circulating, and supplement fatty acids with coronary risk. Ann Intern Med. 2014;161(6):455–6.
Lichtenstein AH, Kennedy E, Barrier P, Danford D, Ernst ND, Grundy SM, et al. Dietary fat consumption and health. Nutr Rev. 1998;56:S3–19.
Lin J, Zhang SM, Cook NR, Lee IM, Buring JE. Dietary fat and fatty acids and risk of colorectal cancer in women. Am J Epidemiol. 2004;160:1011–22.
Laaksonen DE, Nyyssönen K, Niskanen L, Rissanen TH, Salonen JT. Prediction of cardiovascular mortality in middle-aged Men by dietary and serum linoleic and polyunsaturated fatty acids. Arch Intern Med. 2005;165(2):193–9.
Sessler AM, Ntambi JM. Polyunsaturated fatty acid regulation of gene expression. J Nutr. 1998;128(6):923–6.
Bouwens M, van de Rest O, Dellschaft N, Bromhaar MG, de Groot LC, Geleijnse JM, et al. Fish-oil supplementation induces antiinflammatory gene expression profiles in human blood mononuclear cells. Am J Clin Nutr. 2009;90:415–24.
Smith SB, Gill CA, Lunt DK, Brooks MA. Regulation of Fat and fatty acid composition in beef cattle, Asian-Aust. J Anim Sci. 2009;22:1225–33.
MacDougald OA, Lane MD. Adipocyte differentiation. When precursors are also regulators. Curr Biol. 1995;5:618–21.
Poletto AC, Furuya DT, David-Silva A, Ebersbach-Silva P, Santos CL, Corrêa-Giannella ML, et al. Oleic and linoleic fatty acids downregulate Slc2a4/GLUT4 expression via NFKB and SREBP1 in skeletal muscle cells. Mol Cell Endocrinol. 2015;401:65–72.
Ulloth JE, Casiano CA, De Leon M. Palmitic and stearic fatty acids induce caspase-dependent and -independent cell death in nerve growth factor differentiated PC12 cells. J Neurochem. 2003;84(4):655–68.
Kennedy A, Martinez K, Chuang CC, LaPoint K, McIntosh M. Saturated fatty acid-mediated inflammation and insulin resistance in adipose tissue: mechanisms of action and implications. J Nutr. 2009;139(1):1–4.
Benjamini Y, Hochberg Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J Roy Statist Soc Ser B (Methodological). 1995:57(1):289–300.
Furuhashi M, Hotamisligil GS. Fatty acid-binding proteins: role in metabolic diseases and potential as drug targets. Nat Rev Drug Discov. 2008;7(6):489.
Siscovick DS, Raghunathan TE, King I, et al. Dietary intake and cell membrane levels of long-chain n − 3 polyunsaturated fatty acids and the risk of primary cardiac arrest. JAMA. 1995;274:1363–7.
Ayre KJ, Hulbert AJ. Dietary fatty acid profile influences the composition of skeletal muscle phospholipids in rats. J Nutr. 1996;126:653–62.
Furuhashi M, Hotamisligil GS. Fatty acid-binding proteins: role in metabolic diseases and potential as drug targets. Nat Rev Drug Discov. 2008;7(6):489–503.
Marín-García J, Goldenthal MJ. Fatty acid metabolism in cardiac failure: biochemical, genetic and cellular analysis. Cardiovasc Res. 2002;54(3):516–27.
Cesar ASM, Regitano LC, Mourão GB, Tullio RR, Lanna DP, Nassu RT, et al. Genome-wide association study for intramuscular fat deposition and composition in Nellore cattle. BMC Genet. 2014;15:39.
Cesar ASM, Regitano LC, Koltes JE, Fritz-Waters ER, Lanna DP, Gasparin G, et al. Putative regulatory factors associated with intramuscular fat content. PLoS ONE. 2015;10:e0128350.
Gonçalves TM. Differential expression of genes related with meat tenderness in Nellore cattle. 2015. http://www.teses.usp.br/teses/disponiveis/11/11139/tde-12052015-165345/pt-br.php. Accessed 22 Dec 2015.
Liu Y, Beyer A, Aebersold R. On the Dependency of Cellular Protein Levels on mRNA Abundance. Cell. 2016;165(3):535–50.
Ntambi JM. Regulation of stearoyl-CoA desaturase by polyunsaturated fatty acids and cholesterol. J Lipid Res. 1999;40:1549–58.
Dowhan W, Mileykovskaya E, Bogdanov M. Diversity and versatility of lipid-protein interactions revealed by molecular genetic approaches. Biochim Biophys Acta. 2004;1666:19–39.
Schmid A. The role of meat fat in the human diet. Crit Rev Food Sci Nutr. 2011;51:50–66.
Schwartz GJ, Fu J, Astarita G, Li X, Gaetani S, Campolongo P, Cuomo V, Piomelli D. The lipid messenger OEA links dietary fat intake to satiety. Cell Metab. 2008;8:281.
Tsimikas S, Reaven PD. The role of dietary fatty acids in lipoprotein oxidation and atherosclerosis. Curr Opin Lipidol. 1998;9:301–7.
Terés S, Barceló-Coblijn G, Benet M, Alvarez R, Bressani R, Halver JE, et al. Oleic acid content is responsible for the reduction in blood pressure induced by olive oil. Proc Natl Acad Sci U S A. 2008;105:13811–6.
Kritchevsky D. Antimutagenic and some other effects of conjugated linoleic acid. Br J Nutr. 2000;83(5):459–65.
Pauwels EK. The protective effect of the Mediterranean diet: focus on cancer and cardiovascular risk. Med Princ Pract. 2011;20:103–11.
Laaksonen DE, Nyyssönen K, Niskanen L, Rissanen TH, Salonen JT. Prediction of Cardiovascular Mortality in Middle-aged Men by Dietary and Serum Linoleic and Polyunsaturated Fatty Acids. Arch Intern Med. 2005;165:193–9.
Lopez S, Bermudez B, Montserrat-de la Paz S, Jaramillo S, Varela LM, et al. Membrane composition and dynamics: A target of bioactive virgin olive oil constituents. Biochim Biophys Acta. 2014;1838:1638–56.
Smith GI, Atherton P, Reeds DN, Mohammed BS, Rankin D, Rennie MJ, Mittendorfer B. Dietary omega-3 fatty acid supplementation increases the rate of muscle protein synthesis in older adults: a randomized controlled trial. Am J Clin Nutr. 2011;93:402–12.
Serhan CN, Chiang N, Van Dyke TE. Resolving inflammation: dual anti-inflammatory and pro-resolution lipid mediators. Nat Rev Immunol. 2008;8:349–61.
Conquer JA, Tierney MC, Zecevic J, Bettger WJ, Fisher RH. Fatty acid analysis of blood plasma of patients with Alzheimer’s disease, other types of dementia, and cognitive impairment. Lipids. 2000;35(12):1305–12.
Neff LM, Culiner J, Cunningham-Rundles S, Seidman C, Meehan D, Maturi J, et al. Algal docosahexaenoic acid affects plasma lipoprotein particle size distribution in overweight and obese adults. J Nutr. 2011;141:207–13.
Dilzer A, Park Y. Implication of conjugated linoleic acid (CLA) in human health. Crit Rev Food Sci Nutr. 2012;52(6):488–513.
Moloney F, Toomey S, Noone E, Nugent A, Allan B, Loscher CE, Roche HM. Antidiabetic effects of cis-9, trans-11-conjugated linoleic acid may be mediated via anti-inflammatory effects in white adipose tissue. Diabetes. 2007;56(3):574–82.
Schenkel LC, Bakovic M. Palmitic acid and oleic acid differentially regulate choline transporter-like 1 levels and glycerolipid metabolism in skeletal muscle cells. Lipids. 2014;49:731–44.
Swarts HG, Schuurmans Stekhoven FM, De Pont JJ. Binding of unsaturated fatty acids to Na+, K(+)-ATPase leading to inhibition and inactivation. Biochim Biophys Acta. 1990;1024:32–40.
Krey G, Braissant O, L’Horset F, Kalkhoven E, Perroud M, Parker MG, et al. Fatty acids, eicosanoids, and hypolipidemic agents identified as ligands of peroxisome proliferator-activated receptors by coactivator-dependent receptor ligand assay. Mol Endocrinol. 1997;11:779–91.
Calder PC. Functional Roles of Fatty Acids and Their Effects on Human Health. JPEN J Parenter Enteral Nutr. 2015;39:18S–32.
Scaglia N, Igal RA. Stearoyl-CoA desaturase is involved in the control of proliferation, anchorage-independent growth, and survival in human transformed cells. J Biol Chem. 2005;280:25339–49.
Hulver MW, Berggren JR, Carper MJ, Miyazaki M, Ntambi JM, Hoffman EP, et al. Elevated stearoyl-CoA desaturase-1 expression in skeletal muscle contributes to abnormal fatty acid partitioning in obese humans. Cell Metab. 2005;2:251–61.
Schrauwen P, Hardie DG, Roorda B, Clapham JC, Abuin A, Thomason-Hughes M, et al. Improved glucose homeostasis in mice overexpressing human UCP3: a role for AMP-kinase? Int J Obes Relat Metab Disord. 2004;28:824–8.
Clapham JC, Arch JR, Chapman H, Haynes A, Lister C, Moore GB, et al. Mice overexpressing human uncoupling protein-3 in skeletal muscle are hyperphagic and lean. Nature. 2000;406:415–8.
Nabben M, Hoeks J, Briedé JJ, Glatz JF, Moonen-Kornips E, Hesselink MK, et al. The effect of UCP3 overexpression on mitochondrial ROS production in skeletal muscle of young versus aged mice. FEBS Lett. 2008;582:4147–52.
Rasmussen BB, Holmbäck UC, Volpi E, Morio-Liondore B, Paddon-Jones D, Wolfe RR. Malonyl coenzyme A and the regulation of functional carnitine palmitoyltransferase-1 activity and fat oxidation in human skeletal muscle. J Clin Invest. 2002;110:1687–93.
Kimmel AR, Sztalryd C. Perilipin 5, a lipid droplet protein adapted to mitochondrial energy utilization. Curr Opin Lipidol. 2014;25:110–7.
Lim JH, Gerhart-Hines Z, Dominy JE, Lee Y, Kim S, Tabata M, et al. Oleic acid stimulates complete oxidation of fatty acids through protein kinase A-dependent activation of SIRT1-PGC1α complex. J Biol Chem. 2013;288(10):7117–26.
Houten SM, Wanders RJ. A general introduction to the biochemistry of mitochondrial fatty acid β-oxidation. J Inherit Metab Dis. 2010;33:469–77.
Tizioto PC, Coutinho LL, Decker JE, Schnabel RD, Rosa KO, Oliveira PS, et al. Global liver gene expression differences in Nelore steers with divergent residual feed intake phenotypes. BMC Genomics. 2015;16:242.
Zhang Y, Liu C, Zhu L, Jiang X, Chen X, Qi X, et al. PGC-1alpha inhibits oleic acid induced proliferation and migration of rat vascular smooth muscle cells. PLoS ONE. 2007;2:e1137.
Puigserver P, Spiegelman BM. Peroxisome proliferator-activated receptor-γ coactivator 1α (PGC-1α): transcriptional coactivator and metabolic regulator. Endocr Rev. 2003;24:78–90.
Fazio S. Management of mixed dyslipidemia in patients with or at risk for cardiovascular disease: a role for combination fibrate therapy. Clin Ther. 2008;30(2):294–306.
McLeod RS, LeBlanc AM, Langille MA, Mitchell PL, Currie DL. Conjugated linoleic acids, atherosclerosis, and hepatic very-low-density lipoprotein metabolism. Am J Clin Nutr. 2004;79(6 Suppl):1169S–74.
Paez JG, Janne PA, Lee JC, Tracy S, Greulich H, Gabriel S, et al. EGFR mutations in lung cancer: correlation with clinical response to gefitinib therapy. Science. 2004;304(5676):1497–500.
Wong AJ, Ruppert JM, Bigner SH, Grzeschik CH, Humphrey PA, Bigner DS, Vogelstein B. Structural alterations of the epidermal growth factor receptor gene in human gliomas. Proc Natl Acad Sci U S A. 1992;89(7):2965–9.
Smith CM, Steitz JA. Sno storm in the nucleolus: new roles for myriad small RNPs. Cell. 1997;89:669–72.
MacDonald HB. Conjugated linoleic acid and disease prevention: a review of current knowledge. J Am Coll Nutr. 2000;19:S111–8.
Park HS, Cho HY, Ha YL, Park JH. Dietary conjugated linoleic acid increases the mRNA ratio of Bax/Bcl-2 in the colonic mucosa of rats. J Nutr Biochem. 2004;15:229–35.
AMSA - American Meat Science Association. Research guidelines for cookery, sensory evaluation and tenderness measurements of fresh meat. Chicago: National Livestock and Meat Board; 1995.
Hara A, Radin NS. Lipid extraction of tissues with a low-toxicity solvent. Anal Biochem. 1978;90(1):420–6.
Christie WW. The composition and structure of milk lipids. In: Fox PF, editor. Developments in Dairy Chemistry, Volume 2. New York: Applied Science Publishers; 1983. p. 1–35.
Naylor MG, Weiss ST, Lange C. Recommendations for using standardised phenotypes in genetic association studies. Hum Genomics. 2009;3(4):308–19.
Tizioto PC, Decker JE, Taylor JF, Schnabel RD, Mudadu MA, Silva FL, et al. Genome scan for meat quality traits in Nelore beef cattle. Physiol Genomics. 2013;45(21):1012–20.
Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14(4):R36.
Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:R25.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2009;11:R106.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.
Wang C, Gong B, Bushel PR, Thierry-Mieg J, Thierry-Mieg D, Xu J, Fang H, et al. The concordance between RNA-seq and microarray data depends on chemical treatment and transcript abundance. Nat Biotechnol. 2014;32(9):926–32.
Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001;29(9):e45.
Rieu I, Powers SJ. Real-time quantitative RT-PCR: design, calculations, and statistics. Plant Cell. 2009;21(4):1031–3.
This study was conducted with funding from EMBRAPA (Macroprograma 1, 01/2005) and FAPESP (processes number 2014/11871-5, 2014/22884-0 and 2012/23638-8). LCAR and LLC were granted CNPq fellowships. The authors would like to acknowledge the collaborative efforts among EMBRAPA, University of São Paulo, and Iowa State University.
This study was conducted with funding from EMBRAPA (Macroprograma 1, 01/2005) and FAPESP (processes number 2014/11871-5, 2014/22884-0 and 2012/23638-8). LCAR and LLC were granted CNPq fellowships.
ASMC, LCAR, DPDL, RTN, RRT, JR and LLC conceived the idea of this research. ASMC and MDP participated in the experimental design and RNA-sequencing. ASMC and AMF performed the qPCR validation. ASMC, SCSA, JEK and EFW performed data analysis. ASMC drafted the manuscript. MLN, ASC and RRT conducted animal experiments. ASMC, JEK, JMR and PCT and PSNO collaborated with interpretation and discussion of the results interpretation of results and LCAR, JMR, GBM, AZN, and LLC provided the experimental environment and coordination. All authors have read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
All experimental procedures involving steers were approved by the Institutional Animal Care and Use Committee Guidelines from Brazilian Agricultural Research Corporation – EMBRAPA and sanctioned by the president Dr. Rui Machado.
Phenotypic and statistical analysis description for two groups of extreme values of fatty acids content in Longissimus dorsi muscle of Nellore steers. (XLSX 43 kb)
Test of means (t-test) between groups with high and low fatty acid content in Longissimus dorsi muscle in Nellore steers. (XLSX 44 kb)
Summary of total reads before and after filtering, percentage of mapped reads against Bos taurus UMD3.1 reference genome and, information of each group according to specifically fatty acid content from Longissimus dorsi muscle of Nellore steers. (XLSX 53 kb)
Normalized gene expression level and statistical results of differentially expressed genes analysis (FDR 10%) between High and Low groups based on linoleic acid (LA) content from Longissimus dorsi muscle of Nellore steers. (XLSX 15462 kb)
Normalized gene expression level and statistical results of differentially expressed genes analysis (FDR 10%) between High and Low groups based on stearic acid (SA) content from Longissimus dorsi muscle of Nellore steers. (XLSX 15406 kb)
Normalized gene expression level and statistical results of differentially expressed genes analysis (FDR 10%) between High and Low groups based on eicosapentaenoic acid (EPA) content from Longissimus dorsi muscle of Nellore steers. (XLSX 15447 kb)
Normalized gene expression level and statistical results of differentially expressed genes analysis (FDR 10%) between High and Low groups based on docosahexaenoic acid (DHA) content from Longissimus dorsi muscle of Nellore steers. (XLSX 13199 kb)
Normalized gene expression level and statistical results of differentially expressed genes analysis (FDR 10%) between High and Low groups based on palmitic acid (PA) content from Longissimus dorsi muscle of Nellore steers. (XLSX 15040 kb)
Normalized gene expression level and statistical results of differentially expressed genes analysis (FDR 10%) between High and Low groups based on oleic acid (OA) content from Longissimus dorsi muscle of Nellore steers. (XLSX 13046 kb)
Normalized gene expression level and statistical results of differentially expressed genes analysis (FDR 10%) between High and Low groups based on conjugated linoleic acid cis9 trans11 (CLA-c9t11) content from Longissimus dorsi muscle of Nellore steers. (XLSX 15002 kb)
Volcano plot of log2 fold-change (x-axis) versus − log10 FDR-corrected p-value in RNA-Seq data for (A) palmitic acid (PA) content, (B) docosapentaenoic acid (DHA) content, (C) eicosapentaenoic (EPA) content, (D) linoleic acid (LA) content and, (E) stearic acid (SA) content (F) conjugated linoleic acid (CLA) content. (DOCX 887 kb)
Functional enrichment and significant category (BH-adj <10%) are shown from differentially expressed genes (FDR 10%) between High and Low groups based on oleic acid (OA) content from Longissimus dorsi muscle of Nellore steers. (DOCX 127 kb)
Enriched canonical pathways of the differentially expressed genes between High and Low groups based on oleic acid content from Longissimus dorsi muscle of Nellore steers using Ingenuity Pathway Analysis (IPA). (XLS 83 kb)
IPA network analysis of differentially expressed genes analysis (FDR 10%) between High and Low groups based on oleic acid content from Longissimus dorsi muscle of Nellore steers. (XLSX 57 kb)
Uspstrem regulators analysis from DE genes list. (XLS 27 kb)
Biological functions identified by IPA from differentially expressed genes list (FDR 10%) between high and low groups based on oleic acid content in Longissimus dorsi muscle of Nellore steers. (XLSX 57 kb)
Functional enrichment and significant category (BH-adj <10%) are shown from differentially expressed genes (FDR 10%) between High and Low groups based on conjugated linoleic acid (CLA-c9t11) content from Longissimus dorsi muscle of Nellore steers. (DOCX 95 kb)
Summary of functional enrichment analysis from a list of differentially expressed genes between High and Low groups based on conjugated linoleic acid cis9 trans11 content from Longissimus dorsi muscle of Nellore steers using Ingenuity Pathway Analysis. (IPA) (XLSX 12 kb)
Cytokines and transcription regulator identified as upstream regulators by Ingenuity Pathway Analysis (IPA) as inhibited or activated from a list of differential expressed genes between High and Low groups based on conjugated linoleic acid cis9 trans11 (CLA-c9t11) content from Longissimus dorsi muscle of Nellore steers. (DOCX 95 kb)
qPCR results of nine differentially expressed genes by random selection from DataAssist software. (XLSX 49 kb)
qPCR validation of RNA-seq data on a random selection of nine differentially expressed genes. (XLSX 39 kb)