Josep María Folch, Anna Esteve-Codina, Jordi Estellé, Noelia Ibáñez-Escriche, Miguel Pérez-Enciso, Maria Ballester, Yuliaxis Ramayo-Caldas, Jordi Corominas, Anna Castelló, Núria Mach, Ana Isabel Fernández, Consorci CSIC-IRTA-UAB-UB, Center for Research in Agricultural Genomics, Génétique Animale et Biologie Intégrative (GABI), AgroParisTech-Institut National de la Recherche Agronomique (INRA), Genètica i Millora Animal, Instituto de Investigación y Tecnología Agroalimentarias, Departamento de Mejora Genética Animal, Instituto Nacional de Investigaciones Agronomicas, Departament de Ciència Animal i dels Aliments, Facultat de Veterinària, Universitat Autònoma de Barcelona [Barcelona] (UAB), Universitat Autònoma de Barcelona (UAB), and Institut National de la Recherche Agronomique (INRA)-AgroParisTech
Background 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. Results 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. Conclusions 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., 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.