Hypomethylation in the promoter region of ZPBP as a potential litter size indicator in Berkshire pigs.

In pigs, litter size is typically defined as the total number of piglets born (TNB) or the number of piglets born alive (NBA). Increasing pig litter size is of great economic interest as a means to increase productivity. The capacity of the uterus is a critical component of litter size and may play a central role in prolificacy. In this study, we investigated litter-size-related epigenetic markers in uterine tissue from Berkshire pigs with smaller litter size groups (SLGs) and larger litter size groups (LLGs) using genome-wide bisulfite sequencing (GWBS). A total of 3269 differentially methylated regions (DMRs) were identified: 1566 were hypermethylated and 1703 hypomethylated in LLG compared to SLG. The zona pellucida binding protein (ZPBP) gene was significantly hypomethylated in the LLG promoter region, and its expression was significantly upregulated in uterine tissue. Thus, the methylation status of ZPBP gene was identified as a potential indicator of litter size. Furthermore, we verified its negative correlation with litter size traits (TNB and NBA) in whole blood samples from 172 Berkshire sows as a blood-based biomarker by a porcine methylation-specific restriction enzyme polymerase chain reaction (PMP) assay. The results suggest that the methylation status of the ZPBP gene can serve as a valuable epigenetic biomarker for hyperprolific sows.

Identification of a Bromodomain-containing Protein 2 (BRD2) Gene Polymorphic Variant and Its Effects on Pork Quality Traits in Berkshire Pigs.

Bromodomain-containing protein 2 (BRD2) is a nuclear serine/threonine kinase involved in transcriptional regulation. We investigated the expression and association of the BRD2 gene as a candidate gene for meat quality traits in Berkshire pigs. BRD2 mRNA was expressed at relatively high levels in muscle tissue. Statistical analysis revealed that the c.1709G>C polymorphism of the BRD2 gene was significantly associated with carcass weight, meat color (a*, redness), protein content, cooking loss, water-holding capacity, carcass temperatures 4, 12 and 24 h postmortem, and the 24 h postmortem pH in 384 Berkshire pigs. Therefore, this polymorphism in the porcine BRD2 gene may be used as a candidate genetic marker to improve meat quality traits in pigs.

Identification of Association between Supply of Pork and Production of Meat Products in Korea by Canonical Correlation Analysis.

To identify correlation between fresh meat and processed meat products, we performed canonical correlation analysis (CCA) to predict the relationship between pork supply and meat product production in Korea. Results of CCA showed a canonical correlation of 0.8576 in the first canonical pair (p<0.01). The production of meat products showed the highest correlation with pork import but the lowest correlation with the production of domestic pork. Although Korean consumer preferred meat products produced by fresh domestic pork, inexpensively imported pork with high share in meat products was supplied in the market. Therefore, securing domestically produced raw meat is important for expanding consumption of domestic meat products. Results of this study suggest that meat processor and pig producer can achieve the 6th industrialization by combining the production of raw pork materials, meat processing, and sales service.

Association of single-nucleotide polymorphisms in NAT9 and MAP3K3 genes with litter size traits in Berkshire pigs.

Litter size is an economically important trait in the pig industry. We aimed to identify genetic markers associated with litter size, which can be used in breeding programs for improving reproductive traits. Single-nucleotide polymorphisms (SNPs) of Berkshire pigs in the N-acetyltransferase 9 (NAT9) and Mitogen-activated protein kinase kinase kinase 3 (MAP3K3) genes were from RNA sequencing results, and already exist in the databank (NCBI), and were confirmed by polymerase chain reaction and restriction fragment length polymorphism (PCR-RFLP). A total of 272 Berkshire sows were used to examine the genotype, and their association with litter size traits was analyzed. The NAT9 SNP was located in chromosome 12 exon 640 mRNA (A  >  G) and the MAP3K3 SNP was located in chromosome 12 intron 11 (80, C  >  T). Association analysis indicated that the GG genotype of NAT9 and the CT genotype of MAP3K3 had the highest values for litter size traits. The GG genotype expressed higher levels of NAT9 mRNA in the endometrium than the other genotypes did, and a positive correlation was found between litter size traits and NAT9, but not MAP3K3 expression level. These results indicate that the NAT9 and MAP3K3 can be used as candidate genes applicable in breeding program for the improvement of litter size traits in Berkshire pigs.

Effect of Single Nucleotide Polymorphisms in IGFBP2 and IGFBP3 Genes on Litter Size Traits in Berkshire Pigs.

Litter size is among the most important traits in swine breeding. However, information on the genetics of litter size in pigs is lacking. In this study, we identified single nucleotide polymorphisms (SNPs) in the insulin-like growth factor binding protein 2 and 3 (IGFBP2 and IGFBP3) genes in Berkshire pigs and analyzed their association with litter size traits. The IGFBP2 SNP was located on chromosome 15 intron 2 (455, A > T) and the IGFBP3 SNP was on chromosome 18 intron 2 (53, A > G). The AT type of IGFBP2 and the GG type of IGFBP3 had the highest values for all litter size traits including total number born (TNB), number of pigs born alive, and breeding value according to TNB. Homozygous GG pigs expressed higher levels of IGFBP3 mRNA in the endometrium than pigs of other genotypes, and a positive correlation was observed between litter size traits and IGFBP3 but not IGFBP2 expression level. These results suggest that SNPs in the IGFBP2 and the IGFBP3 gene are useful biomarkers for increasing the reproductive productivity of Berkshire pigs.

Squalene epoxidase plays a critical role in determining pig meat quality by regulating adipogenesis, myogenesis, and ROS scavengers.

In mammals, Squalene epoxidase (SQLE) is an enzyme that converts squalene to 2,3-oxidosqualene, in the early stage of cholesterol generation. Here, we identified single nucleotide polymorphisms (SNPs) in the SQLE gene (c.2565 G > T) by RNA Sequencing from the liver tissue of Berkshire pigs. Furthermore, we found that homozygous GG pigs expressed more SQLE mRNA than GT heterozygous and TT homozygous pigs in longissimus dorsi tissue. Next, we showed that the SNP in the SQLE gene was associated with several meat quality traits including backfat thickness, carcass weight, meat colour (yellowness), fat composition, and water-holding capacity. Rates of myogenesis and adipogenesis induced in C2C12 cells and 3T3-L1 cells, respectively, were decreased by Sqle knockdown. Additionally, the expression of myogenic marker genes (Myog, Myod, and Myh4) and adipogenic marker genes (Pparg, Cebpa, and Adipoq) was substantially downregulated in cells transfected with Sqle siRNA. Moreover, mRNA expression levels of ROS scavengers, which affect meat quality by altering protein oxidation processes, were significantly downregulated by Sqle knockdown. Taken together, our results suggest the molecular mechanism by which SNPs in the SQLE gene can affect meat quality.

DNA methylation patterns and gene expression associated with litter size in Berkshire pig placenta.

Increasing litter size is of great interest to the pig industry. DNA methylation is an important epigenetic modification that regulates gene expression, resulting in livestock phenotypes such as disease resistance, milk production, and reproduction. We classified Berkshire pigs into two groups according to litter size and estimated breeding value: smaller (SLG) and larger (LLG) litter size groups. Genome-wide DNA methylation and gene expression were analyzed using placenta genomic DNA and RNA to identify differentially methylated regions (DMRs) and differentially expressed genes (DEGs) associated with litter size. The methylation levels of CpG dinucleotides in different genomic regions were noticeably different between the groups, while global methylation pattern was similar, and excluding intergenic regions they were found the most frequently in gene body regions. Next, we analyzed RNA-Seq data to identify DEGs between the SLG and LLG groups. A total of 1591 DEGs were identified: 567 were downregulated and 1024 were upregulated in LLG compared to SLG. To identify genes that simultaneously exhibited changes in DNA methylation and mRNA expression, we integrated and analyzed the data from bisulfite-Seq and RNA-Seq. Nine DEGs positioned in DMRs were found. The expression of only three of these genes (PRKG2, CLCA4, and PCK1) was verified by RT-qPCR. Furthermore, we observed the same methylation patterns in blood samples as in the placental tissues by PCR-based methylation analysis. Together, these results provide useful data regarding potential epigenetic markers for selecting hyperprolific sows.

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.

The rs196952262 Polymorphism of the AGPAT5 Gene is Associated with Meat Quality in Berkshire Pigs.

High-quality meat is of great economic importance to the pig industry. The 1-acylglycerol-3-phosphate-O-acyltransferase 5 (AGPAT5) enzyme converts lysophosphatidic acid to phosphatidic acid in the mitochondrial membrane. In this study, we found that the porcine AGPAT5 gene was highly expressed in muscle tissue, influencing meat characteristics, and we also identified a non-synonymous single-nucleotide polymorphism (nsSNP) (rs196952262, c.673 A>G) in the gene, associated with a change of isoleucine 225 to valine. The presence of this nsSNP was significantly associated with meat color (lightness), lower cooking loss, and lower carcass temperatures 1, 4, and 12 h after slaughter (items T1, T4, and T12 on the recognized quality scale, respectively), and tended to increase backfat thickness and the water-holding capacity. These results suggest that nsSNP (c.673A>G) of the AGPAT5 gene is a potential genetic marker of high meat quality in pigs.

Pork Quality Traits According to Postmortem pH and Temperature in Berkshire.

This study was performed to investigate the role of pH and temperature postmortem, and to demonstrate the importance of these factors in determining meat quality. Postmortem pH 45min (pH at 45 min postmortem or initial pH) via analysis of Pearson's correlation showed high positive correlation with pH change pH c24 (pH change from pH 45min to pH 24h postmortem). However, postmortem pH after 24 h (pH 24h or ultimate pH) had a high negative correlation with pH change, pH c24 , CIE L*, and protein content. Initial temperature postmortem (T 1h ) was positively associated with a change in temperature from 45 min to 24 h postmortem (T c24 ) and cooking loss, but negatively correlated with water holding capacity. Temperature at 24 h postmortem (T 24h ) was negatively associated with T c24 . Collectively, these results indicate that higher initial pH was associated with higher pH c24 , T 1h , and T c24 . However, higher initial pH was associated with a reduction in carcass weight, backfat thickness, CIE a* and b*, water holding capacity, collagen and fat content, drip loss, and cooking loss as well as decreased shear force. In contrast, CIE a* and b*, drip loss, cooking loss, and shear force in higher ultimate pH was showed by a similar pattern to higher initial pH, whereas pH c24 , carcass weight, backfat thickness, water holding capacity, fat content, moisture content, protein content, T 1h , T 24h , and T c24 were exhibited by completely differential patterns (p<0.05). Therefore, we suggest that initial pH, ultimate pH, and temperatures postmortem are important factors in determining the meat quality of pork.

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.

Proteomic analysis of reproduction proteins involved in litter size from porcine placenta.

A gel-free and label-free quantitative proteomic approach based on a spectral counting strategy was performed to discover prolificacy-related proteins. Soluble proteins of porcine placenta from small litter size group (SLSG) and large litter size group (LLSG) were extracted and subsequently applied to in-solution tryptic digestion followed by liquid chromatography-tandem mass spectrometry analysis. Six and thirteen proteins were highly expressed in SLSG and LLSG, respectively. Of the dominantly expressed proteins, we chose prolificacy-related proteins such as puromycin-sensitive aminopeptidase (PSA) and retinol-binding protein 4 (RBP4). Western blot analysis confirmed that the processed form (70 kDa) of PSA was more expressed and RBP4 (23 kDa) was dominantly expressed in LLSG. These results indicate that PSA and RBP4 are representative proteins involved in porcine fertility traits, and this finding may help to increase litter size of pigs.

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.

RNA-seq analysis and de novo transcriptome assembly of Jerusalem artichoke (Helianthus tuberosus Linne).

Jerusalem artichoke (Helianthus tuberosus L.) has long been cultivated as a vegetable and as a source of fructans (inulin) for pharmaceutical applications in diabetes and obesity prevention. However, transcriptomic and genomic data for Jerusalem artichoke remain scarce. In this study, Illumina RNA sequencing (RNA-Seq) was performed on samples from Jerusalem artichoke leaves, roots, stems and two different tuber tissues (early and late tuber development). Data were used for de novo assembly and characterization of the transcriptome. In total 206,215,632 paired-end reads were generated. These were assembled into 66,322 loci with 272,548 transcripts. Loci were annotated by querying against the NCBI non-redundant, Phytozome and UniProt databases, and 40,215 loci were homologous to existing database sequences. Gene Ontology terms were assigned to 19,848 loci, 15,434 loci were matched to 25 Clusters of Eukaryotic Orthologous Groups classifications, and 11,844 loci were classified into 142 Kyoto Encyclopedia of Genes and Genomes pathways. The assembled loci also contained 10,778 potential simple sequence repeats. The newly assembled transcriptome was used to identify loci with tissue-specific differential expression patterns. In total, 670 loci exhibited tissue-specific expression, and a subset of these were confirmed using RT-PCR and qRT-PCR. Gene expression related to inulin biosynthesis in tuber tissue was also investigated. Exsiting genetic and genomic data for H. tuberosus are scarce. The sequence resources developed in this study will enable the analysis of thousands of transcripts and will thus accelerate marker-assisted breeding studies and studies of inulin biosynthesis in Jerusalem artichoke.

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.