Asep Gunawan

Association study and expression analysis of porcine ESR1 as a candidate gene for boar fertility and sperm quality

Male fertility is impaired through the lack of ESR1 (Estrogen Receptor 1) but little is known about the ESR1 roles in boar spermatogenesis and fertility. Therefore, this research was aimed at investigating the association with sperm quality and boar fertility traits in a total of 300 boars both from purebred Pietrain and Pietrain × Hampshire crosses. A SNP in coding region of ESR1g.672C>T in exon 1 was associated with sperm motility (P<0.05) and plasma droplet rate (P<0.01) while the polymorphism in non-coding region of ESR1g.35756T>C in inton 1 was associated with non-return rate (P<0.05). Furthermore, to analyse the mRNA and protein expression of ESR1 in boar reproductive tissues, a total of six boars were divided into two groups [Group I (G-I) and Group II (G-II)], where G-I had relatively better sperm quality. ESR1 expression was higher in tissues collected from G-I boars than those of collected from G-II boars, and the difference in mRNA expression was significant (P<0.01) in head of epididymis. The ESR1 protein expression results from western blot coincided with the results of qRT-PCR. The ESR1 protein localization observed a strong staining in the cytoplasm of Sertoli cell in the testis, in the epithelial cells in head and tail of epididymis, in smooth muscle in tail of epididymis, and in the post acrosomal region and tail of the spermatozoa. These results will improve the understanding of the functions of the ESR1 in spermatogenesis within the reproductive tract and will shed light on ESR1 as a candidate in the selection of boar with good sperm quality and fertility.

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.

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.

Investigation on association and expression of ESR2 as a candidate gene for boar sperm quality and fertility.

ESR2 is involved in oestrogen-related apoptosis in cell cycle spermatogenesis but their effects have not yet confirmed in pig. Therefore, this study was aimed to investigate the association of ESR2 polymorphism with sperm quality and boar fertility traits and to analyse the ESR2 mRNA and protein expressions in boar reproductive tissues. DNA samples from 203 Pietrain (PI) and 100 Pietrain × Hampshire (PIHA) pigs with records of sperm quality [sperm concentration (SCON), motility (MOT), semen volume (VOL), plasma droplet rate (PDR) and abnormal spermatozoa rate (ASR)] and fertility [non-return rate (NRR) and number of piglet born alive (NBA)] traits were available. A SNP in coding region of ESR2 g.35547A>G in exon 5 was associated with MOT and PDR in the PI and with SCON, VOL, MOT and PDR in PIHA population. For mRNA and protein expression study, a total of six boars were divided into two groups with group I (G-I) and group II (G-II) where G-I characterized for relatively a better sperm quality according to the mean of two groups. mRNA expression was higher in brain and testis than that in all parts of epididymis. Both qRT-PCR and western blot analysis revealed that the ESR2 gene expression and protein expression were significantly higher in testis collected from G-II compared with that of G-I boars. Moreover, ESR2 protein localization in germ cell, Leydig and Sertoli cells, epithelial cells and spermatozoa was remarkable, which indicated the important role of ESR2 in spermatogenesis process. These results might shed new light on the roles of ESR2 in spermatogenesis as candidate for boar fertility, but still the lack of association across populations should be considered.

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.

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

BackgroundBoar 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.ResultsCluster 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.ConclusionsBased 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.

Transcriptome signature of liver tissue with divergent mutton odour and flavour using RNA deep sequencing

Mutton consumption is less popular in many Asian countries including Indonesia, whose consumers often complain about the unpleasant flavour and odour of the meat. The main causes of mutton odour are the two compounds of branched chain fatty acid (BCFA): methylnonanoic (MNA), phenol, 3-methyl (MP), 4-methylnonanoic (MNA) and 4-ethyloctanoic (EOA) present in all the adipose tissue; and the 3-methylindole (MI) or skatole and indole, which are originated from pastoral diets. It is crucial to understand the genetic mechanism of mutton odour and flavour (MOF) to select sheep for lower BCFA and indole thus reduce the unpleasant flavour of meat. The aim of the present study was to investigate transcriptome profiling in liver tissue with divergent MOF using RNA deep sequencing. Liver tissues from higher (n = 3) and lower (n = 3) MOF sheep were analysed using Illumina HiSeq 2500. The total number of reads produced for each liver sample ranged from 21.37 to 25.37 million. Approximately 103 genes were differentially expressed (DEGs) with significance level of p-adjusted value <0.05. Among them, 60 genes were up-regulated, and 43 were down-regulated (p < 0.01, FC > 1.5) in higher MOF group. Differentially regulated genes in high MOF liver samples were enriched in biological processes such as cellular response to chemical stimulus and endogenous stimulus; cellular components such as such as basement membrane and extracellular matrix; and molecular functions such as haeme binding and oxidoreductase activity. Among the DEGs, metabolic phase I related genes belonging to the cytochrome P450 CYP2A6 were dominantly expressed. Additionally, phase II conjugation genes including UDP glucuronosyltransferases UGT2B18, sulfotransferase SULT1C1, and glutathione S-transferase GSTM1 were identified. The dominant candidate genes for SOF could be cytochrome P450, sodium-channel protein, transmembrane protein, glutathione transferase, UDP glucuronosyltransferases and sulfotransferase. Pathway analysis identified steroid hormone biosynthesis and chemical carcinogenesis by cytochrome P450 pathways which may play important roles in MOF-related molecules metabolism. This work highlighted potential genes and gene-networks that may affect meat off flavour and odour in sheep.

Penggunaan Dimethyl Sulfoxide Sebagai Krioprotektan dalam Pembekuan Semen Ayam Kampung (THE USE OF DIMETHYL SULFOXIDE AS CRYOPROTECTIVE AGENT FOR NATIVE CHIKEN FROZEN SEMEN)

This study aim was to determine the best dimethyl sulfoxide (DMSO) concentration for native chicken’s(Gallus gallus domesticus) frozen semen. Semen was collected from three native chickens and then evaluatedmacroscopically and microscopically. Spermatozoa that used in this study had more than 75% spermmotility and more than 25 X 108 mL in sperm concentration. Semen was divided into three equal tubesand each of them diluted with Ringer lactate egg yolk (RLEY) extender consisted of 5, 7 and 9% of DMSO.Diluted semen was added into 0.25 mL of mini straw, then equilibrated at 5 °C for two hours. The semenwas then frozen placed in liquid nitrogen vapors for10 min. Finally, the semen was put into liquid nitrogenwith a temperature of – 196oC. Post thawing quality semen was assessed after 24 hours storage bythawing the straw in warm water (37 °C) for 30 seconds. Data were analyzed using one-way analysis ofvariance. The results demonstrated that the sperm motility and viability diluted with RLEY extenderwith 7% DMSO was better (p<0.05) than those in 5 or 9% DMSO. Semen quality decreased at every stepof freezing, especially between equilibration and thawing. Its concluded that 7% was the best DMSOconcentration for freezing of native chicken semen in Ringer lactate egg yolk extender compared to 5% and9%.

Expression analysis of porcine aromatase (CYP19) as a specific target gene in testis

Purpose: To evaluate storage and temperature effects on extended boar semen characteristics and fertility outcomes. Methods: Extended semen (143 ejaculates/76 boars) were examined within 1 day of collection and after 4 days storage at 20.4 ± 0.1_C (RT) or 14.7 ± 0.2_C (CT). Fertility data (total, live, mummified, stillborn pigs) was derived from 535 inseminations (sow parity range 1–8) using same-batch semen (41/76 boars) stored at 17°C for 0–4 days at stud. Evaluations included CASA (IVOS) sperm motility (total, progressive, rapid, VCL, VAP, VSL, BCF), sperm clump score (0–3), semen pH and temperature. Data were analysed by one-way ANOVA, GLM (stillbirths) and linear regressions (total born and born live) plus chi-square analysis for sow returns and farrowings. Inseminations were performed across four seasons. Results: Initial (receiving) values were: concentration 49.3 ± 1.4 · 106/ml; temperature 18.7 ± 0.2_C; clump score 1.4 ± 0.05; pH 7.45 ± 0.02; motilities: total 57.7 ± 2.4%; progressive 29.3 ± 1.6%; rapid 43.9 ± 2.3%; VCL 154.4 ± 2.8; VAP 75.5 ± 1.6; VSL 43.5 ± 0.9; BCF 33.8 ± 0.3. Both RT and CT storage for 4 days resulted in a rise in pH (7.7 ± 0.02; p < 0.01) and similar (p < 0.01) declines in a number of motility parameters (total, progressive, rapid and BCF) as well as for VSL (p < 0.05). For CT, declines also occurred in VCL and VAP (both p < 0.01). Clump score was not influenced by storage time or temperature. The proportion of sows farrowing vs. returns was similar for semen stored for 1, 2 or 3+ days. However, storage time and season interacted with litter size (p < 0.05), as did storage time, dam line, and parity with % stillborn (both p < 0.01). Total sperm motility positively influenced % live pigs and both total and progressive motility negatively influenced % stillborn (all p < 0.05). Conclusion: Storage of extended semen for 4 days at approximately 15 or 20°C resulted in similar decreases in CASA sperm motility assessments. Overall fertility was not influenced by storage time (0–4 days at 17°C). Despite this, sperm motility measures (total and progressive) influenced % live pigs and % stillborn and storage time interacted with both season and sow factors re. litter size and % stillborn. This research was made possible by financial assistance from Pork CRC Australia.