Gene co-expression network analysis for porcine intramuscular fatty acid composition

In pigs, meat quality depends markedly on the fatty acid (FA) content and composition of the intramuscular fat, which is partly determined by the gene expression in this tissue. The aim of this work was to identify the link between muscle gene expression and its FA composition. In an (Iberian × Duroc) × Duroc backcrossed pig population, we identified modules of co-expressed genes, and correlation analyses were performed for each of them versus the phenotypes, finding four relevant modules. Two of the modules were positively correlated with saturated FAs (SFAs) and monounsaturated FAs (MUFAs), while negatively correlated with polyunsaturated FAs (PUFAs) and the omega-6/omega-3 ratio. The gene-enrichment analysis showed that these modules had over-representation of pathways related with the biosynthesis of unsaturated FAs, the Peroxisome proliferator-activated receptor signalling pathway and FA elongation. The two other relevant modules were positively correlated with PUFA and the n-6/n-3 ratio, but negatively correlated with SFA and MUFA. In this case, they had an over-representation of pathways related with fatty and amino acid degradation, and with oxidative phosphorylation. Using a graphical Gaussian model, we inferred a network of connections between the genes within each module. The first module had 52 genes with 87 connections, and the most connected genes were ADIPOQ, which is related with FA oxidation, and ELOVL6 and FABP4, both involved in FA metabolism. The second module showed 196 genes connected by 263 edges, being FN1 and MAP3K11 the most connected genes. On the other hand, the third module had 161 genes connected by 251 edges and ATG13 was the top neighbouring gene, while the fourth module had 224 genes and 655 connections, and its most connected genes were related with mitochondrial pathways. Overall, this work successfully identified relevant muscle gene networks and modules linked with FA composition, providing further insights on how the physiology of the pigs influences FA composition.