Functional genomics and genetical genomics approaches towards elucidating networks of genes affecting meat performance in pigs

An Integrative Analysis of Transcriptome and GWAS Data to Identify Potential Candidate Genes Influencing Meat Quality Traits in Pigs

Understanding the genetic factors behind meat quality traits is of great significance to animal breeding and production. We previously conducted a genome-wide association study (GWAS) for meat quality traits in a White Duroc × Erhualian F2 pig population using Illumina porcine 60K SNP data. Here, we further investigate the functional candidate genes and their network modules associated with meat quality traits by integrating transcriptomics and GWAS information. Quantitative trait transcript (QTT) analysis, gene expression QTL (eQTL) mapping, and weighted gene co-expression network analysis (WGCNA) were performed using the digital gene expression (DGE) data from 493 F2 pig's muscle and liver samples. Among the quantified 20,108 liver and 23,728 muscle transcripts, 535 liver and 1,014 muscle QTTs corresponding to 416 and 721 genes, respectively, were found to be significantly (p < 5 × 10^-4) correlated with 22 meat quality traits measured on longissiums dorsi muscle (LM) or semimembranosus muscle (SM). Transcripts associated with muscle glycolytic potential (GP) and pH values were enriched for genes involved in metabolic process. There were 42 QTTs (for 32 genes) shared by liver and muscle tissues, of which 10 QTTs represent GP- and/or pH-related genes, such as JUNB, ATF3, and PPP1R3B. Furthermore, a genome-wide eQTL mapping revealed a total of 3,054 eQTLs for all annotated transcripts in muscle (p < 2.08 × 10^-5), including 1,283 cis-eQTLs and 1771 trans-eQTLs. In addition, WGCNA identified five modules relevant to glycogen metabolism pathway and highlighted the connections between variations in meat quality traits and genes involved in energy process. Integrative analysis of GWAS loci, eQTL, and QTT demonstrated GALNT15/GALNTL2 and HTATIP2 as strong candidate genes for drip loss and pH drop from postmortem 45 min to 24 h, respectively. Our findings provide valuable insights into the genetic basis of meat quality traits and greatly expand the number of candidate genes that may be valuable for future functional analysis and genetic improvement of meat quality.

Genome-wide analysis of expression QTL (eQTL) and allele-specific expression (ASE) in pig muscle identifies candidate genes for meat quality traits

BACKGROUND

Genetic analysis of gene expression level is a promising approach for characterizing candidate genes that are involved in complex economic traits such as meat quality. In the present study, we conducted expression quantitative trait loci (eQTL) and allele-specific expression (ASE) analyses based on RNA-sequencing (RNAseq) data from the longissimus muscle of 189 Duroc × Luchuan crossed pigs in order to identify some candidate genes for meat quality traits.

RESULTS

Using a genome-wide association study based on a mixed linear model, we identified 7192 cis-eQTL corresponding to 2098 cis-genes (p ≤ 1.33e-3, FDR ≤ 0.05) and 6400 trans-eQTL corresponding to 863 trans-genes (p ≤ 1.13e-6, FDR ≤ 0.05). ASE analysis using RNAseq SNPs identified 9815 significant ASE-SNPs in 2253 unique genes. Integrative analysis between the cis-eQTL and ASE target genes identified 540 common genes, including 33 genes with expression levels that were correlated with at least one meat quality trait. Among these 540 common genes, 63 have been reported previously as candidate genes for meat quality traits, such as PHKG1 (q-value = 1.67e-6 for the leading SNP in the cis-eQTL analysis), NUDT7 (q-value = 5.67e-13), FADS2 (q-value = 8.44e-5), and DGAT2 (q-value = 1.24e-3).

CONCLUSIONS

The present study confirmed several previously published candidate genes and identified some novel candidate genes for meat quality traits via eQTL and ASE analyses, which will be useful to prioritize candidate genes in further studies.

Transcriptome analysis reveals a molecular understanding of nicotinamide and butyrate sodium on meat quality of broilers under high stocking density

Background

In recent years, increased attention has been focused on breast muscle yield and meat quality in poultry production. Supplementation with nicotinamide and butyrate sodium can improve the meat quality of broilers. However, the potential molecular mechanism is not clear yet. This study was designed to investigate the effects of supplementation with a combination of nicotinamide and butyrate sodium on breast muscle transcriptome of broilers under high stocking density. A total of 300 21-d-old Cobb broilers were randomly allocated into 3 groups based on stocking density: low stocking density control group (L; 14 birds/m²), high stocking density control group (H; 18 birds/m²), and high stocking density group provided with a combination of 50 mg/kg nicotinamide and 500 mg/kg butyrate sodium (COMB; 18 birds/m²), raised to 42 days of age.

Results

The H group significantly increased cooking losses, pH decline and activity of lactate dehydrogenase in breast muscle when compared with the L group. COMB showed a significant decrease in these indices by comparison with the H group (P < 0.05). The transcriptome results showed that key genes involved in glycolysis, proteolysis and immune stress were up-regulated whereas those relating to muscle development, cell adhesion, cell matrix and collagen were down-regulated in the H group as compared to the L group. In contrast, genes related to muscle development, hyaluronic acid, mitochondrial function, and redox pathways were up-regulated while those associated with inflammatory response, acid metabolism, lipid metabolism, and glycolysis pathway were down-regulated in the COMB group when compared with the H group.

Conclusions

The combination of nicotinamide and butyrate sodium may improve muscle quality by enhancing mitochondrial function and antioxidant capacity, inhibiting inflammatory response and glycolysis, and promoting muscle development and hyaluronic acid synthesis.

Biotechnological bases of the development of cloned pig embryos

The term ‘clone’ in animal biotechnology refers to an organism derived from non-sexual reproduction, which is both a direct offspring and a genetic copy of the parent organism. To date, the pig appears to be the most interesting object in cloning research. Somatic cell nuclear transfer in pigs has a wide range of potential applications in various fields of human scientific and economic activities. However, the efficiency of producing cloned embryos in swine is still lower than that of other livestock species, in particular horses and cattle. Somatic cell nuclear transfer is a technically complex multi-stage technology, at each stage of which the pig oocytes, which are more susceptible to changes of surrounding conditions, are affected by various factors (mechanical, physical, chemical). At the stage of oocyte maturation, changes in the cell ultrastructures of the ooplasm occur, which play an important role in the subsequent nuclear reprogramming of the transferred donor cell. Before transfer to the oocyte donor somatic cells are synchronized in the G0/G1 stage of the cell cycle to ensure the normal ploidy of the cloned embryo. When removing the nucleus of pig oocytes maturated in vitro, it is necessary to pay attention to the problem of preserving the viability of cells, which were devoid of their own nuclear material. To perform the reconstruction, a somatic cell is placed, using micro-tools, in the perivitelline space, where the first polar body was previously located, or in the cytoplasm of an enucleated oocyte. The method of manual cloning involves the removal of the oocyte nucleus with subsequent fusion with the donor cell without the use of micromanipulation techniques. The increased sensitivity of oocytes to the environmental conditions causes special requirements for the choice of the system for in vitro culture of cloned pig embryos. In this work, we have reviewed the modern methods used for the production of cloned embryos and identified the technological issues that prevent improving the efficiency of somatic cloning of pigs.

About the existence of common determinants of gene expression in the porcine liver and skeletal muscle

BACKGROUND

The comparison of expression QTL (eQTL) maps obtained in different tissues is an essential step to understand how gene expression is genetically regulated in a context-dependent manner. In the current work, we have compared the transcriptomic and eQTL profiles of two porcine tissues (skeletal muscle and liver) which typically show highly divergent expression profiles, in 103 Duroc pigs genotyped with the Porcine SNP60 BeadChip (Illumina) and with available microarray-based measurements of hepatic and muscle mRNA levels. Since structural variation could have effects on gene expression, we have also investigated the co-localization of cis-eQTLs with copy number variant regions (CNVR) segregating in this Duroc population.

RESULTS

The analysis of differential expresssion revealed the existence of 1204 and 1490 probes that were overexpressed and underexpressed in the gluteus medius muscle when compared to liver, respectively (|fold-change| > 1.5, q-value < 0.05). By performing genome scans in 103 Duroc pigs with available expression and genotypic data, we identified 76 and 28 genome-wide significant cis-eQTLs regulating gene expression in the gluteus medius muscle and liver, respectively. Twelve of these cis-eQTLs were shared by both tissues (i.e. 42.8% of the cis-eQTLs identified in the liver were replicated in the gluteus medius muscle). These results are consistent with previous studies performed in humans, where 50% of eQTLs were shared across tissues. Moreover, we have identified 41 CNVRs in a set of 350 pigs from the same Duroc population, which had been genotyped with the Porcine SNP60 BeadChip by using the PennCNV and GADA softwares, but only a small proportion of these CNVRs co-localized with the cis-eQTL signals.

CONCLUSION

Despite the fact that there are considerable differences in the gene expression patterns of the porcine liver and skeletal muscle, we have identified a substantial proportion of common cis-eQTLs regulating gene expression in both tissues. Several of these cis-eQTLs influence the mRNA levels of genes with important roles in meat (CTSF) and carcass quality (TAPT1), lipid metabolism (TMEM97) and obesity (MARC2), thus evidencing the practical importance of dissecting the genetic mechanisms involved in their expression.

Comparing the mRNA expression profile and the genetic determinism of intramuscular fat traits in the porcine gluteus medius and longissimus dorsi muscles

BACKGROUND

Intramuscular fat (IMF) content and composition have a strong impact on the nutritional and organoleptic properties of porcine meat. The goal of the current work was to compare the patterns of gene expression and the genetic determinism of IMF traits in the porcine gluteus medius (GM) and longissimus dorsi (LD) muscles.

RESULTS

A comparative analysis of the mRNA expression profiles of the pig GM and LD muscles in 16 Duroc pigs with available microarray mRNA expression measurements revealed the existence of 106 differentially expressed probes (fold-change > 1.5 and q-value < 0.05). Amongst the genes displaying the most significant differential expression, several loci belonging to the Hox transcription factor family were either upregulated (HOXA9, HOXA10, HOXB6, HOXB7 and TBX1) or downregulated (ARX) in the GM muscle. Differences in the expression of genes with key roles in carbohydrate and lipid metabolism (e.g. FABP3, ORMDL1 and SLC37A1) were also detected. By performing a GWAS for IMF content and composition traits recorded in the LD and GM muscles of 350 Duroc pigs, we identified the existence of one region on SSC14 (110-114 Mb) displaying significant associations with C18:0, C18:1(n-7), saturated and unsaturated fatty acid contents in both GM and LD muscles. Moreover, we detected several genome-wide significant associations that were not consistently found in both muscles. Further studies should be performed to confirm whether these associations are muscle-specific. Finally, the performance of an eQTL scan for 74 genes, located within GM QTL regions and with available microarray measurements of gene expression, made possible to identify 14 cis-eQTL regulating the expression of 14 loci, and six of them were confirmed by RNA-Seq.

CONCLUSIONS

We have detected significant differences in the mRNA expression patterns of the porcine LD and GM muscles, evidencing that the transcriptomic profile of the skeletal muscle tissue is affected by anatomical, metabolic and functional factors. A highly significant association with IMF composition on SSC14 was replicated in both muscles, highlighting the existence of a common genetic determinism, but we also observed the existence of a few associations whose magnitude and significance varied between LD and GM muscles.

Genetic control of longissimus dorsi muscle gene expression variation and joint analysis with phenotypic quantitative trait loci in pigs

Background

Economically important growth and meat quality traits in pigs are controlled by cascading molecular events occurring during development and continuing throughout the conversion of muscle to meat. However, little is known about the genes and molecular mechanisms involved in this process. Evaluating transcriptomic profiles of skeletal muscle during the initial steps leading to the conversion of muscle to meat can identify key regulators of polygenic phenotypes. In addition, mapping transcript abundance through genome-wide association analysis using high-density marker genotypes allows identification of genomic regions that control gene expression, referred to as expression quantitative trait loci (eQTL). In this study, we perform eQTL analyses to identify potential candidate genes and molecular markers regulating growth and meat quality traits in pigs.

Results

Messenger RNA transcripts obtained with RNA-seq of longissimus dorsi muscle from 168 F2 animals from a Duroc x Pietrain pig resource population were used to estimate gene expression variation subject to genetic control by mapping eQTL. A total of 339 eQTL were mapped (FDR ≤ 0.01) with 191 exhibiting local-acting regulation. Joint analysis of eQTL with phenotypic QTL (pQTL) segregating in our population revealed 16 genes significantly associated with 21 pQTL for meat quality, carcass composition and growth traits. Ten of these pQTL were for meat quality phenotypes that co-localized with one eQTL on SSC2 (8.8-Mb region) and 11 eQTL on SSC15 (121-Mb region). Biological processes identified for co-localized eQTL genes include calcium signaling (FERM, MRLN, PKP2 and CHRNA9), energy metabolism (SUCLG2 and PFKFB3) and redox hemostasis (NQO1 and CEP128), and results support an important role for activation of the PI3K-Akt-mTOR signaling pathway during the initial conversion of muscle to meat.

Conclusion

Co-localization of eQTL with pQTL identified molecular markers significantly associated with both economically important phenotypes and gene transcript abundance. This study reveals candidate genes contributing to variation in pig production traits, and provides new knowledge regarding the genetic architecture of meat quality phenotypes.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5386-2) contains supplementary material, which is available to authorized users.

Integrated Genome-wide association and hypothalamus eQTL studies indicate a link between the circadian rhythm-related gene PER1 and coping behavior

Animal personality and coping styles are basic concepts for evaluating animal welfare. Struggling response of piglets in so-called backtests early in life reflects their coping strategy. Behavioral reactions of piglets in backtests have a moderate heritability, but their genetic basis largely remains unknown. Here, latency, duration and frequency of struggling attempts during one-minute backtests were repeatedly recorded of piglets at days 5, 12, 19, and 26. A genome-wide association study for backtest traits revealed 465 significant SNPs (FDR ≤ 0.05) mostly located in QTL (quantitative trait locus) regions on chromosome 3, 5, 12 and 16. In order to capture genes in these regions, 37 transcripts with significant SNPs were selected for expressionQTL analysis in the hypothalamus. Eight genes (ASGR1, CPAMD8, CTC1, FBXO39, IL19, LOC100511790, RAD51B, UBOX5) had cis- and five (RANGRF, PER1, PDZRN3, SH2D4B, LONP2) had trans-expressionQTL. In particular, for PER1, with known physiological implications for maintenance of circadian rhythms, a role in coping behavior was evidenced by confirmed association in an independent population. For CTC1 a cis-expression QTL and the consistent relationship of gene polymorphism, mRNA expression level and backtest traits promoted its link to coping style. GWAS and eQTL analyses uncovered positional and functional gene candidates for coping behavior.

Identification of quantitative trait transcripts for growth traits in the large scales of liver and muscle samples

Growth-related traits are economically important traits to the pig industry. Identification of causative gene and mutation responsible for growth-related QTL will facilitate the improvement of pig growth through marker-assisted selection. In this study, we applied whole genome gene expression and quantitative trait transcript (QTT) analyses in 497 liver and 586 longissimus dorsi muscle samples to identify candidate genes and dissect the genetic basis of pig growth in a white Duroc × Erhualian F2 resource population. A total of 20,108 transcripts in liver and 23,728 transcripts in muscle with expression values were used for association analysis between gene expression level and phenotypic value. At the significance threshold of P < 0.0005, we identified a total of 169 and 168 QTTs for nine growth-related traits in liver and muscle, respectively. We also found that some QTTs were correlated to more than one trait. The QTTs identified here showed high tissue specificity. We did not identify any QTTs that were associated with one trait in both liver and muscle. Through an integrative genomic approach, we identified SDR16C5 as the important candidate gene in pig growth trait. These findings contribute to further identification of the causative genes for porcine growth traits and facilitate improvement of pig breeding.

Genetic dissection of blood lipid traits by integrating genome-wide association study and gene expression profiling in a porcine model

BACKGROUND

Serum concentrations of total cholesterol (TC), low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C) and triglycerides (TG) are highly heritable traits that are used clinically to evaluate risk for cardiovascular disease in humans. In this study, we applied a genome-wide association study (GWAS) in 1,075 pigs from two populations and gene expression studies on 497 liver samples to dissect the genetic basis of serum lipids in a pig model.

RESULTS

We totally identified 8, 5, 2 and 3 genomic loci harboring 109 SNPs that were significantly associated with LDL-C, TC, TG and the ratio of HDL-C/LDL-C in two experimental populations, respectively. In the F2 population, the most prominent SNP was identified at the SSC3: 124,769,847 bp where APOB is the well-known candidate gene. However, in the Sutai population, the most number of significant SNPs was identified at SSC2: 64.97-82.22 Mb where LDLR was identified as the candidate gene. Furthermore, we firstly reported 4 novel genomic loci in pigs harboring the LDL-C-associated SNPs. We also observed obvious population heterogeneity in the two tested populations. Through whole-genome gene expression analysis, we detected 718 trait-correlated expressions. Many of these transcripts correspond to candidate genes for blood lipids in humans. The GWAS mapped 120 cis-eQTLs and 523 trans-eQTLs for these transcripts. One gene encoding the transcript gnl|UG|Ssc#S35330332 stands out to be an important candidate gene for LDL-C by an integrative analysis of GWAS, eQTL and trait-associated expression.

CONCLUSIONS

We identified the genomic regions or candidate genes associated with blood lipids by an integrative analysis of GWAS, QTT and eQTL mapping in pigs. The findings would benefit the further identification of the causative genes for blood lipid traits in both pigs and humans.

A genetical genomics approach reveals new candidates and confirms known candidate genes for drip loss in a porcine resource population

In this study lean meat water-holding capacity (WHC) of a Duroc × Pietrain (DuPi) resource population with corresponding genotypes and transcriptomes was investigated using genetical genomics. WHC was characterized by drip loss measured in M. longissimus dorsi. The 60K Illumina SNP chips identified genotypes of 169 F2 DuPi animals. Whole-genome transcriptomes of muscle samples were available for 132 F2 animals using the Affymetrix 24K GeneChip® Porcine Genome Array. Performing genome-wide association studies of transcriptional profiles, which are correlated with phenotypes, allows elucidation of cis- and trans-regulation. Expression levels of 1,228 genes were significantly correlated with drip loss and were further analyzed for enrichment of functional annotation groups as defined by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathways. A hypergeometric gene set enrichment test was performed and revealed glycolysis/glyconeogenesis, pentose phosphate pathway, and pyruvate metabolism as the most promising pathways. For 267 selected transcripts, expression quantitative trait loci (eQTL) analysis was performed and revealed a total of 1,541 significant associations. Because of positional accordance of the gene underlying transcript and the eQTL location, it was possible to identify eight eQTL that can be assumed to be cis-regulated. Comparing the results of gene set enrichment and the eQTL detection tests, molecular networks and potential candidate genes, which seemed to play key roles in the expression of WHC, were detected. The α-1-microglobulin/bikunin precursor (AMBP) gene was assumed to be cis-regulated and was part of the glycolysis pathway. This approach supports the identification of trait-associated SNPs and the further biological understanding of complex traits.

Genetically modified pigs for biomedical research

During the last two decades, pigs have been used to develop some of the most important large animal models for biomedical research. Advances in pig genome research, genetic modification (GM) of primary pig cells and pig cloning by nuclear transfer, have facilitated the generation of GM pigs for xenotransplantation and various human diseases. This review summarizes the key technologies used for generating GM pigs, including pronuclear microinjection, sperm-mediated gene transfer, somatic cell nuclear transfer by traditional cloning, and somatic cell nuclear transfer by handmade cloning. Broadly used genetic engineering tools for porcine cells are also discussed. We also summarize the GM pig models that have been generated for xenotransplantation and human disease processes, including neurodegenerative diseases, cardiovascular diseases, eye diseases, bone diseases, cancers and epidermal skin diseases, diabetes mellitus, cystic fibrosis, and inherited metabolic diseases. Thus, this review provides an overview of the progress in GM pig research over the last two decades and perspectives for future development.

Segregation of regulatory polymorphisms with effects on the gluteus medius transcriptome in a purebred pig population

Background

The main goal of the present study was to analyse the genetic architecture of mRNA expression in muscle, a tissue with an outmost economic importance for pig breeders. Previous studies have used F₂ crosses to detect porcine expression QTL (eQTL), so they contributed with data that mostly represents the between-breed component of eQTL variation. Herewith, we have analysed eQTL segregation in an outbred Duroc population using two groups of animals with divergent fatness profiles. This approach is particularly suitable to analyse the within-breed component of eQTL variation, with a special emphasis on loci involved in lipid metabolism.

Methodology/Principal Findings

GeneChip Porcine Genome arrays (Affymetrix) were used to determine the mRNA expression levels of gluteus medius samples from 105 Duroc barrows. A whole-genome eQTL scan was carried out with a panel of 116 microsatellites. Results allowed us to detect 613 genome-wide significant eQTL unevenly distributed across the pig genome. A clear predominance of trans- over cis-eQTL, was observed. Moreover, 11 trans-regulatory hotspots affecting the expression levels of four to 16 genes were identified. A Gene Ontology study showed that regulatory polymorphisms affected the expression of muscle development and lipid metabolism genes. A number of positional concordances between eQTL and lipid trait QTL were also found, whereas limited evidence of a linear relationship between muscle fat deposition and mRNA levels of eQTL regulated genes was obtained.

Conclusions/Significance

Our data provide substantial evidence that there is a remarkable amount of within-breed genetic variation affecting muscle mRNA expression. Most of this variation acts in trans and influences biological processes related with muscle development, lipid deposition and energy balance. The identification of the underlying causal mutations and the ascertainment of their effects on phenotypes would allow gaining a fundamental perspective about how complex traits are built at the molecular level.

Biotechnologies for the management of genetic resources for food and agriculture

In recent years, the land area under agriculture has declined as also has the rate of growth in agricultural productivity while the demand for food continues to escalate. The world population now stands at 7 billion and is expected to reach 9 billion in 2045. A broad range of agricultural genetic diversity needs to be available and utilized in order to feed this growing population. Climate change is an added threat to biodiversity that will significantly impact genetic resources for food and agriculture (GRFA) and food production. There is no simple, all-encompassing solution to the challenges of increasing productivity while conserving genetic diversity. Sustainable management of GRFA requires a multipronged approach, and as outlined in the paper, biotechnologies can provide powerful tools for the management of GRFA. These tools vary in complexity from those that are relatively simple to those that are more sophisticated. Further, advances in biotechnologies are occurring at a rapid pace and provide novel opportunities for more effective and efficient management of GRFA. Biotechnology applications must be integrated with ongoing conventional breeding and development programs in order to succeed. Additionally, the generation, adaptation, and adoption of biotechnologies require a consistent level of financial and human resources and appropriate policies need to be in place. These issues were also recognized by Member States at the FAO international technical conference on Agricultural Biotechnologies for Developing Countries (ABDC-10), which took place in March 2010 in Mexico. At the end of the conference, the Member States reached a number of key conclusions, agreeing, inter alia, that developing countries should significantly increase sustained investments in capacity building and the development and use of biotechnologies to maintain the natural resource base; that effective and enabling national biotechnology policies and science-based regulatory frameworks can facilitate the development and appropriate use of biotechnologies in developing countries; and that FAO and other relevant international organizations and donors should significantly increase their efforts to support the strengthening of national capacities in the development and appropriate use of pro-poor agricultural biotechnologies.

Genetic determinants for intramuscular fat content and water-holding capacity in mice selected for high muscle mass

Intramuscular fat content and water-holding capacity are important traits in livestock as they influence meat quality, nutritive value of the muscle, and animal health. As a model for livestock, two inbred lines of the Berlin Muscle Mouse population, which had been long-term selected for high muscle mass, were used to identify genomic regions affecting intramuscular fat content and water-holding capacity. The intramuscular fat content of the Musculus longissimus was on average 1.4 times higher in BMMI806 than in BMMI816 mice. This was accompanied by a 1.5 times lower water-holding capacity of the Musculus quadriceps in BMMI816 mice. Linkage analyses with 332 G(3) animals of reciprocal crosses between these two lines revealed quantitative trait loci for intramuscular fat content on chromosome 7 and for water-holding capacity on chromosome 2. In part, the identified loci coincide with syntenic regions in pigs in which genetic effects for the same traits were found. Therefore, these muscle-weight-selected mouse lines and the produced intercross populations are valuable genetic resources to identify genes that could also contribute to meat quality in other species.

Integrating expression profiling and whole-genome association for dissection of fat traits in a porcine model

Traits related to fatness, important as economic factors in pork production, are associated with serious diseases in humans. Genetical genomics is a useful approach for studying the effects of genetic variation at the molecular level in biological systems. Here we applied a whole-genome association analysis to hepatic gene expression traits, focusing on transcripts with expression levels that correlated with fatness traits in a porcine model. A total of 150 crossbred pigs [Pietrain × (German Large White × German Landrace)] were studied for transcript levels in the liver. The 24K Affymetrix expression microarrays and 60K Illumina single nucleotide polymorphism (SNP) chips were used for genotyping. A total of 663 genes, whose expression significantly correlated with the trait "fat area," were analyzed for enrichment of functional annotation groups as defined in the Ingenuity Pathways Knowledge Base (IPKB). Genes involved in metabolism of various macromolecules and nutrients as well as functions related to dynamic cellular processes correlated with fatness traits. Regions affecting the transcription levels of these genes were mapped and revealed 4,727 expression quantitative trait loci (eQTL) at P < 10⁻⁵, including 448 cis-eQTL. In this study, genome-wide association analysis of trait-correlated expression was successfully used in a porcine model to display molecular networks and list genes relevant to fatness traits.

Identification of expression QTL (eQTL) of genes expressed in porcine M. longissimus dorsi and associated with meat quality traits

BACKGROUND

Genetic analysis of transcriptional profiles is a promising approach for identifying and dissecting the genetics of complex traits like meat performance. Accordingly, expression levels obtained by microarray analysis were taken as phenotypes in a linkage analysis to map eQTL. Moreover, expression levels were correlated with traits related to meat quality and principle components with high loadings of these traits. By using an up-to-date annotation and localization of the respective probe-sets, the integration of eQTL mapping data and information of trait correlated expression finally served to point to candidate genes for meat quality traits.

RESULTS

Genome-wide transcriptional profiles of M. longissimus dorsi RNAs samples of 74 F2 animals of a pig resource population revealed 11,457 probe-sets representing genes expressed in the muscle. Linkage analysis of expression levels of these probe-sets provided 9,180 eQTL at the suggestive significance threshold of LOD > 2. We mapped 653 eQTL on the same chromosome as the corresponding gene and these were designated as 'putative cis-eQTL'. In order to link eQTL to the traits of interest, probe-sets were addressed with relative transcript abundances that showed correlation with meat quality traits at p ≤ 0.05. Out of the 653 'putative cis-eQTL', 262 transcripts were correlated with at least one meat quality trait. Furthermore, association of expression levels with composite traits with high loadings for meat quality traits generated by principle component analysis were taken into account leading to a list of 85 genes exhibiting cis-eQTL and trait dependent expression.

CONCLUSION

Holistic expression profiling was integrated with QTL analysis for meat quality traits. Correlations between transcript abundance and meat quality traits, combined with genetic positional information of eQTL allowed us to prioritise candidate genes for further study.