Skip to main content
Download PDF

A genome-wide association study on androstenone levels in pigs reveals a cluster of candidate genes on chromosome 6

BMC Genetics volume 11, Article number: 42 (2010) Cite this article

Abstract

Background

In many countries, male piglets are castrated shortly after birth because a proportion of un-castrated male pigs produce meat with an unpleasant flavour and odour. Main compounds of boar taint are androstenone and skatole. The aim of this high-density genome-wide association study was to identify single nucleotide polymorphisms (SNPs) associated with androstenone levels in a commercial sire line of pigs. The identification of major genetic effects causing boar taint would accelerate the reduction of boar taint through breeding to finally eliminate the need for castration.

Results

The Illumina Porcine 60K+SNP Beadchip was genotyped on 987 pigs divergent for androstenone concentration from a commercial Duroc-based sire line. The association analysis with 47,897 SNPs revealed that androstenone levels in fat tissue were significantly affected by 37 SNPs on pig chromosomes SSC1 and SSC6. Among them, the 5 most significant SNPs explained together 13.7% of the genetic variance in androstenone. On SSC6, a larger region of 10 Mb was shown to be associated with androstenone covering several candidate genes potentially involved in the synthesis and metabolism of androgens. Besides known candidate genes, such as cytochrome P450 A19 (CYP2A19), sulfotransferases SULT2A1, and SULT2B1, also new members of the cytochrome P450 CYP2 gene subfamilies and of the hydroxysteroid-dehydrogenases (HSD17B14) were found. In addition, the gene encoding the ß-chain of the luteinizing hormone (LHB) which induces steroid synthesis in the Leydig cells of the testis at onset of puberty maps to this area on SSC6. Interestingly, the gene encoding the α-chain of LH is also located in one of the highly significant areas on SSC1.

Conclusions

This study reveals several areas of the genome at high resolution responsible for variation of androstenone levels in intact boars. Major genetic factors on SSC1 and SSC6 showing moderate to large effects on androstenone concentration were identified in this commercial breeding line of pigs. Known and new candidate genes cluster especially on SSC6. For one of the most significant SNP variants, the difference in the proportion of animals surpassing the threshold of consumer acceptance between the two homozygous genotypes was as much as 15.6%.

Background

In many countries, male piglets are castrated shortly after birth to prevent boar taint, which is an urine-like, unpleasant flavour and odour released at cooking or heating of pork [1]. However, recent discussions on the pain associated with castration of the piglets early in life have led to a ban on castration without anaesthesia in some countries. In addition, studies have shown that un-castrated males grow faster and have an improved feed efficiency due to reduced fat deposition [2–4]. In future, if un-castrated males will be finished, boar taint needs to be prevented. Two of the major components related to the boar taint are androstenone and skatole [5–7]. Androstenone (5α-androst-16-en-3-one) is a male sex pheromone produced by the testes and stored in adipose tissue causing a perspiration-like odour [8, 9]. Androstenone precursors are also transported to the salivary glands which are capable to produce high levels of androstenone during sexual excitement [10, 11]. Skatole possesses strong faecal odour and is produced by the bacterial breakdown of the amino-acid tryptophane in the lower gut [12]. Skatole then diffuses into fat tissue.

There is considerable variation for androstenone and skatole between and within lines of pigs. Especially androstenone has high heritability estimates ranging from 0.25 to 0.88 [13, 14]. Somewhat lower heritabilities have been reported for skatole, between 0.19 and 0.55 [15, 16]. Two linkage studies using microsatellite markers have identified several QTL regions for androstenone and skatole in experimental crosses with 485 and 187 F2 animals, respectively [17, 18] pointing towards several areas in the genome affecting these traits. Also, single candidate genes involved in androstenone synthesis and metabolism of androstenone and skatole have been analyzed at the level of RNA and protein expression and in single SNP association studies (reviewed by Robic et al., 2008 [19]). However, no conclusive results showing functional mutations affecting androstenone and skatole levels in fat tissue have been described until now. Recently, large-scale microarray expression studies have reported hundreds of differentially expressed genes which might be involved in synthesis and degradation of androstenone and skatole in testis and liver [20, 21]. Subsequent analysis of SNPs in 121 differentially expressed genes identified 10 genes associated with one of the two traits [22]. Recently, Markljung et al. (2008) [23] reported 2 QTL for androstenone in 139 animals from a cross between Hamphsire and Landrace animals. Although these studies are of limited size and resolution, they indicate that several genetic factors seem to be involved in determining the levels of these boar taint compounds.

Recently, the first high-density 60K porcine SNP array has been developed [24] that offers a much higher resolution. A genome-wide association study (GWAS) was initiated using the SNP array to identify the chromosomal regions and specific SNPs influencing boar taint levels in a commercial breeding population. However, mean skatole levels (75 ng/g fat) in this population were far below the threshold accepted by consumers of 250 ng/g fat [25]. To reduce genotyping costs, a selective genotyping strategy for androstenone was applied. In this study, we present the results of a GWAS in pigs by genotyping 987 un-castrated male pigs from a commercial breeding population with large phenotypic variability for androstenone levels in fat, using the 60K (64,232) SNP array. The GWA resulted in an increased resolution compared to previous linkage studies. A large cluster of candidate genes within a 10 Mb region on SSC6 was identified. In addition, three new areas on SSC1 were detected that affect androstenone levels in this breeding line.

Methods

Animals and phenotypes

This experiment was conducted strictly in line with the regulations of the Dutch law on the protection of animals. Phenotypic measurements on androstenone were obtained from 1,663 boars slaughtered at a mean hot carcass weight of 95.71 kg. All the boars were purebred animals from a composite Duroc sire-line. Boar taint compounds were measured using fat samples from the neck collected from the left carcass side. The samples were stored under vacuum at -20°C. For androstenone, a fat extraction was done on the fat samples as described by Tuomola et al., 1997 [26]. Thereafter, androstenone concentrations in liquid fat were estimated by time-resolved fluoro-immunoassay [26] at the Hormone laboratory, Oslo.

Androstenone was not normally distributed and therefore log-transformed (ln-androstenone). The log transformed androstenone values were analysed using the following statistical model using ASREML [27]:

Where y = ln-androstenone; hcw = effect of hot carcass weight as covariate; age = effect of age at slaughter as covariate; fat= effect of fat depth at slaughter as covariate; batch= the random effect of the ith batch, pen = the random effect of the jth pen within the ith batch; litter = the random effect of the kth litter within the ith batch; a = additive genetic effect of the lth animal; e = residual effect. Systematic environmental effects were estimated using the full dataset (N = 1663). In the genome-wide association study, androstenone levels adjusted for systematic environmental effects were used and these were calculated as:

Selective genotyping

A simulation study was performed in order to select about 1000 animals from 1663 candidates for genotyping in an optimal way using the existing pedigree [28]. Ten markers and 1 QTL were simulated on 1 chromosome and also 1 chromosome was simulated without a QTL for determining the false-positive rate for a given threshold. Four alternatives for selecting 1000 individuals to be genotyped out of 1663 candidates were compared: 1. random, 2. selecting large half-sib families, 3. selecting high and low phenotypes, 4. selecting high and low phenotypes within full sib families. ANOVA was used to analyze each marker and determine the F-statistic. Selection of high and low phenotypes within full sibs showed the highest power (results not shown). Applying the selection of high and low phenotypes (within full sib families) to our data set consisting of 1663 pigs resulted in 987 pigs selected for genotyping. These pigs originated from 57 sires and 212 dams. Among them, 45 sires and 11 dams were available for genotyping.

Genotyping and quality control

Genotyping was performed using the PorcineSNP60 Beadchip of Illumina (San Diego, CA, USA) [24]. A total of 1043 samples (including sires and dams) were genotyped for 64,232 SNPs at Service XS (Leiden, The Netherlands) and data quality was evaluated. The average call rate for all samples was 98.4% ± 3.4. A total of 63 animals were removed due to pedigree errors (<99% correct genotypes). After quality control, 943 animals were available for the GWAS with 106 singletons and 313 divergent full sib pairs (2 or more full sibs). For the SNPs, a threshold of 30 pedigree errors or more was applied and 190 SNPs were removed. In addition, 10,210 SNPs were removed because of low quality score (GenCall score <0.7). A minor allele frequency of 0.01 was applied removing another 4,925 SNPs of which 980 were monomorphic. In total, 47,897 SNPs remained for the GWAS.

Genome-wide association analysis

Corrected log-transformed androstenone was analyzed as a quantitative trait under an additive model using the QFAM module of PLINK [29]. The more stringent within-sib-ship test within QFAM was performed which is robust for population stratification compared to the total-sib-ship test. Nominal scores were permuted to obtain an empirical p-value while maintaining familial correlation between genotype and phenotype. The permutation procedure employed by QFAM corrects for relatedness within families and was performed 1,000,000 times. Genomic control was used to correct for score inflation introduced by relatedness between family units (sib ships) [30]. False-discovery rate (FDR) was applied to correct for multiple-testing. The R package q-value [31] was used to calculate a FDR-based q value to measure the statistical significance at the genome-wide level for association studies. The cut-off of significant association at the whole genome level was set at q-value ≤ 0.05. The total variance explained by a SNP was calculated using ASReml version 2.0, [27]. For ASReml the full model (as described earlier and including the polygenic effect) was used for the animals genotyped including the SNP as a random effect.

The fraction of the phenotypic variance explained by the .

Linkage disequilibrium (LD) between SNPs was quantified as r² on all the animals of the GWA study using Haploview (V4.2; [32]) and the LD block was defined by the criteria of Gabriel et al. (2002) [33].

Identification of candidate genes

Porcine transcripts and annotation were downloaded from the porcine Ensembl data base (build9) and aligned with the human RefSeq mRNA sequences using BLAT [34]. The human-porcine comparative map was calculated based on the orthologous human-porcine transcripts and for the syntenic regions annotations were downloaded from the NCBI database (build37). Additional candidate genes present in human but not identified in the BLAT search against the human transcriptome were mapped to SSC6 performing a BLAST alignment with the porcine cDNA (SULT2A1) or the human homolog (SULT2B1, HSD17B14) against the porcine genome sequence (build9).

Results

The descriptive statistics of the phenotypic measurements of the boars used for the GWAS are given in Table 1. Animals were slaughtered at a mean age of 179.80 d with an average carcass weight of 95.7 kg. The average androstenone level was 1.88 μg/g melted fat, and the average skatole and indole levels were 91.11 ng/g and 54.15 ng/g, respectively. The change of the distribution of androstenone concentrations after selection of divergent sib pairs is shown in Figure 1 indicating that not only extreme androstenone levels are represented.

Table 1 Descriptive statistics for traits measured.
Full size table
Figure 1
figure1

Distribution of androstenone for the full dataset (N = 1663) and after selective genotyping was applied (N = 987).

Full size image

The GWA analysis using the threshold for FDR of q ≤ 0.05 showed that 37 SNPs were genome-wide significantly associated with log-androstenone (Figure 2 and additional file 1). Among them, thirty-five SNPs are located in regions with multiple significant SNPs. Three regions were identified on SSC1, and one larger region on SSC6. The region between 36.9 Mb and 44.9 Mb on SSC6 encompasses a large cluster of 31 significant SNPs. A single SNP analysis of the most significant SNPs on SSC1 and SSC6 using a mixed model and including a polygenic effect thereby correcting for other genetic factors affecting androstenone (background genome) is shown in additional file 1. The fraction of the phenotypic variance explained by a single SNP varies between 1.5% and 5.8%.

Figure 2
figure2

Association between ln-androstenone and 40,525 mapped SNPs across 18 autosomes using an additive model. Each dot represents one SNP. On the y-axis are -log10 (p-values), and on the x-axis are the physical positions of the SNPs by chromosome. Cut-off value is 4.35 which equals a FDR q-value ≤ 0.05.

Full size image

Figure 3 shows the means for the untransformed androstenone levels of the three genotypes of SNP nr 15 on SSC6. There is a difference of 0.66 μg/g between the mean level of the two homozygous genotypes. Correction for systematic environmental effects hardly affects the differences between the genotypes (data not shown). Moreover, among the animals homozygous for the allele associated with high androstenone levels, 39.6% of the animals surpass the threshold of consumer acceptance (2 μg/g). This proportion is markedly reduced by 15.6% in the homozygous low genotypes (24.0% above 2 μg/g). A more detailed view of SSC6 is shown in additional file 2. The high density of genes presently annotated in EnSembl on SSC6 (n = 351, additional file 3) is even more pronounced in the area of interest with a total of 24 genes between 36,9 Mb and 40 Mb. The homologous region in human on HSA 19q13 between 50 Mb and 52,2 Mb is also very gene-rich with a total of 139 genes and 255 transcripts being annotated until now.

Figure 3
figure3

Box plots of the distribution of the untransformed androstenone concentrations for the SNP MARC0049189 (nr 15). The mean is given in bold.

Full size image

Linkage disequilibrium was calculated between all the SNPs in the region between 33 and 44.9 Mb on SSC6. A large block of strong linkage disequilibrium in this area is observed. A part of this region, the area between 36.9 and 39.7 Mb is shown in Figure 4. All the 29 significant SNPs are present in only three major haplotypes in this population. Two copies of haplotype 1 has an average androstenone level of 2.13 μg/g and two copies of haplotypes 2 and 3 have an average level of 1.44 μg/g and 1.54 μg/g, respectively (Figure 4B). None of the remaining chromosomes show a comparable convincing cluster of closely linked SNPs associated with androstenone levels. Only isolated SNPs approach the significance threshold on SSC6 and SSC16 at 122 Mb and 105 Mb, respectively.

Figure 4
figure4

Linkage disequilibrium plot for the region between 36.9 Mb and 39.7 Mb on SSC6. All 31 significant SNPs (p ≤ 0.05 after FDR) and intervening SNPs for all animals (N = 943) are shown (A). The values in the boxes are pair wise SNP correlations (r²) and the box colour reflects the degree of correlation. B Haplotypes with all SNPs from the LD block are shown. Each line represents a haplotype and the frequency of the haplotype in this population is given at the end of the line. Haplotypes with a frequency below 2% are not included. Two SNPs are tagged and the SNP names are given in C.

Full size image

Discussion

Filtering of SNP data and statistical analyses

Quality control of the SNPs was based on the GenCall score, MAF and pedigree errors. Hardy-Weinberg equilibrium (HWE) was not considered relevant as a quality control tool as HWE is underpowered to detect genotyping errors [35] and only extreme sib pairs have been genotyped. GWA studies are particularly prone to spurious associations because ten thousands of associations are tested inflating the rate of false positives [36, 37]. In this study, FDR was used to control for false-positive associations due to multiple testing. The genomic control approach was used to account for spurious association due to population stratification [30] and because the breeding line is a composite line derived from three different breeds. Correction for the inflation by division reduces the unadjusted p-value to adjusted levels and accounts for relatedness between the sib ships and possible population stratification. However, in this study the deviation from the chi-square distribution under the null-hypothesis (no association) was very low (λGC = 1.06).

QTL areas

Mainly two chromosomes harbour highly significant associations with fat androstenone levels, a rather broad area of 10 Mb on SSC6, and three different regions on SSC1. For each region on SSC1 only 1-4 SNPs pass the significance level, whereas on SSC6 a total of 31 highly significant SNPs are detected. This is for the first time that an association with androstenone or related boar taint traits has been reported on SSC1. However, several studies have described QTL effects for traits related to boar taint on SSC6 (Figure 5). In an experimental F2 cross with Large White x Meishan, Lee et al. (2005) [17] reported a QTL for androstenone level in fat on SSC6 partially overlapping with the area identified here. In the same study, QTLs from a sensory panel for subjective pork flavour and boar flavour in lean meat were described that are also located nearby the area in our study. Finally, Szyda et al. (2003)[38] identified a QTL for smell intensity in a Duroc X Norwegian Landrace cross covering the area of interest. Considering the low resolution of these QTL studies, it is not possible to conclude whether they might be caused by the same genes segregating as in our study. The remaining boar taint QTL previously identified on SSC6 for smell intensity [39], subjective pork odour and skatole measurements [17, 40] are located distal or proximal on SSC6. None of the other QTL studies investigating androstenone or sensory panel traits identified effects on SSC6 [19, 24, 39].

Figure 5
figure5

Location of the QTL from PigQTLdb for boar taint traits on the physical map of Sus scrofa build9 SSC6. The references and traits of the QTLs are given in Table 3. Positions in Mb were deduced from a BLAST alignment with the microsatellite markers. The green bar indicates the region found in this GWA study between 33 Mb and 45 Mb.

Full size image

Candidate Genes

For the two major areas of interest, genes potentially affecting steroid synthesis and metabolism of androstenone are listed in Table 2. The region on SSC 6, which is an extremely gene-dense area, shows several candidate genes located closely together between 37 - 38 Mb. Hydroxysteroid sulfotransferaseA1 (SULT2A1) maps to the homologous region in human and has not been annotated in the pig genome sequence. However, a BLAT search with porcine cDNA identifies 100% homology with two exons and including 3'untranslated sequence (pos.380-508, and 780-999). SULT2A1 catalyses the sulfoconjugation of 16-androstene-steroids in liver [41] and testis [42] and has been analysed earlier as a candidate gene for androstenone. Testicular activity of the enzyme was shown to be negatively correlated with fat androstenone levels in Yorkshire boars [43]. However, Moe et al. (2007) [21] found an increased expression of SULT2A1 mRNA in testis of Duroc and Landrace boars with high androstenone levels. In the same study, SULT2B1 expression was also increased in Landrace animals with high androstenone levels. SULT2B1 is located near SULT2A1 in human, but it has not been annotated in the pig genome sequence and cannot be located by BLAT alignment either. The overexpression of SULT2B1 does not explain the role of the sulfotransferases as inactivating enzymes [21]. SULT2B1 is selective for the sulfation of 3ß-hydroxysteroids, and Falany et al. (2006) [44] suggest a role in regulating the responsiveness of cells to adrenal androgens by reducing their conversion to more potent androgens and estrogens. In human, SULT2B1b is not expressed in the liver, however the different physiological functions of the two isoforms remain to be analyzed. Another conjugating enzyme also located in this area is HSD17B14. Differential expression of different hydroxysteroid dehydrogenases of the HSD17ß family (HSD17B4, HSD17B11) in the testis has been reported by Moe et al. (2007) [21]. HSD17B4 catalyzes the last step of androgen and estrogen synthesis. However, the function of HSD17B14 has only recently been investigated in human, and Lukacik et al. (2007) [45] suggest a role for the local inactivation of steroids in the nervous system and placenta. Northern blot analysis of human tissues shows that the gene is highly expressed in the liver, but not in testis. Adjacent, another new candidate gene, LHB is located, forming the ß-chain precursor of the luteinizing hormone (LH). At onset of puberty, LH is secreted by the pituitary gland and induces steroid synthesis in the Leydig cells of the testis. Interestingly, the gene encoding the α-chain of this glycoprotein hormone (CGA chorionic gonadotropin A) is located on SSC1 in the area of 58 Mb which also shows a significant effect in this study. No known or potentially interesting candidate genes could be pinpointed for the remaining two regions around 146 and 290 Mb on SSC1. Finally, the region on SSC6 extends to a second peak of SNPs nearly reaching the significance threshold around 33 Mb (additional file 2). Several cytochrome P450 genes of the CYP2 family are located there. This CYP2ABFGST cluster contains genes from multiple subfamilies [46]. In human, CYP2A6, CYP2A7, CYP2B6, CYP2A13, CYP2F1 and CYP2S1 cluster together. From these, aromatase (CYP2A19), which is the homolog of CYP2A6 in pigs, is known to catalyze the synthesis of estrogens from androgens. The pig expresses two isoforms in the testis. Moe et al. (2007)[21] have shown an upregulation of both isoforms in testis and liver of high-androstenone boars. Recently, Moe et al. (2009) [22] reported SNPs within candidate genes associated with androstenone levels in a commercial Duroc line. However, none of the candidate genes reported by Moe et al. (2009) [22] overlap with the major regions identified here.

Table 2 Candidate genes derived from porcine Ensembl build9.
Full size table
Table 3 Overview of the identified QTL and flanking microsatellites on SSC6 for traits related to boar taint.
Full size table

Taken together, there is overwhelming evidence from previous QTL studies, candidate genes and differential expression that the region on SSC6 contains genetic elements affecting androstenone levels in boars. In order to disentangle the effects of the regions containing the CYP450 genes and the area around 37 Mb, a mixed-model analysis combining the effects of two SNPs (H3GA0052956 at 33.5 and MARC0049189 at 38.3 Mb) was performed. In this model the fraction of the phenotypic variance explained by both SNPs is 2.1% and 3.6% and together 5.7%. This means that both regions explain a part of the effect of the whole region but due to the high LD between the SNPs they capture the same variation individually (5.76%, additional file 1). Therefore, both areas remain relevant for the determination of androstenone levels in this population. This breed is a composite line which could explain this large extent of LD. More data from other unrelated lines or crossbred animals showing the same effect are needed to further reduce the region of interest.

Effect size and application for breeding

Due to the skewed distribution of androstenone levels, even the use of a single marker would reduce the proportion of animals surpassing the threshold for consumer acceptance of 2 μg/g fat considerably. The difference between the two homozygous genotypes amounts to 15.6% (Figure 3). Sorting all offspring by the estimated androstenone effect of marker 50 and comparing the haplotypes of the 10 highest animals shows that all individuals are homozygous for the first haplotype shown in Figure 4. Furthermore, this haplotype is completely absent in the group of 10 animals with the lowest effects (data not shown). The 5 major SNPs (SNP nr. 1, 5, 6 on SSC1 and SNP nr. 15, 124 on SSC6) on SSC1 and 6 together explain 8.8% of the phenotypic variance, and considering a heritability of 64% [47] they account for 13.7% of the additive genetic variance.

A sustainable breeding scheme takes also into account the correlated effects on other production and reproduction traits. In general, the genetic correlations with growth, fatness and muscle depth are very low and favourable and therefore no serious negative effects on genetic progress due to selection against androstenone are to be expected [47]. Also, the positive genetic correlation with skatole would reduce skatole levels indirectly. However, the genetic correlation with fertility traits needs special attention. Male fertility data are not available on the animals in this study because they were slaughtered as commercial fatteners. Female fertility observations are only available on related animals and therefore estimates of genetic parameters have large standard errors (data not shown). A more extended study is underway to monitor the effects of selection against androstenone on male and female fertility. Furthermore, the effects in other lines that form part of the crossbreeding scheme to produce fattening pigs will be investigated.

Conclusion

This study clearly shows the large increase in resolution of high-density SNP panels compared to earlier linkage studies using microsatellite markers. Several regions in the genome affect androstenone levels in fat in this commercial breeding line of pigs. The genome-wide significant SNPs detected on SSC1 and SSC6 show moderate to large effects explaining a fraction of the phenotypic variance of 2-6%. The candidate genes identified in these areas in the pig genome or via the comparative map in human include genes investigated in earlier reports. In addition, new genes from the pathways of the synthesis and metabolism of androstenone such as LHA, LHB, and HSD17B14 are detected. The rather large LD block seen in this population around 33-45 Mb on SSC6 prevents to disentangle the combined effects of these genes and to pinpoint more specifically the responsible genetic elements. Nevertheless, the most significant SNPs can already be used to accelerate genetic progress in breeding against androstenone in this sire line. However, genetic correlations with production traits and especially possible negative effects on fertility traits will deserve special attention.

References

  1. 1.

    Bonneau M: Boar taint in entire male pigs- a review. Pig News Information. 1997, 18: 15-18.

    Google Scholar 

  2. 2.

    Candek-Potokar M, Zlender B, Lefaucheur L, Bonneau M: Effects of age and/or weight at slaughter on longissimus dorsi muscle: Biochemical traits and sensory quality in pigs. Meat Science. 1998, 48: 287-300. 10.1016/S0309-1740(97)00109-5.

    Article  CAS  PubMed  Google Scholar 

  3. 3.

    Metz C, Claus R: Active immunization of boars against GnRH does not affect growth hormone but lowers IGF-I in plasma. Livestock Production Science. 2003, 81: 129-137. 10.1016/S0301-6226(02)00302-0.

    Article  Google Scholar 

  4. 4.

    Xue J, Dial G, Pettigrew J: Performance, carcass, and meat quality advantages of boars over barrows: A literature review. Swine Health and Production. 1997, 5: 21-28.

    Google Scholar 

  5. 5.

    Patterson R: 5-androst-16-en-3-one, compound responsible for taint in boar fat. Journal of the Science of Food and Agriculture. 1968, 19: 31-38. 10.1002/jsfa.2740190107.

    Article  CAS  Google Scholar 

  6. 6.

    Vold E: Fleischproduktionseigenschaften bei Ebern und Kastraten. IV. Organoleptische und gaschromatografische Untersuchungen wasserdampfflüchtiger Stoffe des Rückenspeckes von Ebern. Meldinger fra Norges Landbrukshøgskole. 1970, 49: 1-25.

    Google Scholar 

  7. 7.

    Malmfors B, Lundström K: Consumer reactions to boar meat-a review. Livestock Production Science. 1983, 10: 187-196. 10.1016/0301-6226(83)90034-9.

    Article  Google Scholar 

  8. 8.

    Bonneau M: Compounds responsible for boar taint, with special emphasis on androstenone: a review. Livestock Production Science. 1982, 9: 687-705. 10.1016/0301-6226(82)90017-3.

    Article  CAS  Google Scholar 

  9. 9.

    Perry G, Patterson R, MacFie H, Stinson C: Pig courtship behaviour: pheromonal property of androstene steroids in male submaximillary secretion. Animal Production. 1980, 331: 191-199.

    Article  Google Scholar 

  10. 10.

    Gower D: 16-Unsaturated C 19 steroids. A review of their chemistry, biochemistry and possible physiological role. Journal of steroid biochemistry. 1972, 3: 45-103. 10.1016/0022-4731(72)90011-8.

    Article  CAS  PubMed  Google Scholar 

  11. 11.

    Claus R: Pheromone bei Säugetieren unter besonderer Berücksichtigung des Ebergeruchsstoffes und seiner Beziehung zu anderen Hodensteroiden. Fortschritte der Tierphysiologie und Tierernährung. 1979, 10: 1-136.

    CAS  Google Scholar 

  12. 12.

    Yokoyama M, Carlson J: Microbial metabolites of tryptophan in the intestinal tract with special reference to skatole. American Journal of Clinical Nutrition. 1979, 32: 173-178.

    CAS  PubMed  Google Scholar 

  13. 13.

    Sellier P: Genetics of meat and carcass traits. The Genetics of the Pig. Edited by: Rothschild MF, Ruvinsky A. 1998, Oxon: CAB International, 463-510.

    Google Scholar 

  14. 14.

    Sellier P, Le Roy P, Fouilloux M, Gruand J, Bonneau M: Responses to restricted index selection and genetic parameters for fat androstenone level and sexual maturity status of young boars. Livestock Production Science. 2000, 63: 265-274. 10.1016/S0301-6226(99)00127-X.

    Article  Google Scholar 

  15. 15.

    Pederson B: Heritability of skatole in back fat. Skatole and boar taint. Edited by: Jensen W. 1998, Roskilde: Danish Meat Research Institute, 129-136.

    Google Scholar 

  16. 16.

    Tajet H, Andresen O, Meuwissen T: Estimation of genetic parameters of boar taint; skatole and androstenone and their correlations with sexual maturation. Acta Veterinaria Scandinavica. 2006, 48 (Suppl 1): S9-10.1186/1751-0147-48-S1-S9.

    Article  Google Scholar 

  17. 17.

    Lee G, Archibald A, Law A, Lloyd S, Wood J, Haley C: Detection of quantitative trait loci for androstenone, skatole and boar taint in a cross between Large White and Meishan pigs. Animal Genetics. 2005, 36: 14-22. 10.1111/j.1365-2052.2004.01214.x.

    Article  CAS  PubMed  Google Scholar 

  18. 18.

    Quintanilla R, Demeure O, Bidanel J, Milan D, Iannuccelli N, Amigues Y, Gruand J, Renard C, Chevalet C, Bonneau M: Detection of quantitative trait loci for fat androstenone levels in pigs 1. Journal of Animal Science. 2003, 81: 385-394.

    CAS  PubMed  Google Scholar 

  19. 19.

    Robic A, Larzul C, Bonneau M: Genetic and metabolic aspects of androstenone and skatole deposition in pig adipose tissue: A review (Open Access publication). Genetics Selection Evolution. 2007, 40: 129-143.

    Google Scholar 

  20. 20.

    Moe M, Lien S, Bendixen C, Hedegaard J, Hornshøj H, Berget I, Meuwissen T, Grindflek E: Gene expression profiles in liver of pigs with extreme high and low levels of androstenone. BMC Veterinary Research. 2008, 4: 29-10.1186/1746-6148-4-29.

    PubMed Central  Article  PubMed  Google Scholar 

  21. 21.

    Moe M, Meuwissen T, Lien S, Bendixen C, Wang X, Conley L, Berget I, Tajet H, Grindflek E: Gene expression profiles in testis of pigs with extreme high and low levels of androstenone. BMC Genomics. 2007, 8: 405-10.1186/1471-2164-8-405.

    PubMed Central  Article  PubMed  Google Scholar 

  22. 22.

    Moe M, Lien S, Aasmundstad T, Meuwissen T, Hansen M, Bendixen C, Grindflek E: Association between SNPs within candidate genes and compounds related to boar taint and reproduction. BMC Genetics. 2009, 10: 32-10.1186/1471-2156-10-32.

    PubMed Central  Article  PubMed  Google Scholar 

  23. 23.

    Markljung E, Braunschweig M, Karlskov-Mortensen P, Bruun C, Sawera M, Cho I, Hedebro-Velander I, Josell Å, Lundström K, Seth G: Genome-wide identification of quantitative trait loci in a cross between Hampshire and Landrace II: Meat quality traits. BMC Genetics. 2008, 9: 22-10.1186/1471-2156-9-22.

    PubMed Central  Article  PubMed  Google Scholar 

  24. 24.

    Ramos A, Crooijmans R, Affara N, Amaral A, Archibald A, Beever J, Bendixen C, Churcher C, Clark R, Dehais P: Design of a High Density SNP Genotyping Assay in the Pig Using SNPs Identified and Characterized by Next Generation Sequencing Technology. PLoS ONE. 2009, 4: e6524-10.1371/journal.pone.0006524.

    PubMed Central  Article  PubMed  Google Scholar 

  25. 25.

    Walstra P, Claudi-Magnussen C, Chevillon P, Von Seth G, Diestre A, Matthews K, Homer D, Bonneau M: An international study on the importance of androstenone and skatole for boar taint: levels of androstenone and skatole by country and season. Livestock Production Science. 1999, 62: 15-28. 10.1016/S0301-6226(99)00054-8.

    Article  Google Scholar 

  26. 26.

    Tuomola M, Harpio R, Knuuttila P, Mikola H, Lovgren T: Time-resolved fluoroimmunoassay for the measurement of androstenone in porcine serum and fat samples. Journal of Agricultural and Food Chemistry. 1997, 45: 3529-3534. 10.1021/jf9702398.

    Article  CAS  Google Scholar 

  27. 27.

    Gilmour A, Gogel B, Cullis B, Welham S, Thompson R: ASReml user guide release 1.0. 2002, VSN International Ltd, Hemel Hempstead, UK, 267-

    Google Scholar 

  28. 28.

    Duijvesteijn N, de Koning D: Selective genotyping in commercial populations. Book of abstract of the 60th Annual meeting of the EAAP, Barcelona, Spain. 4654, (Abstr)

    Google Scholar 

  29. 29.

    Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira M, Bender D, Maller J, Sklar P, de Bakker P, Daly M: PLINK: a tool set for whole-genome association and population-based linkage analyses. The American Journal of Human Genetics. 2007, 81: 559-575. 10.1086/519795.

    Article  CAS  PubMed  Google Scholar 

  30. 30.

    Devlin B, Roeder K: Genomic control for association studies. Biometrics. 1999, 55: 997-1004. 10.1111/j.0006-341X.1999.00997.x.

    Article  CAS  PubMed  Google Scholar 

  31. 31.

    Dabney A, Storey J: R-package: q-value extimation for false discovery rate control. 2009

    Google Scholar 

  32. 32.

    Barrett J, Fry B, Maller J, Daly M: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21: 263-265. 10.1093/bioinformatics/bth457.

    Article  CAS  PubMed  Google Scholar 

  33. 33.

    Kent W: BLAT--the BLAST-like alignment tool. Genome Research. 2002, 12: 656-664.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  34. 34.

    Gabriel S, Schaffner S: The structure of haplotype blocks in the human genome. Science. 2002, 296: 2225-2229. 10.1126/science.1069424.

    Article  CAS  PubMed  Google Scholar 

  35. 35.

    Cox D, Kraft P: Quantification of the power of Hardy-Weinberg equilibrium testing to detect genotyping error. Human Heredity. 2006, 61: 10-14. 10.1159/000091787.

    Article  PubMed  Google Scholar 

  36. 36.

    McCarthy M, Abecasis G, Cardon L, Goldstein D, Little J, Ioannidis J, Hirschhorn J: Genome-wide association studies for complex traits: consensus, uncertainty and challenges. Nature Reviews Genetics. 2008, 9: 356-369. 10.1038/nrg2344.

    Article  CAS  PubMed  Google Scholar 

  37. 37.

    Pearson T, Manolio T: How to interpret a genome-wide association study. Journal of the American Medical Association. 2008, 299: 1335-1344. 10.1001/jama.299.11.1335.

    Article  CAS  PubMed  Google Scholar 

  38. 38.

    Szyda J, Grindflek E, Liu Z, Lien S: Multivariate mixed inheritance models for QTL detection on porcine chromosome 6. Genetics Research. 2003, 81: 65-73. 10.1017/S0016672302006043.

    Article  CAS  Google Scholar 

  39. 39.

    Grindflek E, Szyda J, Liu Z, Lien S: Detection of quantitative trait loci for meat quality in a commercial slaughter pig cross. Mammalian Genome. 2001, 12: 299-304. 10.1007/s003350010278.

    Article  CAS  PubMed  Google Scholar 

  40. 40.

    Varona L, Vidal O, Quintanilla R, Gil M, Sanchez A, Folch J, Hortos M, Rius M, Amills M, Noguera J: Bayesian analysis of quantitative trait loci for boar taint in a Landrace outbred population. Journal of Animal Science. 2005, 83: 301-307.

    CAS  PubMed  Google Scholar 

  41. 41.

    Sinclair P, Squires E: Testicular sulfoconjugation of the 16-androstene steroids by hydroxysteroid sulfotransferase: Its effect on the concentrations of 5 {alpha}-androstenone in plasma and fat of the mature domestic boar. Journal of Animal Science. 2005, 83: 358-365.

    CAS  PubMed  Google Scholar 

  42. 42.

    Sinclair P, Squires E, Raeside J, Renaud R: Synthesis of free and sulphoconjugated 16-androstene steroids by the Leydig cells of the mature domestic boar. Journal of Steroid Biochemistry and Molecular Biology. 2005, 96: 217-228. 10.1016/j.jsbmb.2005.02.017.

    Article  CAS  PubMed  Google Scholar 

  43. 43.

    Sinclair P, Gilmore W, Lin Z, Lou Y, Squires E: Molecular cloning and regulation of porcine SULT2A1: relationship between SULT2A1 expression and sulfoconjugation of androstenone. Journal of Molecular Endocrinology. 2006, 36: 301-311. 10.1677/jme.1.01847.

    Article  CAS  PubMed  Google Scholar 

  44. 44.

    Falany C, He D, Dumas N, Frost A, Falany J: Human cytosolic sulfotransferase 2B1: isoform expression, tissue specificity and subcellular localization. Journal of Steroid Biochemistry and Molecular Biology. 2006, 102: 214-221. 10.1016/j.jsbmb.2006.09.011.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  45. 45.

    Lukacik P, Keller B, Bunkoczi G, Kavanagh K, Lee W, Adamski J, Oppermann U: Structural and biochemical characterization of human orphan DHRS10 reveals a novel cytosolic enzyme with steroid dehydrogenase activity. Biochemical Journal. 2007, 402: 419-427. 10.1042/BJ20061319.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  46. 46.

    Hu S, Wang H, Knisely A, Reddy S, Kovacevic D, Liu Z, Hoffman S: Evolution of the CYP2ABFGST gene cluster in rat, and in fine-scale comparison among rodent and primate species. Genetica. 2008, 133: 215-226. 10.1007/s10709-007-9206-x.

    Article  CAS  PubMed  Google Scholar 

  47. 47.

    Merks J, Hanenberg E, Bloemhof S, Knol E: Genetic opportunities for pork production without castration. Animal Welfare. 2009, 18: 539-544.

    CAS  Google Scholar 

Download references

Acknowledgements

We are grateful to Marcos Ramos, Wageningen University, for quality control of the SNP data and TOPIGS for providing the data of non-castrated boars and tissue sampling.

This research project has been co-financed by the European Commission, within the 6th Framework Programme, contract No. FOOD-CT-2006-016250 ("SABRE"). The text represents the authors' views and does not necessarily represent a position of the Commission who will not be liable for the use made of such information.

Author information

Affiliations

  1. IPG, Institute for Pig Genetics B.V., PO Box 43, 6640AA, Beuningen, the Netherlands

    Naomi Duijvesteijn, Egbert F Knol, Jan WM Merks & Barbara Harlizius

  2. ABG, Animal Breeding and Genomics Centre, Wageningen University, PO Box 338, 67001 AH, Wageningen, the Netherlands

    Richard PMA Crooijmans, Martien AM Groenen & Henk Bovenhuis

Authors
  1. Naomi Duijvesteijn
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Egbert F Knol
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Jan WM Merks
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Richard PMA Crooijmans
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Martien AM Groenen
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Henk Bovenhuis
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Barbara Harlizius
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Naomi Duijvesteijn.

Additional information

Authors' contributions

ND was involved in sample collection, design of selective genotyping, organization of the genotyping experiment, conducted statistical analyses, prepared figures and tables, and took part in writing the paper. EFK was involved in planning the project, statistical supervision and experimental set up. JWMM initiated the study, was involved in planning of the project and organization of the genotyping experiment. RPMAC was involved in DNA collection and the organization of the genotyping experiment. MAMG mapped all SNPs to porcine genome build9, identified human-porcine orthologous genes and was involved in the analysis and discussion of the comparative mapping data. HB was involved in the discussion and evaluation of statistical issues. BH coordinated the study, investigated the comparative map and selection of candidate genes, was involved in discussion on statistical issues, and drafted the paper. All authors have read and approved the final manuscript.

Electronic supplementary material

Additional file 1: Association results for corrected ln-androstenone. If the DbSNP accession number was not available for a SNP, the forward and top genomic sequence of the SNP is given. (XLS 27 KB)

Chromosome 6 close-up with the -log (Pval) from the PLINK analyses

Additional file 2: . Above a -log p-value of 4.35 a SNP is considered significant. The start positions of porcine genes (n = 351) from EnSembl are plotted underneath as vertical grey bar based on the sequence of Sus scrofa build9. Vertical bars indicate the interval chosen for LD analysis in Figure 4. (PDF 169 KB)

The annoted genes on SSC 6

Additional file 3: . The start and stop position of the genes are given in bp and are based on Sus Scrofa build9. BLAT score and human gene name is also provided. (XLS 105 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Duijvesteijn, N., Knol, E.F., Merks, J.W. et al. A genome-wide association study on androstenone levels in pigs reveals a cluster of candidate genes on chromosome 6. BMC Genet 11, 42 (2010). https://doi.org/10.1186/1471-2156-11-42

Download citation

Keywords

  • Significant SNPs
  • Skatole
  • Androstenone
  • Boar Taint
  • Androstenone Level

BMC Genetics

ISSN: 1471-2156

Contact us