Liver transcriptome profile in pigs with extreme phenotypes of intramuscular fatty acid composition
© Ramayo-Caldas et al.; licensee BioMed Central Ltd. 2012
Received: 2 March 2012
Accepted: 8 October 2012
Published: 11 October 2012
New advances in high-throughput technologies have allowed for the massive analysis of genomic data, providing new opportunities for the characterization of the transcriptome architectures. Recent studies in pigs have employed RNA-Seq to explore the transcriptome of different tissues in a reduced number of animals. The main goal of this study was the identification of differentially-expressed genes in the liver of Iberian x Landrace crossbred pigs showing extreme phenotypes for intramuscular fatty acid composition using RNA-Seq.
The liver transcriptomes of two female groups (H and L) with phenotypically extreme intramuscular fatty acid composition were sequenced using RNA-Seq. A total of 146 and 180 unannotated protein-coding genes were identified in intergenic regions for the L and H groups, respectively. In addition, a range of 5.8 to 7.3% of repetitive elements was found, with SINEs being the most abundant elements. The expression in liver of 186 (L) and 270 (H) lncRNAs was also detected. The higher reproducibility of the RNA-Seq data was validated by RT-qPCR and porcine expression microarrays, therefore showing a strong correlation between RT-qPCR and RNA-Seq data (ranking from 0.79 to 0.96), as well as between microarrays and RNA-Seq (r=0.72). A differential expression analysis between H and L animals identified 55 genes differentially-expressed between groups. Pathways analysis revealed that these genes belong to biological functions, canonical pathways and three gene networks related to lipid and fatty acid metabolism. In concordance with the phenotypic classification, the pathways analysis inferred that linolenic and arachidonic acids metabolism was altered between extreme individuals. In addition, a connection was observed among the top three networks, hence suggesting that these genes are interconnected and play an important role in lipid and fatty acid metabolism.
In the present study RNA-Seq was used as a tool to explore the liver transcriptome of pigs with extreme phenotypes for intramuscular fatty acid composition. The differential gene expression analysis showed potential gene networks which affect lipid and fatty acid metabolism. These results may help in the design of selection strategies to improve the sensorial and nutritional quality of pork meat.
Pigs, an important source of human food, accounting for over 40% of the meat produced worldwide. In addition, due to the similarities in anatomy and physiology with humans, they have been used in biomedicine as an important animal model for the study of the genetic basis of metabolic diseases such as obesity, type II diabetes, metabolic syndrome and atherosclerosis. As well it is often mentioned as the preferred animal species for organ xenotransplantation [1, 2].
Over the last decade, a growing awareness of the association between diet and health has led nutritional quality to become a relevant factor in consumers’ food choices. A major development has been the recognition that certain fatty acids (FA), such as oleic acid, and α-linolenic acid (ALA), can improve human health status and prevent disease [3, 4]. Production of meat with a fatty acids profile more in line with public health recommendations has the potential to improve long-term human health without requiring substantial changes in consumer habits. It is well known that the fatty acid meat composition of pigs is largely dependent on genotype, physiological status and environmental factors such as nutrition [5–11].
The liver a highly specialized organ present in vertebrates and other animals regulates a wide variety of metabolic processes, which play a key role in the digestive function, the decomposition of red blood cells, hormone production and detoxification. Together with adipose tissue and skeletal muscle, the liver is crucial in regulating lipid metabolism. In pigs, the liver is the primary site of de novo cholesterol synthesis and fatty acid oxidation, whereas lipogenesis occurs essentially in liver and adipose tissues [12–16].
In the last few years, new high-throughput technologies have been developed for the massive analysis of genomic data. These methodologies yield new opportunities to explore the genetic variability of populations, as well as the characterization of the transcriptome architectures. Until the development of Next-generation sequencing (NGS) technologies, most mRNA expression studies have used microarray or quantitative PCR-based (qPCR) approaches. The development of RNA-Seq, a method based on NGS which consisting of the direct sequencing of RNA molecules present in a given sample, has provided a new tool for both transcriptome characterization and gene expression profiling. In RNA-Seq, the counts corresponding to each transcript can be used for quantification and these sequences can be mapped to the genome for their annotation. In comparison to microarrays, RNA-Seq provides a higher dynamic range, specificity and sensibility . In addition, it provides a picture of the transcriptome, allowing the characterization of alternative splicing, variation in the usage of promoters and polyadenilation sites, non-coding RNAs (ncRNA), single nucleotide variants (SNVs) and transposable elements. Furthermore, RNA-Seq data may allow the discovery of novel transcripts and long intergenic non-coding RNAs (lncRNAs) [18–20].
Recent studies in livestock species have employed RNA-Seq to explore the transcriptome of animal products, such as cow milk , bovine embryos , and tissue as pig gonads , liver, muscle, and abdominal fat , sheep bone  and bovine abomasal . However, most of the RNA-Seq studies in pigs have included analysis of only a few animals, and ignored within group intrinsic variability. For instance, two single animals of different breeds were compared by Esteve-Codina et al., (2011) and three tissues in two phenotypically extreme full-sib F2 females formed the basis of Chen et al., (2011) study.
The main goal of this study was the identification of differentially-expressed genes in the liver of groups of Iberian x Landrace crossbred pigs showing extreme phenotypes for intramuscular fatty acid composition using RNA-Seq. In addition, the porcine hepatic transcriptome was analyzed and transposable elements, new putative protein-coding genes and lncRNAs were identified.
Results and discussion
Phenotypic variation between extreme groups
Mean comparison ± standard deviation between high and low groups for the traits included in the principal components analysis (PCA)
Carcass height (CH)
72.79 ± 9.90
69.57 ± 12.73
Weight of ham (WH)
19.88 ± 2.46
19.66 ± 3.10
Weight of Shoulder (WS)
10.26 ± 1.25
10.93 ± 2.10
Intramuscular fat (IMF)
2.21 ± 0.88
1.49 ± 0.38
Myristic acid (C14:0)
1.22 ± 0.13
1.12 ± 0.12
Palmitic acid (C16:0)
23.78 ± 0.79
21.39 ± 0.69
Heptadecanoic acid (C17:0)
0.20 ± 0.03
0.33 ± 0.09
Stearic acid (C18:0)
14.65 ± 1.16
13.69 ± 0.77
Arachidic acid (C20:0)
0.26 ± 0.05
0.21 ± 0.08
Palmitoleic acid (C16:1 n-7)
2.74 ± 0.24
2.20 ± 0.33
Heptadecenoic acid (C17:1)
0.20 ± 0.05
0.32 ± 0.20
Oleic acid (C18:1 n-9)
42.57 ± 1.34
34.92 ± 2.96
Octadecenoic acid (C18:1 n-9)
4.04 ± 0.27
3.70 ± 0.31
Eicosenoic acid (C20:1 n-9)
0.86 ± 0.08
0.77 ± 0.07
Linoleic acid (C18:2 n-6)
7.16 ± 0.52
15.11 ± 1.65
α-Linolenic acid (C18:2 n-3)
0.46 ± 0.08
1.10 ± 0.56
Eicosadienoic acid (C20:2 n-6)
0.41 ± 0.05
0.63 ± 0.13
Eicosatrienoic acid (C20:3 n-6)
0.16 ± 0.03
0.53 ± 0.15
Arachidonic acid (C20:4 n-6)
0.84 ± 0.17
3.10 ± 0.77
Average Chain Length (ACL)
17.44 ± 0.02
17.49 ± 0.12
Saturated FA (SFA)
40.12 ± 1.51
36.72 ± 1.22
Monounsaturated FA (MUFA)
50.72 ± 1.56
42.38 ± 3.13
Polyunsaturated FA (PUFA)
9.03 ± 0.64
20.46 ± 2.21
Peroxidability index (PI)
12.28 ± 3.80
32.43 ± 4.43
Double-bond index (DBI)
0.65 ± 0.20
0.91 ± 0.03
Unsaturated index (UI)
1.63 ± 0.50
2.48 ± 0.10
Previous studies have reported that, in both backfat and LD muscle, Iberian pigs have higher percentages of palmitic acid, oleic acid, SFA and MUFA, and lower concentrations of LA and ALA acids than commercial breeds [7, 27, 28]. Moreover, Pascual et al., (2007)  reported that Landrace pigs have a higher content of LA and AA acids in their muscle than other commercial breeds. In general, fatter pigs show higher proportions of SFA and MUFA, but less PUFA than lean pigs [6, 29]. The genetic architecture of intramuscular FA composition in the Iberian x Landrace backcross was described in a genome-wide association study (GWAS), showing 43 chromosomal intervals associated with these traits . Since all animals were raised and fed under the same standard management conditions, differences between H and L groups are probably caused by the segregation within the analyzed animals of Iberian and Landrace alleles.
Phenotypic means between groups were compared and significant statistical differences in 73% of the analyzed traits was noted (19/26), mainly relating to intramuscular fatty acid composition (Table 1). The maximum differences between groups were observed for the profiles of essential PUFAs (AA, ETE, LA and ALA acids). Significant differences were also observed for PI and the percentage of palmitic, palmitoleic, heptadecenoic and heptadecanoic acids. From the 20 extreme animals, 10 females were selected for RNA sequencing (five per group). Pedigree information was used to select animals representing the parental genetic diversity. In addition, full-sibs within groups were avoided, animals within groups had different mothers, and four different fathers were selected per group. However, interesting familial relationships between animals of different groups were retained: there were two pairs of full-sibs and two pair of maternal half-sibs belonging to opposite groups. As before, the phenotypic means differed between groups. However, due to the reduced sample size, only sixteen traits showed significant differences (Additional file 1, Table S1).
Mapping and annotation
Number of single-end 100 bp reads obtained and percentages of mapped reads per animal
Proportion of reads mapping to exons, introns or within 1 Kb upstream or downstream of the annotated genes
% _5’ or _3’2
Number of transcripts assembled (TA) with Cufflinks and the percentage they represent in each sample
Gene expression analysis
The total amount of expressed genes in liver was similar between groups (L= 8797 – 10161, H= 8765 – 10083). Taking into account only those genes with a mean FPKM (normalized number of fragments per kilobase of exon per million reads) higher than zero, an aggregate of 10,485 expressed genes in L and 10,626 in H groups was observed. A total of 10,280 common genes were expressed in both groups. The correlation of mean gene expression levels between both groups (H vs L) was very high (r = 0.99), suggesting that the major fraction of the liver transcriptome is conserved between groups. Gene expression distribution reveals that less than 10% of these genes were expressed between 1 – 10 FPKM; around 42% between 10 FPKM - 100 FPKM; 38% among 100 – 1000 FPKM and, approximately, 8% more than 1000 FPKM (Additional file 2, Figure S1).
All 10 individuals were also assayed with the GeneChip® Porcine microarray (Affymetrix, Santa Clara, CA) which allows the expression analysis of 20,201 Sus scrofa genes. After probe normalization, correlation between the expression data of microarrays and RNA-Seq was calculated. In accordance with previous studies [17, 23], a strong Spearman correlation (r=0.72) was observed (Additional file 3, Figure S2). Results from both technologies were, in general, more similar for genes that showed intermediate expression values, whereas major differences were observed for low and high expressed genes in the Affymetrix microarray data. This pattern can be explained by the higher dynamic range of RNA-Seq [17, 33]. Finally, in line with the previous description of liver transcriptome , the top 100 expressed genes showed an overrepresentation in biological gene ontologies related to oxidoreductase activity, transport, proteolysis, translation, signal transduction, cholesterol homeostasis and lipid transport (p < 0.001).
Transposable elements analysis
The percentage of repetitive elements identified in the pig liver transcriptome was around5.8 – 7.3% (Additional file 4, Table S2), similar to that found in male gonad transcriptome (7.3%) . However, it should be noted that the total length of base pairs masked and the total number of transcripts were higher in male gonad transcriptome  than in our study. Two possible explanations may account for these differences: 1) in liver from 7.3 – 12.4 M of single-end reads per individual were obtained (Table 2), whereas in gonads  a total of 20 M paired-end reads were observed and, thus, a better fragment distribution and a higher number of transcripts were analysed; 2) the transcriptome complexity has been reported  to be higher in kidney, testes and brain tissues in comparison to liver and muscle.
Gene orthology and lncRNAs detection
Putative proteins identified in each H and L groups and orthologies detected against Homo sapiens , Bos taurus and Sus scrofa protein databases
Esteve-Codina et al. (2011)
For lncRNAs annotation, the previously reported sequences in pig male gonad transcriptome  were used as a reference database. A total of 186 (L) and 270 (H) of these putative lncRNAs was also expressed in pig liver. Within groups, 101 and 108 lncRNAs were expressed in all L and H animals, respectively, but only 89 lncRNAs were expressed in both groups (Additional file 5, Figure S3).
Differential gene expression analysis
Description of the differentially-expressed genes detected between High and Low groups with fold change ≥ 1.5 and p-value ≤0.005
1.1 x 10-16
3.0 x 10-13
2.9 x 10-10
2.2 x 10-11
1.4 x 10-8
1.3 x 10-8
6.3 x 10-6
5.2 x 10-7
2.0 x 10-4
1.0 x 10-6
3.0 x 10-4
1.4 x 10-6
4.0 x 10-4
7.8 x 10-6
2.0 x 10-3
1.3 x 10-5
3.0 x 10-3
1.6 x 10-5
3.0 x 10-3
1.8 x 10-5
3.0 x 10-3
3.2 x 10-5
1.0 x 10-2
3.5 x 10-5
1.0 x 10-2
3.9 x 10-5
1.0 x 10-2
5.1 x 10-5
1.0 x 10-2
1.4 x 10-4
2.0 x 10-2
1.4 x 10-4
2.0 x 10-2
2.2 x 10-4
2.0 x 10-2
3.0 x 10-4
3.0 x 10-2
3.2 x 10-4
3.0 x 10-2
3.3 x 10-4
3.0 x 10-2
3.7 x 10-4
3.0 x 10-2
5.6 x 10-4
5.0 x 10-2
6.0 x 10-4
5.0 x 10-2
6.5 x 10-4
5.0 x 10-2
6.7 x 10-4
5.0 x 10-2
6.7 x 10-4
5.0 x 10-2
7.9 x 10-4
5.0 x 10-2
8.4 x 10-4
6.0 x 10-2
1.0 x 10-3
7.0 x 10-2
1.0 x 10-3
7.0 x 10-2
1.1 x 10-3
7.0 x 10-2
1.2 x 10-3
7.0 x 10-2
1.2 x 10-3
7.0 x 10-2
1.3 x 10-3
7.0 x 10-2
1.3 x 10-3
7.0 x 10-2
1.3 x 10-3
7.0 x 10-2
1.5 x 10-3
8.0 x 10-2
1.9 x 10-3
8.0 x 10-2
1.9 x 10-3
8.0 x 10-2
1.9 x 10-3
8.0 x 10-2
2.3 x 10-3
9.0 x 10-2
2.7 x 10-3
9.0 x 10-2
2.7 x 10-3
9.0 x 10-2
2.8 x 10-3
9.0 x 10-2
2.8 x 10-3
9.0 x 10-2
3.2 x 10-3
1.3 x 10-1
3.4 x 10-3
1.4 x 10-1
3.5 x 10-3
1.4 x 10-1
3.9 x 10-3
1.6 x 10-1
4.0 x 10-3
1.6 x 10-1
4.3 x 10-3
1.6 x 10-1
4.6 x 10-3
1.7 x 10-1
4.7 x 10-3
1.7 x 10-1
4.8 x 10-3
1.7 x 10-1
5.0 x 10-3
1.7 x 10-1
5.0 x 10-3
1.7 x 10-1
5.0 x 10-3
1.7 x 10-1
In order to validate the expression data obtained by RNA-Seq, five genes (APOA2, LPIN1, ME3, CYP7A1 and CYP2C49) were selected among the differentially-expressed protein-coding genes to perform real time reverse transcription (RT-qPCR) assays. When the pattern of gene expression levels was compared, strong correlations ranking from 0.79 to 0.96 between RT-qPCR and RNA-Seq platforms were observed, confirming the high reproducibility of the data (Additional file 9, Table S3).
Interestingly, one of the studied genes, the CYP2C49 [ENSEMBL_Id: ENSSSCG00000010488] which belongs to the highly diverse superfamily of CYP450  and it is homologue to the human CYP2C9 gene, was located in a genomic region in which copy number variation (CNV) has been previously described in pigs . In order to assess whether observed differences of gene expression were influenced by differences in the CNV between animals, a real time quantitative PCR (qPCR) to determine the number of copies of the CYP2C49 gene was developed. For the first time, CNV affecting the CYP2C49 gene was described with relative quantification values ranging from 1 to 5.2 copies (Additional file 9, Table S3). However, no correlation between the number of copies and gene expression was observed. Therefore, further analysis will be necessary to elucidate the possible role of these structural variants in the fatty acid metabolism.
Differentially-expressed genes previously reported to be associated with the profile of intramuscular fatty acid composition in a genome-wide association study
Ensembl gene ID
C18:1(n-9), C18:2(n-6), MUFA
C16:1(n-7), ratio C16:1(n-7)/c16:0
C16:0, C18:2(n-6), rate MUFA/SFA
Functional clustering of differentially-expressed genes in the liver
From the 55 differentially-expressed protein-coding genes, 26 were up-regulated and 29 were down-regulated in H group in comparison to L (Table 6). To gain insight into the liver tissue processes that differed between groups, the list of the differentially-expressed genes was explored using the core analysis function included in Ingenuity Pathways Analysis (IPA). Initially, the pig gene IDs were converted to human genes but five protein-coding genes did not match with human homologs [Ensembl Ids: ENSSSCG00000007873, ENSSSCG00000003971, ENSSSCG00000014368, ENSSSCG00000013116, ENSSSCG00000001229], and therefore only 50 pig genes were eligible for network construction.
Description of the top seven molecular and cellular biological functions significantly modulated in the liver tissue when comparing H relative to L animals
ABCG8,ALOX15,AQP7,APOB,CYP2C9,THBS1, APOA2,ME3,NR4A3,LPIN1,CYP7A1,MTMR7, CYP2C19,CYP4A11
Small Molecule Biochemistry
ABCG8,ALOX15,AQP7,APOB,CYP2C9,THBS1, APOA2,ME3,GSTO1,FOS,NR4A3,LPIN1,SPTB, CYP7A1,MTMR7,CYP2C19,SDS,SLC11A1, CYP4A11
ABCG8,AQP7,ALOX15,APOB,THBS1,APOA2, GSTO1,FOS,NR4A3,LPIN1,CYP7A1,SLC11A1, FNDC1
Nucleic Acid Metabolism
Vitamin and Mineral Metabolism
Description of the top six canonical pathways significantly modulated in liver tissue when comparing H to L animals
Ingenuity canonical pathways
LPS/IL-1 Mediated Inhibition of RXR Function
Arachidonic Acid Metabolism
Fatty Acid Metabolism
Linoleic Acid Metabolism
In addition, the up-regulation of PPAR-α and RXR were coupled with the increased expression of lipin (LPIN1) and CYP7A1 genes. In mice, it has been reported that LPIN1 selectively activates a subset of coactivator 1α (PGC-1α) target pathways involved in fatty acid oxidation and mitochondrial oxidative phosphorylation, while suppressing the lipogenic program and lowering circulating lipid levels . Lipin activates mitochondrial fatty acid oxidative metabolism by inducing expression of the nuclear receptor PPAR-α, a known PGC-1α target, and via direct physical interactions with PPAR-α and PGC-1α. Furthermore, CYP7A1 has been shown to be a key factor of hepatic cholesterol homeostasis. All together these results suggest that H group may present greater uptake of fatty acids into hepatocytes (mainly LA and ALA acids). It is likely that the higher PUFA bioavailability in liver may affect expression of PPAR-α, RXR and their target genes, inducing a greater stimulation of both peroxisomal and mitochondrial β-oxidation, and leading to reduced triglyceride and cholesterol synthesis, and an enhanced elimination of cholesterol from the liver via bile acid formation. This intriguing possibility remains to be demonstrated, although there is evidence that FA, in particular unsaturated FA, exert many of their biological effects by regulating the activity of numerous transcription factors in liver, such as PPAR-α . Recently,  has demonstrated that FA oxidation is regulated by hepatic MUFA to PUFA ratio through the activation of PPAR-α. In agreement to our results, hepatic expression of PPAR-α was higher in pigs fed with a higher level of PUFA. This is also in line with the lower IMF content in H group than in L animals, and the lower proportion of de novo fatty acids in the IMF. Therefore, these transcriptome changes may reflect counter mechanisms of liver tissue to respond or compensate for changes in IMF fatty acid profile, which depends on possible metabolisation of FAs and the possibility of being synthesized by the pig adipose tissue . However, the question remains how different types of FA control the expression of genes and a direct examination of the effect of each individual FA on porcine muscle fatty acid composition is needed.
In the present study, RNA-Seq was used for the analysis of the pig liver transcriptome in animals of extreme phenotypes for intramuscular fatty acid composition.
The liver plays an important role in lipid metabolism and, thus, the analysis of liver transcriptome in extreme pigs for intramuscular fatty acid composition may be relevant to elucidate its functional complexity. Although the main goal of this study was to find differentially-expressed genes between phenotypically extreme animals, the use of RNA-Seq allowed the identification of transposable elements, lncRNAs and new protein-coding genes in the porcine liver transcriptome.
The first principal component of PCA analysis classified animals in two extreme groups for the fatty acid composition of LD muscle. Group H of animals had a higher content of PUFA, including essential FA such as LA, ALA, ELE and AA acids than group L animals. Conversely, the latter had a higher content of SFA and MUFA, palmitoleic and oleic acids.
The lipid content and fatty acid profile of muscle plays an important role in the tenderness, flavour and juiciness of cooked meat . In swine production, the reduction of intramuscular fat (IMF) in some breeds due to a preference selection for lean pigs, has affected meat quality. From this point of view, PUFA has a negative effect on the oxidative stability of muscle, which, in turn, affects flavour and muscle colour . On the other hand, desirable sensorial characteristics tend to be associated with MUFA and SFA [6, 48, 49]. Lipid and fatty acid compositions of food have an important impact on human health, with a high consumption of SFA associated with obesity, high plasma cholesterol and cardiovascular disease [50, 51]. Conversely, PUFAs, mainly ω−3, have been considered beneficial for human health, by reducing serum low-density lipoprotein-C, total cholesterol concentration and modulating immune functions and inflammatory processes [52–54].
There is increasing awareness of the wide range of health benefits of PUFA in general, and of ω−3 fatty acids in particular. Meat is an important basis of human nutrition, and pork meat is seen to be a major source of human food. The composition of fatty acids stored in adipose tissue in pigs largely reflects that of ingested lipids . Thus, swine meat enriched with ω−3 fatty acids can be achieved by feeding with commercial diets supplemented with this PUFA , and possibly by selective breeding. In fact, there is a genetic basis of PUFA level in pork meat. It is likely that H group presented higher absorption of essential PUFA, increasing their amount reaching the IMF tissue, which in turn could be considered as an important factor in the inhibition of the de novo saturated fatty acid proportion in meat. Furthermore, differences on elongation, desaturation and oxidation of those essential PUFA to longer-chain ω−3 and ω−6 fatty acids cannot be discarded. Therefore, from the human health perspective, increasing H genotypes through breeding programs could be desirable because meat and meat-derived foods are still large contributors to saturated fatty acids intake in humans. These observations, together with several gene expression effects are the major factors leading us to believe that genetic was indeed a significant factor affecting meat IMF PUFA content and composition. However, an inverse relationship exists between nutritional value and eating quality of meat and, as consequence, established selection criteria to all together improve meat quality from the sensorial and nutritional point of view is a complex matter. Therefore, a holistic approach including both nutrigenetic and nutrigenomic disciplines may be required to improve the pork meat quality from both points of view.
We used RNA-Seq as a tool to explore the liver transcriptome of ten female pigs with extreme phenotypes for intramuscular fatty acid composition. Transposable elements, lncRNAs and new putative protein-coding genes were identified. Reproducibility of the data was confirmed by the strong correlation observed between the values of gene expression obtained by RNA-Seq, RT-qPCR and microarrays. A total of 55 genes differentially-expressed between extreme animals were identified. These genes belong to canonical pathways and gene networks related to the lipid and fatty acid metabolism. In concordance with the initial phenotypic classification, pathway analysis inferred that linolenic and arachidonic acid metabolism was altered between extreme animals. The results obtained may help in the design of new selection strategies to improve pork meat quality from both the sensorial and nutritional points of view.
Animal material and phenotypes
The population studied was originated by crossing three Iberian boars (Guadyerbas line) with 31 Landrace sows [55, 56]. Five F1 males were backcrossed with 26 Landrace sows and 144 BC1_LD pigs were obtained. All pigs were raised in a normal intensive system, fed under standard management conditions and were slaughtered at an average age of 179.8 ± 2.6 days following national and institutional guidelines for the ethical use and treatment of animals in experiments.
A total of 48 traits related with growth, carcass quality and intramuscular fatty acid composition were measured. A PCA was performed with the prcomp procedure of R software , including phenotypic information from twenty-six of the total traits. Four of these were related to carcass quality (carcass height, weight of ham, weight of shoulder and intramuscular fat) whereas the rest corresponded to fatty acids composition in muscle and indices of fatty acids metabolism. Animals with extreme phenotypes, according to the first principal component, were selected to generate the High (H) and Low (L) groups with 20 animals per group (Figure 1). Phenotypic mean comparisons between groups were performed using R. Since sex differences in liver transcriptome have been reported in several species, selection was made considering pedigree information representing the parental genetic diversity and only females were retained for RNA sequencing (five per group).
RNA isolation, library preparation and sequencing
From the 10 selected animals, total RNA was isolated from liver using the RiboPureTM Isolation of High Quality Total RNA (Ambion®, Austin, TX) following the manufacturer’s recommendations. RNA was quantified using the NanoDrop ND-1000 spectrophotometer (NanoDrop products, Wilmington, USA) and checked for purity and integrity in a Bioanalyzer-2100 (Agilent Technologies, Inc., Santa Clara CA, USA).
Sequencing libraries were generated using Illumina mRNA-Seq following manufacturer’s instructions and ten index codes were added to attribute sequences to each animal. A total of two channels of an Illumina Hi-Seq 2000 instrument (Fasteris SA, Plan-les-Ouates, Switzerland) were used to sequence two pools of five samples (one pool with five samples per lane with barcoding).
Mapping, assembling and annotation of reads
After removal of sequencing adaptors and low-complexity readsTopHatv1.2.0 software  was employed to map reads using as reference the version 9.61 of pig genome (Sscrofa 9.61) http://www.ensembl.org/info/data/ftp/index.html. Quality control and reads statistics were determined with FASTQC http://www.bioinformatics.bbsrc.ac.uk/projects/fastqc/. Transcript assembly was performed using Cufflinks v0.9.3  , with a minimum alignment count per locus of 10. Finally, S-MART http://urgi.versailles.inra.fr/Tools/S-MART for read annotation was used.
Gene expression quantification and correlation analysis with expression microarrays
Gene expression quantification was performed using the normalized number of fragments per kilobase of exon per million reads (FPKM) as reported in Cufflinks output . Correlations between mean expression values between groups were calculated. All individuals were also assayed with high-density oligonucleotide microarray chips (GeneChip® Porcine) from Affymetrix (Santa Clara, CA) containing a total of 23,937 probe sets (23,256 transcripts), representing 20,201 Sus scrofa genes. Microrarrays were hybridized and scanned at the Institut de Recerca Hospital Universitari Vall d’Hebron (Barcelona, Spain) following Affymetrix standard protocols. Expression data were generated with Gene-Chip Operating Software (GCOS). Probes were adjusted for background noises and normalized using the GCRMA R package . The average probe value per gene was calculated and a total of 6,025 Ensembl gene IDs could be retrieved to estimate the Spearman correlation between the log2 expression values of genes analysed by RNA-Seq and microarrays. Finally, a GO enrichment analysis with the QuickGO browser http://www.ebi.ac.uk/QuickGO/ was performed for the top 100 most expressed genes.
Differential gene expression analysis
Differential expression analysis (DE) between groups was performed using DESeq . This R package uses as input file the unambiguous table of counts per gene obtained from HTseq-count http://www-huber.embl.de/users/anders/HTSeq/doc/overview.html. DESeq models the data using negative binomial distributions assuming that the mean is a good predictor of variance. Therefore, it assumes that genes with similar expression level also have similar variance across replicates . Following the DESeq author’s recommendations, some exploratory diagnostic plots were executed to check the dispersion estimate and data quality. In order to ascertain the base variance the function 'varianceFitDiagnostics' was used and the per-gene estimates of the base variance was plotted against the base levels. The uniformity of the cumulative probabilities estimated by the 'varianceFitDiagnostics' was also verified via the 'residualsEcdfPlot' function.
Differentially-expressed genes were detected through the ‘nbinomTest’ function of DESeq. All the genes with a fold change between H and L groups higher than or equal to 1.5 fold were retained (total of 2051 genes). Then, for this subset of genes, the R package q-value  was employed to calculate the false-discovery rate and genes with a p-value ≤ 0.005 (which is equivalent to a q-value ≤ 0.17) were retained.
Validation of differentially-expressed genes by RT-qPCR and copy number determination by qPCR
In order to evaluate the repeatability and reproducibility of gene expression data obtained by RNA-Seq, a RT-qPCR assay using SYBR Green chemistry (Fast Start Universal Sybr green master, Rox; Roche Applied Science, Mannheim, Germany) and the comparative Ct method  was performed in an ABI PRISM® 7900 HT (Applied Biosystems, Inc., FosterCity, CA).
The isolated RNA of individual samples was reverse-transcribed into cDNA using the High Capacity Reverse cDNA transcription Kit (Applied Biosystems) in a total volume of 20 μl containing 1 μg of total RNA, following the manufacturer’s instructions. PCR primers were designed using Primer Express™ software (Applied Biosystems) and are shown in Additional file 12, Table S4. Two genes: β-2 microglobulin (β2m) and hypoxanthine phosphoribosyltransferase 1 (HPRT1), previously validated as stable expressed control genes in liver tissue by geNorm  were used as endogenous controls (Corominas et al., unpublished data). Due to the comparative Ct method requiring the target and endogenous PCR efficiencies to be nearly equal, validation experiments for each gene were performed. Thus, the log cDNA dilution (1:2, 1:20, 1:200, 1:2,000) versus ΔCt, was plotted to obtain absolute slopes < 0.1 in all cases that allowed the use of the 2-ΔΔCt method. PCR amplifications were performed in a total volume of 20 μl containing 5 μl of cDNA sample diluted 1:125. Depending on the pair primers, various concentrations were utilized (see Additional file 12, Table S4). Each sample was analyzed by triplicate, thermal cycle was: 10 min at 95°C and 40 cycles of 15 sec at 95°C and 1 min at 60°C. A dissociation curve was drawn for each primer pair. Data was analyzed using the SDS v2.4 and DataAssistTM v3.0 software (Applied Biosystems). The sample of lowest expression level was selected as calibrator. Correlation between RNA-seq (Htseq) and RT-qPCR data (2-ΔΔCt) was carried out with R.
Copy number variation was quantified using the assay described above with some modifications. PCR amplification was carried out with 10 ng of genomic DNA isolated from diaphragm samples by the phenol-chloroform method . Primers used to amplify CYP2C49 gene are described in Additional file 12, Table S4. For single copy endogenous control gene amplification, a previously described design on the glucagon (GCG) gene  was used, but a single nucleotide substitution on primer forward was introduced to adapt the primer to the porcine species (Additional file 12, Table S4). A sample with the lowest copy number was selected as a calibrator.
Transposable element analysis
To identify repetitive and transposable elements in pig liver transcriptome RepeatMasker version open-3.3.0 [http://www.repeatmasker.org/] and the ‘quick search’ and ‘pig’ species options with Search Engine: NCBI/RMBLAST and complete Database: 20090604 were used.
Orthology and lncRNA detection
Intergenic expressed regions not yet annotated in the Sscrofa 9.61 pig genome assembly as described in  were extracted. Then, a conservative approach was followed, using only sequences expressed in at least four of the five animals of each group (H and L). To identify which of these transcripts were putative coding transcripts the Augustus software was used , providing exon boundaries and allowing complete protein translation from the forward strand. Finally, BLASTP was employed to check which of these predicted proteins were already annotated in the Homo sapiens, Bos taurus and Sus scrofa protein databases. For lncRNA annotation, the intergenic expressed regions were compared by BLAST with the 2,047 putative porcine lncRNA reported by Esteve-Codina et al. (2011). All transcripts that matched with an expectation value lower than 10-5 were retained.
Gene functional classification, network and canonical pathways analyses
Biological network generation, functional classification and pathways analyses of differentially-expressed genes were carried out using Ingenuity Pathways Analysis software (IPA; Ingenuity Systems, http://www.ingenuity.com). The list of human homologs that correspond to the 50 protein-coding pig genes differentially-expressed was uploaded into the application. Then, each gene identifier was mapped to its corresponding gene object in the Ingenuity Pathways Knowledge Base (IPKB). Networks of these genes were generated based on their connectivity. Network analysis returns a score that ranks networks according to their degree of relevance to the network eligible molecules in the dataset . The network score is based on the hypergeometric distribution and is calculated with the right-tailed Fisher’s exact test. The score is the negative log of this p- value. Only those molecules that demonstrate direct and indirect relationships to other genes, or proteins were integrated into the analysis.
IPA Functional Analysis was employed to identify the most significant biological functions in the comparative dataset of H and L groups. A canonical pathways analysis was generated to identify the pathways from the IPA library that were most significant. Fischer’s exact test was employed to calculate a p-value which determines the probability that each biological functions and/or canonical pathway is due to chance alone. The cut-off for considering a significance association was established by Benjamini & Hochberg (B-H) multiple testing correction of the p-value (FDR < 0.05) .
The full data sets have been submitted to Gene Expression Omnibus (GEO) under accession GSE38588 and at NCBI Sequence Read Archive (SRA) under Accession SRA053452, Bioproject: PRJNA168072.
Ingenuity Pathways Analysis
Retinoid X receptors
Pregnane X receptor
Farnesoid X receptor
Peroxisome proliferator-activated receptors alpha.
This work was funded by MICINN projects AGL2008-04818-C03/GAN and AGL2011-29821-C02 (Ministerio de Economía y Competitividad), and by the Innovation Consolider-Ingenio 2010 Program (CSD2007-00036, Centre for Research in Agrigenomics). Y. Ramayo-Caldas was funded by a FPU PhD grant from the Spanish Ministerio de Educación (AP2008-01450), J. Corominas was funded by a FPI PhD grant from the Spanish Ministerio de Educación (BES-2009-018223), A. Esteve-Codina is recipient of a FPI PhD fellowship from the Ministerio de Educación (BES-2008-005772), Spain. This manuscript has been proofread by Mrs. Valma Ruth Dunstan, BA, BEdSt, CELTA, MA (TESOL), a native English-speaker and instructor of English at the University of Queensland (Brisbane, Australia).
- Cooper D: A brief history of cross-species organ transplantation. Proc (Bayl Univ Med Cent). 2012, 25 (1): 49-57.Google Scholar
- Dodson M, Hausman G, Guan L, Du M, Rasmussen T, Poulos S, Mir P, Bergen W, Fernyhough M, McFarland D, Rhoads R, Soret B, Reecy J, Velleman S, Jiang Z: Lipid metabolism, adipocyte depot physiology and utilization of meat animals as experimental models for metabolic research. Int J Biol Sci. 2010, 6 (7): 691-699.PubMed CentralView ArticlePubMedGoogle Scholar
- Bauman DE, Mather IH, Wall RJ, Lock AL: Major Advances Associated with the Biosynthesis of Milk. J Dairy Sci. 2006, 89 (4): 1235-1243. 10.3168/jds.S0022-0302(06)72192-0.View ArticlePubMedGoogle Scholar
- FAO: Fats and Fatty Acids in Human Nutrition.Report of Expert Consultation. 2008, 308-ISBN: 978-3-8055-9261-1Google Scholar
- Wood JD, Richardson RI, Nute GR, Fisher AV, Campo MM, Kasapidou E, Sheard PR, Enser M: Effects of fatty acids on meat quality: a review. Meat Sci. 2004, 66 (1): 21-32. 10.1016/S0309-1740(03)00022-6.View ArticlePubMedGoogle Scholar
- Wood JD, Enser M, Fisher AV, Nute GR, Sheard PR, Richardson RI, Hughes SI, Whittington FM: Fat deposition, fatty acid composition and meat quality: A review. Meat Sci. 2008, 78 (4): 343-358. 10.1016/j.meatsci.2007.07.019.View ArticlePubMedGoogle Scholar
- Estévez M, Morcuende D, Cava López R: Physico-chemical characteristics of M. Longissimus dorsi from three lines of free-range reared Iberian pigs slaughtered at 90 kg live-weight and commercial pigs: a comparative study. Meat Sci. 2003, 64 (4): 499-506. 10.1016/S0309-1740(02)00228-0.View ArticlePubMedGoogle Scholar
- Guillevic M, Kouba M, Mourot J: Effect of a linseed diet or a sunflower diet on performances, fatty acid composition, lipogenic enzyme activities and stearoyl-CoA-desaturase activity in the pig. Livest Sci. 2009, 124 (1‚Äì3): 288-294.View ArticleGoogle Scholar
- Pascual J, Rafecas M, Canela M, Boatella J, Bou R, Barroeta A, Codony R: Effect of increasing amounts of a linoleic-rich dietary fat on the fat composition of four pig breeds, Part II: Fatty acid composition in muscle and fat tissues. Food Chem. 2007, 100 (4): 1639-1648. 10.1016/j.foodchem.2005.12.045.View ArticleGoogle Scholar
- Casellas J, Noguera JL, Reixach J, Diaz I, Amills M, Quintanilla R: Bayes factor analyses of heritability for serum and muscle lipid traits in Duroc pigs. J Anim Sci. 2010, 88 (7): 2246-2254. 10.2527/jas.2009-2205.View ArticlePubMedGoogle Scholar
- Ntawubizi M, Colman E, Janssens S, Raes K, Buys N, De Smet S: Genetic parameters for intramuscular fatty acid composition and metabolism in pigs. J Anim Sci. 2010, 88: 1286-1294. 10.2527/jas.2009-2355.View ArticlePubMedGoogle Scholar
- Gondret F, Ferre P, Dugail I: ADD-1/SREBP-1 is a major determinant of tissue differential lipogenic capacity in mammalian and avian species. J Lipid Res. 2001, 42 (1): 106-113.PubMedGoogle Scholar
- Nafikov R, Beitz D: Carbohydrate and Lipid Metabolism in Farm Animals. J Nutr. 2007, 137 (3): 702-705.PubMedGoogle Scholar
- O'Hea E, Leveille G: Significance of adipose tissue and liver as sites of fatty acid synthesis in the pig and the efficiency of utilization of various substrates for lipogenesis. J Nutr. 1969, 99 (3): 338-344.PubMedGoogle Scholar
- Mourot J, Kouba M, Peiniau P: Comparative study of in vitro lipogenesis in various adipose tissues in the growing domestic pig (Sus domesticus). Comp Biochem Physiol B Biochem Mol Biol. 1995, 111 (3): 379-384. 10.1016/0305-0491(95)00005-S.View ArticlePubMedGoogle Scholar
- Nguyen P, Leray V, Diez M, Serisier S, Le Bloc'h J, Siliart B, Dumon H: Liver lipid metabolism. J Anim Physiol Anim Nutr (Berl). 2008, 92 (3): 272-283. 10.1111/j.1439-0396.2007.00752.x.View ArticleGoogle Scholar
- Marioni JC, Mason CE, Mane SM, Stephens M, Gilad Y: RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays. Genome Res. 2008, 18 (9): 1509-1517. 10.1101/gr.079558.108.PubMed CentralView ArticlePubMedGoogle Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008, 5 (7): 621-628. 10.1038/nmeth.1226.View ArticlePubMedGoogle Scholar
- Garber M, Grabherr MG, Guttman M, Trapnell C: Computational methods for transcriptome annotation and quantification using RNA-seq. Nat Meth. 2011, 8 (6): 469-477. 10.1038/nmeth.1613.View ArticleGoogle Scholar
- Davey JW, Hohenlohe PA, Etter PD, Boone JQ, Catchen JM, Blaxter ML: Genome-wide genetic marker discovery and genotyping using next-generation sequencing. Nat Rev Genet. 2011, 12 (7): 499-510. 10.1038/nrg3012.View ArticlePubMedGoogle Scholar
- Cánovas A, Rincon G, Islas-Trejo A, Wickramasinghe S, Medrano J: SNP discovery in the bovine milk transcriptome using RNA-Seq technology. Mamm Genome. 2010, 21 (11): 592-598. 10.1007/s00335-010-9297-z.PubMed CentralView ArticlePubMedGoogle Scholar
- Huang W, Khatib H: Comparison of transcriptomic landscapes of bovine embryos using RNA-Seq. BMC Genomics. 2011, 11 (1): 711-View ArticleGoogle Scholar
- Esteve-Codina A, Kofler R, Palmieri N, Bussotti G, Notredame C, Perez-Enciso M: Exploring the gonad transcriptome of two extreme male pigs with RNA-seq. BMC Genomics. 2011, 12 (1): 552-10.1186/1471-2164-12-552.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen C, Ai H, Ren J, Li W, Li P, Qiao R, Ouyang J, Yang M, Ma J, Huang L: A global view of porcine transcriptome in three tissues from a full-sib pair with extreme phenotypes in growth and fat deposition by paired-end RNA sequencing. BMC Genomics. 2011, 12 (1): 448-10.1186/1471-2164-12-448.PubMed CentralView ArticlePubMedGoogle Scholar
- Jager M, Ott C-E, Grunhagen J, Hecht J, Schell H, Mundlos S, Duda G, Robinson P, Lienau J: Composite transcriptome assembly of RNA-seq data in a sheep model for delayed bone healing. BMC Genomics. 2011, 12 (1): 158-10.1186/1471-2164-12-158.PubMed CentralView ArticlePubMedGoogle Scholar
- Li R, Rinaldi M, Capuco A: Characterization of the abomasal transcriptome for mechanisms of resistance to gastrointestinal nematodes in cattle. Vet Res. 2011, 42 (1): 114-10.1186/1297-9716-42-114.PubMed CentralView ArticlePubMedGoogle Scholar
- Ramirez R, Cava R: Carcass composition and meat quality of three different Iberian × Duroc genotype pigs. Meat Sci. 2007, 75 (3): 388-396. 10.1016/j.meatsci.2006.08.003.View ArticlePubMedGoogle Scholar
- Serra X, Gil F, Pérez-Enciso M, Oliver MA, Vázquez JM, Gispert M, Díaz I, Moreno F, Latorre R, Noguera JL: A comparison of carcass, meat quality and histochemical characteristics of Iberian (Guadyerbas line) and Landrace pigs. Livestock Production Science. 1998, 56 (3): 215-223. 10.1016/S0301-6226(98)00151-1.View ArticleGoogle Scholar
- Webb EC, O‚AôNeill HA: The animal fat paradox and meat quality. Meat Sci. 2008, 80 (1): 28-36. 10.1016/j.meatsci.2008.05.029.View ArticlePubMedGoogle Scholar
- Ramayo-Caldas Y, Mercadé A, Castelló A, Yang B, Rodríguez C, Alves E, Díaz I, Ibáñez-Escriche N, Noguera JL, Pérez-Enciso M, Fernández AI, Folch JM: Genome-Wide Association Study for Intramuscular Fatty Acid Composition in an Iberian x Landrace Cross. J Animal Sci. 2012, 90: 1-11. 102527/jas2011-4900.View ArticleGoogle Scholar
- Trapnell C, Pachter L, Salzberg SL: TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009, 25 (9): 1105-1111. 10.1093/bioinformatics/btp120.PubMed CentralView ArticlePubMedGoogle Scholar
- Zytnicki M, Quesneville H: S-MART, A Software Toolbox to Aid RNA-seq Data Analysis. PLoS One. 2011, 6 (10): e25988-10.1371/journal.pone.0025988.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang Z, Gerstein M, Snyder M: RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009, 10 (1): 57-63. 10.1038/nrg2484.PubMed CentralView ArticlePubMedGoogle Scholar
- Yu Y, Ping J, Chen H, Jiao L, Zheng S, Han Z-G, Hao P, Huang J: A comparative analysis of liver transcriptome suggests divergent liver function among human, mouse and rat. Genomics. 2010, 96 (5): 281-289. 10.1016/j.ygeno.2010.08.003.View ArticlePubMedGoogle Scholar
- Ramskold D, Wang ET, Burge CB, Sandberg R: An Abundance of Ubiquitously Expressed Genes Revealed by Tissue Transcriptome Sequence Data. PLoS Comput Biol. 2009, 5 (12): e1000598-10.1371/journal.pcbi.1000598.PubMed CentralView ArticlePubMedGoogle Scholar
- Stanke M, Diekhans M, Baertsch R, Haussler D: Using native and syntenically mapped cDNA alignments to improve de novo gene finding. Bioinformatics. 2008, 24 (5): 637-644. 10.1093/bioinformatics/btn013.View ArticlePubMedGoogle Scholar
- Anders S, Huber W: Differential expression analysis for sequence count data. Genome Biol. 2010, 11 (10): R106-10.1186/gb-2010-11-10-r106.PubMed CentralView ArticlePubMedGoogle Scholar
- Kirischian N, McArthur A, Jesuthasan C, Krattenmacher B, Wilson J: Phylogenetic and Functional Analysis of the Vertebrate Cytochrome P450 2 Family. J Mol Evol. 2011, 72 (1): 56-71. 10.1007/s00239-010-9402-7.View ArticlePubMedGoogle Scholar
- Ramayo-Caldas Y, Castello A, Pena R, Alves E, Mercade A, Souza C, Fernandez A, Perez-Enciso M, Folch J: Copy number variation in the porcine genome inferred from a 60 k SNP BeadChip. BMC Genomics. 2010, 11 (1): 593-10.1186/1471-2164-11-593.PubMed CentralView ArticlePubMedGoogle Scholar
- Wise EM, Ball EG: Malic Enzyme and Lipogenesis. Proc Natl Acad Sci U S A. 1964, 52: 1255-1263. 10.1073/pnas.52.5.1255.PubMed CentralView ArticlePubMedGoogle Scholar
- Cherkaoui-Malki M, Meyer K, Cao WQ, Latruffe N, Yeldandi AV, Rao MS, Bradfield CA, Reddy JK: Identification of novel peroxisome proliferator-activated receptor alpha (PPARalpha) target genes in mouse liver using cDNA microarray analysis. Gene Expr. 2001, 9 (6): 291-304.View ArticlePubMedGoogle Scholar
- Tugwood JD, Holden PR, James NH, Prince RA, Roberts RA: A peroxisome proliferator-activated receptor-alpha (PPARalpha) cDNA cloned from guinea-pig liver encodes a protein with similar properties to the mouse PPARalpha: implications for species differences in responses to peroxisome proliferators. Arch Toxicol. 1998, 72 (3): 169-177. 10.1007/s002040050483.View ArticlePubMedGoogle Scholar
- Cheon Y, Nara TY, Band MR, Beever JE, Wallig MA, Nakamura MT: Induction of overlapping genes by fasting and a peroxisome proliferator in pigs: evidence of functional PPARŒ± in nonproliferating species. Am J Physiol Regul Integr Comp Physiol. 2005, 288 (6): R1525-R1535. 10.1152/ajpregu.00751.2004.View ArticlePubMedGoogle Scholar
- Finck BN, Gropler MC, Chen Z, Leone TC, Croce MA, Harris TE, Lawrence JC, Kelly DP: Lipin 1 is an inducible amplifier of the hepatic PGC-1Œ±/PPARŒ± regulatory pathway. Cell Metab. 2006, 4 (3): 199-210. 10.1016/j.cmet.2006.08.005.View ArticlePubMedGoogle Scholar
- Barnes KM, Winslow NR, Shelton AG, Hlusko KC, Azain MJ: Effect of dietary conjugated linoleic acid on marbling and intramuscular adipocytes in pork. J Anim Sci. 2011, 90 (4): 1142-1149.View ArticlePubMedGoogle Scholar
- Theil PK, Lauridsen C: Interactions between dietary fatty acids and hepatic gene expression in livers of pigs during the weaning period. Livest Sci. 2007, 108 (1‚Äì3): 26-29.View ArticleGoogle Scholar
- Pascual J, Rafecas M, Canela M, Boatella J, Bou R, Baucells M, Codony R: Effect of increasing amounts of a linoleic-rich dietary fat on the fat composition of four pig breeds, Part I: Backfat fatty acid evolution. Food Chem. 2006, 96 (4): 538-548. 10.1016/j.foodchem.2005.02.042.View ArticleGoogle Scholar
- Chizzolini R, Novelli E, Zanardi E: Oxidation in traditional mediterranean meat products. Meat Sci. 1998, 49 (Supplement 1(0)): S87-S99.View ArticleGoogle Scholar
- Carrapiso AI, Bonilla F: Garcƒ±ÃÅa C: Effect of crossbreeding and rearing system on sensory characteristics of Iberian ham. Meat Sci. 2003, 65 (1): 623-629. 10.1016/S0309-1740(02)00256-5.View ArticlePubMedGoogle Scholar
- Katan M, Zock P, Mensink R: Effects of fats and fatty acids on blood lipids in humans: an overview. Am J Clin Nutr. 1994, 60 (6): 1017S-1022S.PubMedGoogle Scholar
- Chizzolini R, Zanardi E, Dorigoni V, Ghidini S: Calorific value and cholesterol content of normal and low-fat meat and meat products. Trends in Food Science & Technology. 1999, 10 (4–5)): 119-128.View ArticleGoogle Scholar
- Rudel LL, Parks JS, Sawyer JK: Compared With Dietary Monounsaturated and Saturated Fat, Polyunsaturated Fat Protects African Green Monkeys From Coronary Artery Atherosclerosis. Arterioscler Thromb Vasc Biol. 1995, 15 (12): 2101-2110. 10.1161/01.ATV.15.12.2101.View ArticlePubMedGoogle 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 (4): 372-387. 10.1016/j.plipres.2011.06.003.View ArticlePubMedGoogle Scholar
- Stables MJ, Gilroy DW: Old and new generation lipid mediators in acute inflammation and resolution. Prog Lipid Res. 2011, 50 (1): 35-51. 10.1016/j.plipres.2010.07.005.View ArticlePubMedGoogle Scholar
- Clop A, Ovilo C, Perez-Enciso M, Cercos A, Tomas A, Fernandez A, Coll A, Folch J, Barragan C, Diaz I, Oliver M, Varona L, Silio L, Sanchez A, Noguera J: Detection of QTL affecting fatty acid composition in the pig. Mamm Genome. 2003, 14 (9): 650-656. 10.1007/s00335-002-2210-7.View ArticlePubMedGoogle Scholar
- Perez-Enciso M, Clop A, Noguera JL, Ovilo C, Coll A, Folch JM, Babot D, Estany J, Oliver MA, Diaz I, Sanchez A: A QTL on pig chromosome 4 affects fatty acid metabolism: evidence from an Iberian by Landrace intercross. J Anim Sci. 2000, 78 (10): 2525-2531.PubMedGoogle Scholar
- Ihaka R, Gentleman R: R: A Language for Data Analysis and Graphics. Journal of Computational and Graphical Statistics. 1996, 5 (3): 299-314.Google Scholar
- Zhang Y, Klein K, Sugathan A, Nassery N, Dombkowski A, Zanger UM, Waxman DJ: Transcriptional Profiling of Human Liver Identifies Sex-Biased Genes Associated with Polygenic Dyslipidemia and Coronary Artery Disease. PLoS One. 2011, 6 (8): e23506-10.1371/journal.pone.0023506.PubMed CentralView ArticlePubMedGoogle Scholar
- Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L: Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotech. 2010, 28 (5): 511-515. 10.1038/nbt.1621.View ArticleGoogle Scholar
- Zhijin W, Irizarry RA, Robert G, Francisco Martinez M, Spencer F: A Model Based Background Adjustment for Oligonucleotide Expression Arrays. 2004, Dept. of Biostatistics, Johns Hopkins University, Available in: http://www.bepress.com/jhubiostat/paper1Google Scholar
- Storey JD, Tibshirani R: Statistical significance for genomewide studies. Proc Natl Acad Sci U S A. 2003, 100 (16): 9440-9445. 10.1073/pnas.1530509100.PubMed CentralView ArticlePubMedGoogle Scholar
- Schmittgen TD, Livak KJ: Analyzing real-time PCR data by the comparative CT method. Nat Protocols. 2008, 3 (6): 1101-1108. 10.1038/nprot.2008.73.View ArticlePubMedGoogle Scholar
- Overgard A-C, Nerland A, Patel S: Evaluation of potential reference genes for real time RT-PCR studies in Atlantic halibut (Hippoglossus Hippoglossus L.); during development, in tissues of healthy and NNV-injected fish, and in anterior kidney leucocytes. BMC Mol Biol. 2010, 11 (1): 36-10.1186/1471-2199-11-36.PubMed CentralView ArticlePubMedGoogle Scholar
- Ausubel F, Brent R, Kingston R, Moore D, Seidman J, Smith J, Struhl K: Current Protocols in Molecular Biology. 1987Google Scholar
- Ballester M, Castelló A, Ibáñez E, Sánchez A, Folch J: Real-time quantitative PCR-based system for determining transgene copy number in transgenic animals. Biotechniques. 2004, 37 (4): 610-613.PubMedGoogle Scholar
- Calvano SE, Xiao W, Richards DR, Felciano RM, Baker HV, Cho RJ, Chen RO, Brownstein BH, Cobb JP, Tschoeke SK, Miller-Graziano C, Moldawer LL, Mindrinos MN, Davis RW, Tompkins RG, Lowry SF: Large Scale Collab. Res P, Host Response to I: A network-based analysis of systemic inflammation in humans. Nature. 2005, 437 (7061): 1032-1037. 10.1038/nature03985.View ArticlePubMedGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society Series B (Methodological). 1995, 57 (1): 289-300.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.