Dissecting the Genetic Basis of Pig Intramuscular Fatty Acid Composition
Total Page:16
File Type:pdf, Size:1020Kb
Facultat de Veterinària Departament de Ciència Animal i dels Aliments Dissecting the Genetic Basis of Pig Intramuscular Fatty Acid Composition Ph. D, Thesis YULIAXIS RAMAYO CALDAS Supervisors Dr. Josep M. Folch Albareda Dr. Miguel Pérez Enciso BELLATERRA, JULY 2013 “I know one thing: that I know nothing” - Socrates “Solo se ve bien con el corazón. Lo esencial es invisible a los ojos” - Antoine de Saint-Exupéry SUMMARY Pigs (Sus scrofa) are relevant to humans, both as a source of food and as an animal model for scientific progress. Technological, nutritional and organoleptic properties of pork meat quality are highly dependent on lipid content and fatty acid (FA) composition. The molecular processes controlling FA composition and metabolism are highly interconnected and not fully understood. Elucidating these molecular processes will aid the technological development towards the improvement of pork meat quality and increase knowledge of FA metabolism underpinning metabolic diseases in humans. This thesis deals with FA metabolism in pigs, analyzed from different perspectives. We have described for the first time a map of Copy Number Variant Regions in the swine genome. It is based on whole genome SNP genotyping chips and some of them may have an effect on FA metabolism. A Genome-Wide Association Study (GWAS) approach was carried out in an Iberian x Landrace backcross across 32 traits related to intramuscular FA composition and 813 associated SNPs, distributed in 43 chromosomal regions, were identified. Besides, 53.5% of these chromosomal intervals had been reported to contain QTL for lipid traits in previous studies, being particular consistent the regions detected on chromosomes SSC4, SSC8 and SSC16. The liver transcriptome of two groups of phenotypically extreme pigs for intramuscular FA composition were sequenced using RNA-Seq. A differential expression analysis identified 55 protein-coding genes differentially-expressed (DE) between groups. Functional analysis revealed that these genes belong to biological functions, canonical pathways and three gene networks related to lipid and FA metabolism. Moreover, a connection was observed among the top three networks, suggesting that the DE genes detected in this study are linked and play an important and coordinate role in lipid and FA metabolism. The complexity and high interconnectivity of the molecular processes controlling FA metabolism identified by our GWAS and RNA-Seq results, motivated us to use a holistic gene network approach based on SNP-to-SNP co-association analysis. Supporting evidence for co-association network predictions were confirmed at tissue- specific manner by gene co-expression analysis in adipose and liver tissues. To the best of our knowledge, we reported for the first time in a vertebrate, the observation that interacting loci could jointly regulate the co-expression patterns of pairs of genes. The analysis of the topological properties of both the co-association and co- expression predicted gene networks, allowed the identification of key transcription factors (TF), candidate genes and metabolic pathways that are likely being determining meat quality and FA composition, as well as controlling energy homeostasis in pigs. Future studies targeting these genes, their pathways and interactions will continue to expand our knowledge of molecular control of FA metabolism and it might lead to the discovery of functional relevant mutations, unfolding new strategies for improve pork meat quality. Keywords: pig, meat quality, intramuscular fatty acid composition, lipid metabolism, candidate genes, copy number variant regions, GWAS, RNA-Seq, gene network, transcription factors, SNPs, co-association, co-expression. RESUMEN El cerdo (Sus scrofa domestica) constituye una de las principales fuentes de carne para la humanidad, y es también un excelente modelo animal para el estudio de enfermedades metabólicas en humanos. Los ácidos grasos (AG) juegan un papel importante actuando como moléculas de señalización celular en diversas rutas metabólicas asociadas a este tipo de patologías y además en la determinación de la calidad organoléptica y nutricional de la carne porcina. Como consecuencia, descifrar la base molecular del metabolismo de los AG podría contribuir al conocimiento de la base genética de enfermedades metabólicas humanas, además de favorecer el desarrollo tecnológico para mejorar la calidad de la carne porcina. El objetivo principal de esta tesis ha sido estudiar la base genética del contenido intramuscular de AG en cerdos. Para ello hemos utilizado diferentes, pero complementarias, aproximaciones analiticas. Inicialmente utilizamos los datos genotipicos del chip porcino de 62.612 SNPs (Illumina) e identificamos 49 regions con variación en el número de copias (CNVRs). Posteriormente, estudios de asociación masivos (GWAS) entre 32 fenotipos relacionados con el perfil intramuscular de AG y los SNPs del chip de Illumina nos permitieron identificar 813 SNPs que distribuidos en 43 regiones están significativamente asociados a 15 de los fenotipos analizados. Cabe destacar que el 53.5% de las regions identificadas contienen QTLs asociados a caracteres lipídicos. Esta observación sugiere que algunos de estos QTLs podrían tener un efecto pleiotropico, especialmente en las regiones identificadas en los cromosomas porcinos SSC4, SCC8 y SSC16. Al comparar el transcriptoma del hígado entre dos grupos de cerdos con fenotipos extremos para la composición intramuscular de AG se han identificado un total de 55 genes diferencialmente expresados. Análisis funcionales revelan que muchos de estos genes pertenecen a funciones biológicas, rutas metabólicas y redes génicas interconectadas relacionadas con el metabolismo de los lípidos y de los AG. Los resultados obtenidos en el estudio de asociación y en el análisis del trasncriptoma hepático confirman la complejidad e interconexión de los procesos moleculares relacionados con el metabolismo de los AG. Por lo tanto, nos propusimos utilizar un enfoque holístico basado en la co-asociación entre genes y/o SNPs a lo largo de varios fenotipos. Este análisis nos permitió identificar redes génicas, rutas metabólicas y factores de transcripción (TF) relevantes para el metabolismo de los AG. Por primera vez, interacciones génicas predichas en base a la co-asociación entre genes fueron válidadas mediante estudios de co-expresión génica en dos tejidos relevantes para el metabolismo de los AG en cerdos. Los resultados obtenidos muestran la utilidad de la biología de sistemas en el estudio de caracteres complejos. Sin embargo, serán necesarios estudios que integren diferentes disciplinas para establecer estrategias más eficientes que permitan mejorar la calidad de la carne de cerdo. Palabras claves: cerdos, calidad de la carne, ácidos grasos, metabolismo lipídico, gen candidato, CNV, GWAS, RNA-Seq, red génica, factor de trasncripción, SNPs, co- asociación, co-expresión. CONTENTS List of Tables……………………………………………………………................. 9 List of Figures……………………………………………………………………… 10 List of Publications……………………………………………………………….. 13 Related publications by the author..…………………………………………….. 14 Abbreviations……………………………………………………………………... 16 Chapter I. General Introduction………..……………………………………….. 18 Pork meat quality…………………………………………………………………. 20 Pig genomics……………………………………………………………………… 23 Mapping complex traits in pigs…………………………………………………… 27 Chapter II. Objectives………….………..………………………………………… 36 Chapter III. General Discussion..……..…………………………………………… 37 Detecting CNVR in the pig Genome……………………………………………….. 38 From GWAS to Gene-Networks through expression and co-expression analysis…………………………….………………………………………………... 44 Future directions…………………………….………………………………………. 50 Chapter IV. Conclusions………..……..……………………………………………. 54 Chapter V. References………..……..………………………………………………. 55 Chapter VI. Papers included in the thesis..………………………………………… 84 Paper I. Copy number variation in the porcine genome inferred from a 60 k SNP BeadChip. Paper II. Genome-wide association study for intramuscular fatty acid composition in an Iberian × Landrace cross. Paper III. Liver transcriptome profile in pigs with extreme phenotypes of intramuscular fatty acid composition. Chapter VII. Annexes. ..…………………………………………………………… 124 From SNP co-association to gene co-expression: Evidence of pleiotropic expression- QTL in lipid metabolism in pigs…………………………………………………….125 LIST OF TABLES General Introduction Table 1. Traits related with the sensorial, technological and nutritional quality of pork meat quality………...…………………………………………………………….. 21 Table 2. Assembly and annotation statistics corresponding to the Sscrofa10.2 genome assembly…………………………………………………………………… 25 Table 3. Functional polymorphisms affecting meat and/or carcass traits in the pig…………………………………………………………………………. 30 General Discussion Table 4. Summary of independent confirmations of 26 CNVRs ....………………. 40 Table 5. Positional concordance among candidate genes for FA metabolism located within CNVRs and QTL deposited in the pig QTL database for fatness related traits………………………………………………………………………… 42 Annexes Table 6. Description of the ten most connected nodes in the co-association network……………………………………………………………………..129 Table 7. Concordance validation between the co-association and the co-expression networks for the best TF trio and in the two tissues……………………….137 LIST OF FIGURES General Introduction Figure 1. Distribution of pork production among the autonomous communities in Spain..………………………………………………………………………. 19 Figure 2. Distribution of reported QTLs