Transcript profiling of candidate genes in testis of pigs exhibiting large differences in androstenone levels

Involvement of porcine and human carbonyl reductases in the metabolism of epiandrosterone, 11-oxygenated steroids, neurosteroids, and corticosteroids

Porcine carbonyl reductases (pCBR1 and pCBR-N1) and aldo-keto reductases (pAKR1C1 and pAKR1C4) exhibit hydroxysteroid dehydrogenase (HSD) activity. However, their roles in the metabolism of porcine-specific androgens (19-nortestosterone and epiandrosterone), 11-oxygenated androgens, neurosteroids, and corticosteroids remain unclear. Here, we compared the steroid specificity of the four recombinant enzymes by kinetic and product analyses. In C18/C19-steroids,11-keto- and 11β-hydroxy-5α-androstane-3,17-diones were reduced by all the enzymes, whereas 5α-dihydronandrolone (19-nortestosterone metabolite) and 11-ketodihydrotestosterone were reduced by pCBR1, pCBR-N1, and pAKR1C1, of which pCBR1 exhibited the lowest (submicromolar) Km values. Product analysis showed that pCBR1 and pCBR-N1 function as 3α/β-HSDs, in contrast to pAKR1C1 and pAKR1C4 (acting as 3β-HSD and 3α-HSD, respectively). Additionally, 17β-HSD activity was observed in pCBR1 and pCBR-N1 (toward epiandrosterone and its 11-oxygenated derivatives) and in pAKR1C1 (toward androsterone, 4-androstene-3,17-dione and their 11-oxygenated derivatives). The four enzymes also showed different substrate specificity for 3-keto-5α/β-dihydro-C21-steroids, including GABAergic neurosteroid precursors and corticosteroid metabolites. 5β-Dihydroprogesterone was reduced by all the enzymes, whereas 5α-dihydroprogesterone was reduced only by pCBR1, and 5α/β-dihydrodeoxycorticosterones by pCBR1 and pCBR-N1. The two pCBRs also reduced the 5α/β-dihydro-metabolites of cortisol, 11-deoxycortisol, cortisone, and corticosterone. pCBR1 exhibited lower Km values (0.3-2.9 μM) for the 3-keto-C21-steroids than pCBR-N1 (Km=10-36 μM). The reduced products of the 3-keto-C21-steroids by pCBR1 and pCBR-N1 were their 3α-hydroxy-metabolites. Finally, we found that human CBR1 has similar substrate specificity for the C18/C19/C21-steroids to pCBR-N1. Based on these results, it was concluded that porcine and human CBRs can be involved in the metabolism of the aforementioned steroids as 3α/β,17β-HSDs.

Genomic Prediction and Genome-Wide Association Study for Boar Taint Compounds

With a perspective future ban on surgical castration in Europe, selecting pigs with reduced ability to accumulate boar taint (BT) compounds (androstenone, indole, skatole) in their tissues seems a promising strategy. BT compound concentrations were quantified in the adipose tissue of 1075 boars genotyped at 29,844 SNPs. Traditional and SNP-based breeding values were estimated using pedigree-based BLUP (PBLUP) and genomic BLUP (GBLUP), respectively. Heritabilities for BT compounds were moderate (0.30-0.52). The accuracies of GBLUP and PBLUP were significantly different for androstenone (0.58 and 0.36, respectively), but comparable for indole and skatole (~0.43 and ~0.47, respectively). Several SNP windows, each explaining a small percentage of the variance of BT compound concentrations, were identified in a genome-wide association study (GWAS). A total of 18 candidate genes previously associated with BT (MX1), reproduction traits (TCF21, NME5, PTGFR, KCNQ1, UMODL1), and fat metabolism (CTSD, SYT8, TNNI2, CD81, EGR1, GIPC2, MIGA1, NEGR1, CCSER1, MTMR2, LPL, ERFE) were identified in the post-GWAS analysis. The large number of genes related to fat metabolism might be explained by the relationship between sexual steroid levels and fat deposition and be partially ascribed to the pig line investigated, which is selected for ham quality and not for lean growth.

Copy Number Variation Analysis Revealed the Evolutionary Difference between Chinese Indigenous Pigs and Asian Wild Boars

Copy number variation (CNV) has been widely used to study the evolution of different species. We first discovered different CNVs in 24 Anqingliubai pigs and 6 Asian wild boars using next-generation sequencing at the whole-genome level with 10× depth to understand the relationship between genetic evolution and production traits in wild boars and domestic pigs. A total of 97,489 CNVs were identified and divided into 10,429 copy number variation regions (CNVRs), occupying 32.06% of the porcine genome. Chromosome 1 had the most CNVRs, and chromosome 18 had the least. Ninety-six CNVRs were selected using VST 1% based on the signatures of all CNVRs, and sixty-five genes were identified in the selected regions. These genes were strongly correlated with traits distinguishing groups by enrichment in Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathways, such as growth (CD36), reproduction (CIT, RLN), detoxification (CYP3A29), and fatty acid metabolism (ELOVL6). The QTL overlapping regions were associated with meat traits, growth, and immunity, which was consistent with CNV analysis. Our findings increase the understanding of evolved genome structural variations between wild boars and domestic pigs, and provide new molecular biomarkers to guide breeding and the efficient use of available genetic resources.

Porcine aldo-keto reductase 1C subfamily members AKR1C1 and AKR1C4: Substrate specificity, inhibitor sensitivity and activators

Most members of the aldo-keto reductase (AKR) 1 C subfamily are hydroxysteroid dehydrogenases (HSDs). Similarly to humans, four genes for AKR1C proteins (AKR1C1-AKR1C4) have been identified in the pig, which is a suitable species for biomedical research model of human diseases and optimal organ donor for xenotransplantation. Previous study suggested that, among the porcine AKR1Cs, AKR1C1 and AKR1C4 play important roles in steroid hormone metabolism in the reproductive tissues; however, their biological functions are still unknown. Herein, we report the biochemical properties of the two recombinant enzymes. Kinetic and product analyses of steroid specificity indicated that AKR1C1 is a multi-specific reductase, which acts as 3α-HSD for 3-keto-5β-dihydro-C₁₉/C₂₁-steroids, 3β-HSD for 3-keto-5α-dihydro-C₁₉-steroids including androstenone, 17β-HSD for 17-keto-C₁₉-steroids including estrone, and 20α-HSD for progesterone, showing K_m values of 0.5-11 µM. By contrast, AKR1C4 exhibited only 3α-HSD activity for 3-keto groups of 5α/β-dihydro-C₁₉-steroids, 5β-dihydro-C₂₁-steroids and bile acids (K_m: 1.0-1.9 µM). AKR1C1 and AKR1C4 also showed broad substrate specificity for nonsteroidal carbonyl compounds including endogenous 4-oxo-2-nonenal, 4-hydroxy-nonenal, acrolein, isocaproaldehyde, farnesal, isatin and methylglyoxal, of which 4-oxo-2-nonenal was reduced with the lowest K_m value of 0.9 µM. Moreover, AKR1C1 had the characteristic of reducing aliphatic ketones and all-trans-retinal. The enzymes were inhibited by flavonoids, synthetic estrogens, nonsteroidal anti-inflammatory drugs, triterpenoids and phenolphthalein, whereas only AKR1C4 was activated by bromosulfophthalein. These results suggest that AKR1C1 and AKR1C4 function as 3α/3β/17β/20α-HSD and 3α-HSD, respectively, in metabolism of steroid hormones and a sex pheromone androstenone, both of which also play roles in metabolism of nonsteroidal carbonyl compounds.

Differences in gene expression and variable splicing events of ovaries between large and small litter size in Chinese Xiang pigs

BACKGROUND

Although lots of quantitative trait loci (QTLs) and genes present roles in litter size of some breeds, the information might not make it clear for the huge diversity of reproductive capability in pig breeds. To elucidate the inherent mechanisms of heterogeneity of reproductive capability in litter size of Xiang pig, we performed transcriptome analysis for the expression profile in ovaries using RNA-seq method.

RESULTS

We identified 1,419 up-regulated and 1,376 down-regulated genes in Xiang pigs with large litter size. Among them, 1,010 differentially expressed genes (DEGs) were differently spliced between two groups with large or small litter sizes. Based on GO and KEGG analysis, numerous members of genes were gathered in ovarian steroidogenesis, steroid biosynthesis, oocyte maturation and reproduction processes.

CONCLUSIONS

Combined with gene biological function, twelve genes were found out that might be related with the reproductive capability of Xiang pig, of which, eleven genes were recognized as hub genes. These genes may play a role in promoting litter size by elevating steroid and peptide hormones supply through the ovary and facilitating the processes of ovulation and in vivo fertilization.

Correlation Networks Provide New Insights into the Architecture of Testicular Steroid Pathways in Pigs

Steroid metabolism is a fundamental process in the porcine testis to provide testosterone but also estrogens and androstenone, which are essential for the physiology of the boar. This study concerns boars at an early stage of puberty. Using a RT-qPCR approach, we showed that the transcriptional activities of several genes providing key enzymes involved in this metabolism (such as CYP11A1) are correlated. Surprisingly, HSD17B3, a key gene for testosterone production, was absent from this group. An additional weighted gene co-expression network analysis was performed on two large sets of mRNA-seq to identify co-expression modules. Of these modules, two containing either CYP11A1 or HSD17B3 were further analyzed. This comprehensive correlation meta-analysis identified a group of 85 genes with CYP11A1 as hub gene, but did not allow the characterization of a robust correlation network around HSD17B3. As the CYP11A1-group includes most of the genes involved in steroid synthesis pathways (including LHCGR encoding for the LH receptor), it may control the synthesis of most of the testicular steroids. The independent expression of HSD17B3 probably allows part of the production of testosterone to escape this control. This CYP11A1-group contained also INSL3 and AGT genes encoding a peptide hormone and an angiotensin peptide precursor, respectively.

Recent genetic advances on boar taint reduction as an alternative to castration: a review

Boar taint is an unpleasant odor in male pig meat, mainly caused by androstenone, skatole, and indole, which are deposited in the fat tissue. Piglet castration is the most common practice to prevent boar taint. However, castration is likely to be banished in a few years due to animal welfare concerns. Alternatives to castration, such as genetic selection, have been assessed. Androstenone and skatole have moderate to high heritability, which makes it feasible to select against these compounds. This review presents the latest results obtained on genetic selection against boar taint, on correlation with other traits, on differences in breeds, and on candidate genes related to boar taint. QTLs for androstenone and skatole have been reported mainly on chromosomes 6, 7, and 14. These chromosomes were reported to contain genes responsible for synthesis and degradation of androstenone and skatole. A myriad of work has been done to find markers or genes that can be used to select animals with lower boar taint. The selection against boar taint could decrease performance of some reproduction traits. However, a favorable response on production traits has been observed by selecting against boar taint. Selection results have shown that it is possible to reduce boar taint in few generations. In addition, modifications in diet and environment conditions could be associated with genetic selection to reduce boar taint. Nevertheless, costs to measure and select against boar taint should be rewarded with incentives from the market; otherwise, it would be difficult to implement genetic selection.

Sub-fertility in crossbred bulls: deciphering testicular level transcriptomic alterations between zebu (Bos indicus) and crossbred (Bos taurus x Bos indicus) bulls

BACKGROUND

The incidence of poor semen quality and sub-fertility/infertility is higher in crossbred as compared to Zebu males. Several attempts have been made to understand the possible reasons for higher incidence of fertility problems in crossbred males, at sperm phenotype, proteome and genome level but with variable results. Since the quality of the ejaculated spermatozoa is determined by the testicular environment, assessing the testicular transcriptome between these breeds would help in identifying the possible mechanisms associated with infertility in crossbred bulls. However, such information is not available. We performed global transcriptomic profiling of testicular tissue from crossbred and Zebu bulls using Agilent Bos taurus GXP 8X60k AMADID: 29411 array. To the best of our knowledge, this is the first study comparing the testicular mRNAs between crossbred and Zebu bulls.

RESULTS

Out of the 14,419 transcripts detected in bovine testis, 1466 were differentially expressed between crossbred and Zebu bulls, in which 1038 were upregulated and 428 were downregulated in crossbred bulls. PI4KB and DPY19L2 genes, reported to be involved in sperm capacitation and acrosome formation respectively, were among the top 10 downregulated transcripts in crossbred testis. Genes involved in ubiquitination and proteolysis were upregulated, while genes involved in cell proliferation, stem cell differentiation, stem cell population maintenance, steroidogenesis, WNT signalling, protein localization to plasma membrane, endocannabinoid signalling, heparin binding, cAMP metabolism and GABA receptor activity were downregulated in crossbred testis. Among the 10 genes validated using qPCR, expression of CCNYL, SOX2, MSMB, SPATA7, TNP1, TNP2 and CRISP2 followed the same trend as observed in microarray analysis with SPATA7 being significantly downregulated and transition proteins (TNP1, TNP2) being significantly upregulated in crossbred bulls.

CONCLUSIONS

Abundant proteolysis by ubiquitination and downregulation of WNT signaling, cell proliferation, differentiation and steroidogenesis might be associated with higher incidence of poor semen quality and/or sub-fertility/infertility in crossbred bulls as compared to Zebu bulls. Downregulation of SPATA7 (Spermatogenesis Associated 7) and upregulation of transition proteins (TNP1 and TNP2) in crossbred bull testis might be associated with impaired spermatogenesis processes including improper chromatin compaction in crossbred bulls.

NRDR inhibits estradiol synthesis and is associated with changes in reproductive traits in pigs

Cumulus cells secreting steroid hormones have important functions in oocyte development. Several members of the short-chain dehydrogenase/reductase (SDR) family are critical to the biosynthesis of steroid hormones. NADPH-dependent retinol dehydrogenase/reductase ( NRDR), a member of the SDR superfamily, is overexpressed in pig breeds that also show high levels of androstenone. However, the potential functions and regulatory mechanisms of NRDR in pig ovaries have not been reported to date. The present study demonstrated that NRDR is highly expressed in pig ovaries and is specifically located in cumulus granulosa cells. Functional studies showed that NRDR inhibition increased estradiol synthesis. Both pregnant mare serum gonadotropin and human chorionic gonadotropin downregulated the expression of NRDR in pig cumulus granulosa cells. When the relationship between reproductive traits and single-nucleotide polymorphisms (SNPs) of the NRDR gene was examined, we found that two SNPs affected reproductive traits. SNP rs701332503 was significantly associated with a decrease in the total number of piglets born during multiparity, and rs326982309 was significantly associated with an increase in the average birth weight during primiparity. Thus, NRDR has an important role in steroid hormone biosynthesis in cumulus granulosa cells, and NRDR SNPs are associated with changes in porcine reproduction traits.

SULFATION PATHWAYS: Expression of SULT2A1, SULT2B1 and HSD3B1 in the porcine testis and epididymis

In the porcine testis, in addition to estrogen sulfates, the formation of numerous sulfonated neutral hydroxysteroids has been observed. However, their functions and the underlying synthetic pathways are still widely unclear. To obtain further information on their formation in postpubertal boars, the expression of sulfotransferases considered relevant for neutral hydroxysteroids (SULT2A1, SULT2B1) was investigated in the testis and defined segments of the epididymis applying real-time RT-qPCR, Western blot and immunohistochemistry (IHC). Sulfotransferase activities were assessed in tissue homogenates or cytosolic preparations applying dehydroepiandrosterone and pregnenolone as substrates. A high SULT2A1 expression was confirmed in the testis and localized in Leydig cells by IHC. In the epididymis, SULT2A1 expression was virtually confined to the body. SULT2B1 expression was absent or low in the testis but increased significantly along the epididymis. Immunohistochemical observations indicate that both enzymes are secreted into the ductal lumen via an apocrine mechanism. The results from the characterization of expression patterns and activity measurements suggest that SULT2A1 is the prevailing enzyme for the sulfonation of hydroxysteroids in the testis, whereas SULT2B1 may catalyze the formation of sterol sulfates in the epididymis. In order to obtain information on the overall steroidogenic capacity of the porcine epididymis, the expression of important steroidogenic enzymes (CYP11A1, CYP17A1, CYP19, HSD3B1, HSD17B3, SRD5A2) was monitored in the defined epididymal segments applying real-time RT-qPCR. Surprisingly, in addition to a high expression of SRD5A2 in the epididymal head, a substantial expression of HSD3B1 was detected, which increased along the organ.

Associations of the Polymorphisms in DHRS4, SERPING1, and APOR Genes with Postmortem pH in Berkshire Pigs

Postmortem pH is a main factor influencing the meat quality in pigs. This study investigated the association of postmortem pH with single-nucleotide polymorphisms (SNPs) in the fourth member of the short-chain dehydrogenase/reductase family (DHRS4), the first member of serpin peptidase inhibitor, clade G (complement inhibitor) (SERPING1), and the apolipoprotein R precursor (APOR) genes in Berkshire pigs. The study included 437 pigs, and genotyping was conducted using the GoldenGate Assay (Illumina, San Diego, CA, USA). DHRS4, SERPING1, and APOR polymorphisms were significantly associated with pH₄₅ or pH₂₄ (p < 0.05). SERPING1 was also statistically significantly associated with water holding capacity (p < 0.05), which is closely associated with postmortem pH. These results suggest that SNPs in the DHRS4, SERPING1, and APOR genes have potential for use as genetic markers for the meat quality in pigs.

A global analysis of CNVs in swine using whole genome sequence data and association analysis with fatty acid composition and growth traits

Copy number variations (CNVs) are important genetic variants complementary to SNPs, and can be considered as biomarkers for some economically important traits in domestic animals. In the present study, a genomic analysis of porcine CNVs based on next-generation sequencing data was carried out to identify CNVs segregating in an Iberian x Landrace backcross population and study their association with fatty acid composition and growth-related traits. A total of 1,279 CNVs, including duplications and deletions, were detected, ranging from 106 to 235 CNVs across samples, with an average of 183 CNVs per sample. Moreover, we detected 540 CNV regions (CNVRs) containing 245 genes. Functional annotation suggested that these genes possess a great variety of molecular functions and may play a role in production traits in commercial breeds. Some of the identified CNVRs contained relevant functional genes (e.g., CLCA4, CYP4X1, GPAT2, MOGAT2, PLA2G2A and PRKG1, among others). The variation in copy number of four of them (CLCA4, GPAT2, MOGAT2 and PRKG1) was validated in 150 BC1_LD (25% Iberian and 75% Landrace) animals by qPCR. Additionally, their contribution regarding backfat and intramuscular fatty acid composition and growth–related traits was analyzed. Statistically significant associations were obtained for CNVR112 (GPAT2) for the C18:2(n-6)/C18:3(n-3) ratio in backfat and carcass length, among others. Notably, GPATs are enzymes that catalyze the first step in the biosynthesis of both triglycerides and glycerophospholipids, suggesting that this CNVR may contribute to genetic variation in fatty acid composition and growth traits. These findings provide useful genomic information to facilitate the further identification of trait-related CNVRs affecting economically important traits in pigs.

Transcript levels of genes implicated in steroidogenesis in the testes and fat tissue in relation to androstenone accumulation in fat of pubertal pigs

The present study was performed to measure messenger RNA levels of steroidogenic enzymes in testes and fat tissue and determine whether they are related to fat androstenone level. Real-time polymerase chain reaction experiments were performed on 26 testes and 12 adipose tissue samples from pubertal boars using 21 genes. The absence of significant correlations between fat androstenone and the transcriptional activity of the SRD5A2 and SRD5A3 genes but the high correlation coefficient with that of the SRD5A1 gene (r = 0.62, P < 0.05) suggests that the enzyme coded by SRD5A1 is mainly responsible for the last step of androstenone synthesis. The testicular transcriptional activities of CYP17, CYP11A1, CYP19A, AKR1C-pig6, SRD5A1, LHCGR, and AR were significantly correlated. Only transcriptional levels of CYP17, CYP11A1, CYP19A, SRD5A1, and AKR1C-pig6 were correlated with the fat concentration of androstenone (0.57 < r < 0.70, P < 0.05) confirming that the amount of androstenone stored in fat is related to the production in testes of androstenone and more generally to all sex steroids. Altogether, our data are in favor of a preponderant role of AKR1C-pig6 instead of HSD17B3 for testicular synthesis of steroids. Concerning fat tissue, our data do not support a significant de novo biosynthesis of steroids in porcine adipose tissues. The presence of transcripts coding for steroid enzymes, especially those of AKR1C-pig6, suggests that steroids can be transformed. None of transcript abundance was related to androstenone accumulation (P > 0.1). Therefore, steroids synthesized elsewhere can be transformed in fat tissue but synthesis of androstenone is unlikely.

Skatole metabolites in urine as a biological marker of pigs with enhanced hepatic metabolism

Boar taint is a quality defect in meat, related to accumulation of skatole and androstenone in male pigs. The levels of skatole and its main metabolites in plasma and urine samples were measured with a validated liquid chromatography-MS method and related to activity of hepatic cytochrome P450 (CYP450) in order to identify 'fast metabolizing' pigs. Urine (n=46), blood (n=12), liver (n=25) and adipose tissue (n=46) were sampled from a total of 46 entire male pigs. Skatole levels in fat were negatively correlated to CYP2E1 activity and positively to 3-hydroxy-3-methyloxindole (HMOI), indole-3-carboxylic acid (ICA) and 2-aminoacetophenone in urine. HMOI and ICA levels in urine were the best predictors of high skatole levels in fat. In summary, the present study provided further evidence for the key role of CYP2E1 in skatole metabolism and suggested that measurement of HMOI and/or ICA in urine might provide information about skatole levels in live pigs.

Cytochrome b5 modulates multiple reactions in steroidogenesis by diverse mechanisms

Cytochrome b5 (cyt-b5) is a relatively small haemoprotein which plays an important role in the regulation of mammalian steroidogenesis. This unique protein has the ability to modulate the activity of key steroidogenic enzymes via a number of diverse reaction mechanisms. Cyt-b5 can augment the 17,20-lyase activity of CYP17A1 by promoting the interaction of CYP17A1 and POR; enhance the 16-ene-synthase activity of CYP17A1 by acting as an electron donor; and enhance the activity of 3βHSD by increasing the affinity of 3βHSD for its cofactor NAD(+). We review the modulation of CYP17A1 and 3βHSD activity by cyt-b5 and discuss the reaction mechanisms associated with each activity. The physiological importance of cyt-b5 in regulating mammalian steroidogenesis is presented and the impact of inactivating cyt-b5 mutations are reviewed. This article is part of a Special Issue entitled 'Steroid/Sterol signaling'.

Preliminary study of FMO1, FMO5, CYP21, ESR1, PLIN2 and SULT2A1 as candidate gene for compounds related to boar taint

An association study between polymorphisms of six genes and boar taint related compounds androstenone, skatole and indole was performed in a boar population (n=370). Significant association (P<0.05) was detected for SNP of FMO5 (g.494A>G) with all boar taint compounds, SNP of CYP21 (g.3911T>C) with skatole and indole, and SNP of ESR1 (g.672C>T) with androstenone and indole. mRNA expression of CYP21 and ESR1 was higher in CAB (castrated boar) compared to non-castrated boars; whereas, the expression of FMO5 and ESR1 was higher in LBT (low boar taint) compared to HBT (high boar taint) in liver tissue. FMO5, CYP21 and ESR1 proteins were less detectable in HBT compared with LBT and CAB in liver tissues. These findings suggest that FMO5, CYP21 and ESR1 gene variants might have effects on the boar taint compounds.

Identification of gene co-expression clusters in liver tissues from multiple porcine populations with high and low backfat androstenone phenotype

Background

Boar taint is principally caused by accumulation of androstenone and skatole in adipose tissues. Studies have shown high heritability estimates for androstenone whereas skatole production is mainly dependent on nutritional factors. Androstenone is a lipophilic steroid mainly metabolized in liver. Majority of the studies on hepatic androstenone metabolism focus only on a single breed and very few studies account for population similarities/differences in gene expression patterns. In this work, we concentrated on population similarities in gene expression to identify the common genes involved in hepatic androstenone metabolism of multiple pig populations. Based on androstenone measurements, publicly available gene expression datasets from three porcine populations were compiled into either low or high androstenone dataset. Gene expression correlation coefficients from these datasets were converted to rank ratios and joint probabilities of these rank ratios were used to generate dataset specific co-expression clusters. Finally, these networks were clustered using a graph clustering technique.

Results

Cluster analysis identified a number of statistically significant co-expression clusters in the dataset. Further enrichment analysis of these clusters showed that one of the clusters from low androstenone dataset was highly enriched for xenobiotic, drug, cholesterol and lipid metabolism and cytochrome P450 associated metabolism of drugs and xenobiotics. Literature references revealed that a number of genes in this cluster were involved in phase I and phase II metabolism. Physical and functional similarity assessment showed that the members of this cluster were dispersed across multiple clusters in high androstenone dataset, possibly indicating a weak co-expression of these genes in high androstenone dataset.

Conclusions

Based on these results we hypothesize that majority of the genes in this cluster forms a signature co-expression cluster in low androstenone dataset in our experiment and that majority of the members of this cluster might be responsible for hepatic androstenone metabolism across all the three populations used in our study. We propose these results as a background work towards understanding breed similarities in hepatic androstenone metabolism. Additional large scale experiments using data from multiple porcine breeds are necessary to validate these findings.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-014-0158-8) contains supplementary material, which is available to authorized users.

Pathways and genes involved in steroid hormone metabolism in male pigs: a review and update

This paper reviews state-of-the-art knowledge on steroid biosynthesis pathways in the pig and provides an updated characterization of the porcine genes involved in these pathways with particular focus on androgens, estrogens, and 16-androstenes. At least 21 different enzymes appear to be involved in these pathways in porcine tissues together with at least five cofactors. Until now, data on several porcine genes were scarce or confusing. We characterized the complete genomic and transcript sequences of the single porcine CYP11B gene. We analyzed the porcine AKR1 gene cluster and identified four AKR1C, one AKR1C like genes and one AKR1E2 gene. We provide evidence that porcine AKR1C genes are not orthologous to human AKR1C. A new nomenclature is thus needed for this gene family in the pig. Thirty-two genes are now described: transcript (30+2 characterized in this study) and genomic (complete: 18+1 and partial: 12+1) sequences are identified. However, despite increasing knowledge on steroid metabolism in the pig, there is still no explanation of why porcine testes can produce androstenone and epiandrosterone, but not dihydrotestosterone (DHT), which is also a reduced steroid.

On the relationship between an Asian haplotype on chromosome 6 that reduces androstenone levels in boars and the differential expression of SULT2A1 in the testis

Background

Androstenone is one of the major compounds responsible for boar taint, a pronounced urine-like odor produced when cooking boar meat. Several studies have identified quantitative trait loci (QTL) for androstenone level on Sus scrofa chromosome (SSC) 6. For one of the candidate genes in the region SULT2A1, a difference in expression levels in the testis has been shown at the protein and RNA level.

Results

Haplotypes were predicted for the QTL region and their effects were estimated showing that haplotype 1 was consistently related with a lower level, and haplotype 2 with a higher level of androstenone. A recombinant haplotype allowed us to narrow down the QTL region from 3.75 Mbp to 1.94 Mbp. An RNA-seq analysis of the liver and testis revealed six genes that were differentially expressed between homozygotes of haplotypes 1 and 2. Genomic sequences of these differentially expressed genes were checked for variations within potential regulatory regions. We identified one variant located within a CpG island that could affect expression of SULT2A1 gene. An allele-specific expression analysis in the testis did not show differential expression between the alleles of SULT2A1 located on the different haplotypes in heterozygous animals. However a synonymous mutation C166T (SSC6: 49,117,861 bp in Sscrofa 10.2; C/T) was identified within the exon 2 of SULT2A1 for which the haplotype 2 only had the C allele which was higher expressed than the T allele, indicating haplotype-independent allelic-imbalanced expression between the two alleles. A phylogenetic analysis for the 1.94 Mbp region revealed that haplotype 1, associated with low androstenone level, originated from Asia.

Conclusions

Differential expression could be observed for six genes by RNA-seq analysis. No difference in the ratio of C:T expression of SULT2A1 for the haplotypes was found by the allele-specific expression analysis, however, a difference in expression between the C over T allele was found for a variation within SULT2A1, showing that the difference in androstenone levels between the haplotypes is not caused by the SNP in exon 2.

RNA deep sequencing reveals novel candidate genes and polymorphisms in boar testis and liver tissues with divergent androstenone levels

Boar taint is an unpleasant smell and taste of pork meat derived from some entire male pigs. The main causes of boar taint are the two compounds androstenone (5α-androst-16-en-3-one) and skatole (3-methylindole). It is crucial to understand the genetic mechanism of boar taint to select pigs for lower androstenone levels and thus reduce boar taint. The aim of the present study was to investigate transcriptome differences in boar testis and liver tissues with divergent androstenone levels using RNA deep sequencing (RNA-Seq). The total number of reads produced for each testis and liver sample ranged from 13,221,550 to 33,206,723 and 12,755,487 to 46,050,468, respectively. In testis samples 46 genes were differentially regulated whereas 25 genes showed differential expression in the liver. The fold change values ranged from −4.68 to 2.90 in testis samples and −2.86 to 3.89 in liver samples. Differentially regulated genes in high androstenone testis and liver samples were enriched in metabolic processes such as lipid metabolism, small molecule biochemistry and molecular transport. This study provides evidence for transcriptome profile and gene polymorphisms of boars with divergent androstenone level using RNA-Seq technology. Digital gene expression analysis identified candidate genes in flavin monooxygenease family, cytochrome P450 family and hydroxysteroid dehydrogenase family. Moreover, polymorphism and association analysis revealed mutation in IRG6, MX1, IFIT2, CYP7A1, FMO5 and KRT18 genes could be potential candidate markers for androstenone levels in boars. Further studies are required for proving the role of candidate genes to be used in genomic selection against boar taint in pig breeding programs.

Effect of dietary chicory on boar taint

Following preliminary screening and feeding trials on farms supplying a commercial abattoir, 360 entire male pigs were used to evaluate the effects of different percentages of chicory (Cichorium intybus L.) on levels of boar taint compounds and sensory aspects in backfat. Pigs were fed 0, 3, 6 or 9% chicory in the diet, 30 pigs being sampled at 3 different times: initially to measure basal levels of skatole and androstenone and after 1 and 2 weeks on the test diets. Cooked samples of backfat were presented to a trained sensory panel for "sniff" tests. Chicory fed at 9% for 2 weeks reduced skatole levels significantly (P<0.001), with 0.55 of pigs below 0.05 μg/g, typical of levels in castrated males. Abnormal odour scores were significantly lower for pigs in this group compared with 0% pigs (P<0.001), however, androstenone concentration was significantly higher in this group after the 2 week feeding period (P<0.005). Thus, feeding 9% chicory for 2 weeks was effective in reducing backfat skatole concentrations and abnormal odour scores of cooked fat but not androstenone concentration.

Genome-wide association scan and phased haplotype construction for quantitative trait loci affecting boar taint in three pig breeds

Background

Boar taint is the undesirable smell and taste of pork meat derived from some entire male pigs. The main causes of boar taint are the two compounds androstenone and skatole (3-methyl-indole). The steroid androstenone is a sex pheromone produced in the testis of the boars. Skatole is produced from tryptophan by bacteria in the intestine of the pigs. In many countries pigs are castrated as piglets to avoid boar taint, however, this is undesirable for animal welfare reasons. Genetic variations affecting the level of boar taint have previously been demonstrated in many breeds. In the study presented in this paper, markers and haplotypes, which can be applied to DNA-based selection schemes in order to reduce or eliminate the boar taint problem, are identified.

Results

Approximately 30,000 SNPs segregating in 923 boars from three Danish breeds; Duroc, Landrace, and Yorkshire, were used to conduct genome wide association studies of boar taint compounds. At 46 suggestive quantitative trait loci (QTL), 25 haplotypes and three single markers with effects were identified. Furthermore, 40% of the haplotypes mapped to previously identified regions. Haplotypes were also analysed for effects of slaughter weight and meat content. The most promising haplotype was identified on Sus scrofa chromosome 1. The gain in fixed effect of having this haplotype on level of androstenone in Landrace was identified to be high (1.279 μg/g). In addition, this haplotype explained 16.8% of the phenotypic variation within the trait. The haplotype was identified around the gene CYB5A which is known to have an indirect impact on the amount of androstenone. In addition to CYB5A, the genes SRD5A2, LOC100518755, and CYP21A2 are candidate genes for other haplotypes affecting androstenone, whereas, candidate genes for the indolic compounds were identified to be SULT1A1 and CYP2E1.

Conclusions

Despite the small sample size, a total of 25 haplotypes and three single markers were identified including genomic regions not previously reported. The haplotypes that were analysed showed large effects on trait level. However, little overlap of QTL between breeds was observed.

Large scale genome-wide association and LDLA mapping study identifies QTLs for boar taint and related sex steroids

Background

Boar taint is observed in a high proportion of uncastrated male pigs and is characterized by an unpleasant odor/flavor in cooked meat, primarily caused by elevated levels of androstenone and skatole. Androstenone is a steroid produced in the testis in parallel with biosynthesis of other sex steroids like testosterone and estrogens. This represents a challenge when performing selection against androstenone in breeding programs, without simultaneously decreasing levels of other steroids. The aim of this study was to use high-density genome wide association (GWA) in combination with linkage disequilibrium-linkage analysis (LDLA) to identify quantitative trait loci (QTL) associated with boar taint compounds and related sex steroids in commercial Landrace (n = 1,251) and Duroc (n = 918) breeds.

Results

Altogether, 14 genome wide significant (GWS) QTL regions for androstenone in subcutaneous fat were obtained from the LDLA study in Landrace and 14 GWS QTL regions in Duroc. LDLA analysis revealed that 7 of these QTL regions, located on SSC 1, 2, 3, 7 and 15, were obtained in both breeds. All 14 GWS androstenone QTLs in Landrace are also affecting the estrogens at chromosome wise significance (CWS) or GWS levels, while in Duroc, 3 of the 14 QTLs affect androstenone without affecting any of the estrogens. For skatole, 10 and 4 QTLs were GWS in the LDLA analysis for Landrace and Duroc respectively, with 4 of these detected in both breeds. The GWS QTLs for skatole obtained by LDLA are located at SSC 1, 5, 6, 7, 10, 11, 13 and 14.

Conclusion

This is the first report applying the Porcine 60 K SNP array for simultaneous analysis of boar taint compounds and related sex hormones, using both GWA and LDLA approaches. Several QTLs are involved in regulation of androstenone and skatole, and most of the QTLs for androstenone are also affecting the levels of estrogens. Seven QTLs for androstenone were detected in one breed and confirmed in the other, i.e. in an independent sample, although the majority of QTLs are breed specific. Most QTLs for skatole do not negatively affect other sex hormones and should be easier to implement into the breeding scheme.

Revealing genetic relationships between compounds affecting boar taint and reproduction in pigs

Boar taint is characterized by an unpleasant taste or odor in intact male pigs and is primarily attributed to increased concentrations of androstenone and skatole and to a lesser extent by increased indole. The boar taint compounds skatole and indole are produced by gut bacteria, metabolized in the liver, and stored in the fat tissue. Androstenone, on the other hand, is synthesized in the testis along with testosterone and estrogens, which are known to be important factors affecting fertility. The main goal of this study was to investigate the relationship between genetic factors involved in the primary boar taint compounds in an attempt to discover ways to reduce boar taint without decreasing fertility-related compounds. Heritabilities and genetic correlations between traits were estimated for compounds related to boar taint (androstenone, skatole, indole) and reproduction (testosterone, 17β-estradiol, and estrone sulfate). Heritabilities in the range of 0.47 to 0.67 were detected for androstenone concentrations in both fat and plasma, whereas those for skatole and indole were slightly less (0.27 to 0.41). The genetic correlations between androstenone in plasma and fat were extremely high (0.91 to 0.98) in Duroc and Landrace. In addition, genetic correlations between androstenone (both plasma and fat) and the other sex steroids (estrone sulfate, 17β-estradiol, and testosterone) were very high, in the range of 0.80 to 0.95. Furthermore, a genome-wide association study (GWA) and a combined linkage disequilibrium and linkage analysis (LDLA) were conducted on 1,533 purebred Landrace and 1,027 purebred Duroc to find genome regions involved in genetic control of the boar taint compounds androstenone, skatole, and indole, and sex hormones related to fertility traits. Up to 3,297 informative SNP markers were included for both breeds, including SNP from several boar taint candidate genes. From the GWA study, we found that altogether 27 regions were significant at a genome-wide level (P < 0.05) and an additional 7 regions were significant at a chromosomal level. From the LDLA study, 7 regions were significant on a genome-wide level and an additional 7 regions were significant at a chromosomal level. The most convincing associations were obtained in 6 regions affecting skatole and indole in fat on chromosomes 1, 2, 3, 7, 13, and 14, 1 region on chromosome 6 affecting androstenone in plasma only, and 5 regions on chromosomes 3, 4, 13, and 15 affecting androstenone, testosterone, and estrogens.

New investigations around CYP11A1 and its possible involvement in an androstenone QTL characterised in Large White pigs

BACKGROUND

Previously, in boars with extreme androstenone levels, differential expression of the CYP11A1 gene in the testes has been characterised. CYP11A1 is located in a region where a QTL influencing boar fat androstenone levels has been detected in a Large White pig population. Clarifying the role of CYP11A1 in boar taint is important because it catalyses the initial step of androstenone synthesis and also of steroid synthesis.

RESULTS

A genome-wide association study located CYP11A1 at approximately 1300 kb upstream from SNP H3GA0021967, defining the centre of the region containing the QTL for androstenone variation. In this study, we partially sequenced the CYP11A1 gene and identified several new single nucleotide polymorphisms (SNP) within it. Characterisation of one animal, heterozygous for CYP11A1 testicular expression but homozygous for a haplotype of a large region containing CYP11A1, revealed that variation of CYP11A1 expression is probably regulated by a mutation located downstream from the SNP H3GA0021967. We analysed CYP11A1 expression in LW families according to haplotypes of the QTL region's centre. Effects of haplotypes on CYP11A1 expression and on androstenone accumulation were not concordant.

CONCLUSION

This study shows that testicular expression of CYP11A1 is not solely responsible for the QTL influencing boar fat androstenone levels. As a conclusion, we propose to refute the hypothesis that a single mutation located near the centre of the QTL region could control androstenone accumulation in fat by regulating the CYP11A1 expression.

Expression levels of 25 genes in liver and testis located in a QTL region for androstenone on SSC7q1.2

A quantitative trait locus (QTL) for boar fat androstenone levels has been identified near the SSC7 centromere in a Large White × Meishan cross. Backcrosses were produced to isolate the Chinese haplotype in a European genetic background. The expression of 25 genes from the QTL region was studied in the testes and livers of 5-month-old backcross boars, with the aim of identifying the causal gene. Using Fluidigm, a new high-throughput technology, the expression of 25 genes was measured in a single real-time PCR experiment. This study found six significantly down-regulated genes (C6ORF106, C6ORF81, CLPS, SLC26A8, SRPK1 and MAPK14) in the testes of MS-LW backcross boars. However, according to current knowledge, none of the genes appear to be related to androstenone metabolism. In the livers, none of the genes were significantly up- or down-regulated, including TEAD3, which was previously designated as a possible candidate to explain this QTL.

Copy number variation in the porcine genome inferred from a 60 k SNP BeadChip

Background

Recent studies in pigs have detected copy number variants (CNVs) using the Comparative Genomic Hybridization technique in arrays designed to cover specific porcine chromosomes. The goal of this study was to identify CNV regions (CNVRs) in swine species based on whole genome SNP genotyping chips.

Results

We used predictions from three different programs (cnvPartition, PennCNV and GADA) to analyze data from the Porcine SNP60 BeadChip. A total of 49 CNVRs were identified in 55 animals from an Iberian x Landrace cross (IBMAP) according to three criteria: detected in at least two animals, contained three or more consecutive SNPs and recalled by at least two programs. Mendelian inheritance of CNVRs was confirmed in animals belonging to several generations of the IBMAP cross. Subsequently, a segregation analysis of these CNVRs was performed in 372 additional animals from the IBMAP cross and its distribution was studied in 133 unrelated pig samples from different geographical origins. Five out of seven analyzed CNVRs were validated by real time quantitative PCR, some of which coincide with well known examples of CNVs conserved across mammalian species.

Conclusions

Our results illustrate the usefulness of Porcine SNP60 BeadChip to detect CNVRs and show that structural variants can not be neglected when studying the genetic variability in this species.