<title>Background</title> <p>Boar taint is an off-odour and off-flavour in pig carcasses that is primarily caused by high levels of 16-androstene steroids [##UREF##0##1##] and/or skatole [##UREF##1##2##] in adipose tissue. Male pigs used for meat production are normally castrated very early in life to prevent boar taint in the meat. However, castration also removes the source of natural anabolic androgens that stimulate lean growth and, as a result, uncastrated males have improved feed efficiency and greater lean yield of the carcass compared to barrows [##UREF##2##3##]. Detection of genetic factors influencing boar taint may facilitate implementation of selective breeding practices to produce pigs with little or no taint. Because of unfavourable correlations between androstenone and other sex steroids [##REF##11465375##4##], a direct selection against high levels of androstenone would result in decreased production of testosterone and estrogens, with associated negative effects on performance and sexual maturation. Therefore, a comprehensive understanding of the complex genetic system controlling boar taint is required before genetic improvement can be achieved.</p> <p>Androstenone (5α-androst-16-en-3-one) is produced in the testis and is transported by blood to the salivary gland where it functions as a pheromone to stimulate mating stance in gilts [##REF##3516959##5##]. It is metabolised in the liver, producing α-androstenol and β-androstenol [##REF##15013816##6##,##REF##6648034##7##], and deficient degradation may lead to the accumulation of androstenone in adipose tissue. Skatole (3-methylindole) is produced by the metabolism of tryptophan from the gut and is also catabolised in liver. Diaz <italic>et al</italic>. [##REF##10497141##8##] identified seven metabolites of skatole in pig liver microsomes. An increase in levels of skatole in adipose tissue occurs in boars around puberty, but not in barrows or sows, indicating that skatole metabolism is regulated by testicular steroids [##REF##10064031##9##,##REF##12044562##10##].</p> <p>Enzymes responsible for metabolism of androstenone and skatole in liver have been identified but relatively few genes have been investigated. Several studies have shown that 3β-hydroxysteroid dehydrogenase (3β-HSD) is involved in androstenone metabolism, with low mRNA, protein and enzyme expression levels correlated to high androstenone levels in adipose tissue [##REF##15013816##6##,##REF##16971583##11##,##REF##17686889##12##]. Cytochrome P450 2E1 (CYP2E1) is involved in metabolism of skatole [##REF##9535343##13##,##REF##9303469##14##] and deficient CYP2E1 induction is suggested to be a major cause of high levels of skatole in adipose tissue [##REF##12044562##10##]. A relationship between metabolism of androstenone and skatole in liver has been observed [##REF##10064031##9##], with evidence to suggest that androstenone blocks CYP2E1 induction by skatole [##REF##12044562##10##]. Sulfotransferase enzymes including SULT2A1, SULT2B1 and SULT1A1 are associated with metabolic clearance of androstenone and skatole since they are known to conjugate steroid hormones and drugs into more water soluble compounds that facilitates excretion [##REF##17686889##12##,##REF##16595701##15##, ####REF##15644508##16##, ##REF##17210083##17##, ##REF##15014971##18##, ##REF##9535344##19##, ##REF##12746105##20####12746105##20##].</p> <p>Although both androstenone and skatole are responsible for boar taint, androstenone seems to be the biggest problem in the Duroc and Norwegian Landrace populations [##UREF##3##21##]. Levels of androstenone are also predominantly genetically determined whereas levels of skatole are more heavily influenced by feeding and environmental factors. Levels of androstenone will be a function of both production in testes and elimination in liver. We have previously studied differentially expressed genes related to the production of androstenone in boar testis [##REF##17210083##17##], and the objective of this study was to investigate global gene expression profiles related to the degradation of androstenone in pig livers. Microarray technology was used to obtain gene expression profiles from liver samples of Duroc (D) and Norwegian Landrace (NL) boars with extreme high or low levels of androstenone in adipose tissue. These samples were hybridised together with a common reference sample to cDNA microarrays representing approximately 20,000 porcine gene transcripts. A total of 29 animals were selected from each group, resulting in a total of 116 arrays. A subset of the differentially expressed genes was verified using a quantitative PCR-based method.</p>
<title>Materials and methods</title> <title>Animals</title> <p>Animals used in this study were Duroc (D) and Norwegian Landrace (NL) boars from NORSVIN's three boar testing stations. The D boars were on average 156 days old at 100 kg live weight compared to the NL boars that were on average 143 days old at 100 kg live weight. The boars were slaughtered on average 14 days later. Tissue samples from liver were frozen in liquid N<sub>2 </sub>immediately after slaughter and stored at -80°C until used for RNA extraction as described below. Samples from adipose tissue were collected from the neck at slaughter and stored at -20°C until used for androstenone measurements. Androstenone levels were measured by a modified time-resolved fluoroimmunoassay at the hormone laboratory, Norwegian School of Veterinary Science (NVH) [##UREF##8##61##] using an antiserum produced at NVH [##UREF##9##62##]. The androstenone measurements were performed on more than 2500 boars, and statistical power calculations showed that selecting animals from each tail of the androstenone distribution would yield sufficient power to detect differentially expressed genes with a limited number of arrays. The 30 most extreme boars from each tail of the androstenone level distribution were therefore selected. Due to poor RNA quality for two of the samples, 29 samples were subsequently used from each group. 42 of the animals used were the same individuals as those used in our previous study examining gene expression in boar testis [##REF##17210083##17##], with 16 being new animals as liver samples were not available for all of the previously studied animals. Average androstenone levels for the selected boars were 1.17 ppm and 3.22 ppm for NL and D, respectively. (See additional file ##SUPPL##4##5##: Androstenone levels). Average values for the groups were 5.95 ± 2.04 ppm for NL high (NLH), 0.14 ± 0.04 ppm for NL low (NLL), 11.57 ± 3.2 ppm for D high (DH) and 0.37 ± 0.17 ppm for D low (DL). In order to reduce family effects, a maximum of two and three half sibs were chosen from NL and D, respectively. The selected animals were used for expression profiling by microarrays and for the following verification of selected genes by rcPCR.</p> <title>Expression profiling using microarrays</title> <p>The present work utilises and extends methods described in our previous microarray experiment [##REF##17210083##17##]. The porcine cDNA microarrays were produced at the Faculty of Agricultural Sciences, University of Aarhus and contained 27,774 features printed in duplicates. 26,877 features were PCR products amplified from cDNA clones produced by the Sino-Danish Porcine Genome Sequencing project [##REF##17407547##63##,##REF##15656902##64##], and 867 were control features. The 26,877 features represent approximately 20K gene transcripts. Additional information about the porcine cDNA microarray can be found at NCBIs Gene Expression Omnibus (GEO, [##UREF##10##65##]) using the platform accession number GPL3585. This is a different batch of microarrays compared to the one we used in the testis experiment.</p> <p>Total RNA was extracted from liver tissue using Qiagen's RNeasy midi kit according to manufacturer's instruction (Qiagen, CA, USA). RNA quantity was measured on a NanoDrop ND-1000 Spectrophotometer (NanoDrop Technologies, DE, USA) and RNA quality was evaluated by the 28S:18S rRNA ratio using a RNA 6000 Nano LabChip<sup>® </sup>Kit on 2100 Bioanalyzer (Agilent Technologies, CA, USA). The microarray experiment was conducted as a common reference design using RNA purified from liver tissue sampled from an unrelated Danish Landrace × Hampshire pig as a reference. Aminoallyl-cDNA was synthesised from 15 μg of total RNA using the SuperScript indirect cDNA labelling system (Invitrogen Corporation, CA, USA) and labelled using ARES Alexa Fluor labelling kits (Molecular Probes, OR, USA). Amino-modified and fluorescently labelled cDNA was purified using NucleoSpin 96 Extract II PCR Clean-up kits (Macherey-Nagel, Düren, Germany). The individual samples were labelled with Alexa Fluor 647 and the reference was labelled with Alexa Fluor 555. "Green" spike-in RNA from the Lucidea Universal ScoreCard (Amersham Biosciences) was added to the individuals and "red" spike-in RNA was added to the reference. Hybridisation was performed in a Discovery XT hybridisation station (Ventana Discovery Systems, AZ, USA), followed by manual washing and drying by centrifugation. The microarrays were scanned using a ScanArray Express scanner (Perkin Elmer Inc., MA, USA) and image analyses were conducted using GenePix Pro 6 software (version 6.0.1.26, Molecular Devices Corp., CA, USA). Statistical analyses were carried out in R version 2.3.1 [##UREF##11##66##] using the software package Linear Models for Microarray Analysis (Limma version 2.7.2) [##REF##14597310##67##, ####UREF##12##68##, ##UREF##13##69####13##69##] which is part of the Bioconductor package [##REF##15461798##70##]. Log-transformed ratios of mean foreground intensities (not background corrected) were print tip loess normalised. The duplicate correlation function in Limma was used to consider duplicate printing of each feature. To evaluate the analyses, MA-plots (M = log<sub>2</sub>594/log<sub>2</sub>488, A = (log<sub>2</sub>594 + log<sub>2</sub>488)/2), image plots and box plots were constructed using the Limma tools for visualisation both before and after normalisation (see additional file ##SUPPL##5##6##: Boxplots). To assess differential expression, Limma uses linear models in combination with an empirical Bayes method to moderate the standard errors of the estimated log-fold changes [##UREF##12##68##]. The nominal p-values were corrected for multiple testing by false discovery rates (FDR) using Benjamini and Hochberg approach [##UREF##14##71##]. Each of the groups DH, DL, NLH and NLL animals were hybridised in individual batches, causing some confounding between androstenone levels and hybridisation batch. However, as neither boxplots nor MA-plots were found to differ between hybridisation batches, it was assumed that the impact of hybridisation batch on the obtained data can be ignored. The top 1% of the genes was considered for further analyses. The array features were mapped to a LocusLink identifier and an annotation package was built using the Bioconductor package AnnBuilder (version 1.9.14). GO terms (p &lt; 0.01 and more than 10 significant genes) were analysed for overrepresentation using the GOHyperG function of the Bioconductor package GOstats (version 1.6.0) [##REF##17098774##72##,##REF##18030337##73##]. More detailed descriptions of the microarray experiments are available at the GEO database through the series accession number GSE 11073.</p> <title>Quantitative real competitive PCR analysis</title> <p>A real competitive PCR (rcPCR) gene expression analysis was used to verify a subset of the results from the microarray study. Quantitative Gene Expression (QGE) was performed using MassARRAY methodology and the iPLEX protocol (Sequenom, CA, USA). Total RNA was isolated from liver using an automatic DNA/RNA extractor (BioRobot M48 workstation; Qiagen; CA, USA). Total RNA was treated with TURBO DNA-free™ (Ambion, Huntingdon, UK) for removal of contaminating DNA and first strand cDNA synthesis was conducted on 0.5 μg total RNA using SuperScript™-II Rnase H<sup>- </sup>Reverse Transcriptase (Invitrogen, Carlsbad, CA, USA). Assays for the genes included in this study (See additional file ##SUPPL##6##7##: Gene transcripts included in the rcPCR analyses) were designed and multiplexed into a single reaction using MassARRAY QGE Assay Design software (Sequenom, CA, USA). The competitor, a synthetic DNA molecule matching the targeted cDNA sequence at all positions except for one single base, served as an internal standard for each transcript. A 10-fold competitor dilution was initially used over a wide range of concentrations to determine an approximate equivalence point. Following this, a 7-fold dilution of competitor from 4.04 × 10<sup>-11 </sup>to 1.43 × 10<sup>-19 </sup>was used to achieve accurate quantification. The cDNA and competitor were co-amplified in the same PCR reaction with the following conditions: 95°C for 15 minutes, 45 cycles each of 95°C for 20 second, 56°C for 30 seconds and 72°C for 1 minute, and 72°C for 3 minutes. A clean-up step was performed to remove unincorporated nucleotides. The iPLEX reaction cocktail mix and PCR conditions were according to manufacturer's instructions [##UREF##15##74##]. Parallel PCR reactions were performed for all samples and the products were printed with two replicates on a SpectroCHIP. The primer extension reaction uses PCR products as templates and generates distinct mass signals for competitor and cDNA-derived products. Mass spectrometric analysis generated signals from which peak areas were calculated. Gene expression levels were analysed using TITAN software version 1.0–13 [##REF##15797563##75##,##UREF##16##76##] that runs in the R statistical environment. Raw data from the Genotype Analyzer Software (Sequenom) were imported into TITAN and analysed using the default values of linear least squares polynomial regression and 4000 bootstrap replicates. cDNA concentrations were corrected with respect to the housekeeping gene (<italic>TFRC</italic>), and p-values and confidence intervals for fold changes were calculated.</p>
<title>Results</title> <p>Porcine cDNA microarrays were used to obtain global gene expression profiles in liver tissues sampled from 58 D and 58 NL boars with extreme high and low levels of androstenone. Average levels of androstenone were 11.57 ± 3.2 ppm for D high (DH) and 0.37 ± 0.17 ppm for D low (DL) boars. Average levels of androstenone were 5.95 ± 2.04 ppm for NL high (NLH) and 0.14 ± 0.04 ppm for NL low (NLL) boars. Linear models were used to identify significantly differentially expressed genes and the top 1% (269) most significant genes were subsequently inspected as possible candidate genes affecting androstenone levels (See additional file ##SUPPL##0##1##: Microarray results for Duroc, additional file ##SUPPL##1##2##: Microarray results for Norwegian Landrace). The top 100 genes of D and NL are presented in table ##TAB##0##1## and ##TAB##1##2##, respectively. Among the 269 affected genes, 25 were found to be in common for the two breeds and expression differences were generally more significant in D compared to NL.</p> <p>Testing for overrepresentation of gene ontology (GO) terms among the 269 affected genes was performed relative to the global representation of GO terms on the microarray. A majority of the significant genes were related to oxidoreductase activities, glucuronosyltransferase activity and the binding of iron ion, haem, DNA and oxygen (Figure ##FIG##0##1##). Additionally, cellular metabolic process pathways such as amino acid metabolism and catabolism, electron transport and inflammatory response were overrepresented (Figure ##FIG##1##2##). Furthermore, the genes were classified according to their cellular component ontology (See additional file ##SUPPL##2##3##: Gene ontology (GO) results for the cellular component ontology in Duroc, additional file ##SUPPL##3##4##: Gene ontology (GO) results for the cellular component ontology in Norwegian Landrace). For D, cytoplasm, microsomes and endoplasmatic reticulum were among the important cellular components, none of the cellular components reached the level of significance, i.e. had more than 10 significant genes related to the CC terms.</p> <p>To validate the microarray results, a real competitive PCR (rcPCR) approach was applied to five selected genes: Flavin-containing monooxygenase 1 (<italic>FMO1</italic>), transcription factor GATA-4 (<italic>GATA4</italic>), 17β-hydroxysteroid dehydrogenase 13 (<italic>HSD17B13</italic>), 17β-hydroxysteroid dehydrogenase 2 (<italic>HSD17B2</italic>) and N-acetyltransferase 12 (<italic>NAT12</italic>). Their expression levels were normalised to the housekeeping gene transferrin receptor (<italic>TFRC</italic>). Differential expression of <italic>FMO1</italic>, <italic>HSD17B13</italic>, <italic>HSD17B2 </italic>and <italic>NAT12 </italic>was confirmed in both D and NL breeds (p &lt; 0.05) while altered expression of <italic>GATA4 </italic>was not found in either of the two breeds (Figure ##FIG##2##3##).</p>
<title>Discussion</title> <p>Elevated levels of androstenone in adipose tissue can result from both increased biosynthesis in the testes and deficient metabolism in liver. We have previously conducted a microarray experiment investigating gene expression profiles in the testes of boars with extreme levels of androstenone, and in this study we have examined gene expression profiles in the liver of some representatives from that study along with new individuals. Liver metabolism can be divided into phase I and phase II reactions [##REF##8354165##22##]. Phase I metabolism involves oxidation and hydroxylation reactions that make the substrate more water soluble. Phase II conjugation reactions further increase hydrophilicity by adding polar groups. As a result of these metabolic reactions, compounds including endogenous steroids, fatty acids and drugs are inactivated and eliminated. We have identified a number of differentially expressed genes that function in pathways affecting both phase I and phase II reactions involved in metabolism of androstenone in the liver.</p> <title>Phase I metabolism</title> <p>The most significantly differentially expressed gene involved in phase I oxidation reactions was cytochrome P450 2C49 (<italic>CYP2C49</italic>). This gene is a member of the CYP2C family [##UREF##4##23##] and genes of this family encode monooxygenases that catalyse several reactions involved in the metabolism of drugs and endogenous compounds. Different CYP2C isoforms show some cross reactivity towards substrates which makes it difficult to differentiate CYP2C activities, and substrate specificity for CYP2C49 is not known [##REF##16611125##24##]. In our study, <italic>CYP2C49 </italic>transcripts were shown to have a significant up-regulation in boars with high levels of androstenone in both D and NL lines (DH and NLH) but differences were clearly most prominent in the D breed.</p> <p>Another member of the cytochrome P450 superfamily, that was differentially expressed in our study, is cytochrome P450 2E1 (<italic>CYP2E1</italic>). In contrast to expression of the <italic>CYP2C49 </italic>gene, <italic>CYP2E1 </italic>showed a down-regulation in DH and NLH boars. <italic>CYP2E1 </italic>encodes an enzyme that metabolises many endogenous and exogenous substrates such as alcohols, ketones and drugs [##REF##16611125##24##]. Androstenone has been shown to block skatole induced expression of CYP2E1 in pig hepatocytes and to inhibit activity of CYP2E1 in porcine liver microsomes [##REF##12044562##10##,##REF##17156906##25##]. Liver metabolism of skatole also involves this gene and enzyme [##REF##9535343##13##,##REF##9303469##14##,##UREF##5##26##, ####REF##16167985##27##, ##REF##16501006##28####16501006##28##], and the observed pubertal increase in skatole levels has been attributed to the inhibition of CYP2E1 by the sex steroids androstenone and 17β-estradiol [##REF##17156906##25##]. This is the first study, however, to show gene expression differences for <italic>CYP2E1 </italic>in relation to androstenone.</p> <p>Cytochrome P450 member 2A19 (<italic>CYP2A19</italic>), which is the pig ortholog of human <italic>CYP2A6</italic>, catalyses 7-hydroxylation of coumarin in pigs [##UREF##5##26##], in addition to being involved in skatole metabolism [##REF##15998359##29##]. In this study, <italic>CYP2A19 </italic>was significantly down-regulated in DH boars. Gene expression differences in animals with high and low androstenone levels have not previously been reported for porcine <italic>CYP2A19</italic>, however its protein expression has been shown to be inhibited by androstenone [##REF##17908921##30##]. Diaz and Squires [##UREF##5##26##] found that both CYP2A6 protein content and its enzymatic activity were negatively correlated with skatole levels in adipose tissue. In contrast, results presented by Terner <italic>et al</italic>. [##REF##16501006##28##] indicate that the CYP2A6 enzyme is not important for metabolism of skatole in primary cultured porcine hepatocytes.</p> <p>Two other members of the cytochrome P450 family, <italic>CYP27A1 </italic>and <italic>CYP2C33</italic>, were also up-regulated in DH. CYP27A1 is a multifunctional enzyme that oxygenates cholesterol, bile acids and vitamin D. It has not previously been identified as a candidate gene associated with androstenone, but it is regulated by other sex steroids in human cells [##REF##17482558##31##]. <italic>CYP2C33 </italic>belongs to the same sub-family as <italic>CYP2C49</italic>, but its function has yet to be described in the literature. Like the <italic>CYP2C49 </italic>gene, <italic>CYP2C33 </italic>was found to be up-regulated in DH and NLH boars. The differential expression of <italic>CYP2A19, CYP27A1 </italic>and <italic>CYP2C33 </italic>in D and not NL boars might suggest that more genes are involved in phase I reactions in D. This is supported by a high number of oxidoreductase pathways significant in D (figure ##FIG##0##1##). Monooxygenase reactions are, however, clearly important in both breeds according to our gene ontology results (figure ##FIG##0##1##). Down-regulation of <italic>CYP2E1 </italic>and <italic>CYP2A19 </italic>might suggest a different role for these genes compared to the up-regulated <italic>CYP2C18</italic>, <italic>CYP2C33 </italic>and <italic>CYP27A1</italic>.</p> <p>Our study implicates another class of genes involved in phase I oxidation reactions, namely the flavin-containing monooxygenases (FMOs). The FMO family of enzymes converts lipophilic compounds into more polar metabolites and decreases activity of the compounds, a similar activity to that of the cytochrome P450s [##REF##16112078##32##]. Microarray expression results show that flavin-containing monooxygenase 1 (<italic>FMO1</italic>) was significantly up-regulated in DH in the microarray results. The expression pattern in DH boars was confirmed by rcPCR and the rcPCR analysis also revealed that <italic>FMO1 </italic>was up-regulated in NLH boars. As shown in figure ##FIG##1##2##, the fold change in NLH was only half of DH, which might explain why this gene was not found to be significantly differentially expressed in the microarray results. Regulation of FMO involves sex steroids. In male mice, castration was reported to increase FMO1 expression [##REF##9186481##33##], while in rats the opposite effect has been shown, with positive regulation of FMO by testosterone and negative by estradiol [##REF##3733730##34##]. The monooxygenase activity of FMO1 in addition to its regulation by steroid hormones makes it an interesting candidate gene for boar taint.</p> <title>Phase II metabolism</title> <p>Phase I oxidation reactions of androstenone by enzymes including cytochrome P450s and FMOs are often followed by phase II conjugation reactions that are catalysed by glucuronosyltransferases, sulfotransferases, acetyltransferases and glutathione S-transferases [##REF##8354165##22##]. Through the addition of polar moieties, these enzymes increase substrate solubility. Conjugation reactions are an important means of excreting steroid hormones and other compounds, these reactions have also been proposed to keep inactive steroids easily available in cells. The hydroxysteroid sulfotransferases SULT2A1 and SULT2B1 have been associated with androstenone levels in previous studies [##REF##17686889##12##,##REF##16595701##15##, ####REF##15644508##16##, ##REF##17210083##17####17210083##17##]. These genes were not found to be significantly differentially expressed in our study, however the estrogen sulfotransferase <italic>STE </italic>(<italic>SULT2E1</italic>) was found to be up-regulated in DH boars and down-regulated in NLH boars. A previous microarray study performed by our group showed an up-regulation of <italic>CYP19 </italic>gene expression in the testis of D and NL boars with high androstenone levels [##REF##17210083##17##]. These genes encode enzymes biosynthesising estrogens and might explain why an estrogen conjugation enzyme is differentially expressed in livers of the same pigs. Furthermore, boars secrete large quantities of estrogens from the testis [##REF##7191612##35##], and both fat and plasma estrogen levels have previously been shown to be highly correlated with levels of androstenone in adipose tissue [##REF##16324073##36##,##REF##17609472##37##].</p> <p>Phase II conjugation by glucuronidation is another major pathway of liver elimination of endogenous and exogenous compounds, and in D this molecular function was significantly overrepresented. It is catalysed by uridine diphospho-glucuronosyltransferases (UGTs) which transfer glucuronic acid to substrates to increase solubility [##REF##10419020##38##]. The UGTs have been divided into two subfamilies, the UGT1s and the UGT2s, and in this study the family member <italic>UGT1A5 </italic>was found to be up-regulated in DH boars. Furthermore, a transcript similar to <italic>UGT2B15 </italic>was up-regulated in D and a transcript similar to <italic>UGT2A1 </italic>was down-regulated in NLH. UGT2B15 has been shown to conjugate several androgens in humans [##REF##10419020##38##], whereas UGT1A5 has been found to be catalytically active against some exogenous compounds [##REF##16120810##39##]. UGT2A1 contributes to glucuronidation of steroids and phenolic compounds in olfactory tissue and has been found to eliminate odourants [##REF##10359671##40##]. The UGT family of conjugation enzymes have previously been associated with the androstenone metabolites α-androstenol and β-androstenol [##UREF##6##41##,##REF##4462588##42##].</p> <p>N-acetyltransferases (NATs) are another family of conjugation enzymes and act by transferring the acetyl group of acetyl coenzyme A to aromatic amines to increase their water solubility [##REF##8354165##22##]. NAT type 12 (<italic>NAT12</italic>) was found to be up-regulated in NLH based on the microarray results and this was verified by rcPCR. The rcPCR results also revealed an up-regulation in DH boars. Androgens have been shown to increase expression of NAT type 1 in human cancer cells [##REF##17210686##43##], but the enzyme family has not previously been associated with levels of androstenone. Tryptophan, the precursor of skatole, contains an aromatic, suggesting a role for this family of conjugation enzymes in the control of skatole levels, possibly through regulation by steroid hormones.</p> <p>Glutathione S-transferases (GSTs) are functionally diverse enzymes mostly known to catalyse conjugation reactions of endogenous substances, haem, fatty acids, xenobiotics and products of oxidative processes [##REF##3068032##44##]. They have also been implicated in the intracellular transport of steroids to their site of action [##REF##3068032##44##]. In this study, a glutathione S-transferase gene was found to be down-regulated in NLH. The gene is described in the database as <italic>LOC396850</italic>. We have previously described up-regulation of two GST genes, <italic>GSTO1 </italic>and <italic>MGST1</italic>, in association with testicular androstenone levels [##REF##17210083##17##], and here we additionally propose a role for glucuronidation in liver metabolism. Expression of phase II metabolic genes in D and NL pigs suggest breed specific mechanisms. We found differential expression of a GST variant (<italic>LOC396850</italic>) in NL but not D boars, up-regulation of <italic>STE </italic>in DH but down-regulation in NLH, and up-regulation of two <italic>UGT </italic>genes in DH but down-regulation of another family member in NLH. Breed specific phase II mechanisms for these breeds are in agreement with our previous finding for the protein <italic>SULT2B1 </italic>[##REF##17686889##12##].</p> <title>Regulation of steroid availability</title> <p>17β-hydroxysteroid dehydrogenases (17β-HSDs) regulate the availability of androgens and estrogens in tissues by catalysing interconvertion of active and inactive forms of steroids [##REF##11165032##45##]. They do so by regulating the occupancy of steroid nuclear hormone receptors like androgen receptor, estrogen receptor and progesterone receptor [##REF##15645014##46##]. Gene expression of <italic>17</italic>β-HSD in porcine liver has been shown to be negatively associated with levels of androstenone in adipose tissue [##REF##17609472##37##]. In a previous study, we found that the <italic>HSD17B4 </italic>gene was significantly up-regulated in testicle samples from DH and NLH [##REF##17210083##17##], and in this study <italic>HSD17B4 </italic>was up-regulated in the liver of DH boars. Additionally, we found that <italic>HSD17B2 </italic>was down-regulated in DH and that the isoform <italic>HSD17B13 </italic>was down-regulated in both DH and NLH. Differential expression of both <italic>HSD17B2 </italic>and <italic>HSD17B13 </italic>was verified by rcPCR, which also showed significant down-regulation of <italic>HSD17B2 </italic>in NLH boars. The HSD17B2 enzyme catalyses the interconversion of testosterone and androstenedione, as well as estradiol and estrone [##REF##8099587##47##]. The function of isoform <italic>HSD17B13</italic>, also known as short-chain dehydrogenase/reductase 9, has yet to be described but its expression has been characterised in human liver [##REF##17311113##48##]. The short-chain dehydrogenase/reductase 8 gene (<italic>DHRS8</italic>), also known as HSD17B11, was found to be down-regulated in NLH boars. HSD17B11 is involved in androgen metabolism [##REF##12697717##49##] and we previously found this gene to be down-regulated in the testis of DH boars [##REF##17210083##17##]. The HSDs belong to two superfamilies: the short-chain dehydrogenase/reductases (SDRs) and the aldo-keto reductases (AKRs) [##REF##15645014##46##]. The four HSDs described above all belong to the SDR superfamily but a gene from the AKRs was also found to be differentially expressed. The aldo-keto reductase family member 1D1 (<italic>AKR1D1</italic>) was down-regulated in NLH animals. AKR1D1 is a liver specific enzyme that regulates the hormone levels of several steroids [##REF##15645014##46##]. The enzyme has a 5β-reductase activity and is necessary for the catabolism of aldosterone, cortisol and androgens [##REF##18243262##50##]. Aldosterone, cortisol and androgens are three of four metabolic compounds formed from pregnenolone, with the fourth compound being androstenone [##REF##3516959##5##]. A role for AKR1D1 in androstenone catabolism might therefore be possible. We previously reported up-regulation of the family member <italic>AKR1C4 </italic>in the testis of DH and NLH boars [##REF##17210083##17##], but <italic>AKR1D1 </italic>has not formerly been associated with porcine androstenone levels.</p> <p>Another class of enzymes that regulate the availability of steroids is the plasma proteins. Binding of plasma protein is a reversible reaction that has no physiological effect, but it controls the amount of free drugs and hormones available to tissues by protecting them from metabolism [##REF##11058758##51##]. Alpha-1-acid glycoprotein (AGP) is an acute phase serum protein synthesised in the liver and secreted to plasma where it binds and carries drugs and steroid hormones [##REF##11058758##51##]. The protein has been shown to interact with CYP3A and CYP2C19 in humans and to inhibit cytochrome P450 activity [##REF##12368055##52##,##REF##12597994##53##]. It is suggested that serum protein interaction is an important factor for cytochrome P450 mediated metabolism, and that the interaction is isoform specific [##REF##12597994##53##]. AGP has also been studied in pigs as a binding factor of pheromones, but has not been found to bind progesterone in pig nasal mucosa [##REF##14578123##54##]. AGP has two variants in humans: orosomucoid 1 and 2 (ORM1 and ORM2) [##REF##11058758##51##]. In this study we found that porcine transcript homologous to bovine <italic>AGP </italic>and human <italic>ORM1 </italic>were significantly down-regulated in DH and NLH boars. Significant expression of these genes might explain the inflammatory response ontology term in figure ##FIG##1##2## as these genes have mostly been described in relation to inflammatory response. AGP has not been investigated as a binding protein and transporter for androstenone, but different levels of androstenone could be related to differences in its availability through plasma binding by AGP.</p> <title>Regulatory factors</title> <p>One of the most significant molecular function terms in NL was nucleic acid binding. This term includes transcription factor binding and several transcription factors were found to be differentially expressed in this study. The GATA factor 4 (<italic>GATA4</italic>) was found to be down-regulated based on the microarray data, but we were not able to confirm this result with rcPCR. Other transcription factors identified in this study should be further investigated, including RAR-related orphan receptor A (<italic>RORA</italic>), transcription elongation factor B3 (<italic>TCEB3</italic>), nuclear factor I/X (<italic>NFIX</italic>), transcription factor 8 (<italic>TCF8</italic>) and heat shock transcription factor 1 (<italic>HSF1</italic>). We have previously identified iron ion binding, ferric ion binding and electron transport as being associated with levels of androstenone [##REF##17210083##17##], possibly through interaction with the haem-containing cytochrome P450s. In this study we also found iron ion binding and electron transport, in addition to haem binding, to be important pathways, supporting our previous findings. Cytochrome b5 (CYB5) is involved in electron transfer to cytochrome P450s and has been proposed as a candidate gene for androstenone through its interaction with cytochrome P450 c17 (CYP17) [##REF##8240978##55##]. In this study we found <italic>CYB5 </italic>to be down-regulated in NLH boars, however we previously identified <italic>CYB5 </italic>as being up-regulated in the testis of both DH and NLH [##REF##17210083##17##]. Both <italic>CYP2E1 </italic>and <italic>CYP2A6 </italic>have been found to be activated by <italic>CYB5 </italic>in humans [##REF##15728263##56##], supporting the down-regulation of <italic>CYB5 </italic>together with <italic>CYP2E1 </italic>and <italic>CYP2A6 </italic>in this study. Up-regulation in testis and down-regulation in liver might also support a key regulatory role of <italic>CYB5 </italic>in both tissues.</p> <title>Breed differences</title> <p>Breed differences in boar taint candidate genes have previously been shown between D and NL [##REF##17686889##12##,##REF##17210083##17##]. Consistent with these findings, we found breed differences in the expression profiles of genes involved in both phase I metabolism and phase II metabolism. However, rcPCR results might suggest that genes appearing in one breed can be differentially expressed also in the other breed. Additional and more comprehensive rcPCR studies are planned to clarify this. In general, the D breed showed higher levels of significance compared to the NL breed, potentially due to the larger contrasts between high and low androstenone groups in this breed. A number of new candidate genes have been identified in this study, however, we do not know their effect on performance and sexual maturation in pig. Additional studies are needed to investigate potential roles of the genes identified for phenotypes related to fertility.</p> <title>Relationship with skatole</title> <p>Pigs in this study were selected based on extreme androstenone values. However, levels of androstenone and skatole have previously been found to be highly correlated, showing genetic correlations of 0.62 for D and 0.36 for NL [##UREF##3##21##]. Some of the genes found differentially expressed in this study have previously been found to be associated with tryptophan or skatole and this could be explained by high correlations with androstenone. This might also be the reason for the significance of GO terms such as aromatic amino acid metabolism and catabolism, since the skatole precursor tryptophan belongs to this family of amino acids. The tryptophan 2, 3-dioxygenase gene (<italic>TDO2</italic>), which oxidises tryptophan, was differentially expressed in DH. A gene of the FMO family, kynurenine 3-monooxygenase (<italic>KMO</italic>), was found to be down-regulated in NLH. KMO is involved in tryptophan degradation and might therefore be interesting in regards to skatole [##REF##16973376##57##]. The aldehyde oxidase (<italic>AOX1</italic>) gene was significantly down-regulated in high androstenone animals of both breeds and this gene has not previously been associated with androstenone levels. It has, however, been shown to play an important role in skatole metabolism in several species including pigs [##UREF##7##58##, ####REF##9671536##59##, ##REF##8825199##60####8825199##60##].</p>
<title>Conclusion</title> <p>In this study we compared global gene expression profiles from the livers of boars with extreme high and low levels of androstenone from two breeds, Duroc and Norwegian Landrace. Breed differences are evident for molecular functions and biological processes involved in metabolism of androstenone, however many of the same genes are differentially expressed in the two breeds as well. Genes encoding different oxidising enzymes including the cytochrome P450 family (<italic>CYP2E1</italic>, <italic>CYP2A19</italic>, <italic>CYP2C49</italic>, <italic>CYP27A1 </italic>and <italic>CYP2C33</italic>) and the flavin-containing monooxygenase family (<italic>FMO1 </italic>and <italic>KMO</italic>) were significantly differentially expressed. Furthermore, genes involved in conjugation reactions, including the UDP-glucuronosyltransferases (<italic>UGT1A5</italic>, <italic>UGT2A1 </italic>and <italic>UGT2B15</italic>), sulfotransferases (<italic>STE</italic>), N-acetyltransferases (<italic>NAT12</italic>) and glutathione S-transferase were significant, in addition to genes of the 17β-hydroxysteroid dehydrogenase family (<italic>HSD17B2</italic>, <italic>HSD17B4</italic>, <italic>HSD17B11 </italic>and <italic>HSD17B13</italic>), which are known to regulate availability of active steroids. We suggest a novel role for plasma proteins including <italic>AGP </italic>and <italic>ORM1 </italic>in regulating availability of androstenone in pigs. This is the first published microarray experiment describing liver metabolism of androstenone. A number of new candidate genes have been identified, both from phase I and phase II metabolism as well as pathways regulating steroid availability.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>Boar taint is the unpleasant odour and flavour of the meat of uncastrated male pigs that is primarily caused by high levels of androstenone and skatole in adipose tissue. Androstenone is a steroid and its levels are mainly genetically determined. Studies on androstenone metabolism have, however, focused on a limited number of genes. Identification of additional genes influencing levels of androstenone may facilitate implementation of marker assisted breeding practices. In this study, microarrays were used to identify differentially expressed genes and pathways related to androstenone metabolism in the liver from boars with extreme levels of androstenone in adipose tissue.</p> <title>Results</title> <p>Liver tissue samples from 58 boars of the two breeds Duroc and Norwegian Landrace, 29 with extreme high and 29 with extreme low levels of androstenone, were selected from more than 2500 individuals. The samples were hybridised to porcine cDNA microarrays and the 1% most significant differentially expressed genes were considered significant. Among the differentially expressed genes were metabolic phase I related genes belonging to the cytochrome P450 family and the flavin-containing monooxygenase <italic>FMO1</italic>. Additionally, phase II conjugation genes including UDP-glucuronosyltransferases <italic>UGT1A5</italic>, <italic>UGT2A1 </italic>and <italic>UGT2B15</italic>, sulfotransferase <italic>STE</italic>, N-acetyltransferase <italic>NAT12 </italic>and glutathione S-transferase were identified. Phase I and phase II metabolic reactions increase the water solubility of steroids and play a key role in their elimination. Differential expression was also found for genes encoding 17beta-hydroxysteroid dehydrogenases (<italic>HSD17B2</italic>, <italic>HSD17B4</italic>, <italic>HSD17B11 </italic>and <italic>HSD17B13</italic>) and plasma proteins alpha-1-acid glycoprotein (<italic>AGP</italic>) and orosomucoid (<italic>ORM1</italic>). 17beta-hydroxysteroid dehydrogenases and plasma proteins regulate the availability of steroids by controlling the amount of active steroids accessible to receptors and available for metabolism. Differences in the expression of <italic>FMO1</italic>, <italic>NAT12</italic>, <italic>HSD17B2 </italic>and <italic>HSD17B13 </italic>were verified by quantitative real competitive PCR.</p> <title>Conclusion</title> <p>A number of genes and pathways related to metabolism of androstenone in liver were identified, including new candidate genes involved in phase I oxidation metabolism, phase II conjugation metabolism, and regulation of steroid availability. The study is a first step towards a deeper understanding of enzymes and regulators involved in pathways of androstenone metabolism and may ultimately lead to the discovery of markers to reduce boar taint.</p>
<title>List of abbreviations</title> <p>D: Duroc; NL: Norwegian Landrace; DH: Duroc high androstenone; NLH: Norwegian Landrace high androstenone; DL: Duroc low androstenone; NLL: Norwegian Landrace high androstenone; 17β-HSD: 17β-hydroxysteroid dehydrogenase; 3β-HSD: 3β-hydroxysteroid dehydrogenase; AGP: alpha-1-acid glycoprotein; AKR: aldo-keto reductase; AKR1C4: aldo-keto reductase family member 1C4; AKR1D1: aldo-keto reductase family member 1D1; AOX1: aldehyde oxidase; CYB5: cytochrome b5; CYP17: cytochrome P450 family 17; CYP19: cytochrome P450 family 19; CYP27A1: cytochrome P450 family member 27A1; CYP2A: cytochrome P450 family 2A; CYP2A19: cytochrome P450 family member 2A19; CYP2A6: cytochrome P450 family member 2A6; CYP2C19: cytochrome P450 family member 2C19; CYP2C33: cytochrome P450 family member 2C33; CYP2C49: cytochrome P450 family member 2C49; CYP2E1: cytochrome P450 family member 2E1; CYP3A: cytochrome P450 family member 3A; DHRS8: short-chain dehydrogenase/reductase 8; FMO: flavin-containing monooxygenase; FMO1: flavin-containing monooxygenase; GATA4: transcription factor GATA-4; GO: gene ontology; GST: glutathione S-transferase; GSTO1: glutathione S-transferase omega; HSDs: hydroxysteroid dehydrogenases; HSD17B2: hydroxysteroid 17β dehydrogenase 2; HSD17B4: hydroxysteroid 17β dehydrogenase 4; HSD17B11: hydroxysteroid 17β dehydrogenase 11; HSD17B13: hydroxysteroid 17β dehydrogenase 13; HSF1: heat shock transcription factor 1; KMO: kynurenine 3-monooxygenase; MGST1: glutathione S-transferase; NAT: N-acetyltransferase; NAT12: N-acetyltransferase family member 12; NFIX: nuclear factor I/X; ORM1: orosomucoid 1; ORM2: orosomucoid 2; QGE: quantitative gene expression; RORA: RAR-related orphan receptor A; rcPCR: real competitive PCR; SDR: short-chain dehydrogenase/reductase; STE: estrogen sulfotransferase; SULT1A1: sulfotransferase family member 1A1; SULT2E1: sulfotransferase family member 1E1; SULT2A1: sulfotransferase family member 2A1; SULT2B1: sulfotransferase family member 2B1; TCEB3: transcription elongation factor B3; TCF8: transcription factor 8; TDO2: tryptophan 2,3-dioxygenase; TFRC: transferrin receptor; UGT: uridine diphospho-glucuronosyltransferase; UGT1A5: UDP-glucuronosyltransferase family member 1A5; UGT2A1: UDP-glucuronosyltransferase family member 2A1; UGT2B15: UDP-glucuronosyltransferase family member 2B15</p> <title>Authors' contributions</title> <p>MM carried out the microarray experimental work, performed statistical analysis and drafted the paper. SL was involved in planning the project, provided laboratory facilities for rcPCR work and took part in writing the paper. CB was involved in planning the project and was in charge of lab facilities for microarray studies. JH was involved in microarray experimental work and took part in writing the paper. HH carried out bioinformatics work. IB was involved in statistical analysis. TM was involved in power calculations and statistical supervision. EG coordinated the study, was involved in planning the project, carried out rcPCR molecular work and took part in writing the paper.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>This study was financed by the Norwegian pig breeders association (NORSVIN) and The Research Council of Norway. The work conducted at University of Aarhus was supported by the EU-project SABRE. We want to thank Ellen Dahl and Øystein Andresen for being in charge of the androstenone analyses performed at the hormone laboratory at the Norwegian School of Veterinary Science (NVH). We also thank Elin Bergseth at NORSVIN for collecting samples, Dr. Peter Sørensen and Helle Jensen at the microarray platform at University of Aarhus for kind assistance, Dr. Paul Oeth at Sequenom for carrying out the rcPCR primer design and Dr. Matthew Baranski for valuable help with English grammar.</p>
<fig position="float" id="F1"><label>Figure 1</label></fig> <fig position="float" id="F2"><label>Figure 2</label></fig> <fig position="float" id="F3"><label>Figure 3</label></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Top 100 genes identified in Duroc.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">ID</td><td align="left">Name</td><td align="left">gene_id</td><td align="left">FoldChange</td><td align="left">adj.P.Val</td></tr></thead><tbody><tr><td align="right">103725</td><td align="left">Cytochrome P450 2C49 (CYP2C49)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_214420.1">NM_214420.1</ext-link></td><td align="left">1.227</td><td align="left">8.26e-24</td></tr><tr><td align="right">104194</td><td align="left">Cytochrome P450 2C49 (CYP2C49)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_214420.1">NM_214420.1</ext-link></td><td align="left">0.853</td><td align="left">4.98e-20</td></tr><tr><td align="right">102872</td><td align="left">Glycine-N-acyltransferase (GLYAT)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_201648.1">NM_201648.1</ext-link></td><td align="left">1.113</td><td align="left">1.27e-19</td></tr><tr><td align="right">103265</td><td align="left">Cytochrome P450 2C49 (CYP2C49)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_214420.1">NM_214420.1</ext-link></td><td align="left">1.480</td><td align="left">4.76e-18</td></tr><tr><td align="right">103234</td><td align="left">Hepatic flavin-containing monooxygenase (FMO) (FMO1)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_214064.1">NM_214064.1</ext-link></td><td align="left">1.179</td><td align="left">2.28e-17</td></tr><tr><td align="right">206700</td><td align="left">Cytochrome P450 2C49 (CYP2C49)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_214420.1">NM_214420.1</ext-link></td><td align="left">1.099</td><td align="left">2.71e-16</td></tr><tr><td align="right">102682</td><td align="left">Serum amyloid P component (APCS)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_213887.1">NM_213887.1</ext-link></td><td align="left">1.294</td><td align="left">2.71e-16</td></tr><tr><td align="right">102284</td><td align="left">Sterile alpha motif domain containing 8 (Samd8)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_026283.2">NM_026283.2</ext-link></td><td align="left">0.969</td><td align="left">3.44e-16</td></tr><tr><td align="right">102668</td><td align="left">Hepatic flavin-containing monooxygenase (FMO) (FMO1)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_214064.1">NM_214064.1</ext-link></td><td align="left">1.271</td><td align="left">8.01e-16</td></tr><tr><td align="right">100748</td><td align="left">Regulatory factor X domain containing 2 (RFXDC2)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_022841.3">NM_022841.3</ext-link></td><td align="left">0.594</td><td align="left">7.47e-14</td></tr><tr><td align="right">206864</td><td align="left">Similar to Alcohol dehydrogenase 6 (LOC523510)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="XM_601810.2">XM_601810.2</ext-link></td><td align="left">-0.925</td><td align="left">3.49e-13</td></tr><tr><td align="right">102390</td><td align="left">Cytochrome P450 2C49 (CYP2C49)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_214420.1">NM_214420.1</ext-link></td><td align="left">1.152</td><td align="left">5.05e-13</td></tr><tr><td align="right">102098</td><td align="left">Carnitine O-octanoyltransferase (CROT)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_021151.2">NM_021151.2</ext-link></td><td align="left">0.886</td><td align="left">5.30e-13</td></tr><tr><td align="right">102991</td><td align="left">Cytochrome P450. family 2. subfamily E. polypeptide 1 (CYP2E1)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_214421.1">NM_214421.1</ext-link></td><td align="left">-1.040</td><td align="left">2.17e-12</td></tr><tr><td align="right">216026</td><td align="right">216026</td><td/><td align="left">-0.728</td><td align="left">3.72e-12</td></tr><tr><td align="right">102413</td><td align="left">Small nuclear RNA activating complex. polypeptide 1. 43 kDa (SNAPC1)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_003082.2">NM_003082.2</ext-link></td><td align="left">0.994</td><td align="left">4.84e-12</td></tr><tr><td align="right">103432</td><td align="right">103432</td><td/><td align="left">0.527</td><td align="left">2.07e-11</td></tr><tr><td align="right">210250</td><td align="left">Metallothionein (MT1A)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001001266.1">NM_001001266.1</ext-link></td><td align="left">-0.965</td><td align="left">3.09e-11</td></tr><tr><td align="right">103062</td><td align="left">FKBP1A-like (LOC654323)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001038000.1">NM_001038000.1</ext-link></td><td align="left">-1.078</td><td align="left">3.09e-11</td></tr><tr><td align="right">103279</td><td align="left">Orosomucoid 1 (ORM1)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_000607.1">NM_000607.1</ext-link></td><td align="left">-0.947</td><td align="left">3.99e-11</td></tr><tr><td align="right">103021</td><td align="left">Hydroxysteroid (17-beta) dehydrogenase 13 (HSD17B13)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_178135.2">NM_178135.2</ext-link></td><td align="left">-0.869</td><td align="left">6.10e-11</td></tr><tr><td align="right">201488</td><td align="right">201488</td><td/><td align="left">0.462</td><td align="left">1.08e-10</td></tr><tr><td align="right">100575</td><td align="left">17beta-estradiol dehydrogenase (HSD17B4)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_214306.1">NM_214306.1</ext-link></td><td align="left">0.452</td><td align="left">1.24e-10</td></tr><tr><td align="right">100723</td><td align="left">Aspartate aminotransferase (GOT2)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_213928.1">NM_213928.1</ext-link></td><td align="left">0.493</td><td align="left">3.71e-10</td></tr><tr><td align="right">103194</td><td align="left">Alpha-1 acid glycoprotein (AGP)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="XM_585370.2">XM_585370.2</ext-link></td><td align="left">-0.932</td><td align="left">4.52e-10</td></tr><tr><td align="right">202082</td><td align="left">Family with sequence similarity 10. member A4 (FAM10A4)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NR_002183.1">NR_002183.1</ext-link></td><td align="left">0.424</td><td align="left">6.45e-10</td></tr><tr><td align="right">100724</td><td align="left">Phytanoyl-CoA hydroxylase (PHYH)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_006214.3">NM_006214.3</ext-link></td><td align="left">0.518</td><td align="left">8.69e-10</td></tr><tr><td align="right">210077</td><td align="left">Cytochrome P450 2C49 (CYP2C49)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_214420.1">NM_214420.1</ext-link></td><td align="left">0.872</td><td align="left">9.95e-10</td></tr><tr><td align="right">101277</td><td align="left">Estrogen sulfotransferase (STE)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_213992.1">NM_213992.1</ext-link></td><td align="left">0.557</td><td align="left">1.45e-09</td></tr><tr><td align="right">216031</td><td align="left">Insulin-like-growth factor 2 (IGF2)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_213883.1">NM_213883.1</ext-link></td><td align="left">-0.687</td><td align="left">1.88e-09</td></tr><tr><td align="right">210204</td><td align="left">Similar to HLA class I histocompatibility antigen, A-11 alpha chain precursor (LOC642049)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="XM_936199.1">XM_936199.1</ext-link></td><td align="left">-0.591</td><td align="left">1.88e-09</td></tr><tr><td align="right">103031</td><td align="left">Type I iodothyronine deiodinase (DIO1)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001001627.1">NM_001001627.1</ext-link></td><td align="left">0.499</td><td align="left">2.85e-09</td></tr><tr><td align="right">210252</td><td align="left">Sarcosine dehydrogenase (SARDH)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_007101.2">NM_007101.2</ext-link></td><td align="left">-0.273</td><td align="left">3.55e-09</td></tr><tr><td align="right">102976</td><td align="left">Similar to UDP-glucuronosyltransferase 2B15 precursor (UDPGT) (HLUG4) (LOC653180)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="XM_931558.1">XM_931558.1</ext-link></td><td align="left">0.623</td><td align="left">4.69e-09</td></tr><tr><td align="right">201274</td><td align="left">Sorbitol dehydrogenase (SORD). mRNA</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_003104.3">NM_003104.3</ext-link></td><td align="left">0.523</td><td align="left">5.70e-09</td></tr><tr><td align="right">103801</td><td align="left">Similar to Sorbitol dehydrogenase (L-iditol 2-dehydrogenase) (LOC650043)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="XM_939131.1">XM_939131.1</ext-link></td><td align="left">0.452</td><td align="left">1.83e-08</td></tr><tr><td align="right">204461</td><td align="left">Glycine-N-acyltransferase (GLYAT). nuclear gene encoding mitochondrial protein</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_005838.2">NM_005838.2</ext-link></td><td align="left">0.370</td><td align="left">2.47e-08</td></tr><tr><td align="right">102538</td><td align="right">102538</td><td/><td align="left">0.746</td><td align="left">2.65e-08</td></tr><tr><td align="right">204863</td><td align="left">Hydroxysteroid (17-beta) dehydrogenase 13 (HSD17B13)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_178135.2">NM_178135.2</ext-link></td><td align="left">-0.630</td><td align="left">3.89e-08</td></tr><tr><td align="right">101843</td><td align="left">Abhydrolase domain containing 3 (ABHD3)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_138340.3">NM_138340.3</ext-link></td><td align="left">0.370</td><td align="left">9.87e-08</td></tr><tr><td align="right">215880</td><td align="left">Transcription elongation factor B (SIII). polypeptide 3 (110 kDa. elongin A) (TCEB3)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_003198.1">NM_003198.1</ext-link></td><td align="left">-0.612</td><td align="left">9.87e-08</td></tr><tr><td align="right">200322</td><td align="left">Taxilin beta (TXLNB)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_153235.2">NM_153235.2</ext-link></td><td align="left">0.385</td><td align="left">9.87e-08</td></tr><tr><td align="right">210206</td><td align="left">Metallothionein (MT1A)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001001266.1">NM_001001266.1</ext-link></td><td align="left">-0.536</td><td align="left">1.25e-07</td></tr><tr><td align="right">201503</td><td align="left">Carboxylesterase (CES3)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_214246.1">NM_214246.1</ext-link></td><td align="left">0.472</td><td align="left">1.44e-07</td></tr><tr><td align="right">202070</td><td align="right">202070</td><td/><td align="left">0.485</td><td align="left">1.57e-07</td></tr><tr><td align="right">102106</td><td align="right">102106</td><td/><td align="left">0.335</td><td align="left">1.86e-07</td></tr><tr><td align="right">206300</td><td align="left">Protein tyrosine phosphatase type IVA. member 1 (PTP4A1)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_003463.3">NM_003463.3</ext-link></td><td align="left">0.337</td><td align="left">2.01e-07</td></tr><tr><td align="right">103047</td><td align="left">Similar to mouse 2310016A09Rik gene (LOC134147)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_138809.3">NM_138809.3</ext-link></td><td align="left">0.456</td><td align="left">2.01e-07</td></tr><tr><td align="right">204653</td><td align="left">FKBP1A-like (LOC654323)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001038000.1">NM_001038000.1</ext-link></td><td align="left">0.256</td><td align="left">2.06e-07</td></tr><tr><td align="right">103327</td><td align="left">UDP glucuronosyltransferase 1 family, polypeptide A5 (UGT1A5)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_019078.1">NM_019078.1</ext-link></td><td align="left">0.264</td><td align="left">2.42e-07</td></tr><tr><td align="right">102469</td><td align="left">Chromosome 1 open reading frame 128 (C1orf128)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_020362.2">NM_020362.2</ext-link></td><td align="left">0.190</td><td align="left">2.42e-07</td></tr><tr><td align="right">201919</td><td align="left">FKBP1A-like (LOC654323)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001038000.1">NM_001038000.1</ext-link></td><td align="left">0.303</td><td align="left">2.53e-07</td></tr><tr><td align="right">103349</td><td align="left">Kininogen 1 (KNG1)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_000893.2">NM_000893.2</ext-link></td><td align="left">0.428</td><td align="left">3.03e-07</td></tr><tr><td align="right">103940</td><td align="left">UDP glucuronosyltransferase 1 family, polypeptide A5 (UGT1A5)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_019078.1">NM_019078.1</ext-link></td><td align="left">0.191</td><td align="left">4.03e-07</td></tr><tr><td align="right">204455</td><td align="left">Ornithine carbamoyltransferase (OTC)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_000531.3">NM_000531.3</ext-link></td><td align="left">0.357</td><td align="left">4.61e-07</td></tr><tr><td align="right">212598</td><td align="left">L1 cell adhesion molecule (L1CAM)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_024003.1">NM_024003.1</ext-link></td><td align="left">-0.167</td><td align="left">4.72e-07</td></tr><tr><td align="right">207639</td><td align="right">207639</td><td/><td align="left">0.383</td><td align="left">4.83e-07</td></tr><tr><td align="right">200534</td><td align="right">200534</td><td/><td align="left">0.388</td><td align="left">4.89e-07</td></tr><tr><td align="right">102183</td><td align="left">Protein tyrosine phosphatase type IVA. member 1 (PTP4A1)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_003463.3">NM_003463.3</ext-link></td><td align="left">0.360</td><td align="left">5.72e-07</td></tr><tr><td align="right">206424</td><td align="left">Similar to transcription elongation factor B polypeptide 3 binding protein 1 isoform 1 (LOC533427)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="XM_612827.2">XM_612827.2</ext-link></td><td align="left">0.265</td><td align="left">6.41e-07</td></tr><tr><td align="right">103306</td><td align="left">Tryptophan 2.3-dioxygenase (TDO2)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_005651.1">NM_005651.1</ext-link></td><td align="left">0.443</td><td align="left">6.41e-07</td></tr><tr><td align="right">104125</td><td align="left">Sorbitol dehydrogenase (SORD)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_003104.3">NM_003104.3</ext-link></td><td align="left">0.428</td><td align="left">7.73e-07</td></tr><tr><td align="right">102709</td><td align="right">102709</td><td/><td align="left">0.444</td><td align="left">7.76e-07</td></tr><tr><td align="right">100472</td><td align="left">S-adenosylhomocysteine hydrolase (AHCY)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001011727.1">NM_001011727.1</ext-link></td><td align="left">0.321</td><td align="left">8.64e-07</td></tr><tr><td align="right">206129</td><td align="right">206129</td><td/><td align="left">-0.484</td><td align="left">8.86e-07</td></tr><tr><td align="right">102624</td><td align="left">Ornithine carbamoyltransferase (OTC)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_000531.3">NM_000531.3</ext-link></td><td align="left">0.393</td><td align="left">8.86e-07</td></tr><tr><td align="right">213735</td><td align="left">Metallothionein (MT1A)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001001266.1">NM_001001266.1</ext-link></td><td align="left">0.422</td><td align="left">1.01e-06</td></tr><tr><td align="right">102836</td><td align="left">Carboxylesterase (CES3)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_214246.1">NM_214246.1</ext-link></td><td align="left">0.563</td><td align="left">1.10e-06</td></tr><tr><td align="right">201162</td><td align="left">FKBP1A-like (LOC654323)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001038000.1">NM_001038000.1</ext-link></td><td align="left">0.268</td><td align="left">1.42e-06</td></tr><tr><td align="right">201961</td><td align="right">201961</td><td/><td align="left">0.293</td><td align="left">1.42e-06</td></tr><tr><td align="right">216074</td><td align="left">Metallothionein (MT1A)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001001266.1">NM_001001266.1</ext-link></td><td align="left">-0.772</td><td align="left">1.64e-06</td></tr><tr><td align="right">207278</td><td align="right">207278</td><td/><td align="left">0.282</td><td align="left">1.67e-06</td></tr><tr><td align="right">101716</td><td align="right">101716</td><td/><td align="left">0.328</td><td align="left">1.75e-06</td></tr><tr><td align="right">204421</td><td align="left">Meiosis-specific nuclear structural 1 (MNS1)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_018365.1">NM_018365.1</ext-link></td><td align="left">0.364</td><td align="left">1.99e-06</td></tr><tr><td align="right">103156</td><td align="right">103156</td><td/><td align="left">0.333</td><td align="left">2.05e-06</td></tr><tr><td align="right">103376</td><td align="left">Hydroxysteroid (17-beta) dehydrogenase 13 (HSD17B13)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_178135.2">NM_178135.2</ext-link></td><td align="left">-0.497</td><td align="left">2.38e-06</td></tr><tr><td align="right">100064</td><td align="left">Similar to UDP-galactoseN-acetylgalactosamine-alpha-R beta 1,3-galactosyltransferase (LOC645551)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="XM_928571.1">XM_928571.1</ext-link></td><td align="left">0.159</td><td align="left">2.43e-06</td></tr><tr><td align="right">204950</td><td align="left">Abhydrolase domain containing 3 (ABHD3)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_138340.3">NM_138340.3</ext-link></td><td align="left">0.236</td><td align="left">2.43e-06</td></tr><tr><td align="right">204048</td><td align="left">Suppression of tumorigenicity 13 (colon carcinoma) (Hsp70 interacting protein) (ST13)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_003932.3">NM_003932.3</ext-link></td><td align="left">0.379</td><td align="left">2.74e-06</td></tr><tr><td align="right">204867</td><td align="right">204867</td><td/><td align="left">0.308</td><td align="left">2.95e-06</td></tr><tr><td align="right">207963</td><td align="right">207963</td><td/><td align="left">0.254</td><td align="left">3.08e-06</td></tr><tr><td align="right">200952</td><td align="right">200952</td><td/><td align="left">0.236</td><td align="left">3.08e-06</td></tr><tr><td align="right">103294</td><td align="left">Similar to Mob4B protein (MGC124888)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001033891.1">NM_001033891.1</ext-link></td><td align="left">0.369</td><td align="left">3.79e-06</td></tr><tr><td align="right">208427</td><td align="left">RAR-related orphan receptor A (RORA)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_134262.1">NM_134262.1</ext-link></td><td align="left">0.246</td><td align="left">4.25e-06</td></tr><tr><td align="right">208058</td><td align="left">FKBP1A-like (LOC654323)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001038000.1">NM_001038000.1</ext-link></td><td align="left">0.247</td><td align="left">4.48e-06</td></tr><tr><td align="right">205360</td><td align="left">FLJ20105 protein (FLJ20105)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_017669.2">NM_017669.2</ext-link></td><td align="left">0.278</td><td align="left">4.73e-06</td></tr><tr><td align="right">207838</td><td align="left">3-hydroxymethyl-3-methylglutaryl-Coenzyme A lyase (hydroxymethylglutaricaciduria) (HMGCL)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_000191.2">NM_000191.2</ext-link></td><td align="left">0.156</td><td align="left">4.78e-06</td></tr><tr><td align="right">101413</td><td align="left">Replication protein A3. 14 kDa (RPA3)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_002947.3">NM_002947.3</ext-link></td><td align="left">0.253</td><td align="left">5.25e-06</td></tr><tr><td align="right">103180</td><td align="left">T-complex 1 (TCP1)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001008897.1">NM_001008897.1</ext-link></td><td align="left">0.310</td><td align="left">5.69e-06</td></tr><tr><td align="right">202298</td><td align="right">202298</td><td/><td align="left">0.262</td><td align="left">5.73e-06</td></tr><tr><td align="right">100384</td><td align="left">Alpha-methylacyl-CoA racemase (AMACR)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_014324.4">NM_014324.4</ext-link></td><td align="left">0.394</td><td align="left">6.68e-06</td></tr><tr><td align="right">103434</td><td align="left">Acyl-Coenzyme A oxidase 1, palmitoyl (ACOX1)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_004035.5">NM_004035.5</ext-link></td><td align="left">0.368</td><td align="left">7.11e-06</td></tr><tr><td align="right">202477</td><td align="left">Hypothetical LOC505518 (LOC505518)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="XM_581816.2">XM_581816.2</ext-link></td><td align="left">0.240</td><td align="left">7.42e-06</td></tr><tr><td align="right">102858</td><td align="right">102858</td><td/><td align="left">0.587</td><td align="left">7.52e-06</td></tr><tr><td align="right">103058</td><td align="right">103058</td><td/><td align="left">0.381</td><td align="left">7.67e-06</td></tr><tr><td align="right">211612</td><td align="left">Gamma-aminobutyric acid (GABA) A receptor. alpha 3 (GABRA3)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_000808.2">NM_000808.2</ext-link></td><td align="left">-0.135</td><td align="left">8.97e-06</td></tr><tr><td align="right">201347</td><td align="left">HRAS-like suppressor 3 (HRASLS3)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_007069.1">NM_007069.1</ext-link></td><td align="left">0.352</td><td align="left">8.97e-06</td></tr><tr><td align="right">102840</td><td align="left">Sorbitol dehydrogenase (SORD)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_003104.3">NM_003104.3</ext-link></td><td align="left">0.270</td><td align="left">8.97e-06</td></tr><tr><td align="right">201105</td><td align="left">COP9 constitutive photomorphogenic homolog subunit 2 (Arabidopsis) (COPS2)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_004236.1">NM_004236.1</ext-link></td><td align="left">0.242</td><td align="left">9.06e-06</td></tr><tr><td align="right">201848</td><td align="left">Kinesin family member C1 (KIFC1)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_002263.2">NM_002263.2</ext-link></td><td align="left">0.274</td><td align="left">9.06e-06</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Top 100 genes identified in Norwegian Landrace.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">ID</td><td align="left">Name</td><td align="left">gene_id</td><td align="left">FoldChange</td><td align="left">adj.P.Val</td></tr></thead><tbody><tr><td align="right">206864</td><td align="left">Similar to Alcohol dehydrogenase 6 (LOC523510)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="XM_601810.2">XM_601810.2</ext-link></td><td align="left">-0.730</td><td align="left">1.07e-13</td></tr><tr><td align="right">204853</td><td align="right">204853</td><td/><td align="left">-0.441</td><td align="left">2.39e-13</td></tr><tr><td align="right">209551</td><td align="left">Isochorismatase domain containing 1 (ISOC1)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_016048.1">NM_016048.1</ext-link></td><td align="left">-0.403</td><td align="left">3.28e-13</td></tr><tr><td align="right">103021</td><td align="left">Hydroxysteroid (17-beta) dehydrogenase 13 (HSD17B13)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_178135.2">NM_178135.2</ext-link></td><td align="left">-0.900</td><td align="left">3.28e-13</td></tr><tr><td align="right">204863</td><td align="left">Hydroxysteroid (17-beta) dehydrogenase 13 (HSD17B13)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_178135.2">NM_178135.2</ext-link></td><td align="left">-0.730</td><td align="left">3.78e-13</td></tr><tr><td align="right">209525</td><td align="left">Zinc finger protein 471 (ZNF471)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_020813.1">NM_020813.1</ext-link></td><td align="left">-0.335</td><td align="left">1.69e-12</td></tr><tr><td align="right">102789</td><td align="right">102789</td><td/><td align="left">-0.364</td><td align="left">1.69e-12</td></tr><tr><td align="right">204725</td><td align="right">204725</td><td/><td align="left">-0.367</td><td align="left">1.69e-12</td></tr><tr><td align="right">103376</td><td align="left">Hydroxysteroid (17-beta) dehydrogenase 13 (HSD17B13)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_178135.2">NM_178135.2</ext-link></td><td align="left">-0.589</td><td align="left">2.41e-12</td></tr><tr><td align="right">205190</td><td align="left">Similar to heterogeneous nuclear ribonucleoprotein C isoform b (LOC654074)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="XM_945342.1">XM_945342.1</ext-link></td><td align="left">0.249</td><td align="left">7.01e-12</td></tr><tr><td align="right">204302</td><td align="left">Aldehyde oxidase 1 (AOX1)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001159.3">NM_001159.3</ext-link></td><td align="left">-0.481</td><td align="left">9.60e-12</td></tr><tr><td align="right">204466</td><td align="right">204466</td><td/><td align="left">-0.345</td><td align="left">1.20e-11</td></tr><tr><td align="right">203752</td><td align="right">203752</td><td/><td align="left">-0.274</td><td align="left">1.52e-11</td></tr><tr><td align="right">206901</td><td align="left">Non-histone protein HMG1 (LOC445521)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001004034.1">NM_001004034.1</ext-link></td><td align="left">-0.286</td><td align="left">5.34e-11</td></tr><tr><td align="right">218781</td><td align="left">KIAA1344 (KIAA1344)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_020784.1">NM_020784.1</ext-link></td><td align="left">0.280</td><td align="left">5.34e-11</td></tr><tr><td align="right">216031</td><td align="left">Insulin-like-growth factor 2 (IGF2)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_213883.1">NM_213883.1</ext-link></td><td align="left">-0.847</td><td align="left">5.34e-11</td></tr><tr><td align="right">213651</td><td align="right">213651</td><td/><td align="left">-0.437</td><td align="left">9.01e-11</td></tr><tr><td align="right">215880</td><td align="left">Transcription elongation factor B (SIII). polypeptide 3 (110 kDa. elongin A) (TCEB3)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_003198.1">NM_003198.1</ext-link></td><td align="left">-0.797</td><td align="left">9.45e-11</td></tr><tr><td align="right">204619</td><td align="left">Cytochrome b-5 (CYB5)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001001770.1">NM_001001770.1</ext-link></td><td align="left">-0.458</td><td align="left">9.45e-11</td></tr><tr><td align="right">204520</td><td align="right">204520</td><td/><td align="left">-0.290</td><td align="left">1.45e-10</td></tr><tr><td align="right">103177</td><td align="left">Cytochrome b-5 (CYB5)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001001770.1">NM_001001770.1</ext-link></td><td align="left">-0.489</td><td align="left">2.06e-10</td></tr><tr><td align="right">102774</td><td align="left">Non-histone protein HMG1 (LOC445521)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001004034.1">NM_001004034.1</ext-link></td><td align="left">-0.302</td><td align="left">3.48e-10</td></tr><tr><td align="right">204081</td><td align="right">204081</td><td/><td align="left">-0.251</td><td align="left">3.48e-10</td></tr><tr><td align="right">100528</td><td align="left">ATP synthase. H+ transporting. mitochondrial F0 complex. subunit F6 (ATP5J)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001003703.1">NM_001003703.1</ext-link></td><td align="left">0.250</td><td align="left">3.55e-10</td></tr><tr><td align="right">211561</td><td align="left">KIAA0427 (KIAA0427)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_014772.1">NM_014772.1</ext-link></td><td align="left">0.263</td><td align="left">4.63e-10</td></tr><tr><td align="right">209674</td><td align="right">209674</td><td/><td align="left">0.269</td><td align="left">7.12e-10</td></tr><tr><td align="right">208386</td><td align="left">Dual specificity phosphatase 12 (DUSP12)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_007240.1">NM_007240.1</ext-link></td><td align="left">-0.318</td><td align="left">7.12e-10</td></tr><tr><td align="right">207228</td><td align="left">Small nuclear ribonucleoprotein polypeptide E (SNRPE)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_003094.2">NM_003094.2</ext-link></td><td align="left">0.252</td><td align="left">7.37e-10</td></tr><tr><td align="right">219297</td><td align="left">Purine-rich element binding protein A (PURA)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_005859.3">NM_005859.3</ext-link></td><td align="left">0.499</td><td align="left">8.70e-10</td></tr><tr><td align="right">218205</td><td align="left">SMC1 structural maintenance of chromosomes 1-like 1 (yeast) (SMC1L1)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_006306.2">NM_006306.2</ext-link></td><td align="left">0.300</td><td align="left">1.11e-09</td></tr><tr><td align="right">103194</td><td align="left">Alpha-1 acid glycoprotein (AGP)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="XM_585370.2">XM_585370.2</ext-link></td><td align="left">-1.071</td><td align="left">1.16e-09</td></tr><tr><td align="right">219538</td><td align="right">219538</td><td/><td align="left">0.352</td><td align="left">1.67e-09</td></tr><tr><td align="right">201622</td><td align="left">Hypothetical protein LOC285016 (LOC285016)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001002919.1">NM_001002919.1</ext-link></td><td align="left">-0.417</td><td align="left">1.99e-09</td></tr><tr><td align="right">217117</td><td align="left">Amyloid beta (A4) precursor protein-binding. family B. member 3 (APBB3)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_006051.2">NM_006051.2</ext-link></td><td align="left">0.294</td><td align="left">2.38e-09</td></tr><tr><td align="right">102677</td><td align="left">4-aminobutyrate aminotransferase (ABAT)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_000663.2">NM_000663.2</ext-link></td><td align="left">0.288</td><td align="left">2.38e-09</td></tr><tr><td align="right">216026</td><td align="right">216026</td><td/><td align="left">-0.840</td><td align="left">2.51e-09</td></tr><tr><td align="right">216643</td><td align="left">Leupaxin (Lpxn)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_134152.1">NM_134152.1</ext-link></td><td align="left">0.271</td><td align="left">2.93e-09</td></tr><tr><td align="right">205721</td><td align="right">205721</td><td/><td align="left">-0.256</td><td align="left">3.02e-09</td></tr><tr><td align="right">202919</td><td align="right">202919</td><td/><td align="left">0.272</td><td align="left">3.02e-09</td></tr><tr><td align="right">220257</td><td align="left">Similar to 60S ribosomal protein L39 (LOC651724)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="XM_940942.1">XM_940942.1</ext-link></td><td align="left">0.273</td><td align="left">3.10e-09</td></tr><tr><td align="right">103136</td><td align="left">Mitochondrial translational initiation factor 2 (MTIF2)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001005369.1">NM_001005369.1</ext-link></td><td align="left">-0.250</td><td align="left">3.35e-09</td></tr><tr><td align="right">201942</td><td align="right">201942</td><td/><td align="left">-0.286</td><td align="left">4.06e-09</td></tr><tr><td align="right">216638</td><td align="right">216638</td><td/><td align="left">0.304</td><td align="left">4.06e-09</td></tr><tr><td align="right">204605</td><td align="left">Similar to Maltase-glucoamylase, intestinal (LOC642103)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="XM_936233.1">XM_936233.1</ext-link></td><td align="left">-0.234</td><td align="left">4.12e-09</td></tr><tr><td align="right">208567</td><td align="right">208567</td><td/><td align="left">-0.307</td><td align="left">4.71e-09</td></tr><tr><td align="right">103361</td><td align="left">Phenylalanine hydroxylase (PAH)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_000277.1">NM_000277.1</ext-link></td><td align="left">-0.627</td><td align="left">5.15e-09</td></tr><tr><td align="right">202659</td><td align="left">FKBP1A-like (LOC654323)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001038000.1">NM_001038000.1</ext-link></td><td align="left">-0.312</td><td align="left">5.49e-09</td></tr><tr><td align="right">207982</td><td align="left">La ribonucleoprotein domain family. member 4 (LARP4)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_052879.3">NM_052879.3</ext-link></td><td align="left">-0.247</td><td align="left">5.63e-09</td></tr><tr><td align="right">204513</td><td align="left">Deoxyribonuclease II (DNASE2)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_214196.1">NM_214196.1</ext-link></td><td align="left">-0.255</td><td align="left">5.73e-09</td></tr><tr><td align="right">206370</td><td align="left">Syndecan binding protein (syntenin) (SDCBP)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001007068.1">NM_001007068.1</ext-link></td><td align="left">-0.208</td><td align="left">6.33e-09</td></tr><tr><td align="right">214304</td><td align="left">Heat shock transcription factor 1 (HSF1)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_005526.1">NM_005526.1</ext-link></td><td align="left">0.331</td><td align="left">7.28e-09</td></tr><tr><td align="right">103062</td><td align="left">FKBP1A-like (LOC654323)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001038000.1">NM_001038000.1</ext-link></td><td align="left">-1.046</td><td align="left">7.37e-09</td></tr><tr><td align="right">200969</td><td align="left">Family with sequence similarity 70. member A (FAM70A)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_017938.2">NM_017938.2</ext-link></td><td align="left">0.424</td><td align="left">7.37e-09</td></tr><tr><td align="right">206063</td><td align="left">Chromosome 14 open reading frame 43 (C14orf43)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_194278.2">NM_194278.2</ext-link></td><td align="left">-0.230</td><td align="left">8.05e-09</td></tr><tr><td align="right">206129</td><td align="right">206129</td><td/><td align="left">-0.447</td><td align="left">1.10e-08</td></tr><tr><td align="right">206577</td><td align="left">Leukemia inhibitory factor receptor (LIFR)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_002310.3">NM_002310.3</ext-link></td><td align="left">-0.265</td><td align="left">1.39e-08</td></tr><tr><td align="right">103378</td><td align="left">Similar to Putative steroid dehydrogenase KIK-I (LOC508455)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="XM_585231.2">XM_585231.2</ext-link></td><td align="left">-0.314</td><td align="left">1.49e-08</td></tr><tr><td align="right">204838</td><td align="left">FKBP1A-like (LOC654323)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001038000.1">NM_001038000.1</ext-link></td><td align="left">-0.309</td><td align="left">1.49e-08</td></tr><tr><td align="right">104413</td><td align="left">Pyruvate carboxylase (PC)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_214349.1">NM_214349.1</ext-link></td><td align="left">0.246</td><td align="left">1.49e-08</td></tr><tr><td align="right">219405</td><td align="left">Coiled-coil domain containing 42 (CCDC42)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_144681.1">NM_144681.1</ext-link></td><td align="left">0.447</td><td align="left">1.49e-08</td></tr><tr><td align="right">209512</td><td align="right">209512</td><td/><td align="left">-0.310</td><td align="left">1.72e-09</td></tr><tr><td align="right">221025</td><td align="right">221025</td><td/><td align="left">0.446</td><td align="left">1.74e-08</td></tr><tr><td align="right">200362</td><td align="left">Chromosome 3 open reading frame 43, transcript variant 1 (C3orf43)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="XM_173087.5">XM_173087.5</ext-link></td><td align="left">-0.199</td><td align="left">1.88e-08</td></tr><tr><td align="right">206386</td><td align="right">206386</td><td/><td align="left">-0.262</td><td align="left">2.65e-08</td></tr><tr><td align="right">211845</td><td align="left">Alanyl-tRNA synthetase domain containing 1 (AARSD1)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_025267.2">NM_025267.2</ext-link></td><td align="left">0.250</td><td align="left">3.31e-08</td></tr><tr><td align="right">208761</td><td align="left">Ring finger protein 148 (RNF148)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_198085.1">NM_198085.1</ext-link></td><td align="left">-0.189</td><td align="left">3.31e-08</td></tr><tr><td align="right">209449</td><td align="left">Similar to Zinc finger CCCH-type domain containing protein 11A transcript variant 3 (LOC441155)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="XM_930970.1">XM_930970.1</ext-link></td><td align="left">-0.237</td><td align="left">3.31e-08</td></tr><tr><td align="right">205142</td><td align="left">Bone morphogenetic protein receptor, type IA (BMPR1A)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_004329.2">NM_004329.2</ext-link></td><td align="left">-0.340</td><td align="left">3.31e-08</td></tr><tr><td align="right">203722</td><td align="right">203722</td><td/><td align="left">-0.217</td><td align="left">3.32e-08</td></tr><tr><td align="right">103279</td><td align="left">Orosomucoid 1 (ORM1)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_000607.1">NM_000607.1</ext-link></td><td align="left">-0.999</td><td align="left">3.32e-08</td></tr><tr><td align="right">101898</td><td align="left">Bifunctional apoptosis regulator (BFAR)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_016561.1">NM_016561.1</ext-link></td><td align="left">-0.189</td><td align="left">3.32e-08</td></tr><tr><td align="right">104972</td><td align="left">Cyclin-dependent kinase inhibitor 3 (CDKN3)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_214320.1">NM_214320.1</ext-link></td><td align="left">0.364</td><td align="left">3.45e-08</td></tr><tr><td align="right">202997</td><td align="right">202997</td><td/><td align="left">-0.270</td><td align="left">3.45e-08</td></tr><tr><td align="right">102991</td><td align="left">Cytochrome P450. family 2. subfamily E. polypeptide 1 (CYP2E1)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_214421.1">NM_214421.1</ext-link></td><td align="left">-1.079</td><td align="left">3.68e-08</td></tr><tr><td align="right">208187</td><td align="left">Early growth response 3 (EGR3)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_004430.2">NM_004430.2</ext-link></td><td align="left">0.378</td><td align="left">3.77e-08</td></tr><tr><td align="right">205347</td><td align="left">Histone H3.3A (H3F3A)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_213930.1">NM_213930.1</ext-link></td><td align="left">-0.237</td><td align="left">3.77e-08</td></tr><tr><td align="right">203762</td><td align="left">Zinc finger protein 7 (KOX 4. clone HF.16) (ZNF7)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_003416.1">NM_003416.1</ext-link></td><td align="left">-0.234</td><td align="left">4.37e-08</td></tr><tr><td align="right">102955</td><td align="left">L-kynurenine 3-monooxygenase Fpk (KMO)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_214076.1">NM_214076.1</ext-link></td><td align="left">-0.294</td><td align="left">4.81e-08</td></tr><tr><td align="right">212532</td><td align="left">TBC1 domain family. member 10B (TBC1D10B)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_015527.2">NM_015527.2</ext-link></td><td align="left">0.231</td><td align="left">4.99e-08</td></tr><tr><td align="right">206466</td><td align="left">Similar to glycosyltransferase-like domain containing 1 isoform a (LOC534677)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="XM_871633.1">XM_871633.1</ext-link></td><td align="left">-0.299</td><td align="left">5.15e-08</td></tr><tr><td align="right">104900</td><td align="left">Cullin 4A (CUL4A). transcript variant 1</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001008895.1">NM_001008895.1</ext-link></td><td align="left">0.349</td><td align="left">5.59e-08</td></tr><tr><td align="right">208650</td><td align="left">FKBP1A-like (LOC654323)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001038000.1">NM_001038000.1</ext-link></td><td align="left">-0.214</td><td align="left">5.83e-08</td></tr><tr><td align="right">205233</td><td align="left">5'-3' exoribonuclease 2 (XRN2)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_012255.3">NM_012255.3</ext-link></td><td align="left">0.217</td><td align="left">5.99e-08</td></tr><tr><td align="right">102672</td><td align="left">Cytochrome b-5 (CYB5)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001001770.1">NM_001001770.1</ext-link></td><td align="left">-0.528</td><td align="left">6.21e-08</td></tr><tr><td align="right">101034</td><td align="left">Programmed cell death 6 interacting protein (PDCD6IP)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_013374.3">NM_013374.3</ext-link></td><td align="left">-0.247</td><td align="left">6.21e-08</td></tr><tr><td align="right">206253</td><td align="left">Cyclin-dependent kinase 6 (CDK6)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001259.5">NM_001259.5</ext-link></td><td align="left">-0.260</td><td align="left">6.21e-08</td></tr><tr><td align="right">209516</td><td align="right">209516</td><td/><td align="left">-0.208</td><td align="left">6.21e-08</td></tr><tr><td align="right">100909</td><td align="left">Myeloid-associated differentiation marker (MYADM)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001020819.1">NM_001020819.1</ext-link></td><td align="left">-0.239</td><td align="left">7.42e-08</td></tr><tr><td align="right">205225</td><td align="left">Family with sequence similarity 80, member B (FAM80B)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_020734.1">NM_020734.1</ext-link></td><td align="left">0.247</td><td align="left">7.42e-08</td></tr><tr><td align="right">200683</td><td align="right">200683</td><td/><td align="left">-0.302</td><td align="left">7.64e-08</td></tr><tr><td align="right">208724</td><td align="left">Dickkopf-like 1 (soggy) (DKKL1)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_014419.3">NM_014419.3</ext-link></td><td align="left">0.275</td><td align="left">7.64e-08</td></tr><tr><td align="right">218253</td><td align="left">Mesoderm induction early response 1, family member 2 (MIER2)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_017550.1">NM_017550.1</ext-link></td><td align="left">0.265</td><td align="left">7.97e-08</td></tr><tr><td align="right">209231</td><td align="left">Heparin binding protein (HBP15/L22)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_213987.1">NM_213987.1</ext-link></td><td align="left">-0.343</td><td align="left">8.00e-08</td></tr><tr><td align="right">104000</td><td align="left">Ribosomal protein L39 (RPL39)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001000.2">NM_001000.2</ext-link></td><td align="left">0.264</td><td align="left">8.37e-08</td></tr><tr><td align="right">205725</td><td align="right">205725</td><td/><td align="left">-0.208</td><td align="left">8.73e-08</td></tr><tr><td align="right">102726</td><td align="left">Complement factor I (CFI)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_000204.1">NM_000204.1</ext-link></td><td align="left">-0.380</td><td align="left">8.77e-08</td></tr><tr><td align="right">100609</td><td align="left">Proline rich 13 (PRR13)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_001005355.1">NM_001005355.1</ext-link></td><td align="left">0.267</td><td align="left">9.03e-08</td></tr><tr><td align="right">104957</td><td align="left">Succinyl-CoA:alpha-ketoacid coenzyme A transferase (OXCT)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="NM_213938.1">NM_213938.1</ext-link></td><td align="left">0.392</td><td align="left">9.03e-08</td></tr><tr><td align="right">216727</td><td align="right">216727</td><td/><td align="left">-0.287</td><td align="left">9.12e-08</td></tr><tr><td align="right">200320</td><td align="left">Ubiquitin specific peptidase 32 (USP32)</td><td align="left"><ext-link ext-link-type="gen" xlink:href="XM_939750.1">XM_939750.1</ext-link></td><td align="left">0.439</td><td align="left">9.28e-08</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p>Microarray results for Duroc. Gene expression profiling was performed using 58 arrays and the 1% most differentially expressed genes were considered significant. The clone names are sequences with a hit to pig, human, mouse or bovine genes. Some genes are represented by several different clones on the array and may therefore show up more than once in the table, while some have no hits to the abovementioned species. The ID column gives the feature ID on the microarray, M value indicates fold change, t gives the t-statistics, P.value is the nominal p-value and adj.P.value is the FDR corrected p-value.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title><p>Microarray results for Norwegian Landrace. Gene expression profiling was performed using 58 arrays and the 1% most differentially expressed genes were considered significant. The clone names are sequences with a hit to pig, human, mouse or bovine genes. Some genes are represented by several different clones on the array and may therefore show up more than once in the table, while some have no hits to the abovementioned species. The ID column gives the feature ID on the microarray, M value indicates fold change, t gives the t-statistics, P.value is the nominal p-value and adj.P.value is the FDR corrected p-value.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional file 3</title><p>Gene ontology (GO) results for the cellular component ontology in Duroc. The top 1% differentially expressed genes in Duroc were analysed for over-represented (p &lt; 0.01) GO terms in the cellular component ontology.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S4"><caption><title>Additional file 4</title><p>Gene ontology (GO) results for the cellular component ontology in Norwegian Landrace. The top 1% differentially expressed genes in Norwegian Landrace were analysed for over-represented (p &lt; 0.01) GO terms in the cellular component ontology.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S5"><caption><title>Additional file 5</title><p>Androstenone values. The androstenone values (ppm) in Duroc high (DH), Duroc low (DL), Landrace high (NLH) and Landrace low (LL) animals used in this study.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S6"><caption><title>Additional file 6</title><p><bold>Boxplots</bold>. Boxplots of normalised arrays for Duroc (D) and Norwegian Landrace (L).</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S7"><caption><title>Additional file 7</title><p>Gene transcripts included in the rcPCR analyses.</p></caption></supplementary-material>
<table-wrap-foot><p>Gene expression profiles were identified in Limma using empirical Bayes moderated t-statistics. Fold change and statistical significance is shown (FDR-adjusted p-values). The clone names refer to hits to pig, human, mouse or bovine sequences. Some genes are represented by several different clones on the array and may therefore show up more than once in the table, while others have no hits to the abovementioned species.</p></table-wrap-foot> <table-wrap-foot><p>Gene expression profiles were identified in Limma using empirical Bayes moderated t-statistics. Fold change and statistical significance is shown (FDR-adjusted p-values). The clone names refer to hits to pig, human, mouse or bovine sequences. Some genes are represented by several different clones on the array and may therefore show up more than once in the table, while others have no hits to the abovementioned species.</p></table-wrap-foot>
<graphic xlink:href="1746-6148-4-29-1"/> <graphic xlink:href="1746-6148-4-29-2"/> <graphic xlink:href="1746-6148-4-29-3"/>
<media xlink:href="1746-6148-4-29-S1.txt" mimetype="text" mime-subtype="plain"><caption><p>Click here for file</p></caption></media> <media xlink:href="1746-6148-4-29-S2.txt" mimetype="text" mime-subtype="plain"><caption><p>Click here for file</p></caption></media> <media xlink:href="1746-6148-4-29-S3.txt" mimetype="text" mime-subtype="plain"><caption><p>Click here for file</p></caption></media> <media xlink:href="1746-6148-4-29-S4.txt" mimetype="text" mime-subtype="plain"><caption><p>Click here for file</p></caption></media> <media xlink:href="1746-6148-4-29-S5.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1746-6148-4-29-S6.png" mimetype="image" mime-subtype="png"><caption><p>Click here for file</p></caption></media> <media xlink:href="1746-6148-4-29-S7.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Patterson"], "given-names": ["RLS"], "article-title": ["5alpha-androst-16-ene-3-one: \u2013 Compound responsible for taint in boar fat"], "source": ["Journal of the science of food and agriculture"], "year": ["1968"], "fpage": ["19"]}, {"surname": ["Vold"], "given-names": ["E"], "article-title": ["Fleischproduktionseigenschaften bei ebern und kastraten"], "source": ["Meldinger fra Norges Landbruksh\u00f8gskole"], "year": ["1970"], "volume": ["49"], "fpage": ["1"], "lpage": ["25"]}, {"surname": ["Babol", "Squires"], "given-names": ["J", "EJ"], "article-title": ["Quality of meat from entire male pigs"], "source": ["Food research international"], "year": ["1995"], "volume": ["28"], "fpage": ["201"], "lpage": ["212"]}, {"surname": ["Tajet", "Andresen", "Meuwissen"], "given-names": ["H", "\u00d8", "THE"], "article-title": ["Prevention of boar taint in pig production: the 19th symposium of the nordic committee for veterinary scientific cooperation gardermoen, norway. 21\u201322 november 2005. Abstracts"], "source": ["Acta Veterinaria Scandinavica Supplement"], "year": ["2005"], "volume": ["48"], "fpage": ["22"], "lpage": ["23"]}, {"surname": ["Goldstein"], "given-names": ["JA"], "article-title": ["Clinical relevance of genetic polymorphisms in the human CYP2C subfamily"], "source": ["Journal of clinical pharmacology"], "year": ["2001"], "volume": ["52"], "fpage": ["349"], "lpage": ["355"]}, {"surname": ["Diaz", "Squires"], "given-names": ["GJ", "EJ"], "article-title": ["Metabolism of 3-methylindole by porcine liver microsomes: responsible cytochrome P450 enzymes"], "source": ["Toxicology science"], "year": ["2000"], "volume": ["55"], "fpage": ["284"], "lpage": ["292"]}, {"surname": ["Sinclair", "Hancock", "Gilmore", "Squires"], "given-names": ["PA", "S", "WJ", "EJ"], "article-title": ["Metabolism of the 16-androstene steroids in primary cultured porcine hepatocytes"], "source": ["J Steroid Biochem Mol Biol"], "year": ["2005"]}, {"surname": ["Diaz", "Squires"], "given-names": ["GJ", "EJ"], "article-title": ["Role of aldehyde oxidase in the hepatic in vitro metabolism of 3-methylindole in pigs"], "source": ["Journal of agricultur food chemistry"], "year": ["2000"], "volume": ["48"], "fpage": ["833"], "lpage": ["837"]}, {"surname": ["Tuomola", "Harpio", "Knuuttila", "Mikola", "L\u00f8vgren"], "given-names": ["M", "R", "P", "H", "T"], "article-title": ["Time-resolved fluoroimmunoassay for the measurement of androstenone in porcine serum and fat samples"], "source": ["Journal of agricultural food chemistry"], "year": ["1997"], "volume": ["45"], "fpage": ["3529"], "lpage": ["3534"]}, {"surname": ["Andresen"], "given-names": ["\u00d8"], "article-title": ["Development of radioimunoassay for 5alpha-adrost-16-en-3-one in pig peripheral plasma"], "source": ["Acta endochrinologia"], "year": ["1974"], "volume": ["76"], "fpage": ["377"], "lpage": ["387"]}, {"year": ["2007"], "comment": ["Type: Electronic Citation"]}, {"surname": ["Ihaka", "Gentleman"], "given-names": ["R", "R"], "article-title": ["R: A language for data analysis and graphics"], "source": ["Journal of computational and graphical statistics"], "year": ["1996"], "volume": ["5"], "fpage": ["299"], "lpage": ["314"]}, {"surname": ["Smyth"], "given-names": ["GK"], "article-title": ["Linear models and empirical bayes methods for assessing differential expression in microarray experiments"], "source": ["Statistical applications in genetics and molecular biology"], "year": ["2004"], "fpage": ["3"]}, {"surname": ["Smyth", "Gentleman R, Carey V, Dudoit S, Irizarry R, Huber W"], "given-names": ["GK"], "article-title": ["Limma: linear models for microarray data"], "source": ["Bioinformatics and Computational Biology Solutions using R and Bioconductor"], "year": ["2005"], "publisher-name": ["New York: Springer"], "fpage": ["397"], "lpage": ["420"]}, {"surname": ["Benjamini", "Hochberg"], "given-names": ["Y", "Y"], "article-title": ["Controlling the false discovery rate: A practical and powerful approach to multiple testing"], "source": ["Journal Of The Royal Statistical Society Series B"], "year": ["1995"], "volume": ["57"], "fpage": ["289"], "lpage": ["300"]}, {"year": ["2007"]}, {"year": ["2007"]}]
BMC Vet Res. 2008 Aug 6; 4:29