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

Differential Regulation of Male-Hormones-Related Enhancers Revealed by Chromatin Accessibility and Transcriptional Profiles in Pig Liver

Surgical castration can effectively avoid boar taint and improve pork quality by removing the synthesis of androstenone in the testis, thereby reducing its deposition in adipose tissue. The expression of genes involved in testis-derived hormone metabolism was altered following surgical castration, but the upstream regulatory factors and underlying mechanism remain unclear. In this study, we systematically profiled chromatin accessibility and transcriptional dynamics in liver tissue of castrated and intact full-sibling Yorkshire pigs. First, we identified 897 differentially expressed genes and 6864 differential accessible regions (DARs) using RNA- and ATAC-seq. By integrating the RNA- and ATAC-seq results, 227 genes were identified, and a significant positive correlation was revealed between differential gene expression and the ATAC-seq signal. We constructed a transcription factor regulatory network after motif analysis of DARs and identified a candidate transcription factor (TF) SP1 that targeted the HSD3B1 gene, which was responsible for the metabolism of androstenone. Subsequently, we annotated DARs by incorporating H3K27ac ChIP-seq data, marking 2234 typical enhancers and 245 super enhancers involved in the regulation of all testis-derived hormones. Among these, four typical enhancers associated with HSD3B1 were identified. Furthermore, an in-depth investigation was conducted on the androstenone-related enhancers, and an androstenone-related mutation was identified in a newfound candidatetypical enhancer (andEN) with dual-luciferase assays. These findings provide further insights into how enhancers function as links between phenotypic and non-coding area variations. The discovery of upstream TF and enhancers of HSD3B1 contributes to understanding the regulatory networks of androstenone metabolism and provides an important foundation for improving pork quality.

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.

A Simple, Affordable, and Rapid Visual CRISPR-Based Field Test for Sex Determination of Earlier Porcine Embryos and Pork Products

Sex selection technologies have immensely impacted swine production globally. Conventional earlier embryo sex identification methods require professional technicians and sophisticated laboratory instruments. Rapid on-site gender identification of porcine embryos and pork products remains challenging. In this study, we developed a CRISPR/Cas12a-based fluorescence visualization point-of-care sex determination test that is rapid, accurate and easy to implement on-site. The CRISPR/Cas12a assay coupled with either the polymerase chain reaction (PCR) or loop-mediated isothermal amplification (LAMP) employs precisely designed primers and single-guide RNAs targeting the sex-determining region Y (SRY) and the zinc finger protein X-linked (ZFX) genes. PCR and LAMP amplicons were cleaved with the subsequent generation of fluorescing products detectable with portable blue and ultraviolet light transilluminators. Approximately two copies per microliter of the ZFX and SRY genes were detected using the RApid VIsual CRISPR (RAVI-CRISPR) assay. This method is a sensitive, inexpensive, versatile, and point-of-care test. The technology has other potential applications like determining the sex of diverse livestock species, detecting livestock disease-causing pathogens and evaluating the quality of meat products.

Hepatic transcriptome analysis identifies genes, polymorphisms and pathways involved in the fatty acids metabolism in sheep

Fatty acids (FA) in ruminants, especially unsaturated FA (USFA) have important impact in meat quality, nutritional value, and flavour quality of meat, and on consumer’s health. Identification of the genetic factors controlling the FA composition and metabolism is pivotal to select sheep that produce higher USFA and lower saturated (SFA) for the benefit of sheep industry and consumers. Therefore, this study was aimed to investigate the transcriptome profiling in the liver tissues collected from sheep with divergent USFA content in longissimus muscle using RNA deep-sequencing. From sheep (n = 100) population, liver tissues with higher (n = 3) and lower (n = 3) USFA content were analysed using Illumina HiSeq 2500. The total number of reads produced for each liver sample were ranged from 21.28 to 28.51 million with a median of 23.90 million. Approximately, 198 genes were differentially regulated with significance level of p-adjusted value <0.05. Among them, 100 genes were up-regulated, and 98 were down-regulated (p<0.01, FC>1.5) in the higher USFA group. A large proportion of key genes involved in FA biosynthesis, adipogenesis, fat deposition, and lipid metabolism were identified, such as APOA5, SLC25A30, GFPT1, LEPR, TGFBR2, FABP7, GSTCD, and CYP17A. Pathway analysis revealed that glycosaminoglycan biosynthesis- keratan sulfate, adipokine signaling, galactose metabolism, endocrine and other factors-regulating calcium metabolism, mineral metabolism, and PPAR signaling pathway were playing important regulatory roles in FA metabolism. Importantly, polymorphism and association analyses showed that mutation in APOA5, CFHR5, TGFBR2 and LEPR genes could be potential markers for the FA composition in sheep. These polymorphisms and transcriptome networks controlling the FA variation could be used as genetic markers for FA composition-related traits improvement. However, functional validation is required to confirm the effect of these SNPs in other sheep population in order to incorporate them in the sheep breeding program.

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.

Transcriptomic comparison of liver tissue between Anqing six-end-white pigs and Yorkshire pigs based on RNA sequencing

Chinese indigenous pig and Western commercial pig breeds show different patterns of lipid metabolism, fat deposition, and fatty acid composition; for these reasons, they have become vitally important models of energy metabolism and obesity in humans. To compare the mechanisms underlying lipid metabolism between Yorkshire pigs (lean type) and Anqing six-end-white pigs (obese type), the liver transcriptomes of six castrated boars with a body weight of approximately 100 kg (three Yorkshire and three Anqing) were analyzed by RNA-seq. The total number of reads produced for each liver sample ranged from 47.05 to 62.6 million. Among 362 differentially expressed genes, 142 were up-regulated and 220 were down-regulated in Anqing six-end-white pigs. Based on these data, 79 GO terms were significantly enriched. The top 10 (the 10 with lowest corrected P-value) significantly enriched GO terms were identified, including lipid metabolic process and carboxylic acid metabolic process. Pathway analysis revealed three significantly enriched KEGG pathways including PPAR signaling pathway, steroid hormone biosynthesis, and retinol metabolism. Based on protein–protein interaction networks, multiple genes responsible for lipid metabolism were identified, such as PCK1, PPARA, and CYP7A1, and these were considered promising candidate genes that could affect porcine liver lipid metabolism and fat deposition. Our results provide abundant transcriptomic information that will be useful for animal breeding and biomedical research.

Comparative analysis of differentially expressed genes between the ovaries from pregnant and nonpregnant goats using RNA-Seq

Background

A multitude of genes tightly regulate ovarian follicular development and hormone secretion. These complex and coordinated biological processes are altered during pregnancy. In order to further understand the regulatory role of these genes during pregnancy, it is important to screen the differentially expressed genes (DEGs) in the ovaries of pregnant and nonpregnant mammals. To detect the genes associated with the development of pregnancy in goats, DEGs from the ovaries from pregnant and nonpregnant Anhui white goats (pAWGs and nAWGs, respectively) were analyzed using RNA sequencing technology (RNA-Seq).

Results

In this study, 13,676,394 and 13,549,560 clean reads were generated from pAWGs and nAWGs, respectively, and 1724 DEGs were identified between the two libraries. Compared with nAWGs, 1033 genes were upregulated and 691 genes were downregulated in pAWGs, including PGR, PRLR, STAR and CYP19A1, which play important roles in goat reproduction. Gene Ontology analysis showed that the DEGs were enriched for 49 functional GO terms. Kyoto Encyclopedia of Genes and Genomes analysis revealed that 397 DEGs were significantly enriched in 13 pathways, including “cell cycle”, “cytokine–cytokine receptor interaction” and “steroid biosynthesis”, suggesting that the genes may be associated with cell cycle regulation, follicular development and hormone secretion to regulate the reproduction process. Additionally, quantitative real-time PCR was used to verify the reliability of the RNA-Seq data.

Conclusions

The data obtained in this work enrich the genetic resources of goat and provide a further understanding of the complex molecular regulatory mechanisms occurring during the development of pregnancy and reproduction in goats. The DEGs screened in this study may play an important role in follicular development and hormone secretion and they would provide scientific basis for related research in the future.

Electronic supplementary material

The online version of this article (10.1186/s40709-019-0095-9) contains supplementary material, which is available to authorized users.

The first comprehensive description of the expression profile of genes involved in differential body growth and the immune system of the Jeju Native Pig and miniature pig

Sus scrofa provides a major source of animal protein for humans as well as being an excellent biomedical model. This study was carried out to understand, in detail, the genetic and functional variants of Jeju Native Pigs and miniature pigs through differential expression profiling of the genes controlling their immune response, growth performance, and meat quality. The Illumina HiSeq 2000 platform was used for generating 1.3 billion 90 bp paired-end reads, which were mapped to the S. scrofa genome using TopHat2. A total of 2481 and 2768 genes were differentially expressed with 8-log changes in muscle and liver samples, respectively. Five hundred forty-eight genes in muscle and 642 genes in liver samples had BLAST matches within the non-redundant database. GO process and pathway analyses showed enhanced biological processes related to the extracellular structural organization and skeletal muscle cell differentiation in muscle tissue, whereas the liver tissue shares functions related to the inflammatory response. Herein, we identify inflammatory regulatory genes in miniature pigs and growth response genes in Jeju Native Pigs, information which can provide a stronger base for the selection of breeding stock and facilitate further in vitro and in vivo studies for therapeutic purposes.

Systems genomics study reveals expression quantitative trait loci, regulator genes and pathways associated with boar taint in pigs

Boar taint is an offensive odour and/or taste from a proportion of non-castrated male pigs caused by skatole and androstenone accumulation during sexual maturity. Castration is widely used to avoid boar taint but is currently under debate because of animal welfare concerns. This study aimed to identify expression quantitative trait loci (eQTLs) with potential effects on boar taint compounds to improve breeding possibilities for reduced boar taint. Danish Landrace male boars with low, medium and high genetic merit for skatole and human nose score (HNS) were slaughtered at ~100 kg. Gene expression profiles were obtained by RNA-Seq, and genotype data were obtained by an Illumina 60K Porcine SNP chip. Following quality control and filtering, 10,545 and 12,731 genes from liver and testis were included in the eQTL analysis, together with 20,827 SNP variants. A total of 205 and 109 single-tissue eQTLs associated with 102 and 58 unique genes were identified in liver and testis, respectively. By employing a multivariate Bayesian hierarchical model, 26 eQTLs were identified as significant multi-tissue eQTLs. The highest densities of eQTLs were found on pig chromosomes SSC12, SSC1, SSC13, SSC9 and SSC14. Functional characterisation of eQTLs revealed functions within regulation of androgen and the intracellular steroid hormone receptor signalling pathway and of xenobiotic metabolism by cytochrome P450 system and cellular response to oestradiol. A QTL enrichment test revealed 89 QTL traits curated by the Animal Genome PigQTL database to be significantly overlapped by the genomic coordinates of cis-acting eQTLs. Finally, a subset of 35 cis-acting eQTLs overlapped with known boar taint QTL traits. These eQTLs could be useful in the development of a DNA test for boar taint but careful monitoring of other overlapping QTL traits should be performed to avoid any negative consequences of selection.

Differential expression and co-expression gene networks reveal candidate biomarkers of boar taint in non-castrated pigs

Boar taint (BT) is an offensive odour or taste observed in pork from a proportion of non-castrated male pigs. Surgical castration is effective in avoiding BT, but animal welfare issues have created an incentive for alternatives such as genomic selection. In order to find candidate biomarkers, gene expression profiles were analysed from tissues of non-castrated pigs grouped by their genetic merit of BT. Differential expression analysis revealed substantial changes with log-transformed fold changes of liver and testis from -3.39 to 2.96 and -7.51 to 3.53, respectively. Co-expression network analysis revealed one module with a correlation of -0.27 in liver and three modules with correlations of 0.31, -0.44 and -0.49 in testis. Differential expression and co-expression analysis revealed candidate biomarkers with varying biological functions: phase I (COQ3, COX6C, CYP2J2, CYP2B6, ACOX2) and phase II metabolism (GSTO1, GSR, FMO3) of skatole and androstenone in liver to steroidgenesis (HSD17B7, HSD17B8, CYP27A1), regulation of steroidgenesis (STARD10, CYB5R3) and GnRH signalling (MAPK3, MAP2K2, MAP3K2) in testis. Overrepresented pathways included "Ribosome", "Protein export" and "Oxidative phosphorylation" in liver and "Steroid hormone biosynthesis" and "Gap junction" in testis. Future work should evaluate the biomarkers in large populations to ensure their usefulness in genomic selection programs.

Long noncoding RNA expression profile changes associated with dietary energy in the sheep testis during sexual maturation

Spermatogenesis can be affected by nutrition, which operates through normal physiological processes by changing the testicular mass and hormone levels profoundly. However, little is known regarding how testis development is regulated by long noncoding RNA (lncRNA). In this study, we investigated the effects of high-grain (HG) feeding on testis development during sexual maturation mediated by lncRNA. The HG diet group showed an increase in growth hormone (GH), insulin-like growth factor-1 (IGF-1) and testosterone (T) levels, and in the number of sperm in the seminiferous tubules compared with the hay-fed group (p  < 0.05). Moreover, we found 59 differentially expressed (DE) lncRNAs and 229 DE mRNAs in sheep testis between the two groups. qRT-PCR results of 20 randomly selected DE lncRNAs and mRNAs were also consistent with the RNA-seq data. Through functional enrichment analysis and lncRNA-mRNA interaction network analysis, we screened several lncRNAs that may be enriched for male reproduction such as spermatogenesis, sperm motility, steroid hormones, MAPK and ErbB signaling pathways. This study provides a first insight into the development of the testis with HG feeding in sheep and shows that these changes are associated with alterations in lncRNA expression.

Seasonal differences in the testicular transcriptome profile of free-living European beavers (Castor fiber L.) determined by the RNA-Seq method

The European beaver (Castor fiber L.) is an important free-living rodent that inhabits Eurasian temperate forests. Beavers are often referred to as ecosystem engineers because they create or change existing habitats, enhance biodiversity and prepare the environment for diverse plant and animal species. Beavers are protected in most European Union countries, but their genomic background remains unknown. In this study, gene expression patterns in beaver testes and the variations in genetic expression in breeding and non-breeding seasons were determined by high-throughput transcriptome sequencing. Paired-end sequencing in the Illumina HiSeq 2000 sequencer produced a total of 373.06 million of high-quality reads. De novo assembly of contigs yielded 130,741 unigenes with an average length of 1,369.3 nt, N50 value of 1,734, and average GC content of 46.51%. A comprehensive analysis of the testicular transcriptome revealed more than 26,000 highly expressed unigenes which exhibited the highest homology with Rattus norvegicus and Ictidomys tridecemlineatus genomes. More than 8,000 highly expressed genes were found to be involved in fundamental biological processes, cellular components or molecular pathways. The study also revealed 42 genes whose regulation differed between breeding and non-breeding seasons. During the non-breeding period, the expression of 37 genes was up-regulated, and the expression of 5 genes was down-regulated relative to the breeding season. The identified genes encode molecules which are involved in signaling transduction, DNA repair, stress responses, inflammatory processes, metabolism and steroidogenesis. Our results pave the way for further research into season-dependent variations in beaver testes.

Genome-Wide Analysis and Functional Characterization of the Polyadenylation Site in Pigs Using RNAseq Data

Polyadenylation, a critical step in the production of mature mRNA for translation in most eukaryotes, involves cleavage and poly(A) tail addition at the 3′ end of mRNAs at the polyadenylation site (PAS). Sometimes, one gene can have more than one PAS, which can produce the alternative polyadenylation (APA) phenomenon and affect the stability, localization and translation of the mRNA. In this study, we discovered 28,363 PASs using pig RNAseq data, with 13,033 located in 7,403 genes. Among the genes, 41% were identified to have more than one PAS. PAS distribution analysis indicated that the PAS position was highly variable in genes. Additionally, the analysis of RNAseq data from the liver and testis showed a difference in their PAS number and usage. RT-PCR and qRT-PCR were performed to confirm our findings by detecting the expression of 3′UTR isoforms for five candidate genes. The analysis of RNAseq data under a different androstenone level and salmonella inoculation indicated that the functional usage of PAS might participate in the immune response and may be related to the androstenone level in pigs. This study provides new insights into pig PAS and facilitates further functional research of PAS.

Multi-omic data integration and analysis using systems genomics approaches: methods and applications in animal production, health and welfare

In the past years, there has been a remarkable development of high-throughput omics (HTO) technologies such as genomics, epigenomics, transcriptomics, proteomics and metabolomics across all facets of biology. This has spearheaded the progress of the systems biology era, including applications on animal production and health traits. However, notwithstanding these new HTO technologies, there remains an emerging challenge in data analysis. On the one hand, different HTO technologies judged on their own merit are appropriate for the identification of disease-causing genes, biomarkers for prevention and drug targets for the treatment of diseases and for individualized genomic predictions of performance or disease risks. On the other hand, integration of multi-omic data and joint modelling and analyses are very powerful and accurate to understand the systems biology of healthy and sustainable production of animals. We present an overview of current and emerging HTO technologies each with a focus on their applications in animal and veterinary sciences before introducing an integrative systems genomics framework for analysing and integrating multi-omic data towards improved animal production, health and welfare. We conclude that there are big challenges in multi-omic data integration, modelling and systems-level analyses, particularly with the fast emerging HTO technologies. We highlight existing and emerging systems genomics approaches and discuss how they contribute to our understanding of the biology of complex traits or diseases and holistic improvement of production performance, disease resistance and welfare.

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.

Gonadal Transcriptome Analysis in Sterile Double Haploid Japanese Flounder

Sterility is a serious problem that can affect all bionts. In teleosts, double haploids (DHs) induced by mitogynogenesis are often sterile. This sterility severely restricts the further application of DHs for production of clones, genetic analysis, and breeding. However, sterile DH individuals are good source materials for investigation of the molecular mechanisms of gonad development, especially for studies into the role of genes that are indispensable for fish reproduction. Here, we used the Illumina sequencing platform to analyze the transcriptome of sterile female DH Japanese flounder in order to identify major genes that cause sterility and to provide a molecular basis for an intensive study of gonadal development in teleosts. Through sequencing, assembly, and annotation, we obtained 52,474 contigs and found that 60.7% of these shared homologies with existing sequences. A total of 1225 differentially expressed unigenes were found, including 492 upregulated and 733 downregulated genes. Gene Ontology and KEGG analyses showed that genes showing significant upregulation, such as CYP11A1, CYP11B2, CYP17, CYP21, HSD3β, bcl2l1, and PRLR, principally correlated with sterol metabolic process, steroid biosynthetic process, and the Jak-stat signaling pathway. The significantly downregulated genes were primarily associated with immune response, antigen processing and presentation, cytokine–cytokine receptor interaction, and protein digestion and absorption. Using a co-expression network analysis, we conducted a comprehensive comparison of gene expression in the gonads of fertile and sterile female DH Japanese flounder. Identification of genes showing significantly different expression will provide further insights into DH reproductive dysfunction and oocyte maturation processes in teleosts.

Identification of Genes Related to Growth and Lipid Deposition from Transcriptome Profiles of Pig Muscle Tissue

Transcriptome profiles established using high-throughput sequencing can be effectively used for screening genome-wide differentially expressed genes (DEGs). RNA sequences (from RNA-seq) and microRNA sequences (from miRNA-seq) from the tissues of longissimus dorsi muscle of two indigenous Chinese pig breeds (Diannan Small-ear pig [DSP] and Tibetan pig [TP]) and two introduced pig breeds (Landrace [LL] and Yorkshire [YY]) were examined using HiSeq 2000 to identify and compare the differential expression of functional genes related to muscle growth and lipid deposition. We obtained 27.18 G clean data through the RNA-seq and detected that 18,208 genes were positively expressed and 14,633 of them were co-expressed in the muscle tissues of the four samples. In all, 315 DEGs were found between the Chinese pig group and the introduced pig group, 240 of which were enriched with functional annotations from the David database and significantly enriched in 27 Gene Ontology (GO) terms that were mainly associated with muscle fiber contraction, cadmium ion binding, response to organic substance and contractile fiber part. Based on functional annotation, we identified 85 DEGs related to growth traits that were mainly involved in muscle tissue development, muscle system process, regulation of cell development, and growth factor binding, and 27 DEGs related to lipid deposition that were mainly involved in lipid metabolic process and fatty acid biosynthetic process. With miRNA-seq, we obtained 23.78 M reads and 320 positively expressed miRNAs from muscle tissues, including 271 known pig miRNAs and 49 novel miRNAs. In those 271 known miRNAs, 20 were higher and 10 lower expressed in DSP-TP than in LL-YY. The target genes of the 30 miRNAs were mainly participated in MAPK, GnRH, insulin and Calcium signaling pathway and others involved cell development, growth and proliferation, etc. Combining the DEGs and the differentially expressed (DE) miRNAs, we drafted a network of 46 genes and 18 miRNAs for regulating muscle growth and a network of 15 genes and 16 miRNAs for regulating lipid deposition. We identified that CAV2, MYOZ2, FRZB, miR-29b, miR-122, miR-145-5p and miR-let-7c, etc, were key genes or miRNAs regulating muscle growth, and FASN, SCD, ADORA1, miR-4332, miR-182, miR-92b-3p, miR-let-7a and miR-let-7e, etc, were key genes or miRNAs regulating lipid deposition. The quantitative expressions of eight DEGs and seven DE miRNAs measured with real-time PCR certified that the results of differential expression genes or miRNAs were reliable. Thus, 18,208 genes and 320 miRNAs were positively expressed in porcine longissimus dorsi muscle. We obtained 85 genes and 18 miRNAs related to muscle growth and 27 genes and 16 miRNAs related to lipid deposition, which provided new insights into molecular mechanism of the economical traits in pig.

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.

Review on porcine endogenous retrovirus detection assays--impact on quality and safety of xenotransplants

Xenotransplantation of porcine organs, tissues, and cells inherits a risk for xenozoonotic infections. Viable tissues and cells intended for transplantation have to be considered as potentially contaminated non-sterile products. The demands on microbial testing, based on the regulatory requirements, are often challenging due to a restricted shelf life or the complexity of the product itself. In Europe, the regulatory framework for xenogeneic cell therapy is based on the advanced therapy medicinal products (ATMP) regulation (2007), the EMA CHMP Guideline on xenogeneic cell-based medicinal products (2009), as well as the WHO and Council of Europe recommendations. In the USA, FDA guidance for industry (2003) regulates the use of xenotransplants. To comply with the regulations, validated test methods need to be established that reveal the microbial status of a transplant within its given shelf life, complemented by strictly defined action alert limits and supported by breeding in specific pathogen-free (SPF) facilities. In this review, we focus on assays for the detection of the porcine endogenous retroviruses PERV-A/-B/-C, which exhibit highly polymorphic proviral loci in pig genomes. PERVs are transmitted vertically and cannot be completely eliminated by breeding or gene knock out technology. PERVs entail a public health concern that will persist even if no evidence of PERV infection of xenotransplant recipients in vivo has been revealed yet. Nevertheless, infectious risks must be minimized by full assessment of pigs as donors by combining different molecular screening assays for sensitive and specific detection as well as a functional analysis of the infectivity of PERV including an adequate monitoring of recipients.

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.

Fatty acid composition of subcutaneous adipose tissue from entire male pigs with extremely divergent levels of boar taint compounds--an exploratory study

This exploratory study investigated the variability of fatty acid composition in entire male pigs with extremely divergent levels of boar taint compounds. Fatty acids were quantified in back fat samples from 20 selected carcasses of Pietrain*F1 sired boars (average carcass weight 84 kg) with extremely low (LL) or extremely high (HH) levels of androstenone, skatole, and indole. Concentrations of polyunsaturated fatty acids (PUFA) were significantly (p<0.05) increased in LL boars (23.4%) compared to HH boars (19.7%). This was mainly due to increased levels of linoleic acid (C18:2 n-6) and α-linolenic acid (C18:3 n-3). Correspondingly, unsaturated fatty acids (SFA) were significantly lower (p<0.05) in LL boars (35.2%) compared to HH boars (37.7%). The findings are discussed with respect to potential effects on flavor formation in boar fat and meat. Further research is needed to study the gender specificity and the interplay of the synthesis and the metabolism of steroids, lipids, and the clearance of skatole in pigs.

Pathway based analysis of genes and interactions influencing porcine testis samples from boars with divergent androstenone content in back fat

One of the primary factors contributing to boar taint is the level of androstenone in porcine adipose tissues. A majority of the studies performed to identify candidate biomarkers for the synthesis of androstenone in testis tissues follow a reductionist approach, identifying and studying the effect of biomarkers individually. Although these studies provide detailed information on individual biomarkers, a global picture of changes in metabolic pathways that lead to the difference in androstenone synthesis is still missing. The aim of this work was to identify major pathways and interactions influencing steroid hormone synthesis and androstenone biosynthesis using an integrative approach to provide a bird's eye view of the factors causing difference in steroidogenesis and androstenone biosynthesis. For this purpose, we followed an analysis procedure merging together gene expression data from boars with divergent levels of androstenone and pathway mapping and interaction network retrieved from KEGG database. The interaction networks were weighted with Pearson correlation coefficients calculated from gene expression data and significant interactions and enriched pathways were identified based on these networks. The results show that 1,023 interactions were significant for high and low androstenone animals and that a total of 92 pathways were enriched for significant interactions. Although published articles show that a number of these enriched pathways were activated as a result of downstream signaling of steroid hormones, we speculate that the significant interactions in pathways such as glutathione metabolism, sphingolipid metabolism, fatty acid metabolism and significant interactions in cAMP-PKA/PKC signaling might be the key factors determining the difference in steroidogenesis and androstenone biosynthesis between boars with divergent androstenone levels in our study. The results and assumptions presented in this study are from an in-silico analysis done at the gene expression level and further laboratory experiments at genomic, proteomic or metabolomic level are necessary to validate these findings.

Identification of genes regulating growth and fatness traits in pig through hypothalamic transcriptome analysis

Previous studies on Iberian × Landrace (IBMAP) pig intercrosses have enabled the identification of several quantitative trait locus (QTL) regions related to growth and fatness traits; however, the genetic variation underlying those QTLs are still unknown. These traits are not only relevant because of their impact on economically important production traits, but also because pig constitutes a widely studied animal model for human obesity and obesity-related diseases. The hypothalamus is the main gland regulating growth, food intake, and fat accumulation. Therefore, the aim of this work was to identify genes and/or gene transcripts involved in the determination of growth and fatness in pig by a comparison of the whole hypothalamic transcriptome (RNA-Seq) in two groups of phenotypically divergent IBMAP pigs. Around 16,000 of the ∼25.010 annotated genes were expressed in these hypothalamic samples, with most of them showing intermediate expression levels. Functional analyses supported the key role of the hypothalamus in the regulation of growth, fat accumulation, and energy expenditure. Moreover, 58,927 potentially new isoforms were detected. More than 250 differentially expressed genes and novel transcript isoforms were identified between the two groups of pigs. Twenty-one DE genes/transcripts that colocalized in previously identified QTL regions and/or whose biological functions are related to the traits of interest were explored in more detail. Additionally, the transcription factors potentially regulating these genes and the subjacent networks and pathways were also analyzed. This study allows us to propose strong candidate genes for growth and fatness based on expression patterns, genomic location, and network interactions.

Identification of the novel candidate genes and variants in boar liver tissues with divergent skatole levels using RNA deep sequencing

Boar taint is the unpleasant odour of meat derived from non-castrated male pigs, caused by the accumulation of androstenone and skatole in fat. Skatole is a tryptophan metabolite produced by intestinal bacteria in gut and catabolised in liver. Since boar taint affects consumer's preference, the aim of this study was to perform transcriptome profiling in liver of boars with divergent skatole levels in backfat by using RNA-Seq. The total number of reads produced for each liver sample ranged from 11.8 to 39.0 million. Approximately 448 genes were differentially regulated (p-adjusted <0.05). Among them, 383 genes were up-regulated in higher skatole group and 65 were down-regulated (p<0.01, FC>1.5). Differentially regulated genes in the high skatole liver samples were enriched in metabolic processes such as small molecule biochemistry, protein synthesis, lipid and amino acid metabolism. Pathway analysis identified the remodeling of epithelial adherens junction and TCA cycle as the most dominant pathways which may play important roles in skatole metabolism. Differential gene expression analysis identified candidate genes in ATP synthesis, cytochrome P450, keratin, phosphoglucomutase, isocitrate dehydrogenase and solute carrier family. Additionally, polymorphism and association analysis revealed that mutations in ATP5B, KRT8, PGM1, SLC22A7 and IDH1 genes could be potential markers for skatole levels in boars. Furthermore, expression analysis of exon usage of three genes (ATP5B, KRT8 and PGM1) revealed significant differential expression of exons of these genes in different skatole levels. These polymorphisms and exon expression differences may have impacts on the gene activity ultimately leading to skatole variation and could be used as genetic marker for boar taint related traits. However, further validation is required to confirm the effect of these genetic markers in other pig populations in order to be used in genomic selection against boar taint in pig breeding programs.