Next Article in Journal
Volatile Fingerprinting and Sensory Profiles of Coffee Cascara Teas Produced in Latin American Countries
Next Article in Special Issue
1H NMR Metabolic Profiling and Meat Quality in Three Beef Cattle Breeds from Northeastern Thailand
Previous Article in Journal
Investigation of Chemical Compounds and Evaluation of Toxicity, Antibacterial, and Anti-Inflammatory Activities of Three Selected Essential Oils and Their Mixtures with Moroccan Thyme Honey
Previous Article in Special Issue
Metabolomic Analysis on the Mechanism of Nanoselenium Biofortification Improving the Siraitia grosvenorii Nutritional and Health Value
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Genetic Architecture of Meat Quality Traits in a Crossbred Commercial Pig Population

1
National Engineering Research Center for Breeding Swine Industry, College of Animal Science, South China Agricultural University, Guangzhou 510642, China
2
Guangdong Zhongxin Breeding Technology Co., Ltd., Guangzhou 511466, China
3
College of Animal Science and Technology, Zhongkai University of Agriculture and Engineering, Guangzhou 510225, China
4
Guangdong Provincial Key Laboratory of Agro-Animal Genomics and Molecular Breeding, South China Agricultural University, Guangzhou 510642, China
5
Yunfu Subcenter of Guangdong Laboratory for Lingnan Modern Agriculture, Yunfu 527400, China
*
Author to whom correspondence should be addressed.
Foods 2022, 11(19), 3143; https://doi.org/10.3390/foods11193143
Submission received: 21 September 2022 / Revised: 7 October 2022 / Accepted: 8 October 2022 / Published: 9 October 2022
(This article belongs to the Special Issue Omics Technologies in Food Science 2.0)

Abstract

:
Meat quality is of importance in consumer acceptance and purchasing tendency of pork. However, the genetic architecture of pork meat quality traits remains elusive. Herein, we conducted genome-wide association studies to detect single nucleotide polymorphisms (SNPs) and genes affecting meat pH and meat color (L*, lightness; a*, redness; b*, yellowness) in 1518 three-way crossbred pigs. All individuals were genotyped using the GeneSeek Porcine 50K BeadChip. In sum, 30 SNPs and 20 genes are found to be associated with eight meat quality traits. Notably, we detect one significant quantitative trait locus (QTL) on SSC15 with a 143 kb interval for meat pH (pH_12h), together with the most promising candidate TNS1. Interestingly, two newly identified SNPs located in the TTLL4 gene demonstrate the highest phenotypic variance of pH_12h in this QTL, at 2.67%. The identified SNPs are useful for the genetic improvement of meat quality traits in pigs by assigning higher weights to associated SNPs in genomic selection.

1. Introduction

Pork contributes to a large share of the world’s meat provision, especially in the Chinese market. People pay more attention to meat quality with the improvement in living standards. Pork quality can increase both customer purchasing tendency and pork consumption, and then enable the pig industry to breed pigs with better meat quality. Studies demonstrate that meat pH and meat color are closely connected with the quality of meat, affecting freshness, shearing force, and drip loss of pork [1,2,3]. Therefore, applying molecular genetic approaches to meat pH and color traits are useful strategies to improve the meat quality [4,5,6], and this can accelerate the breeding cycle in pigs. Meat pH and meat color (L*, lightness; a*, redness; b*, yellowness) of pork are low-heritability traits and varies among different populations [7,8], implying that meat quality is difficult to improve by traditional selection [9]. Due to high phenotyping costs of meat quality traits (because these phenotypes usually are measured post-mortem) and genotyping costs of large-scale pigs, genetic processes in improving pork meat quality were slow. The effort to identify quantitative trait loci (QTLs) and genetic markers has been ongoing for many years using genome-wide association study (GWAS), which made it possible to discern the genetic basis of meat pH and meat color traits. Genome-wide association studies are widely used to detect associations of genotypes with phenotypes by testing for differences in the allele frequency of single nucleotide polymorphisms (SNPs) between individuals [10]. To date, a large number of QTLs (>900) have been detected as being associated with pork meat pH and meat color, as reported in the pig QTL database [11]. Previous studies also report that the PRKAG3 gene plays an important role in affecting pork meat pH and color [12,13,14]. Ma et al. [15] demonstrated that a causal loss-of-function mutation in the PHKG1 gene causes high glycogen content and low meat quality in pigs. Although numerous QTLs and genes affecting meat pH and meat color traits have been identified [16,17,18], the genetic architecture of these polygenic quantitative traits remains elusive.
The aim of this study was to detect QTLs and genes associated with meat quality traits in pigs. To this end, we conducted GWAS for eight meat pH and meat color traits in a three-way crossbred commercial Duroc × (Landrace × Yorkshire) (DLY) pig population, numbering 1518 individuals. We detect numerous significant SNPs and promising candidate genes affecting meat quality and further reveal the complexity of the genetic architecture of meat quality. In terms of improving the genetic process of meat quality, the identified significant SNPs can be used for the genetic improvement of meat pH and color in pigs by assigning higher weights to associated SNPs in genomic selection. The expansion of the basis of meat quality can potentially provide novel insights into future molecular breeding of pork and achieve better meat quality.

2. Materials and Methods

2.1. Ethics Statement

All animals used in this study met the guidelines for the care and use of experimental animals established by the Ministry of Agriculture of China. Tissue samples from pigs were collected with the approval of the ethics committee of South China Agricultural University (Guangzhou, China) under 2018F098.

2.2. Animal Samples and Meat Quality Traits Phenotyping

The experiment animals were from a three-way crossbred commercial pig population. Briefly, 84 Duroc boars were mated to 397 Landrace × Yorkshire sows to produce a large-scale offspring. All pigs sustained uniform feeding conditions and were raised on four farms of Wen’s Foodstuffs Group Co., Ltd. (Guangdong, China). After fattening, 1518 individuals born from 2018 to 2019 were slaughtered for phenotype recording at an average body weight of 115 kg in 13 batches. For each animal, longissimus thoracis (LT) muscle was removed from the left side of each carcass and was used to determine meat quality traits. Meat pH and meat color traits measurements were performed on LT muscle as described previously [19]. Briefly, a portable pH meter (equipped with an insertion glass electrode) was used to determine the pH of LT muscle samples between the 11th and 12th ribs at 45 min (pH_45min) and 12h (pH_12h) after slaughter, respectively. The pH meter was calibrated before measurements using standard phosphate buffers (pH = 4.01 and 7.00) and adjusted to the actual temperature of sample measurement following the instrumental user’s manual. Three meat color parameters (L*, a*, and b*) were measured on the exposed cut surface after blooming for 30 min of the LT at 45 min and 12 h post-mortem, respectively, using a CM-2600d/2500d Minolta Chromameter (Tokyo, Japan) with an 8 mm measuring port, D65 illuminant, and one permanent observer to avoid human factor. Meat samples of LT were kept at 4 °C inside the refrigerator until the pH, L*, a*, and b* were measured at 12 h post-mortem.

2.3. Genotyping and Quality Control

Genomic DNA of each pig was extracted from ear tissue via a standard phenol/chloroform method and was diluted to 50 ng/μL [20] for genotyping procedure. The 1518 DLY pigs were genotyped using a GeneSeek Porcine 50K BeadChip (Neogen, Lincoln, NE, USA), which contained 50,703 SNPs. After genotyping, the genotype dataset was converted from Sus scrofa genome (Sscrofa) 10.2 to build Sscrofa 11.1. Quality control (QC) procedures were conducted using PLINK v1.07 software [21] with the following criteria: individual call rate > 95%, SNP call rate > 99%, minor allele frequency > 1%, and p > 10−6 for the Hardy–Weinberg equilibrium test. Only SNPs located on the autosome chromosomes were retained in this study. After QC, we removed 203 SNPs because of missing genotype data, 18,004 SNPs due to Hardy–Weinberg exact test, and 45 SNPs due to minor allele threshold. Moreover, 4245 SNPs not located on autosome chromosomes were discarded. Finally, 1518 pigs and 28,206 SNPs remained for subsequent analyses.

2.4. Population Structure and Linkage Disequilibrium Estimation

The qualified SNPs were used to conduct principal component analysis (PCA) and linkage disequilibrium (LD) analysis to investigate the population structure of the DLY pigs. PCA was conducted with GCTA software [22]. The SNPs distributed in the pig genomes were used to calculate the average LD decay distance, which was estimated as squared correlation of allele frequencies (r2). The window size of LD block was set as 1000 using PLINK v1.07 software.

2.5. Genome-Wide Association Study

We conducted GWAS with GEMMA software [23] for eight meat pH and meat color traits including pH_45min, L*_45min, a*_45min, b*_45min, pH_12h, L*_12h, a*_12h, and b*_12h. The mixed linear model used in this study was as follows:
y = W α + x β + g + ε
where y is the vector of each meat quality trait in DLY pig population; W is a incidence matrix of covariates (fixed effects) including farms (four levels), sex (two levels), slaughter batches (13 levels), and the top five eigenvectors of PCs and a column of 1s; α is a vector of corresponding coefficients that includes the intercept; x is the vector of all SNP marker genotypes; β refers to the corresponding effect of the SNP; g refers to an n × 1 vector of random effects, with g ~ MVN (0, K σ g 2 ), and ε is the vector of random residuals, with ε ~ MVN (0, I σ e 2 ); K is genomic relatedness matrix and σ g 2 is the additive genetic variance; I is the identity matrix and σ e 2 is the residual variance; n refers to the number of analyzed DLY pigs; and MVN denotes multivariate normal distribution.
In the current study, we used Bonferroni correction to acquire the genome-wide significant and chromosome-wide significant thresholds. A SNP was considered to have genome-wide significance at p < 0.05/N and chromosome-wide significance at p < 1/N, where N refers to the total number of qualified SNPs. The quantile–quantile (Q–Q) plots of the eight meat quality traits were constructed for the assessment of population stratification effects on GWAS results.
We used Haploview v4.2 software [24] to perform haplotype block analysis within the identified significant QTLs.

2.6. Estimation of Heritability and Phenotypic Variance

In this study, we conducted SNP-based heritability estimation of each meat quality trait and calculated trait-associated SNP’s contribution to phenotypic variance via GCTA software [22]. The SNP-based heritability and the proportion of phenotypic variance explained by significant SNPs were calculated using:
y = X β + g + ε   with   var ( y ) = A g σ g 2 + Ι σ e 2
where y is the phenotypic value of each meat quality trait; β is the vector including covariates as described above; X is an incidence matrix; g is an n × 1 vector of the total genetic effects of the analyzed pigs; Ι is the identity matrix; A g refers to the genomic relatedness matrix between animals; σ g 2 is the additive genetic variance estimated by the restricted maximum likelihood approach; and σ e 2 is the residual variance.

2.7. Functional Candidate Genes Search

The functional genes were searched for based on Sscrofa 11.1 genome version (http://asia.ensembl.org/Sus_scrofa/Info/Index, accessed on 9 July 2022). Genes nearest the significant SNPs are list in the Tables. Furthermore, we manually queried PubMed and the literature for information about the association between all candidate genes nearest peak SNPs and the analyzed meat quality traits.

3. Results and Discussion

3.1. Phenotype Statistic and Heritability Estimation

The descriptive statistics of eight meat quality traits in DLY pigs are listed in Table 1. For meat pH and meat color values, they are important parameters affecting sensory quality of pork. In this study, differences in phenotypic values are observed between the meat quality traits measured at 45 min and 12 h post-mortem. A similar pattern is also observed in three Chinese indigenous breeds [19]. Genetic factors that contribute to the phenotype variations in pH, L*, a*, and b* are explored (see below). The heritability of meat pH ranges from 0.12 to 0.17, while the heritability of meat color traits ranges from 0.03 to 0.12. As previous studies report [7,8], meat pH and meat color (L*, a*, b*) of pork are low-heritability traits. Miar et al. [8] report the heritability of meat color varies from 0.1 to 0.4, and pH is 0.15 in commercial crossbred pigs. Analyses performed in the present study indicate that the heritability of meat color traits in DLY pigs is lower than that reported in a previous study. For low-heritability traits in pigs, traditional breeding strategies have weak power in terms of improving the genetic process. Therefore, meat quality may be improved using genetic methods, such as genomic selection, to accelerate the breeding cycle in pigs [8].

3.2. Assessment of Population Structure and Linkage Disequilibrium Decay

The average r2 of 0.2 is about 200 kb apart (Figure 1). Therefore, our three-way crossbreed DLY pig population has a low LD pattern [25,26], which is useful for mining QTLs affecting meat quality. In this study, we added the first five principal components into the association model as covariates to correct the potential population structure. As shown in the Q–Q plots (together with the Manhattan plots: Figure 2, Figure 3 and Figure 4), the genomic inflation factors (lambda) of the GWAS results for all traits are close to 1, which indicates no evidence of population stratification.

3.3. Genome-Wide Association Studies for Meat pH

The GWAS for pH_45min identifies two SNPs that are located on SSC6 and 16 (Figure 2A; Table 2). Figure 2B shows that the lambda is 1.01. These two SNPs surpass the chromosome-wide threshold (p < 3.56 × 10−5) and contribute to 2.36% and 1.27% of the phenotypic variance, respectively. The SNP (rs81274518) on SSC6 is located within the GRIK5 gene and the SNP (rs81324442) on SSC16 is located within the MTMR12 gene. GRIK5 is a protein-coding gene that belongs to the glutamate-gated ionic channel family. This gene encodes glutamate ionotropic receptor kainite type subunit 5 in diverse ophthalmologic and vascular disorders when its expression levels are reduced [27]. MTMR12 is a protein-coding gene and the protein encoded by this gene functions as an adaptor subunit in a complex with an active PtdIns(3)P 3-phosphatase. There is no study that highlights the MTMR12 gene as a candidate gene related to meat quality traits in pigs. Significant SNPs detected by GWAS for meat pH_12h of the DLY pigs are shown in Table 3. In sum, 16 SNPs are found to be associated with pH_12h (Figure 2C). Figure 2D shows that the lambda is 0.993. Of these, 6 surpass the genome-wide threshold (p < 1.78 × 10−6) and 10 surpass the chromosome-wide threshold (p < 3.56 × 10−5). Notably, we detect one consistent QTL on SSC15 for meat quality, and the top SNP demonstrate the highest phenotypic variance of pH_12h in this QTL, at 2.67% (see below). We also identify three new SNPs associated with pH_12h, which are located on SSC2 and SSC9, respectively. The rs81303631 SNP and rs81295472 SNP on SSC2, which explain 0.72% and 0.31% of the phenotypic variance, respectively, are close to the protein-coding genes FAM170A and PRR16. The FAM170A gene acts as a nuclear transcription factor that positively regulates the expression of heat shock genes. The PRR16 is a protein-coding gene involving in the positive regulation of cell size and positive regulation of translation [28]. PRR16 is highly expressed in human placental endothelial cells. Although it is shown to regulate cell size, PRR16’s role in meat pH is unclear. The rs81316230 SNP on SSC9, which explains 0.31% of the phenotypic variance, is located inside the PHOX2A gene, which encodes a protein that contains a paired-like homeodomain most similar to that of the Drosophila aristaless gene product. It is shown that the PHOX2A gene is associated with respiratory rhythm and autonomic nervous system development [29].

3.4. Genome-Wide Association Studies for Meat Color Traits

In this study, we measured meat color parameters at two different time points (45 min and 12 h) after slaughter, which represent the meat color of fresh meat and chilled meat with different economic values.
In total, six SNPs (Table 2) in five chromosomes (SSC4, 8, 11, 13, and 14) surpass the chromosome-wide threshold (p < 3.56 × 10−5) and are found to be associated with three meat color traits (L*_45min, a*_45min, and b*_45min). For L*_45min, the SNP rs80971313 on SSC4 (Figure 3A), which explains 0.90% of the phenotypic variance, is located within the KCNB2 gene. Figure 3B shows that the lambda is 1.04. The KCNB2 gene is an important regulator of neuron excitability in mammal brains. It is reported that the KCNB2 gene is involved in the muscular growth pathway that plays an important role in calcium and potassium transport, and in meat tenderization through its involvement in the proteolytic system responsible for post-mortem tenderization and muscle contraction [30,31]. However, post-mortem tenderization has an effect on pH decline during early post-mortem, which can affect other proteins and enzymatic systems that control meat quality characteristics such as meat color [32]. Therefore, we proposed the KCNB2 gene as a promising candidate affecting meat color such as L*_45min. For a*_45min, three SNPs are detected (Figure 3C) and the lambda is 1.03 (Figure 3D). The strongest GWAS signal occurs on SSC14, and the peak SNP (rs80944912) is located within the CNNM2 gene, but only explains 0.01% of the phenotypic variance, which implies a minor effect. Two SNPs are found to be significantly associated with b*_45min (Figure 3E) and the lambda is 1.04 (Figure 3F). The rs342146775 SNP on SSC11 is located within the NALCN gene and rs343103580 SNP on SSC13 is located within the RFTN1 gene. NALCN is a protein-coding gene that plays an important role in insulin release [33]. Its related pathways are transport of glucose and other sugars, metal ions, and amine compounds. The RFTN1 gene is involved in the B-cell receptor signaling pathway, membrane raft assembly, and positive regulation of growth rate. This gene is important in the formation or maintenance of membrane lipid rafts [34].
We performed GWAS in DLY pigs to detect SNPs and genes affecting meat color traits (12 h), including L*_12h, a*_12h, and b*_12h. In sum, five SNPs (Table 3) surpass the chromosome-wide threshold (p < 3.56 × 10−5) and are found to be associated with L*_12h and b*_12h. Three SNPs are associated with L*_12h (Figure 4A) and the lambda is 1.01 (Figure 4B). The rs81316230 SNP on SSC9, which explains 0.40% of the phenotypic variance, is close to the PHOX2A gene. This common result implies that rs81316230 SNP, together with the nearest gene PHOX2A, may have an effect on both meat color and pH. For a*_12h, no GWAS signal occurs on the autosomes (Figure 4C) and the lambda is 1.02 (Figure 4D). However, the rs81303631 SNP does not surpassed the chromosome-wide threshold, but is also highlighted as a candidate genetic marker for a*_12h due to its effect for pH_12h, implying that this SNP may have pleiotropic effects on meat quality in pigs. As for b*_12h, we identify two significant SNPs (Figure 4E) that exceed the genome-wide threshold (p < 1.78 × 10−6) and the lambda is 1.02 (Figure 4F). The rs343786555 SNP is located at 17.67 Mb near the SHAS2 gene on SSC4. The rs346116771 SNP is close to the protein-encoding gene EPYC. However, these two SNPs have minor effect for b*_12h, which only contributes less than 0.09% of the phenotypic variance. In addition, it is reported that the expression of the SHAS2 gene is significantly higher in baicalin-treated than control oocytes in the study related to pig oocytes and subsequent embryo development following parthenogenetic activation [35]. The EPYC gene is a member of the small leucine-rich repeat proteoglycan family. It regulates fibrillogenesis by interacting with collagen fibrils and other extracellular matrix proteins [36]. The potential role of the two genes in meat quality needs further investigation.
In this study, there are no overlapped GWAS signals among meat color traits measured in 45 min and 12 h. Our results demonstrate that the genetic factors affecting meat color in the two stages may be different. Meat color is an important economic trait and the genes causing the different meat color are shown in different pig breeds [37], implying that meat color is regulated by complex genetic networks.

3.5. Constant QTL for Meat pH on SSC15 Detected by GWAS

In this study, 13 SNPs associated with pH_12h are located in a QTL region on SSC15 between 117.25–122.99 Mb (Sscofa 11.1). Figure 5A is a region plot of this QTL and shows the LD pattern between the GWAS peak (rs335443100) and other significant SNPs. Notably, we detect one significant QTL on SSC15 with a 143 kb interval (Figure 5B) for pH_12h trait, together with the most promising candidate gene TNS1. The GWAS peak is approximately 733 kb away from the PRKAG3 gene (also located on SSC15), a gene affecting glycogen in Hampshire and meat quality traits in pigs [12,13,14]. Based on the GWAS results in this study, the region around the PRKAG3 gene gives a small p-value. Moreover, low LD (r2) is observed between the SNPs in the TNS1 gene and PRKAG3 region (Figure 5A). Beyond the PRKAG3 gene, functional variants may exist in the identified QTL regions in the current study [38]. The top SNP (rs335443100) explains 1.94% phenotypic variance of the pH_12h trait. Previous studies report the potential influence of the TNS1 gene on meat quality [39]. For instance, the TNS1 gene is identified as a candidate in GWAS for pork meat pH measured 24 h after slaughter in Finnish Yorkshire pigs [38] and in Chinese Laiwu pigs [17]. The TNS1 gene plays an important role in TGF-beta-induced myofibroblast differentiation [40], which indicates the involvement of this gene in the muscle development of pigs in the early period. These results imply that the TNS1 gene on SSC15 of this QTL is likely to be responsible for meat pH. However, two SNPs located in the TTLL4 gene explain the highest level of phenotypic variance of pH_12h in this QTL region, at 2.67%. We then evaluated the phenotype distribution pattern of the peak GWAS signal (Figure 5C) and the SNPs within the TTLL4 gene (Figure 5D) in the DLY pigs. For the two SNPs, they have similar a decrease effect between GG vs. AA (GWAS peak: 5.79 vs. 5.65, respectively) and AA vs. CC (TTLL4: 5.79 vs. 5.67, respectively) genotype. However, the LD between the two SNPs is 0.48, implying an independent effect on meat pH by the TTLL4 genes. Therefore, the TTLL4 gene is proposed as another candidate gene responsible for meat pH in this QTL. TTLL4 is a protein-coding gene and Gene Ontology annotations related to this gene include ligase activity and tubulin binding. It is involved in KLF4 glutamylation, which impedes its ubiquitination, thereby leading to somatic cell reprogramming, pluripotency maintenance, and embryogenesis [41]. The potential role of the TTLL4 gene in meat pH needs further investigation.

4. Conclusions

This study conducted a GWAS for eight meat quality traits in a 1518 DLY pig population, and provided valuable insights into elucidating the genetic architecture of meat pH and meat color. We identify one significant QTL on SSC15 with a 143 kb interval affecting pork pH, and two related genes (TNS1, TTLL4) are likely to be responsible for the meat pH difference among individuals according to their functions. From this, we detect several genetic markers affecting meat pH and color. Results from this study are useful for the genetic improvement of pork meat quality in swine by assigning higher weights to associated SNPs in genomic selection.

Author Contributions

Conceptualization, J.Y. and Z.W.; data curation, Z.Z., J.W., C.X., D.R., Y.Q., S.Z., R.D., J.Q., M.Y., and E.Z.; investigation, J.W., C.X., D.R., Y.Q., S.Z., and E.Z.; software and visualization, Z.Z., J.W., R.D., and J.Q.; formal analysis, Z.Z. and J.Y.; writing—original draft, Z.Z. and J.Y.; writing—review and editing, J.Y. and Z.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by a Local Innovative and Research Teams Project of Guangdong Province (2019BT02N630), a Natural Science Foundation of China project (31972540), and a Natural Science Foundation of Guangdong Province project (2018B030313011). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Institutional Review Board Statement

All animals used in this study met the guidelines for the care and use of experimental animals established by the Ministry of Agriculture of China. Tissue samples from pigs were collected with the approval of the ethics committee of South China Agricultural University (Guangzhou, China) under 2018F098.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article.

Conflicts of Interest

R.D. is from the Guangdong Zhongxin Breeding Technology Co., Ltd. (Guangzhou, China). The company had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript and did not affect the results and the conclusions. There are no potential conflicts of interest.

References

  1. Shen, L.; Gan, M.; Chen, L.; Zhao, Y.; Niu, L.; Tang, G.; Jiang, Y.; Zhang, T.; Zhang, S.; Zhu, L. miR-152 targets pyruvate kinase to regulate the glycolytic activity of pig skeletal muscles and affects pork quality. Meat Sci. 2022, 185, 108707. [Google Scholar] [CrossRef]
  2. Hamoen, J.R.; Vollebregt, H.M.; van der Sman, R.G. Prediction of the time evolution of pH in meat. Food Chem. 2013, 141, 2363–2372. [Google Scholar] [CrossRef]
  3. Mancini, R.A.; Hunt, M.C. Current research in meat color. Meat Sci. 2005, 71, 100–121. [Google Scholar] [CrossRef]
  4. Zhang, F.; Zhu, F.; Yang, F.X.; Hao, J.P.; Hou, Z.C. Genomic selection for meat quality traits in Pekin duck. Anim. Genet. 2022, 53, 94–100. [Google Scholar] [CrossRef]
  5. Lopes, F.B.; Baldi, F.; Passafaro, T.L.; Brunes, L.C.; Costa, M.F.O.; Eifert, E.C.; Narciso, M.G.; Rosa, G.J.M.; Lobo, R.B.; Magnabosco, C.U. Genome-enabled prediction of meat and carcass traits using Bayesian regression, single-step genomic best linear unbiased prediction and blending methods in Nelore cattle. Animal 2021, 15, 100006. [Google Scholar] [CrossRef]
  6. Magalhaes, A.F.B.; Schenkel, F.S.; Garcia, D.A.; Gordo, D.G.M.; Tonussi, R.L.; Espigolan, R.; Silva, R.M.O.; Braz, C.U.; Fernandes Junior, G.A.; Baldi, F.; et al. Genomic selection for meat quality traits in Nelore cattle. Meat Sci. 2019, 148, 32–37. [Google Scholar] [CrossRef]
  7. Khanal, P.; Maltecca, C.; Schwab, C.; Gray, K.; Tiezzi, F. Genetic parameters of meat quality, carcass composition, and growth traits in commercial swine. J. Anim. Sci. 2019, 97, 3669–3683. [Google Scholar] [CrossRef]
  8. Miar, Y.; Plastow, G.S.; Moore, S.S.; Manafiazar, G.; Charagu, P.; Kemp, R.A.; Van Haandel, B.; Huisman, A.E.; Zhang, C.Y.; McKay, R.M.; et al. Genetic and phenotypic parameters for carcass and meat quality traits in commercial crossbred pigs. J. Anim. Sci. 2014, 92, 2869–2884. [Google Scholar] [CrossRef] [Green Version]
  9. Davoli, R.; Braglia, S. Molecular approaches in pig breeding to improve meat quality. Brief. Funct. Genom. Proteom. 2007, 6, 313–321. [Google Scholar] [CrossRef]
  10. Uffelmann, E.; Huang, Q.Q.; Munung, N.S.; de Vries, J.; Okada, Y.; Martin, A.R.; Martin, H.C.; Lappalainen, T.; Posthuma, D. Genome-wide association studies. Nat. Rev. Methods Prim. 2021, 1, 59. [Google Scholar] [CrossRef]
  11. Hu, Z.L.; Park, C.A.; Reecy, J.M. Bringing the Animal QTLdb and CorrDB into the future: Meeting new challenges and providing updated services. Nucleic Acids Res. 2022, 50, D956–D961. [Google Scholar] [CrossRef]
  12. Otto, G.; Roehe, R.; Looft, H.; Thoelking, L.; Knap, P.W.; Rothschild, M.F.; Plastow, G.S.; Kalm, E. Associations of DNA markers with meat quality traits in pigs with emphasis on drip loss. Meat Sci. 2007, 75, 185–195. [Google Scholar] [CrossRef]
  13. Ciobanu, D.; Bastiaansen, J.; Malek, M.; Helm, J.; Woollard, J.; Plastow, G.; Rothschild, M. Evidence for new alleles in the protein kinase adenosine monophosphate-activated gamma(3)-subunit gene associated with low glycogen content in pig skeletal muscle and improved meat quality. Genetics 2001, 159, 1151–1162. [Google Scholar] [CrossRef]
  14. Milan, D.; Jeon, J.T.; Looft, C.; Amarger, V.; Robic, A.; Thelander, M.; Rogel-Gaillard, C.; Paul, S.; Iannuccelli, N.; Rask, L.; et al. A mutation in PRKAG3 associated with excess glycogen content in pig skeletal muscle. Science 2000, 288, 1248–1251. [Google Scholar] [CrossRef]
  15. Ma, J.; Yang, J.; Zhou, L.; Ren, J.; Liu, X.; Zhang, H.; Yang, B.; Zhang, Z.; Ma, H.; Xie, X.; et al. A splice mutation in the PHKG1 gene causes high glycogen content and low meat quality in pig skeletal muscle. PLoS Genet. 2014, 10, e1004710. [Google Scholar] [CrossRef] [Green Version]
  16. Liu, X.; Zhou, L.; Xie, X.; Wu, Z.; Xiong, X.; Zhang, Z.; Yang, J.; Xiao, S.; Zhou, M.; Ma, J.; et al. Muscle glycogen level and occurrence of acid meat in commercial hybrid pigs are regulated by two low-frequency causal variants with large effects and multiple common variants with small effects. Genet. Sel. Evol. GSE 2019, 51, 46. [Google Scholar] [CrossRef] [Green Version]
  17. Xiong, X.; Liu, X.; Zhou, L.; Yang, J.; Yang, B.; Ma, H.; Xie, X.; Huang, Y.; Fang, S.; Xiao, S.; et al. Genome-wide association analysis reveals genetic loci and candidate genes for meat quality traits in Chinese Laiwu pigs. Mamm. Genome Off. J. Int. Mamm. Genome Soc. 2015, 26, 181–190. [Google Scholar] [CrossRef]
  18. Cho, I.C.; Yoo, C.K.; Lee, J.B.; Jung, E.J.; Han, S.H.; Lee, S.S.; Ko, M.S.; Lim, H.T.; Park, H.B. Genome-wide QTL analysis of meat quality-related traits in a large F2 intercross between Landrace and Korean native pigs. Genet. Sel. Evol. GSE 2015, 47, 7. [Google Scholar] [CrossRef] [Green Version]
  19. Huang, Y.; Zhou, L.; Zhang, J.; Liu, X.; Zhang, Y.; Cai, L.; Zhang, W.; Cui, L.; Yang, J.; Ji, J.; et al. A large-scale comparison of meat quality and intramuscular fatty acid composition among three Chinese indigenous pig breeds. Meat Sci. 2020, 168, 108182. [Google Scholar] [CrossRef]
  20. Ozsensoy, Y.; Sahin, S. Comparison of different DNA isolation methods and use of dodecyle trimethyl ammonium bromide (DTAB) for the isolation of DNA from meat products. J. Adv. Vet. Anim. Res. 2016, 3, 368–374. [Google Scholar] [CrossRef]
  21. Purcell, S.; Neale, B.; Todd-Brown, K.; Thomas, L.; Ferreira, M.A.; Bender, D.; Maller, J.; Sklar, P.; de Bakker, P.I.; Daly, M.J.; et al. PLINK: A tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 2007, 81, 559–575. [Google Scholar] [CrossRef] [Green Version]
  22. Yang, J.; Lee, S.H.; Goddard, M.E.; Visscher, P.M. GCTA: A tool for genome-wide complex trait analysis. Am. J. Hum. Genet. 2011, 88, 76–82. [Google Scholar] [CrossRef] [Green Version]
  23. Zhou, X.; Stephens, M. Genome-wide efficient mixed-model analysis for association studies. Nat. Genet. 2012, 44, 821–824. [Google Scholar] [CrossRef] [Green Version]
  24. Barrett, J.C.; Fry, B.; Maller, J.; Daly, M.J. Haploview: Analysis and visualization of LD and haplotype maps. Bioinformatics 2005, 21, 263–265. [Google Scholar] [CrossRef] [Green Version]
  25. Ai, H.; Huang, L.; Ren, J. Genetic diversity, linkage disequilibrium and selection signatures in chinese and Western pigs revealed by genome-wide SNP markers. PLoS ONE 2013, 8, e56001. [Google Scholar] [CrossRef] [Green Version]
  26. Ji, J.; Zhou, L.; Huang, Y.; Zheng, M.; Liu, X.; Zhang, Y.; Huang, C.; Peng, S.; Zeng, Q.; Zhong, L.; et al. A whole-genome sequence based association study on pork eating quality traits and cooking loss in a specially designed heterogeneous F6 pig population. Meat Sci. 2018, 146, 160–167. [Google Scholar] [CrossRef]
  27. Unlu, G.; Gamazon, E.R.; Qi, X.; Levic, D.S.; Bastarache, L.; Denny, J.C.; Roden, D.M.; Mayzus, I.; Breyer, M.; Zhong, X.; et al. GRIK5 Genetically Regulated Expression Associated with Eye and Vascular Phenomes: Discovery through Iteration among Biobanks, Electronic Health Records, and Zebrafish. Am. J. Hum. Genet. 2019, 104, 503–519. [Google Scholar] [CrossRef] [Green Version]
  28. Yamamoto, K.; Gandin, V.; Sasaki, M.; McCracken, S.; Li, W.; Silvester, J.L.; Elia, A.J.; Wang, F.; Wakutani, Y.; Alexandrova, R.; et al. Largen: A molecular regulator of mammalian cell size control. Mol. Cell 2014, 53, 904–915. [Google Scholar] [CrossRef] [Green Version]
  29. Carbon, S.; Ireland, A.; Mungall, C.J.; Shu, S.; Marshall, B.; Lewis, S.; Ami, G.O.H.; Web Presence Working Group. AmiGO: Online access to ontology and annotation data. Bioinformatics 2009, 25, 288–289. [Google Scholar] [CrossRef] [Green Version]
  30. Srivastava, S.; Srikanth, K.; Won, S.; Son, J.H.; Park, J.E.; Park, W.; Chai, H.H.; Lim, D. Haplotype-Based Genome-Wide Association Study and Identification of Candidate Genes Associated with Carcass Traits in Hanwoo Cattle. Genes 2020, 11, 551. [Google Scholar] [CrossRef]
  31. Mateescu, R.G.; Garrick, D.J.; Reecy, J.M. Network Analysis Reveals Putative Genes Affecting Meat Quality in Angus Cattle. Front. Genet. 2017, 8, 171. [Google Scholar] [CrossRef] [Green Version]
  32. Bhat, Z.F.; Morton, J.D.; Mason, S.L.; Bekhit, A.E.A. Applied and Emerging Methods for Meat Tenderization: A Comparative Perspective. Compr. Rev. Food Sci. Food Saf. 2018, 17, 841–859. [Google Scholar] [CrossRef] [Green Version]
  33. Jonnalagadda, V.G.; Ram Raju, A.V.; Pittala, S.; Shaik, A.; Selkar, N.A. The prelude on novel receptor and ligand targets involved in the treatment of diabetes mellitus. Adv. Pharm. Bull. 2014, 4, 209–217. [Google Scholar]
  34. Schaeren-Wiemers, N.; Bonnet, A.; Erb, M.; Erne, B.; Bartsch, U.; Kern, F.; Mantei, N.; Sherman, D.; Suter, U. The raft-associated protein MAL is required for maintenance of proper axon--glia interactions in the central nervous system. J. Cell Biol. 2004, 166, 731–742. [Google Scholar] [CrossRef] [Green Version]
  35. Guo, Q.; Xuan, M.F.; Luo, Z.B.; Wang, J.X.; Jin, S.S.; Yin, X.J.; Kang, J.D. Baicalin improves IVM of pig oocytes and subsequent preimplantation embryo development by inhibiting apoptosis. Reprod. Fertil. Dev. 2019, 31, 983–992. [Google Scholar] [CrossRef]
  36. Knudson, C.B.; Knudson, W. Cartilage proteoglycans. Semin. Cell Dev. Biol. 2001, 12, 69–78. [Google Scholar] [CrossRef] [Green Version]
  37. Liu, H.; Hou, L.; Zhou, W.; Wang, B.; Han, P.; Gao, C.; Niu, P.; Zhang, Z.; Li, Q.; Huang, R.; et al. Genome-Wide Association Study and FST Analysis Reveal Four Quantitative Trait Loci and Six Candidate Genes for Meat Color in Pigs. Front. Genet. 2022, 13, 768710. [Google Scholar] [CrossRef]
  38. Verardo, L.L.; Sevon-Aimonen, M.L.; Serenius, T.; Hietakangas, V.; Uimari, P. Whole-genome association analysis of pork meat pH revealed three significant regions and several potential genes in Finnish Yorkshire pigs. BMC Genet. 2017, 18, 13. [Google Scholar] [CrossRef] [Green Version]
  39. Piórkowska, K.; Tyra, M.; Żukowski, K.; Ropka-Molik, K. The SSC15 QTL-Rich Region Mutations Affecting Intramuscular Fat and Production Traits in Pigs. Ann. Anim. Sci. 2020, 20, 425–444. [Google Scholar] [CrossRef]
  40. Bernau, K.; Torr, E.E.; Evans, M.D.; Aoki, J.K.; Ngam, C.R.; Sandbo, N. Tensin 1 Is Essential for Myofibroblast Differentiation and Extracellular Matrix Formation. Am. J. Respir. Cell Mol. Biol. 2017, 56, 465–476. [Google Scholar] [CrossRef] [Green Version]
  41. Ye, B.; Liu, B.; Hao, L.; Zhu, X.; Yang, L.; Wang, S.; Xia, P.; Du, Y.; Meng, S.; Huang, G.; et al. Klf4 glutamylation is required for cell reprogramming and early embryonic development in mice. Nat. Commun. 2018, 9, 1261. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Linkage disequilibrium decay across the whole genome of the association panel. The red line represents the linkage disequilibrium threshold for the association panel (r2 = 0.2).
Figure 1. Linkage disequilibrium decay across the whole genome of the association panel. The red line represents the linkage disequilibrium threshold for the association panel (r2 = 0.2).
Foods 11 03143 g001
Figure 2. Manhattan plots of the GWAS and Q–Q plots for meat pH traits. (A) GWAS for pH_45min; (B) Q–Q plot for pH_45min; (C) GWAS for pH_12h; (D) Q–Q plot for pH_12h. The x-axis represents the chromosomes, and the y-axis represents the −log10 (p-value).
Figure 2. Manhattan plots of the GWAS and Q–Q plots for meat pH traits. (A) GWAS for pH_45min; (B) Q–Q plot for pH_45min; (C) GWAS for pH_12h; (D) Q–Q plot for pH_12h. The x-axis represents the chromosomes, and the y-axis represents the −log10 (p-value).
Foods 11 03143 g002
Figure 3. Manhattan plots of the GWAS and Q–Q plots for meat color traits (45 min). (A) GWAS for L*_45min (lightness_45min); (B) Q–Q plot for L*_45min (lightness_45min); (C) GWAS for a*_45min (redness_45min); (D) Q–Q plot for a*_45min (redness_45min); (E) GWAS for b*_45min (yellowness_45min); (F) Q–Q plot for b*_45min (yellowness_45min). The x-axis represents the chromosomes, and the y-axis represents the −log10 (p-value).
Figure 3. Manhattan plots of the GWAS and Q–Q plots for meat color traits (45 min). (A) GWAS for L*_45min (lightness_45min); (B) Q–Q plot for L*_45min (lightness_45min); (C) GWAS for a*_45min (redness_45min); (D) Q–Q plot for a*_45min (redness_45min); (E) GWAS for b*_45min (yellowness_45min); (F) Q–Q plot for b*_45min (yellowness_45min). The x-axis represents the chromosomes, and the y-axis represents the −log10 (p-value).
Foods 11 03143 g003
Figure 4. Manhattan plots of the GWAS and Q–Q plots for meat color traits (12 h). (A) GWAS for L*_12h (lightness_12h); (B) Q–Q plot for L*_12h (lightness_12h); (C) GWAS for a*_12h (redness_12h); (D) Q–Q plot for a*_12h (redness_12h; (E) GWAS for b*_12h (yellowness_12h); (F) Q–Q plot for b*_12h (yellowness_12h). The x-axis represents the chromosomes, and the y-axis represents the -log10 (p-value).
Figure 4. Manhattan plots of the GWAS and Q–Q plots for meat color traits (12 h). (A) GWAS for L*_12h (lightness_12h); (B) Q–Q plot for L*_12h (lightness_12h); (C) GWAS for a*_12h (redness_12h); (D) Q–Q plot for a*_12h (redness_12h; (E) GWAS for b*_12h (yellowness_12h); (F) Q–Q plot for b*_12h (yellowness_12h). The x-axis represents the chromosomes, and the y-axis represents the -log10 (p-value).
Foods 11 03143 g004
Figure 5. (A) Regional plots of rs335443100 at 117.25–122.99 Mb on SSC15 for meat pH_12h in DLY pigs. (B) represents the 143 kb linkage disequilibrium block in the significant region on SSC15. (C) Violin plot showing the differences in meat pH among three genotypes of the GWAS top SNP (rs335443100). (D) Violin plot showing the differences in meat pH among three genotypes of the SNP explaining the most phenotypic variance within TTLL4 gene.
Figure 5. (A) Regional plots of rs335443100 at 117.25–122.99 Mb on SSC15 for meat pH_12h in DLY pigs. (B) represents the 143 kb linkage disequilibrium block in the significant region on SSC15. (C) Violin plot showing the differences in meat pH among three genotypes of the GWAS top SNP (rs335443100). (D) Violin plot showing the differences in meat pH among three genotypes of the SNP explaining the most phenotypic variance within TTLL4 gene.
Foods 11 03143 g005
Table 1. Summary statistics of meat quality traits in DLY pigs.
Table 1. Summary statistics of meat quality traits in DLY pigs.
TraitsNMean (±SD)h2 (±SE)
pH_45min14806.33 ± 0.320.17 ± 0.03
L*_45min151843.86 ± 2.420.12 ± 0.04
a*_45min1516−0.71 ± 1.020.03 ± 0.03
b*_45min15177.31 ± 1.220.05 ± 0.03
pH_12h14985.73 ± 0.280.12 ± 0.04
L*_12h151548.4 ± 4.050.12 ± 0.04
a*_12h1515−0.46 ± 1.100.06 ± 0.03
b*_12h15127.83 ± 1.450.03 ± 0.03
Table 2. Significant SNPs associated with meat quality traits (45 min) in DLY pigs.
Table 2. Significant SNPs associated with meat quality traits (45 min) in DLY pigs.
TraitsChrSNP IDPosition (bp)MAFp-ValueR2 (%) aNearest GeneDistance (bp)
pH_45min6rs81274518498833730.36 5.56 × 10−62.36 GRIK5within
16rs81324442184783680.46 3.46 × 10−61.27 MTMR12within
L*_45min4rs80971313630760560.40 1.70 × 10−50.90 KCNB2within
a*_45min8rs33272607954038890.22 3.09 × 10−60.10 STK32B14440
8rs81400902643885310.28 7.81 × 10−60.47 //
14rs809449121139384950.11 2.81 × 10−60.01 CNNM2within
b*_45min11rs342146775698989910.42 2.59 × 10−61.01 NALCNwithin
13rs34310358035463650.12 3.95 × 10−60.37 RFTN1within
a Proportion of total phenotypic variation explained by each SNP.
Table 3. Significant SNPs associated with meat quality traits (12 h) in DLY pigs.
Table 3. Significant SNPs associated with meat quality traits (12 h) in DLY pigs.
TraitsChrSNP IDPosition (bp)MAFp-ValueR2 (%)aNearest GeneDistance (bp)
pH_12h2rs813036311236002920.26 2.40 × 10−70.72 FAM170A133583
2rs812954721239983740.12 1.86 × 10−50.31 PRR16201134
9rs8131623068152820.23 6.68 × 10−60.31 PHOX2A2903
15rs814546721199952030.30 2.91 × 10−51.33 TNS1within
15rs814547301200833970.29 1.59 × 10−51.61 TNS1within
15rs808186101201060660.48 6.13 × 10−61.68 TNS1within
15rs3354431001201218910.30 4.62 × 10−71.94 TNS1within
15rs809173551202136660.35 1.43 × 10−51.52 RUFY4204
15rs3382386421203378150.33 1.63 × 10−51.61 GPBAR1within
15rs808167881206963510.33 1.10 × 10−62.52 //
15rs3201303591206991440.33 8.96 × 10−72.52 //
15/1207705900.32 7.44 × 10−72.67 TTLL42595
15/1208012380.32 8.91 × 10−72.67 TTLL4within
15rs3453185431209386020.35 5.27 × 10−62.04 WNT10A172
15/1209824520.32 2.93 × 10−51.89 CDK5R22793
15rs812186481210143410.28 6.87 × 10−62.46 CRYBA2within
L*_12h9rs8131623068152820.23 1.38 × 10−50.40 PHOX2A2903
11rs80993821716196770.11 2.95 × 10−50.62 SLC10A2270898
14rs80985792273422140.20 3.32 × 10−50.98 //
a*_12h2rs81303631*1236002920.26 9.18 × 10−50.11 FAM170A133583
b*_12h4rs343786555176786260.08 3.84 × 10−70.01 SHAS23040
5rs346116771920031970.05 9.56 × 10−70.09 EPYC125750
a Proportion of total phenotypic variation explained by each SNP.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhuang, Z.; Wu, J.; Xu, C.; Ruan, D.; Qiu, Y.; Zhou, S.; Ding, R.; Quan, J.; Yang, M.; Zheng, E.; et al. The Genetic Architecture of Meat Quality Traits in a Crossbred Commercial Pig Population. Foods 2022, 11, 3143. https://doi.org/10.3390/foods11193143

AMA Style

Zhuang Z, Wu J, Xu C, Ruan D, Qiu Y, Zhou S, Ding R, Quan J, Yang M, Zheng E, et al. The Genetic Architecture of Meat Quality Traits in a Crossbred Commercial Pig Population. Foods. 2022; 11(19):3143. https://doi.org/10.3390/foods11193143

Chicago/Turabian Style

Zhuang, Zhanwei, Jie Wu, Cineng Xu, Donglin Ruan, Yibin Qiu, Shenping Zhou, Rongrong Ding, Jianping Quan, Ming Yang, Enqin Zheng, and et al. 2022. "The Genetic Architecture of Meat Quality Traits in a Crossbred Commercial Pig Population" Foods 11, no. 19: 3143. https://doi.org/10.3390/foods11193143

APA Style

Zhuang, Z., Wu, J., Xu, C., Ruan, D., Qiu, Y., Zhou, S., Ding, R., Quan, J., Yang, M., Zheng, E., Wu, Z., & Yang, J. (2022). The Genetic Architecture of Meat Quality Traits in a Crossbred Commercial Pig Population. Foods, 11(19), 3143. https://doi.org/10.3390/foods11193143

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop