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.

were excluded if they had low minor allele frequency (MAF) <5%. The deviation from Hardy-Weinberg equilibrium was not considered because a three-way crossbreed pig population was used, where deviation from Hardy Weinberg equilibrium can be expected due to discordant allele frequencies in the parental breeds. The average call rate for all samples was 99.8% ± 0.2. In total, 11,505 SNPs were removed because of low call rates (GenCall score < 0.4) or MAF.

Whole-genome expression profi ling
Gene expression profi ling of the liver was conducted with the same 150 animals. In brief, total RNA of the liver was isolated using TRI Reagent (Sigma, Taufkirchen, Germany) and used for target preparation for microarray hybridization. According to Affymetrix protocols, 500 ng of total RNA were reversely transcribed into cDNA, transcribed into cRNA, and labeled using the Affymetrix One cycle synthesis and labeling kit (Affymetrix, UK) to prepare antisense biotinylated RNA targets. Liver expression patterns were produced using 150 GeneChip Porcine Genome Arrays (Affymetrix). The quality of hybridization was assessed in all samples following the manufacturer's recommendations. Data were analyzed with the Affymetrix GCOS 1.1.1 software, using global scaling to a target signal of 500 for background correction and probe summarization. Data were then imported into the Expression Console software (Affymetrix) for subsequent analysis. First the data were processed with the MAS5.0 algorithm to generate probe cell intensity values and to evaluate presence and absence of transcripts. Using default settings with detection of P < 0.04 for "present"', P у 0.04 and P р 0.06 for "marginal", and P < 0.06 for "absent," only "present" calls were used. Then quantitative expression levels of the "present" transcripts were estimated, and normalization was performed using probe logarithmic intensity error (PLIER). PLIER data were used for further statistical analysis of expression levels. The microarray data related to all samples were deposited in the Gene Expression Omnibus public repository (GEO accession number GSE25445).

Correlation between traits and expression levels
Organismal phenotypes and expression levels were adjusted for systematic effects by ANOVA performed with the "Mixed" procedure of SAS software (SAS version 9.1, SAS Institute, Cary, NC) before analyzing their correlation. "Gender" was used as a fi xed effect, "sire" and "slaughter day" as random effects, and "carcass weight" as a covariate. Subsequently, Pearson correlation coeffi cients were calculated between the residuals of log2-transformed expression intensities and the trait "fat area" (FA) ( i.e., the area of the subcutaneous back fat on M. longissimus dorsi between the 13th and 14th rib. Genes that showed correlation at P р 0.05 were analyzed further. siderable proportion of these genes and their proteins, it is not yet clear what their downstream effects and functions are. Genetical genomics is a useful approach for studying the effect of genetic variation at the molecular level in biological systems. Genome-wide studies of gene expression are starting to provide information on genetic variation that affects gene expression levels ( 5 ). Recent studies in a variety of organisms have shown that levels of gene expression are often highly heritable (6)(7)(8)(9) and that for many genes, it is possible to map cis -and trans -acting factors using linkage (8)(9)(10)(11)(12) or association mapping ( 7,(13)(14)(15)(16).
We aimed to identify candidate genes for the fatness trait "fat area" by combining (i) genome-wide expression profi ling and genotyping, (ii) analyzing correlation of transcript abundance with organismal phenotypic characteristics (fatness/obesity), and (iii) studying subsequent association of genomic variation with variation of expression of transcripts showing trait-correlated expression to determine expression quantitative trait loci (eQTL). Performance-tested animals of commercial crossbred herds were analyzed using 24K expression and 60K genotyping microarrays.

Animals and tissue collection
Animal care and tissue collection processes followed the guidelines of the German Law of Animal Protection, and the experimental protocol was approved by the Institutional Animal Care Committee. This study was based on trait measurements, genotyping records, expression profi les, and genome-wide association analyses done with 150 performance-tested pigs from commercial herds of the crossbreed Pietrain×(German Large White × German Landrace) [Pi×(DE×DL)]. These animals represented 70 full-sibling families derived from 39 sires and 67 dams. The carcass and fat traits data were collected according to the guidelines of the Zentral Verband der Deutschen Schweineproduktion e.V. (ZDS) ( 17 ). Fat area was measured between the 13th and 14th rib on the M. longissimus dorsi (mean and standard error, 15.2 ± 3.2 cm 2 ).

Whole-genome scan and quality control
Illumina bead array technology was used to carry out all genotyping reactions in accordance with the manufacturer's protocol for the single nucleotide polymorphism (SNP) Infi nium HD assay (http://www.illumina.com). Genotyping was performed using the PorcineSNP60 BeadChip of Illumina (Illumina, Inc., San Diego, CA). In brief, 200 ng of DNA were used for genome-wide amplifi cation and subsequent fragmentation. DNA was hybridized to the 62,163 locus-specifi c 50mers covalently linked to the beads distributed on the surface of the microarray. Single-base extension of the oligos on the BeadChip was then performed using the captured DNA as a template, incorporating detectable labels on the Bead-Chip. The signals of each wavelength were determined using an Illumina iScan that converted the images to intensity data. The intensity data for each SNP were normalized and assigned a cluster position and genotype with the GenomeStudio software (Illumina, Inc.), and a quality score for each genotype was generated. A total of 150 samples were genotyped for 62,163 SNPs. Quality of the data was evaluated, and samples with call rates <95% were removed. Markers Fischer's exact test was used to calculate a P value for the networks. A score was determined according to the IPA 7 Feature Manual that refl ected the relevance of the network based on the number of regulated genes within a network and the size of the network.

Genome-wide association analysis
For the association analysis, log2-transformed expression levels of 150 individuals were used as traits for whole-genome association analysis with SNPs genotyped by a mixed model ANOVA using JMP Genomics (SAS Institute). "Genotype" and "gender" were used as fi xed effects, "sire" and "slaughter day" as random effects, and "carcass weight" as a covariate. The sequences fl anking the SNPs represented on the Illumina SNP chip, which were significantly associated with the level of transcription of any probe set of the Affymetrix expression microarrays, were assigned to the porcine genome sequences (Ensembl_Sscrofa_9, released April 2009). Annotation and localization of SNP sites and probe sets allowed discriminating cis -and trans -regulation. We defi ned an eQTL as " cis " if an associated SNP was located within an area less than 10 Mb from the probe set/gene. All other eQTL were considered as " trans ." Results are reported at thresholds of P < 1.4 × 10 Ϫ 9 (corresponding to a Bonferroni-corrected level of signifi cances of P р 0.05) and P < 10 Ϫ 5 . A signifi cance threshold of P < 10 Ϫ 5 was used,

Pathway analysis and network generation
On the basis of BLAST comparison of the Affymetrix porcine target sequences with the porcine genome sequence (Ensembl_ Sscrofa_9, released April 2009), a total of 20,689 of the 24,123 probe sets on the Affymetrix Porcine GeneChip were annotated ( 18 ) for use in this study. The list of signifi cant trait-correlated transcripts was analyzed, referring to predefi ned pathways and functional categories of the Ingenuity Pathways Knowledge Base (IPKB) using Ingenuity Pathways Analysis (IPA). Canonical pathways were also identifi ed from the IPA library, which were most signifi cant to the input data. The signifi cance of the association between the dataset and the predefi ned pathways and functional categories was measured by Fischer's exact test, providing a P value to determine the probability that the association between the genes in the dataset and the pathway was explained by chance alone. Furthermore, identifi ers of genes with correlated expression were assigned to the corresponding gene objects of the IPKB. These genes, called "focus genes," were overlaid onto a global molecular network developed from information contained in the IPKB. Networks of these focus genes were then algorithmically generated based on their connectivity. The functional analysis of members of a network identifi ed the biological functions (and/or diseases) that were most signifi cant to the genes in the network. 24K expression microarray. Out of 62,163 SNPs of the 60K SNP-Chip, a total of 52,727 SNPs with genotyping call rates greater than 95% and MAF greater than 5% were selected for genome-wide association analysis. The GeneChip Porcine Genome Arrays (Affymetrix) contained 24,123 probe sets representing 20,689 known genes ( 18 ), of which 10,945 probe sets were used for further analysis that showed consistent expression according to MAS5 analysis. Analysis of correlation between residuals derived from the mixed model analysis of the trait FA and the expression levels of 10,945 probe sets, respectively, revealed 663 probe sets that showed signifi cant correlation with FA (supple-representing a tradeoff between avoidance of false-positive associations while taking into account the likely higher dependence among the expression traits (correlated expression) and among the SNP genotypes (linkage disequilibrium).

Expression profi ling and genotyping
To characterize the genetic architecture of hepatic gene expression relevant to fatness traits, 150 performancetested pigs were genotyped using the PorcineSNP60K BeadChip (Illumina) and expression-profi led using the

Whole-genome association analyses for abundance of transcripts with FA-correlated expression (eQTL)
Association mapping of gene expression traits allowed us to identify eQTL (i.e., genomic regions that have a regulatory effect on the expression traits). Two types of eQTL could be distinguished: those that mapped near (less than 10 Mb from) the gene that encodes the transcript (" cis ") and those that mapped elsewhere in the genome (" trans "). To link eQTL to an organismal phenotype, we focused on the transcripts that showed expression levels correlated with the classical trait FA, and only these were subjected to association analysis using 52,727 SNPs. The association P values were adjusted to control for multiple testing using the stringent, conservative Bonferroni correction method to constrain the study-wise signifi cance level. In total, 52,727 SNPs were tested for association to each of the 663 expression traits. Therefore, the Bonferroni-adjusted P threshold was computed as 0.05/(52,727 × 663) = 1.43 × 10 Ϫ 9 . At this threshold, 555 eQTL were identifi ed, which were composed of 96 cis -and 364 trans -eQTL. However, the Bonferroni adjustment method is likely too stringent, given the dependence among expression traits and SNP genotypes. For example, the expression level of ALDH18A1 is associated with a number of SNPs in the surrounding sequences of this gene where strong linkage disequilibrium (LD) exists ( Fig. 2 ). Because the Bonferroni adjustment might be overly conservative, we set a more relaxed threshold at nominal values of P < 1.0 × 10 Ϫ 5 , which is used in many studies (http://www.genome.gov/gwastudies) ( 19,20 ). At this threshold, 663 probe sets revealed 4,727 eQTL with a mean negative log of the P value of 6.98 (supplementary Table II). For 3,745 of these, there was information on the genomic localization of both the SNP marker and the transcript represented by the probe set. This allowed assignment of either cis or trans regulation. In fact, 448 cis -eQTL and 3,297 trans -eQTL were detected with mean negative log of the P values of 7.81 and 6.90, respectively. For the rest of the eQTL (982), either the SNP marker or the corresponding probe set could not be localized. The 448 cis -eQTL corresponded to 71 genes, and the 3,297 trans -eQTL were related to 408 genes ( Table 2 ). Functional candidate genes with either a positive or negative correlation with FA and with either a cis -or trans -eQTL are shown in Tables 3-6 .

Biologically relevant pathways and network generation
To determine the biologically relevant networks beyond canonical pathways, pathway analysis was performed for gene sets with cis -(71 genes) and trans -(408 genes) eQTL using the IPKB. The networks describe functional relationships between gene products based on known interactions reported in the literature.  Table I). The correlation coeffi cients at P р 0.05 ranged |0.16-0.33| . FA was positively correlated with 213 genes and negatively correlated with 450 genes. These genes were functional candidates for traits related to fatness.

Detection of transcripts with fat area-correlated expression
We tested the list of signifi cant positively and negatively correlated functional candidate genes for the enrichment in functional annotation groups as defi ned in the IPKB. Genes with expression levels that correlated positively with FA mainly belonged to categories related to metabolism of various macromolecules and nutrients, whereas genes with negative correlations belonged to functions related to dynamic cellular processes ( Table 1 ).
To further refi ne the functional annotation of these gene sets with FA-correlated expression, their assignment to canonical pathways was explored using the IPKB. The positively and negatively correlated expression levels of functional candidate genes were uploaded as gene datasets for analysis. Signifi cantly enriched canonical pathways of the transcription profi les with positive correlation to FA belonged to pregnane X receptor (PXR) / retinoid X receptor (RXR) activation (SCD, GSTM2, CYP1A2, FOXO1, CYP3A4, and CYP7A1); xenobiotic metabolism signaling (PNPLA7, ALDH4A1, GSTT1, MAOB, GSTM2, CYP1A2, CYP3A4, FMO1, CES1, including EG:1066, and CITED2); lipopolysaccharide (LPS) / interleukin (IL)-1-mediated inhibition of RXR function (ALDH4A1, GSTT1, MAOB, GSTM2, SLC10A1, CYP3A4, CYP7A1, and FMO1); arginine and proline metabolism (CKB, ALDH4A1, MAOB, ASS1, and SAT2); bile acid biosynthesis (ALDH4A1, CYP3A4, CYP7A1, and ADHFE1); arachidonic acid metabolism (GSTT1, GPX3, CYP4F2, CYP1A2, CYP3A4, and PRDX6); glutathione metabolism (GSTT1, GPX3, GSTM2, and PRDX6); and serotonin receptor signaling (MAOB, PCBD1, and QDPR) ( Fig.1A ).  regions ( 21 ). The success of these studies prompted us to investigate the potential use of a commercial pig herd (outbred) for eQTL studies compared with the F2 population from our previous studies (22)(23)(24). Here, a wholegenome association study for hepatic gene expression profiles was carried out in a commercial herd of Pi×(DE×DL). The application of global gene expression analysis to studies of cell and tissues has provided a wealth of data relevant to complex traits. Accounting for 10,945 probe sets representing genes found expressed in liver, more than 70,000 eQTL can be expected at P < 10 Ϫ 5 . However, only 6% of liver transcripts out of the 10,945 probe sets were found to correlate with the trait "fat area." Among the carcass traits related to fatness that are measured during performance testing in commercial breeding schemes, "fat area" is highly heritable and can be considered an indicator of obesity and body composition ( 25 ). FA shows high genetic correlation with plasma leptin levels, which in turn are correlated with adipositas and habitus (25)(26)(27). Finding statistically signifi cant correlations between a trait and particular genes suggest a biologic relationship between them ( 22,28,29 ). IPA essentially evaluates the enrichment of particular biological processes and molecular functions of gene sets by examining information collected by databases such as Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and IPKB. In this study, many genes involved in lipid metabolism were highest ranking networks deduced from the list of functional candidate genes having either cis -eQTL or trans -eQTL. The network of cis -regulated genes comprised 19 focus genes belonging to lipid metabolism, small molecule biochemistry, and cell death ( Fig. 3A ). The second network of genes with trans -regulation contained 21 focus genes associated with lipid metabolism, small molecule biochemistry, and genetic disorder ( Fig. 3B ).

DISCUSSION
To get more insight into the genetic control of traits related to obesity and to address candidate genes for these organismal phenotypes, the genetics of expression was studied in porcine liver, a metabolically active tissue that is critical to a number of core biological processes and that plays a role not only in carcass traits in pigs but also in body composition, obesity and a number of common metabolic diseases in humans. We achieved detecting functional candidate genes for these traits by exploiting trait-correlated expression and assigning eQTL to these genes through genome-wide SNP-array based association analysis. Ghazalpour et al. ( 21 ) found that the mapping resolution for eQTL was signifi cantly greater in the outbred population than in F2 intercrosses and showed how this precise mapping could resolve previously identifi ed loci (in F2 intercross studies), which affect many different transcript levels (known as eQTL 'hotspots), into distinct positively correlated with FA, whereas genes with negative correlation were associated with cellular development, cellular growth, and proliferation. These fi ndings were in agreement with our previous conclusion that the genetic differences between obese and lean pig breeds were associated with upregulation of lipid metabolic pathways in obese pigs and upregulation of cell growth and/or maintenance and cell proliferation pathways in lean pigs ( 30 ).

Canonical pathway analysis
Canonical pathway analysis highlighted pathways involving signaling via PXR and RXR. These pathways play an important regulatory role in metabolic pathways (glucose, fatty acid, and cholesterol metabolism) ( 31 ). Individual genes from the lists that were detected here and assigned to the pathways of "PXR/RXR activation," "xenobiotic metabolism signaling," and "LPS/IL-1-mediated inhibition of RXR function" point to lipid metabolism instead of any function related to detoxifi cation. Several of the genes in these pathways have been shown to infl uence adipose tissue mass in other studies. For example, forkhead box protein O1 (FOXO1) plays an important role in modulating metabolic functions. It induces the expression of gluconeogenic genes and contributes to hyperglycemia during hepatic insulin resistance ( 32,33 ). The other two pathways (bile acid biosynthesis and arachidonic acid metabolism) are involved in cholesterol metabolism. Similar pathways were also found positively correlated with fat mass in mice ( 34,35 ). The conversion of cholesterol to bile acid is the major pathway for cholesterol catabolism. In bile acid pathways, cholesterol serves as a substrate that is converted into bile acid through the activation of cholesterol 7 ␣ hydroxylase (CYP7A1), the transcript for which is signifi cantly increased in pigs with a high fat area. This protein is the fi rst rate-limiting enzyme in bile acid synthesis in the liver. FoxO1 has been shown to stimulate rat Cyp7a1 gene transcription, and mutations in the FoxO1 binding site attenuated rat Cyp7a1 promoter activity ( 36 ). In our study, both FOXO1 and CYP7A1 transcripts were signifi cantly increased in pigs with a high fat area. The aldehyde dehydrogenase (ALDH )4A1 gene family, which is involved in arginine and proline metabolism, was also found to be signifi cantly correlated with FA. The increased expression of aldehyde dehydrogenase in mice with higher subcutaneous fat-pad mass has previously been reported ( 37 ). Another pathway with a signifi cant correlation with high FA was glutathione metabolism. Glutathione is a major endogenous cellular anti-oxidant. We found that four of the genes in this pathway were signifi cantly overexpressed in pigs with high FA. Glutathione peroxidase (GPX)3 is involved in cellular protection against oxidative damage through the reduction of peroxides ( 38 ). The cytosilic isoform of GPX3 has been associated with obesity ( 39 ), and there is evidence for causal relationship of GPX3 to obesity and for the idea that GPX3 overexpression modifi es insulin resistance ( 40 ). Recently, transgenic and   knockout mouse models of GPX3 resulted in signifi cant changes in obesity-related traits ( 41 ). PRDX6 is an antioxidant enzyme with calcium-independent phospholipase A2 activity ( 42 ). PRDX6 was identifi ed in a recent study as a possible candidate gene underlying a novel obesity locus on chromosome 1q24 in an isolated population of Cilento, Italy ( 43 ). Recently, a study of 21,000 individuals identifi ed PRDX6 associated with body mass index (BMI) as an indicator of obesity ( 44 ). Using the porcine model providing reliable phenotypes and samples, we were able to associate canonical pathways rather than only single genes with the organismal phenotype related to obesity. For some elements of these canonical pathways there are supporting data from functional studies in rodent models and epidemiological studies in human.
In this study, genes whose expression was found to be negatively correlated with fatness were involved in acute phase response signaling, protein ubiquitination, IL-10 signaling, and the antigen presentation pathway. These pathways, indicating a higher cellular turnover, activity, and reactivity, were also found to correlate with leptinmediated weight loss in ob/ob mice after leptin administration ( 45 ).

Whole-genome association analyses for abundance of transcripts with "fat area"-correlated expression (eQTL)
We present here the fi rst eQTL study based on wholegenome association analysis in a commercial crossbred population with focus on the FA trait. In the fi rst step, we analyzed the correlation between gene expression and the phenotype related to fatness, which revealed biologically meaningful relationships. In the second step, eQTL were identifi ed for transcripts that showed trait-correlated expression, which informed us about the genomic location of putative regulatory loci. This strategy reduced the number of several thousand eQTL that were not associated with FA. A genome-wide association screen for expression SNPs (eSNP) was performed on 150 pigs from commercial herds. The eSNPs were assessed using a mixed model analysis with "genotype" and "sex" as fi xed effects, "sire" and "slaughter days" as random effects, and "carcass weight" as a covariate. For the identifi cation of sequence variants having cis -and trans -regulatory effects on expression traits, SNPs located within a 10 Mb window at the location of the probe set were defi ned as cis -eQTL. The trans -acting eQTL represent transcripts whose abundance is regulated by loci remote from the genomic locus of each of these genes. In our study, the proportion of trans -eQTL was higher (70%) compared with cis -eQTL (10%) at P < 1.0 × 10 Ϫ 5 ; the remaining eQTL could not be assigned cis or trans . Most studies in model organisms have shown that the number of traits with signifi cant trans -eQTL is fewer than the number of expression traits with signifi cant cis -eQTL (9)(10)(11)46 ). Here, eQTL analysis was focused only on biologically meaningful relationship with the traits.

Biologically relevant pathways and network generation of cis -and trans -regulated genes
A network of cis -regulated genes comprised 19 focus genes with a score of 40 belonging to lipid metabolism, small molecule biochemistry, and cell death ( Fig. 2A ). In this network, 10 probe sets representing cis -regulated genes had hepatocyte nuclear factor 4 alpha (HNF4A), a common transcription factor. HNF4A, the protein encoded by this gene, is a nuclear transcription factor that binds DNA as a homodimer to the promoters of several genes involved in hepatic pathways. The 10 commonly cisregulated target genes of HNF4A are likely to exhibit polymorphisms affecting their own expression. Moreover, these polymorphic sites may be located in the target region of the common transcription factor HNF4A. Thus the coregulation of the 10 probe sets representing cis -regulated genes is likely due to differences in the affi nity of allelic variants of the HNF4A binding sites in their promoters. Another central component of the network was FOXO1. The transcript of this molecule not only significantly correlated with FA, as described above, but was also cis -and trans -regulated. Another node within the network is the tumor necrosis factor (TNF) interacting with GSTT1, PLIN4, ANXA5, SOD2, VTI1B, and SLC12A1, some of which were previously associated with fat traits in human ( 47,48 ). Perilipins (PLIN) are proteins localized at the surface of the lipid droplets in adipocytes, steroid-producing cells, and play a key role in the cellular regulation of triglyceride deposition and mobilization ( 49 ). PLIN polymorphisms have been associated with obesity-related phenotypes ( 47,50 ). Mitochondrial superoxide dismutase-2 (SOD-2) is an endogenous anti-oxidant enzyme. Gene variants of SOD-2 have been shown to be associated with coronary heart disease in diabetes mellitus ( 48 ). For many genes mentioned here, functional or epidemiological association to fatness traits has been shown in rodents and humans. Using the porcine model, we were able to provide functional and genetic links between these genes and to derive networks.
For the network of genes with trans -eQTL, central but not differentially regulated molecules were the retinoid receptor, estrogen receptor, thyroid hormone receptor, and growth hormone. This fi nding is consistent with published eQTL studies, which have indicated that trans -eQTL intervals encode many different gene classes, including receptors and their ligands, which can affect their target transcript via a context-dependent mechanism (51)(52)(53). Together, the central components of the trans -networks typically represented hierarchically high positions in signaling cascades, whereas downstream components of signaling pathways were represented in the cis -networks.
The correlation between expression levels and the organismal phenotype tended to be higher for trans -regulated genes. This correlation comprised genetic causal effects that the transcript may have on the phenotype as well as common dependence on environmental factors. Marker association with expression levels tended to be more signifi cant for cases with cis -regulation, where genetic effects were expected to be prevailing. As we noticed in our previous study, eQTL mapping with its potential to categorize cis -and trans -effects provided a way to discriminate between effect and cause of trait-associated differential expression ( 22-24, 54, 55 ). This fi nding was also supported by other studies in model organisms; environmentally induced changes in gene expression seem to act through prominent trans -effects ( 56,57 ).
The number of GWAS in humans and mice continues to rise rapidly. In the GWA studies to date, 10-15% of the top hits have affected a known eQTL in a public data set ( 58 ). Recently, a study has shown that trait-associated SNPs are more likely to be eQTLs ( 20 ). These fi ndings are robust across a range of thresholds for establishing eQTL ( P range from 10 Ϫ 4 to 10 Ϫ 8 ).
To our knowledge, this is the fi rst report of a comprehensive scan for eQTL associated with obesity in pigs. Many of the candidate genes previously identifi ed in humans and mice have been confi rmed here in the porcine model. Moreover, lists of novel candidate genes for fat traits were identifi ed for the fi rst time here in the pig model. These candidate genes for fat traits are valuable for human research due to the many similarities in physiology and genome between pigs and humans. The application of this information can facilitate the detection of trait-associated SNPs and a better understanding of the biology of complex traits.