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 pH45 or pH24 (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.

Identification of Differentially Expressed Genes Associated with Litter Size in Berkshire Pig Placenta.

Improvement in litter size has become of great interest in the pig industry because fecundity is directly related to sow reproductive life. Improved reproduction has thus been achieved by elucidating the molecular functions of genes associated with fecundity. In the present study, we identified differentially expressed genes (DEGs) via transcriptomic analysis using RNA-sequencing (RNA-Seq) in Berkshire pig placentas from larger (LLG, mean litter size >12) and smaller (SLG, mean litter size < 6.5) litter size groups. In total 588 DEGs were identified (p < 0.05, > 1.5-fold change), of which 98 were upregulated, while 490 were downregulated in the LLG compared with the SLG. Gene Ontology (GO) enrichment was also performed. We concluded that 129 of the 588 DEGs were closely related to litter size according to reproduction related genes selected based on previous reports, as 110 genes were downregulated and 19 upregulated in the LLG compared with the SLG. RT-qPCR utilizing specific primers targeting the early growth response 2 (EGR2), pheromaxein c subunit (PHEROC) and endothelial lipase (LIPG) genes showed high accordance with RNA-Seq results. Furthermore, we investigated the upstream regulators of these three genes in the placenta. We found that WNT9B, a Wnt signaling pathway molecule, and IL-6, known inducers of EGR2 and LIPG, respectively, were significantly increased in LLG compared with SLG. We believe that the induction of IL-6 and LIPG may play an important role in increasing nutrition supply through the placenta from the sow to the piglet during gestation. These results provide novel molecular insights into pig reproduction.

Single-Nucleotide Polymorphisms in Pig EPHX1 Gene are Associated with Pork Quality Traits.

Epoxide hydrolase 1 (EPHX1) plays an important role in both the activation and detoxification of exogenous chemicals. Reverse transcription-polymerase chain reaction (RT-PCR) analysis showed that the highest level of EPHX1 expression occurred in Berkshire liver, which is an organ that plays a key role in detoxification. We examined EPHX1 SNPs to analyze effect on increased expression of EPHX1 gene in Berkshire liver by total of 192 pigs of a pure Berkshire line (males = 97; females = 95). As a result, two nonsynonymous SNPs (nsSNPs) of EPHX1 were found from c.685T>G and c.776C > T, and located in 5th and 6th exons, respectively, which constitute the A/b hydrolase 1 domain of epoxide hydrolase. The nsSNP c.685T > G was significant differences in meat color, protein content, collagen content, and pH24 hr. Especially, T and G alleles of the nsSNP c.685T > G were significantly associated with CIE a*/CIE b* and protein content/pH24 hr, respectively. The nsSNP c.776C > T was significant differences in drip loss and protein content. Among meat quality traits to associate with SNPs, the protein content was only significantly associated with sex. Therefore, it is suggested that nsSNP c.685T > G in EPHX1 gene is a potential to apply as appropriate DNA markers for improvement of porcine economic traits.

Selection of appropriate reference genes for RT-qPCR analysis in Berkshire, Duroc, Landrace, and Yorkshire pigs.

Reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) is the most reliable molecular biology technique for assessment of mRNA expression levels. However, to obtain the accurate RT-qPCR results, the expression levels of genes of interest should be normalized with appropriate reference genes and optimal numbers of reference genes. In this study, we assessed the expression stability of 15 well-known candidate reference genes (ACTB, ALDOA, B2M, GAPDH, HPAR1, HSPCB, PGK1, POLR2G, PPIA, RPL4, RPS18, SDHA, TBP, TOP2B, and YWHAZ) in seven body tissues (liver, lung, kidney, spleen, stomach, small intestine, and large intestine) of Berkshire, Landrace, Duroc, and Yorkshire pigs using three excel-based programs, geNorm, NormFinder, and BestKeeper. Combination analysis of these three programs showed that the stable and appropriate reference genes are PPIA, TBP, and HSPCB in Berkshire pigs; PPIA, TBP, RPL4, and RPS18 in Landrace pigs; PPIA and TBP in Duroc pigs; and PPIA, TOP2B, RPL4, and RPS18 in Yorkshire pigs. Because the four pig breeds had different suitable reference genes, the selection of appropriate reference genes is essential in RT-qPCR analyses. Taken together, our data could help to select reliable reference genes for the normalization of expression levels of various target genes in pigs.

Structural importance of the C-terminal region in pig aldo-keto reductase family 1 member C1 and their effects on enzymatic activity.

BACKGROUND: Pig aldo-keto reductase family 1 member C1 (AKR1C1) belongs to AKR superfamily which catalyzes the NAD(P)H-dependent reduction of various substrates including steroid hormones. Previously we have reported two paralogous pig AKR1C1s, wild-type AKR1C1 (C-type) and C-terminal-truncated AKR1C1 (T-type). Also, the C-terminal region significantly contributes to the NADPH-dependent reductase activity for 5α-DHT reduction. Molecular modeling studies combined with kinetic experiments were performed to investigate structural and enzymatic differences between wild-type AKR1C1 C-type and T-type. RESULTS: The results of the enzyme kinetics revealed that Vmax and kcat values of the T-type were 2.9 and 1.6 folds higher than those of the C-type. Moreover, catalytic efficiency was also 1.9 fold higher in T-type compared to C-type. Since x-ray crystal structures of pig AKR1C1 were not available, three dimensional structures of the both types of the protein were predicted using homology modeling methodology and they were used for molecular dynamics simulations. The structural comparisons between C-type and T-type showed that 5α-DHT formed strong hydrogen bonds with catalytic residues such as Tyr55 and His117 in T-type. In particular, C3 ketone group of the substrate was close to Tyr55 and NADPH in T-type. CONCLUSIONS: Our results showed that 5α-DHT binding in T-type was more favorable for catalytic reaction to facilitate hydride transfer from the cofactor, and were consistent with experimental results. We believe that our study provides valuable information to understand important role of C-terminal region that affects enzymatic properties for 5α-DHT, and further molecular mechanism for the enzyme kinetics of AKR1C1 proteins.

Functional mechanism of C-terminal tail in the enzymatic role of porcine testicular carbonyl reductase: a combined experiment and molecular dynamics simulation study of the C-terminal tail in the enzymatic role of PTCR.

Porcine testicular carbonyl reductase, PTCR which is one of the short chain dehydrogenases/reductases (SDR) superfamily catalyzes the NADPH-dependent reduction of carbonyl compounds including steroids and prostaglandins. Previously we reported C-terminal tail of PTCR was deleted due to a nonsynonymous single nucleotide variation (nsSNV). Here we identified from kinetic studies that the enzymatic properties for 5α-dihydrotestosterone (5α-DHT) were different between wild-type and C-terminal-deleted PTCRs. Compared to wild-type PTCR, C-terminal-deleted PTCR has much higher reduction rate. To investigate structural difference between wild-type and C-terminal-deleted PTCRs upon 5α-DHT binding, we performed molecular dynamics simulations for two complexes. Using trajectories, molecular interactions including hydrogen bonding patterns, distance between 5α-DHT and catalytic Tyr193, and interaction energies are analyzed and compared. During the MD simulation time, the dynamic behavior of C-terminal tail in wild-type PTCR is also examined using essential dynamics analysis. The results of our simulations reveal that the binding conformation of 5α-DHT in C-terminal-deleted PTCR is more favorable for reduction reaction in PTCR, which shows strong agreement with kinetic data. These structural findings provide valuable information to understand substrate specificity of PTCR and further kinetic properties of enzymes belonging to the SDR superfamily.

Important role of the C-terminal region of pig aldo-keto reductase family 1 member C1 in the NADPH-dependent reduction of steroid hormones.

The NADPH-dependent reduction activities of two paralogous pig AKR1C1s with and without 19 additional amino acid residues in C-terminus were evaluated against steroid hormones including 5alpha-dihydrotestosterone, testosterone, progesterone, androstenedione and 5alpha-androstane-3.17-dione, which act as substrates of the AKR1C1S. Among the hormones, the AKR1C1s exhibited the highest activity against 5alpha-dihydrotestosterone and the lowest activity against testosterone and progesterone. Furthermore, the AKR1C1s showed the largest differential activities against; 5alpha-dihydrotestosterone, but no such change of activities was found against progestrone and testosterone. These results suggest that the C-terminal region of AKR1C1 plays an important effect in the reduction activities of pig AKR1C1. Thus, the differential activities of two AKR1C1 paralogs observed in the present study provide important insights in understanding the molecular evolution.

Swine liver methylomes of Berkshire, Duroc and Landrace breeds by MeDIPS.

Using a methyl-DNA immunoprecipitation technique in combination with next-generation deep sequencing, we conducted comprehensive DNA methylation profiling of liver genomes from three pig breeds: Berkshire, Duroc and Landrace. The profiles revealed that the distribution patterns of methylation signals along the genome are conserved among the three pig breeds. Specifically, many signals in coding genes were found in introns, and most signals in the repetitive elements were identified in non-long terminal repeat (LTR) retrotransposons such as long and short interspersed repetitive elements, implying a significant association with alternative splicing and expression of retrotransposable elements respectively. Differentially methylated regions among the three pig breeds were identified in the non-LTR retrotransposons, suggesting that they may lead to differential retrotransposable element activity. Altogether, this study provides advanced swine methylome data and valuable resources for understanding the function of DNA methylation in the evolutionary divergence of different pig breeds.

RNA-Seq approach for genetic improvement of meat quality in pig and evolutionary insight into the substrate specificity of animal carbonyl reductases.

Changes in meat quality traits are strongly associated with alterations in postmortem metabolism which depend on genetic variations, especially nonsynonymous single nucleotide variations (nsSNVs) having critical effects on protein structure and function. To selectively identify metabolism-related nsSNVs, next-generation transcriptome sequencing (RNA-Seq) was carried out using RNAs from porcine liver, which contains a diverse range of metabolic enzymes. The multiplex SNV genotyping analysis showed that various metabolism-related genes had different nsSNV alleles. Moreover, many nsSNVs were significantly associated with multiple meat quality traits. Particularly, ch7:g.22112616A>G SNV was identified to create a single amino acid change (Thr/Ala) at the 145th residue of H1.3-like protein, very close to the putative 147th threonine phosphorylation site, suggesting that the nsSNV may affect multiple meat quality traits by affecting the epigenetic regulation of postmortem metabolism-related gene expression. Besides, one nonsynonymous variation, probably generated by gene duplication, led to a stop signal in porcine testicular carbonyl reductase (PTCR), resulting in a C-terminal (E281-A288) deletion. Molecular docking and energy minimization calculations indicated that the binding affinity of wild-type PTCR to 5α-DHT, a C(21)-steroid, was superior to that of C-terminal-deleted PTCR or human carbonyl reductase, which was very consistent with experimental data, reported previously. Furthermore, P284 was identified as an important residue mediating the specific interaction between PTCR and 5α-DHT, and phylogenetic analysis showed that P284 is an evolutionarily conserved residue among animal carbonyl reductases, which suggests that the C-terminal tails of these reductases may have evolved under evolutionary pressure to increase the substrate specificity for C(21)-steroids and facilitate metabolic adaptation. Altogether, our RNA-Seq revealed that selective nsSNVs were associated with meat quality traits that could be useful for successful marker-assisted selection in pigs and also represents a useful resource to enhance understanding of protein folding, substrate specificity, and the evolution of enzymes such as carbonyl reductase.