A porcine gluteus medius muscle genome-wide transcriptome analysis: dietary effects of omega-6 and omega-3 fatty acids on biological mechanisms
© The Author(s) 2017
Received: 12 October 2016
Accepted: 9 January 2017
Published: 31 January 2017
The level of omega-6 and omega-3 polyunsaturated fatty acids can affect many cellular systems and function via nuclear receptors or the bioactive lipid regulation of gene expression. The objective of this study was to investigate changes in the muscle transcriptome and the biological functions regulated by increased consumption of omega-3 and omega-6 fatty acids in the pig gluteus medius muscle.
The transcriptome of the gluteus medius muscle was studied for pigs subjected to either a control diet or a diet supplemented with linseed and rapeseed oil to increase polyunsaturated fatty acid content. Next-generation sequencing (NGS) was used to generate the muscle tissue transcriptome database pointing differentially expressed genes (DEG). Comparative expression analyses identified 749 genes significantly differing at least in the twofold of change between two groups of animals fed with divergent level of omega-3 and omega-6 fatty acids. The expression of 219 genes was upregulated, and the expression of 530 genes was downregulated in the group of pigs supplemented with omega-3 and omega-6 fatty acids in relation to control group pigs. Results of RNA-seq indicated a role of fatty acid in the regulation of the expression of genes which are essential for muscle tissue development and functioning. Functional analysis revealed that the identified genes were important for a number of biological processes including inflammatory response, signaling, lipid metabolism, and homeostasis.
Summarizing, obtained results provide strong evidence that omega-6 and omega-3 fatty acids regulate fundamental metabolic processes in muscle tissue development and functioning.
KeywordsGluteus medius muscle Omega-3 fatty acid Omega-6 fatty acid Pig Transcriptome
The expression of animal genes can be modified in response to environmental factors. One such factor is feed, especially the composition of its bioactive dietary components . The newly emerging field of nutrigenomics investigates the relationship between dietary components and genomic expression profiles . Recently, nutrigenomic research has focused much attention on omega-6 and omega-3 fatty acids.
Omega-6 and omega-3 fatty acids are polyunsaturated fatty acids (PUFA), which belong to the essential fatty acid family and must be obtained from food sources . Mammals can, however, metabolize PUFA by utilizing the omega-6 and omega-3 fatty acid conversion pathways to synthesize other omega fatty acids (i.e., arachidonic acid from linoleic acid (LA), as well as eicosapentaenoic acid and docosahexaenoic acid from alpha-linolenic acid (ALA)) . Omega-6 and omega-3 fatty acids modulate gene expression and provide substrates for the production of signaling molecules or functioning mediators . Additionally, a correct balance of these fatty acids is important for the proper function and development of various cell types . Presence of omega-6 and omega-3 fatty acids and maintaining their optimal ratio in the diet helps to prevent the development of metabolic diseases, such as heart attacks, atherosclerosis, thrombosis, arrhythmia, stroke, immune-inflammatory disorders, asthma, arthritis, cancer proliferation, type II diabetes mellitus, obesity, and psychiatric disorders . Although PUFA are known to have significant biological roles, their mechanisms are not well understood.
Fatty acids are an important cellular component and serve as a source of energy for animals . A major site of fatty acid catabolism in mammals is the skeletal muscle . In fact, it has been suggested that PUFA content can modulate both the function and metabolism of the skeletal muscle. The effect of PUFA on the skeletal muscle also has the potential to substantially impact whole body functioning . PUFA content is important for muscle tissue as well, regulating immunological, hormonal, and metabolic processes . Gluteus medius muscle was selected as the target muscle tissue for this study due to its dynamic energy metabolism and role in independent gait . The gluteus medius muscle contains IIb and IIx fibers as well as IIa (fast oxidative-glycolytic) fibers , which can switch between glycolytic and oxidative metabolism in reaction to accessibility of the fatty acids . Although the immune-stimulating effects of PUFA have been studied for decades , knowledge of their integrative effects is limited. Therefore, the objective of this study was to investigate changes in the pig gluteus medius muscle transcriptome, including the regulated biological functions as a result of increased omega-3 and omega-6 fatty acid concentration. In our studies, we used pigs, which are one of the best models for studies on effects of omega-6 and omega-3 fatty acids supplementation. Moreover, they are a proven and recognized large animal model due to their anatomical and physiological similarities to humans. We hypothesize that supplementation of omega-6 and omega-3 fatty acids in the diet will modulate the expression of genes essential for muscle homeostasis.
Composition of pigs diets
Gas chromatography determination of the ratio of omega-6 to omega-3 fatty acid
Quantitative data for the ratio of omega-6 to omega-3 fatty acid in all gluteus medius muscle samples were obtained by gas chromatography using flame ionization detection (GC–FID) analysis. For each sample, 1 g of muscle tissue was homogenized and extracted with chloroform-methanol 2:1(v/v) according to Folch et al. . Fatty acids were analyzed using a GC-7890 Agilent gas chromatograph (Agilent Technologies, Inc., Santa Clara, CA,) equipped with a 60-m capillary column (Hewlett-Packard-88, Agilent J&W GC Columns, Santa Clara, CA) with 0.25-mm inner diameter and coating thickness of 0.20 μm. Helium was used as a carrier gas at a flow rate of 50 mL min−1. The injector and detector were both maintained at 260 °C. Column oven temperature was programmed to increase from 140 °C (held for 5 min) at a rate of 4 °C min−1 to 190 °C to 215 °C at a rate of 0.8 °C min−1. Individual fatty acid peaks were identified by comparison with known reference methyl esters (Supelco 37 Component FAME Mix, 47885-U, Sigma-Aldrich Co. Warsaw, Poland) and expressed as a percentage of total fatty acid concentration .
To prepare samples for next-generation sequencing , 20 mg of frozen tissue were pulverized with a mortar and pestle precooled with liquid nitrogen. The RNA was isolated using an mRNA extraction kit (RNeasy Fibrous Tissue Mini Kit—Qiagen, Crawley, UK) according to the manufacturer’s recommendations. For real-time PCR, 20 mg of frozen tissue was homogenized in 800 μl of TRIzol reagent (Invitrogen, Carlsbad, CA) using tubes with ceramic beads on Magna Lyser (Roche Applied Science, Penzberg, Germany). The extracts of total cellular RNA were isolated from homogenate according to Chomczynski and Sacchi . The concentration and purity of total RNA were checked using a UV-Vis spectrophotometer (Actgene, Piscataway, NJ) at 260 nm. The integrity of ribosomal RNA was evaluated electrophoretically on the Bioanalyser (Agilent Technologies Inc.). Only high-quality RNA isolates, characterized by an RNA integrity number greater than 7, were used for further analysis.
Equal amounts of total RNA of six randomly selected individuals per group were used to construct libraries for sequencing. The mRNA-Seq library was constructed and sequenced by Genomed (http://www.genomed.pl/, Poland)  in accordance with the manufacturer’s instructions for the Illumina MiSeq2000 system (Illumina, Inc., San Diego, CA, USA). The TruSeq RNA Sample Pre Kit (Illumina Inc.) was used to construct the libraries. To remove rRNA, bead-based depletion was used. Libraries were synthetized twice from freshly prepared RNA pools (each pool contained material from six pigs) and were sequenced independently on the MiSeq high-throughput sequencing instrument with 150 paired-end sequencing. Obtained sequence reads were aligned to the Sus scrofa 10.2 genome (NCBI GenBank ID GCA_000003025.4). Preprocessing and mapping was performed using the CLC Genomics workbench (CLC bio, Denmark) and DNASTAR Lasergene Genomic Suite (DNASTAR, Inc., Madison, WI). The read mapping was conducted using the RNA-Seq analysis package default settings (mismatch cost 2, insertion cost 3, length fraction 0.8, similarity fraction 0.8, and maximum number of reads 1) within the CLC software. The number of mapped reads was normalized and measured in terms of reads per kilobase per million (RPKM) to determine transcript abundance. For this study, pooled sample was used in order to remove animal-specific differences, allowing for the general reaction mechanism to be highlighted (for comparison, see te Pas et al.  and Szostak et al. ). As there was a single sample comprised of pooled individuals, in each condition, genes with the highest expression variance were estimated by using the variance across the two conditions. This is known as a “blind” method of dispersion estimate. The analysis assumes that for most genes, there is no true differential expression and that a valid mean-variance relationship can be estimated from treating the two samples as if they were replicates what will lead to emergence of a minority of differentially abundant genes acting as outliers. This method results in overly strict statistical significance testing of differential gene expression . Genes that showed at least a twofold change in expression between control and experimental group and had a t test P value ≤0.05 adjusted for false discovery rate (FDR) were considered for detailed functionally grouped network analysis using DAVID software and Cytoscape 3.2 ClueGO plugin.
Reverse transcription PCR
DNase digestion was performed with RQ1 RNase-free DNase (Promega Corp., Madison, WI) by incubating samples with DNase for 15 min at 20 °C. One microgram of total RNA was reverse transcribed into cDNA using oligo-dT primers and the Transcriptor First Strand cDNA Synthesis Kit (Roche Applied Science) following the manufacturer’s recommendations.
Real-time PCR primer
Product length (bp)
Retinol-binding protein 4
ADAM metallopeptidase domain 19
Fatty acid desaturase 1
Fatty acid desaturase 2
Adenosine A1 receptor
Mitogen-activated protein kinase 1
Fas (TNFRSF6)-associated via death domain
C-C motif chemokine ligand 4
Topoisomerase (DNA) II beta
Functional bioinformatics analysis
Differentially expressed genes (DEG) that showed at least a twofold change in expression between control and experimental group and had a t test P value ≤0.05 adjusted for FDR were selected for DAVID analysis (The Database for Annotation, Visualization, and Integrated Discovery: https://david.ncifcrf.gov/ ). DAVID software calculates P values for the enrichment of the number of genes in biological mechanisms. These P values indicate the importance of the general biological mechanisms active in all animals but do not indicate the inter-individual differences. All DEG were analyzed. The gene list was transferred to NCBI’s ENTREZ search engine of cross-database, which recognized 493 unique gene names. Of these 493 genes, 481 of the gene IDs were transferred to human gene IDs, to compensate for poor pig gene annotation. Further analyses were performed using the total human genome as a background. An error rate of 0.05 was assumed significant according to the Benjamini and Hochberg multiple test correction method, and a FDR less than 5% was applied.
To further elucidate the biological relevance of the DEG, a functional enrichment analysis was performed using ClueGO software, a Cytoscape 3.1.0 plug-in (Institute for Genomics and Bioinformatics, Graz University of Technology, Graz, Austria) . ClueGO facilitates the visualization of functionally related genes by displaying the genes as a clustered network and chart. The statistical test used to determine the enrichment score for terms and groups was based on a right-sided hypergeometric distribution option with a Benjamini-Hochberg correction. A kappa statistic was calculated to determine the strength of the connection between the terms in the network, based on similarity of associated genes . To obtain complex, clear, and readable networks, various kappa scores were applied. The kappa score was preset to 0.9 for analyses covering a complete list of genes. A minimum level of 7 and a maximum level of 15 were set as the GO level interval with a minimum of two genes per category.
There were no significant differences in carcass nor meat quality traits between pigs of the two dietary groups. For animals from the experimental group, we observed a tendency in the reduction of backfat thickness and fodder intake.
Fatty acid profile [g/100 g FAME] of the muscle of pigs from two examined groups (mean ± SD)
1.39 ± 0.18
1.51 ± 0.32
23.77 ± 1.35
25.17 ± 1.67
13.96 ± 0.54
13.75 ± 0.57
0.27 ± 0.13
0.11 ± 0.07
2.78 ± 0.34
3.34 ± 0.42
36.23 ± 1.93
38.92 ± 2.02
3.34 ± 0.32
3.92 ± 0.31
9.44 ± 1.49
8.14 ± 2.46
4.45 ± 1.09
0.57 ± 0.52
1.24 ± 0.24
1.36 ± 0.53
Ʃ saturated FA
41.85 ± 1.30
42.87 ± 1.46
Ʃ monounsaturated FA
42.77 ± 2.41
46.48 ± 2.54
Ʃ polyunsaturated FA
15.27 ± 2.66
10.19 ± 2.20
Ʃ n-6 FA
10.69 ± 1.67
9.51 ± 2.79
Ʃ n-3 FA
4.47 ± 1.14
0.89 ± 0.55
2.57 ± 1.09
11.57 ± 8.81
NGS was performed to identify changes induced by omega-6 and omega-3 supplementation in the pig muscle transcriptome. The number of reads per sample obtained from NGS was 25,244,670. Among these reads, 16,040,530 (63.54%) were mapped in pairs and 792,066 (3.14%) in broken pairs and 8,412,074 (33.32%) were unmapped. The number of counted fragments was 8,020,265 (63.54%); among these fragments, 7,475,220 (59.22%) were unique and 545,04 (4,32%) were not specific. The number of uncounted fragments was 4,602,070 (36.46%). The sequence reads have been submitted to the NCBI Gene Expression Omnibus under accession series GSE83729. The expression of 14246 genes was registered in the muscle tissue of pigs subjected to NGS. An Additional file 1 shows the heat maps of the transcription levels between control and experimental group. Between genes complying the condition of significance (at least a twofold change in expression between control and experimental group and a t test P value ≤0.05 adjusted for FDR), expression of 219 genes was upregulated and the expression of 530 genes was downregulated.
Real-time PCR verification of NGS data
Functional analysis by DAVID
Results of DAVID functional clustering analysis of differentially expressed genes (DEG)
Number of DEG
Response to wounding/inflammation
Response to endogenous stimulus and organic substance
Regulation of response to external stimulus
Complement and coagulation cascades
Retinol binding/vitamin A
Negative regulation of response to stimulus
Acute inflammatory response
Functional analysis by Cytoscape
The objective of this study was to determine the impact of omega-6 and omega-3 fatty acid concentration on the pig muscle transcriptome. Results provide new insight regarding the possible underlying regulatory mechanisms of the muscle tissue’s biological response to increased concentration of omega-6 and omega-3 fatty acids.
RNA-seq of muscle tissue identified 759 genes with at least a twofold change in gene expression. The number of downregulated genes in pigs characterized by decreased omega-6 to omega-3 ratio in the muscle (for the supplemented diet) was more than two times greater than the number of upregulated genes in pigs fed with the control diet. Duan et al.  reported a dietary ratio of omega-6/omega-3 fatty acids of 1:1 and 5:1, which suppressed gene expression levels of lipid metabolism-related genes and proteins, as well as inflammatory cytokines. The authors concluded that their reported ratio of omega-6 to omega-3 fatty acids is beneficial for lipid metabolism and function of the inflammatory system. However, until now, a few studies have investigated the effect of fatty acid consumption on the muscle transcriptome using high-throughput technology.
Omega-6 and omega-3 fatty acids regulate the inflammatory response by nuclear receptors, which regulate gene expression  and synthesize pro- and anti-inflammatory biomarkers. Indeed, the most affected process pointed out both by DAVID and Cytoscape-ClueGo was the response to wounding and inflammation. We observed decreased expression of a number of chemokines: CCL2, CCL4, CCL5, CCL16, CCL19, and CCL21, recruiting macrophages and monocytes to the site of inflammation. That may indicate a diminished activity of macrophages and monocytes themselves. Similar results were obtained by Liang et al. , who registered decreased expression of CCL2 after the administration of omega-3 fatty acid to mice. Further, we observed diminished expression of several complement components, C4BPA, C5, C7, C8G, and C9, that are components of the complement system, a part of the innate immune system. The usefulness of polyunsaturated fatty acid immunomodulatory activities is well known. Their possible impact on a group of chemokines and complement components explore mechanisms of its acting. That supports their utility in the management of inflammatory and autoimmune diseases.
Among genes with divergent expression that are involved in coagulation, two coagulation factors—F2 and F5—showed reduced expression. These genes encode essential enzymes of the blood coagulation cascade. There are a number of studies on the impact of polyunsaturated fatty acids on coagulation process; however, their results are inconsistent. It is suggested that increased consumption of polyunsaturated fatty acids decreased the blood levels of several coagulation factors and may prevent thrombosis [29–31]. Here, we suggest that the possible mechanism of this phenomenon could be related to diminished levels of F2 and F5 due to increased consumption of omega-6 and omega-3 fatty acids.
The significant number of affected genes was related to signal transduction. Long-chain fatty acids can function as lipid second messengers of signaling cascade systems or as modulators and regulators of intracellular signaling pathways. The specific functions depend on the bioavailability of fatty acids (governed by their interaction with soluble proteins capable of binding fatty acids) and the lipid and protein components of biological membranes, including membrane-associated fatty acid-binding proteins . In our studies, affected genes were related to membrane signaling components, such as receptors, ion channels, and membrane glycoproteins. Among them, genes belonging to several families deserve to be highlighted. The differentially expressed genes (DEG) of the ADAM family play an important role in signaling, cell adhesion, and migration. This finding was consistent with results reported by Reiss et al. , suggesting that free fatty acids are critical regulators of ADAM family gene function. Another family of genes is signaling membrane proteins. Our results showed that encoded by genes differentially expressed between experimental groups are G protein-coupled receptors and solute carrier proteins. Furthermore, a significant number of DEG were involved in ion transport, especially responsible for calcium and iron which showed different expression between examined groups. The regulation of the free ion concentration is not only important for signaling pathways but is also fundamental in the proper functioning of cells and the whole body. This observation confirmed the role of omega-6 and omega-3 fatty acids in ion and overall homeostasis. Likewise, affected genes were involved in the regulation of Cdc42 protein that is involved in a variety of signaling events and a number of fundamental for proper muscle functioning cellular pathways, determining cell growth, differentiation, and apoptosis. Previous genome-wide studies on pig muscle transcriptome, where compared pigs were divergent in fatty acid composition in the muscle, also pointed out the pathway dedicated to the signaling as one of the most affected . The biological functions of polyunsaturated fatty acids are inter alia, due to their ability to act as secondary messengers or modulators of activities of functionally important proteins. Here, we present that adaptation to omega-6 and omega-3 supplementation is accomplished through signal transduction. Taken together, this finding indicates that the regulatory mechanisms underlying the biological effects of omega-6 and omega-3 fatty acids seem to be strictly connected to signaling pathways.
We have detected differential expression of genes related to the lipid binding between the control and experimental groups. Noteworthy are our observation of the downregulation of many apolipoprotein genes, including APOA1, APOA2, APOA4, APOA5, APOC3, APOE, and APON. Previous studies proved that ALA decreased the level of triacylglycerol, total cholesterol, and low-density lipoprotein (LDL) cholesterol, as well as lipoprotein-A, APOA1, and apolipoprotein B (APOB) . Likewise, Fox et al.  perceived that the supplementation with oils containing higher content of PUFA resulted in longer term effect on decreasing the level of APOA1 transcript in comparison to oils containing lower levels of PUFA.
Many other observations made by scientists lead to the conclusion that the control of lipoprotein metabolism through proper supplementation with bioactive components is the essential part of metabolic diseases prevention. A number of studies have shown that polyunsaturated fatty acids have a marked hypolipidemic effect. In our studies, we have detected differential expression of genes related to the lipid-binding protein between the groups of pigs. This is consistent with other research made by Canovas et al. , where animals showing higher levels of PUFA in the muscle were characterized by decreased level of serum lipids: total cholesterol, high-density lipoprotein (HDL) cholesterol, LDL cholesterol, and triglycerides. The diminished concentration of mentioned lipids could be related with observed in our studies downregulation of apolipoproteins.
APOA1 and APOA2 are the major proteins constituent of HDLs. While APOA4 participates in intestinal lipid absorption, APOA5 regulates VLDLs synthesis and secretion . Changes in their expression may suggest regulation of this particles metabolism. Results of previous studies on the impact of PUFA on the expression of APOA family members are consistent with our findings [35, 39]. It was previously suggested that APOA level can be regulated by PUFA through molecular mechanism or indirect mechanism via changes in plasma lipid or lipoprotein concentration . Nevertheless, additional studies are needed to clarify the mechanism by which omega-6 and omega-3 fatty acids affect APOA expression and further the lipoprotein metabolism.
APOC3 is a component of chylomicrons, both very low and HDLs; inhibits the activity of lipoprotein lipase; and delays very LDLs clearance. The reduction in APOC3 expression observed in the group supplemented with omega-3 and omega-6 fatty acid may have been a result of a direct effect of omega-3 fatty acids on a reducing APOC3 production or as an indirect result of synthesis and secretion of fewer very LDL particles. Increased concentrations of APOC3 are an indicator of a risk factor for coronary heart disease and atherosclerotic lesions . Therefore, it could be hypothesized that reduction in APOC3 that we observed after omega-6 and omega-3 fatty acids supplementation could have implications for coronary heart disease risk.
In accordance with our results, Canovas et al.  observed decreased expression of APOE in pigs containing increased levels of PUFA in the muscle. APOE removes chylomicron’s remnants from circulation through interaction with LDL receptors. It was proven that APOE knockout mice exhibit less body fat stores and smaller adipocytes and are more resistant to diet-induced obesity . Possibly, decreased APOE expression that is thought to facilitate obesity decreased risk of this disorder.
It became evident that omega-6 and omega-3 fatty acids through regulation of lipoproteins play the essential role in postprandial metabolism. Therefore, they are associated with a lower risk of atherogenesis, heart disease, obesity, and diabetes. Here, we reported that omega-6 and omega-3 fatty acids affect a number of genes related to lipid transport. An important group of regulated genes belongs to apolipoprotein family. Our results highlight the need of further studies that would broaden knowledge about interactions between PUFA, lipoproteins, and metabolic disorders.
Finally, among the pathways overrepresented between both groups of animals was homeostasis. PUFA are powerful endogenous bioregulators, affecting whole body functioning and therefore homeostasis. In fact, they strongly regulate the muscle, an important organ that adjusts systemic metabolic homeostasis. First and foremost, we have observed significant changes in the expression of genes related to ion homeostasis. Homeostatic regulation of ionic gradients is critical for most functions. Here, we reported that omega-6 and omega-3 fatty acids impact expression of a number of genes involved in ion transport, including ion channels, solute carriers, G protein-coupled receptors, and tyrosine kinases. Possibly, through regulation of genes that control calcium, iron sodium and potassium ion concentration which examined fatty acids affect fundamental regulatory system of cellular processes.
Among the strengths of the present study is the scale of NGS analysis that allowed for whole genome sequencing and validation of NGS results with real-time PCR. In addition, pigs as a choice of the model organism are one of the best models for studies on effect of omega-6 and omega-3 fatty acids, which improves studies. However, limitation of lack of biological replicates should be considered in the interpretation of our results.
In our study, we compared the muscle transcriptome of pigs fed with a control diet to that of pigs fed with a diet supplemented with additional fatty acids, causing different muscular fatty acid profiles. Our results showed that the omega-3 and omega-6 fatty acids regulate cellular muscle metabolism at the transcriptome level. This showed that these fatty acids are regulators of fundamental metabolic processes in muscle tissue development and functioning. A total of 749 unique genes were differently expressed between the two diet groups. The expression of 219 genes was upregulated, and the expression of 530 genes was downregulated in the group of pigs supplemented with omega-3 and omega-6 fatty acids in relation to expression in the control group pigs. The roles of these genes were indicated based on their involvement in particular biological functions. Among these biological processes were basic physiological processes such as inflammatory response, signaling pathways and lipid metabolism, leading to a finding that influence of omega-6 and omega-3 fatty on the skeletal muscle substantially impacts the whole body functioning. Obtained results reveal additional information regarding transcriptome-level regulation of biological processes underlying differences in omega-6 and omega-3 fatty acid levels.
This study was funded by the project “BIOFOOD – innovative, functional products of animal origin” No. POIG.01.01.02-014-090/09 and was co-financed by the European Union from the European Regional Development Fund under Innovative Economy Operational Program 2007–2013.
Availability of data and materials
The datasets generated during the current study are available in the Gene Expression Omnibus data repository under accession series GSE83729.
MO participated in the laboratory analysis, analyzed and interpreted the data, and prepared the draft of the manuscript. AS participated in the laboratory analysis and data analysis. MFWTP participated in the data analysis and coordinated the manuscript writing. EP made substantial contributions to the conception and preparation of the manuscript. PU conducted the GC analysis. TB and JOH participated in designing the experiment. CSP and EJK gave comments to the manuscript. JRD participated in the manuscript creation. MP designed and coordinated the experiment and participated in the bioinformatic analyses. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
The Local Ethics Commission for Experimentation on Animals in Warsaw (Poland) approved all experiments involving animals (approval number: 27/2009). All applicable institutional and national guidelines for the care and use of animals were followed.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Pegorier JP, Le May C, Girard J. Control of gene expression by fatty acids. J Nutr. 2004;134:2444S–9.PubMedGoogle Scholar
- Ouhtit A. Nutrigenomics: from promise to practice. Sultan Qaboos Univ Med J. 2014;14:e1–3.View ArticlePubMedPubMed CentralGoogle Scholar
- Chondrogianni N, Sakellari M, Lefaki M, Papaevgeniou N, Gonos ES. Proteasome activation delays aging in vitro and in vivo. Free Radic Biol Med. 2014;71:303–20.View ArticlePubMedGoogle Scholar
- Liu JJ, Green P, John Mann J, Rapoport SI, Sublette ME. Pathways of polyunsaturated fatty acid utilization: implications for brain function in neuropsychiatric health and disease. Brain Res. 2015;1597:220–46.View ArticlePubMedGoogle Scholar
- Kang JX, Liu A. The role of the tissue omega-6/omega-3 fatty acid ratio in regulating tumor angiogenesis. Cancer Metastasis Rev. 2013;32:201–10.View ArticlePubMedGoogle Scholar
- Duan Y, Li F, Li L, Fan J, Sun X, Yin Y. n-6:n-3 PUFA ratio is involved in regulating lipid metabolism and inflammation in pigs. Br J Nutr. 2014;111:445–51.View ArticlePubMedGoogle Scholar
- Lands B. Consequences of essential fatty acids. Nutrients. 2012;4:1338–57.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang W, Zhang J, Cui L, Ma J, Chen C, Ai H, et al. Genetic architecture of fatty acid composition in the longissimus dorsi muscle revealed by genome-wide association studies on diverse pig populations. Genet Sel Evol. 2016;48:5.View ArticlePubMedPubMed CentralGoogle Scholar
- Jeukendrup AE. Regulation of fat metabolism in skeletal muscle. Ann N Y Acad Sci. 2002;967:217–35.View ArticlePubMedGoogle Scholar
- Jeromson S, Gallagher IJ, Galloway SD, Hamilton DL. Omega-3 fatty acids and skeletal muscle health. Mar Drugs. 2015;13:6977–7004.View ArticlePubMedPubMed CentralGoogle Scholar
- Kohlstadt I. Advancing medicine with food and nutrients. 2nd ed. Boca Raton: CRC Press, Taylor & Francis Group; 2013.Google Scholar
- Lee SK, Lee SY, Jung JM. Muscle activity of the gluteus medius at different gait speeds. J Phys Ther Sci. 2014;26:1915–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Bee G, Solomon MB, Czerwinski SM, Long C, Pursel VG. Correlation between histochemically assessed fiber type distribution and isomyosin and myosin heavy chain content in porcine skeletal muscles. J Anim Sci. 1999;77:2104–11.View ArticlePubMedGoogle Scholar
- Borer KT. Advanced exercise endocrinology. Champaign: Human Kinetics; 2013.Google Scholar
- Weylandt KH, Serini S, Chen YQ, Su HM, Lim K, Cittadini A, et al. Omega-3 polyunsaturated fatty acids: the way forward in times of mixed evidence. Biomed Res Int. 2015;2015:143109.PubMedPubMed CentralGoogle Scholar
- Szostak A, Ogluszka M, te Pas MFW, Polawska E, Urbanski P, Juszczuk-Kubiak E, et al. Effect of a diet enriched with omega-6 and omega-3 fatty acids on the pig liver transcriptome. Genes Nutr. 2016;11:9.Google Scholar
- Folch J, Lees M, Sloane Stanley GH. A simple method for the isolation and purification of total lipides from animal tissues. J Biol Chem. 1957;226:497–509.PubMedGoogle Scholar
- Polawska E, Horbanczuk JO, Pierzchala M, Strzalkowska N, Jozwik A, Wojcik A, et al. Effect of dietary linseed and rapeseed supplementation on fatty acid profiles in the ostrich. Part 1. Muscles. Anim Sci Paper Rep. 2013;31:239–48.Google Scholar
- Hicks AA, Pramstaller PP, Johansson A, Vitart V, Rudan I, Ugocsai P, et al. Genetic determinants of circulating sphingolipid concentrations in European populations. PLoS Genet. 2009;5:e1000672.View ArticlePubMedPubMed CentralGoogle Scholar
- Chomczynski P, Sacchi N. Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. Anal Biochem. 1987;162:156–9.View ArticlePubMedGoogle Scholar
- Genomed. http://www.genomed.pl/. Accessed 10 Oct 2016.
- te Pas MFW, Wijnberg ID, Hoekman AJW, de Graaf-Roelfsema E, Keizer HA, van Breda E, et al. Skeletal muscle transcriptome profiles related to different training intensities and detraining in Standardbred horses: a search for overtraining biomarkers. Vet J. 2013;197:717–23.View ArticleGoogle Scholar
- Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11:R106.View ArticlePubMedPubMed CentralGoogle Scholar
- Ramakers C, Ruijter JM, Deprez RH, Moorman AF. Assumption-free analysis of quantitative real-time polymerase chain reaction (PCR) data. Neurosci Lett. 2003;339:62–6.View ArticlePubMedGoogle Scholar
- da Huang W, Sherman BT, Stephens R, Baseler MW, Lane HC, Lempicki RA. DAVID gene ID conversion tool. Bioinformation. 2008;2:428–30.View ArticlePubMed CentralGoogle Scholar
- Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, Kirilovsky A, et al. ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics. 2009;25:1091–3.View ArticlePubMedPubMed CentralGoogle Scholar
- Siriwardhana N, Kalupahana NS, Fletcher S, Xin WT, Claycombe KJ, Quignard-Boulange A, et al. n-3 and n-6 polyunsaturated fatty acids differentially regulate adipose angiotensinogen and other inflammatory adipokines in part via NF-kappa B-dependent mechanisms. J Nutr Biochem. 2012;23:1661–7.View ArticlePubMedGoogle Scholar
- Liang P, Henning SM, Schokrpur S, Wu L, Doan N, Said J, et al. Effect of dietary omega-3 fatty acids on tumor-associated macrophages and prostate cancer progression. Prostate. 2016;76:1293–302.View ArticlePubMedGoogle Scholar
- Kristensen SD, Iversen AM, Schmidt EB. n-3 polyunsaturated fatty acids and coronary thrombosis. Lipids. 2001;36 Suppl:S79–82.View ArticlePubMedGoogle Scholar
- Shahar E, Folsom AR, Wu KK, Dennis BH, Shimakawa T, Conlan MG, et al. Associations of fish intake and dietary n-3 polyunsaturated fatty acids with a hypocoagulable profile. The Atherosclerosis Risk in Communities (ARIC) Study. Arterioscler Thromb. 1993;13:1205–12.View ArticlePubMedGoogle Scholar
- Miller GJ. Effects of diet composition on coagulation pathways. Am J Clin Nutr. 1998;67:542S–5.PubMedGoogle Scholar
- Glatz JFC, Luiken JJFP. Fatty acids in cell signaling: historical perspective and future outlook. Prostaglandins Leukot Essent Fat Acids. 2015;92:57–62.View ArticleGoogle Scholar
- Reiss K, Cornelsen I, Husmann M, Gimpl G, Bhakdi S. Unsaturated fatty acids drive disintegrin and metalloproteinase (ADAM)-dependent cell adhesion, proliferation, and migration by modulating membrane fluidity. J Biol Chem. 2011;286:26931–42.View ArticlePubMedPubMed CentralGoogle Scholar
- Corominas J, Ramayo-Caldas Y, Puig-Oliveras A, Estelle J, Castello A, Alves E, et al. Analysis of porcine adipose tissue transcriptome reveals differences in de novo fatty acid synthesis in pigs with divergent muscle fatty acid composition. BMC Genomics. 2013;14:843.View ArticlePubMedPubMed CentralGoogle Scholar
- Poudyal H, Panchal SK, Diwan V, Brown L. Omega-3 fatty acids and metabolic syndrome: effects and emerging mechanisms of action. Prog Lipid Res. 2011;50:372–87.View ArticlePubMedGoogle Scholar
- Fox JC, McGill Jr HC, Carey KD, Getz GS. In vivo regulation of hepatic LDL receptor mRNA in the baboon. Differential effects of saturated and unsaturated fat. J Biol Chem. 1987;262:7014–20.PubMedGoogle Scholar
- Canovas A, Quintanilla R, Amills M, Pena RN. Muscle transcriptomic profiles in pigs with divergent phenotypes for fatness traits. BMC Genomics. 2010;11:372.View ArticlePubMedPubMed CentralGoogle Scholar
- Dominiczak MH, Caslake MJ. Apolipoproteins: metabolic role and clinical biochemistry applications. Ann Clin Biochem. 2011;48:498–515.View ArticlePubMedGoogle Scholar
- Bloedon LT, Balikai S, Chittams J, Cunnane SC, Berlin JA, Rader DJ, et al. Flaxseed and cardiovascular risk factors: results from a double blind, randomized, controlled clinical trial. J Am Coll Nutr. 2008;27:65–74.View ArticlePubMedGoogle Scholar
- Berthou L, Saladin R, Yaqoob P, Branellec D, Calder P, Fruchart JC, et al. Regulation of rat liver apolipoprotein A-I, apolipoprotein A-II and acyl-coenzyme A oxidase gene expression by fibrates and dietary fatty acids. Eur J Biochem. 1995;232:179–87.View ArticlePubMedGoogle Scholar
- Davidson MH, Maki KC, Bays H, Carter R, Ballantyne CM. Effects of prescription omega-3-acid ethyl esters on lipoprotein particle concentrations, apolipoproteins AI and CIII, and lipoprotein-associated phospholipase A(2) mass in statin-treated subjects with hypertriglyceridemia. J Clin Lipidol. 2009;3:332–40.View ArticlePubMedGoogle Scholar