Quantitative trait loci analysis of swine meat quality traits

Transcriptome analysis in comparing carcass and meat quality traits of Jiaxing Black Pig and Duroc × Duroc × Berkshire × Jiaxing Black Pig crosses

This study compares two typical strains: Chinese local excellent meat quality of Jiaxing Black (JXB) Pig and quadratic crossbred pig strain Duroc × Duroc × Berkshire × Jiaxing Black (DDBJ). It was found that between the two pig strains, carcass traits and meat quality traits differed significantly. This is exemplified by the leanness and dressing out percent of DDBJ that were significantly higher than JXB pigs of the same age (P < 0.05) and the better growth rate of DDBJ pigs as to JXB pigs was shown by quantifying muscle proliferation and differentiation of longissimus dorsi muscle employing Hematoxylin and Eosin staining of longissimus dorsi muscle. Nutrients such as inosinic acid, intramuscular fat, and free amino acids in the longissimus dorsi muscle were significantly higher in JXB pigs than DDBJ pigs (p < 0.0001); saturated fatty acids were higher in JXB than in DDBJ pigs (p = 0.0097); essential amino acids and fresh taste amino acids (serine, glutamic acid, proline, glycine, alanine) of JXB pigs was higher than that of DDBJ pigs (p < 0.0001) and amino acids in longissimus dorsi muscle of JXB pigs surpasses the amino acid concentration of DDBJ pigs (p < 0.0001), thus showing the superiority of JXB in terms of meat quality. However, the content of polyunsaturated fatty acids, which is responsible for poor meat quality, was significantly higher in the longissimus dorsi muscle of DDBJ pig than JXB pigs (p < 0.0001); RNA-seq analysis of 5 biological replicates from two of the strains was performed. The screening of 164 up-regulated genes and 183 down-regulated genes found in longissimus dorsi muscle of DDBJ was done and the results identified differentially expressed genes related to muscle development, adipogenesis, amino acid metabolism, fatty acid metabolism and inosine synthesis. In conclusion, the study identified functional genes, elucidated the mechanisms associated with carcass quality traits, meat quality traits and other related traits, and provided means of genetic enhancement to improve meat quality traits and carcass traits in Chinese commercial pigs.

Combined GWAS and LDLA approaches to improve genome-wide quantitative trait loci detection affecting carcass and meat quality traits in pig

Many QTL affecting meat quality and carcass traits have been reported. However, in most of the cases these QTL have been detected in non-commercial populations. Therefore, a family structured population of 457 F2 pigs issued from an inter-cross between 2 commercial sire lines was used to detect QTL affecting meat quality and carcass traits. All animals were genotyped using the Illumina PorcineSNP60 BeadChip platform. Genome-wide association studies were used in combination with linkage disequilibrium-linkage analysis to identify QTL. A total of 32 QTL were detected. Nine of these QTL exceeded the genome-wide 5% significance threshold. We detected 18 QTL affecting carcass composition traits and 16 QTL affecting meat quality traits. Using post-QTL bioinformatics analysis we highlighted 26 functional candidate genes related to fatness, muscle development, meat color and meat pH. Finally, our results shed light on the advantage of using different QTL detection methodologies to get a global overview of the QTL present in the studied population.

Genome-wide association study in an F2 Duroc x Pietrain resource population for economically important meat quality and carcass traits

Meat quality is essential for consumer acceptance, it ultimately impacts pork production profitability and it is subject to genetic control. The objective of this study was to map genomic regions associated with economically important meat quality and carcass traits. We performed a genome-wide association (GWA) analysis to map regions associated with 38 meat quality and carcass traits recorded for 948 F2 pigs from the Michigan State University Duroc × Pietrain resource population. The F0, F1, and 336 F2 pigs were genotyped with the Illumina Porcine SNP60 BeadChip, while the remaining F2 pigs were genotyped with the GeneSeek Genomic Profiler for Porcine Low Desnisty (LD) chip, and imputed with high accuracy ( = 0.97). Altogether the genomic dataset comprised 1,019 animals and 44,911 SNP. A Gaussian linear mixed model was fitted to estimate the breeding values and the variance components. A linear transformation was performed to estimate the marker effects and variances. Type I error rate was controlled at a False Discovery Rate of 5%. Seven putative QTL found in this study were previously reported in other studies. Two novel QTL associated with tenderness (TEN) were located on SSC3 [135.6:137.5Mb; False Discovery rate (FDR) < 0.03] and SSC5 (67.3:69.1Mb; FDR < 0.02). The QTL region identified on SSC15 includes Protein Kinase AMP-activated ɣ 3-subunit gene (), which has been associated with 24-h pH (pH24), drip loss (DL) and cook yield (CY). Also, novel candidate genes were identified for TEN in the region on SSC5 [A Kinase (PRKA) Anchor Protein 3 (], and for tenth rib backfat thickness (BF10) [Carnitine O-Acetyltransferase ()] on SSC1. The association of gene polymorphisms with pork quality traits has been reported for several pig populations. However, there are no SNP for this gene on the chip used, thus we genotyped the animals for 2 non-synonymous variants ( and ). We then performed a GWA conditioning on the genotype of both SNP and was associated with pH24, DL, protein content (PRO) and CY ( < 0.004) and T30N with Juiciness, TEN, shear force, pH24, PRO, and CY < 0.04). Finally, we performed a GWA conditioning on the genotype of the SNP peak detected in this study, and T30N remained associated only with PRO ( < 0.02). Therefore, in this study we identified 2 novel QTL regions, suggest 2 novel candidate genes, and conclude that other SNP in PRKAG3 or nearby gene(s) explain the observed associations on SSC15 in this population.

Identification of single nucleotide polymorphisms in carnosine-related genes and effects of genotypes on pork meat quality attributes

Carnosine has pH-buffering and antioxidant properties that may bring advantages in terms of meat quality attributes. This study aimed at identifying polymorphisms in carnosine-related genes (CARNS1, SLC6A6, SLC15A3, SLC15A4) that might associate with muscle carnosine content and meat quality traits in pigs (Duroc, Landrace, Yorkshire). Twenty seven SNPs were identified and association analyses performed for SLC15A3 c.*35C>T and c.*52C>T (3' UTR region), and SLC15A4 c.658A>G (Ile220Val) and c.818G>A (Ser273Asn) SNPs. Associations were observed for SNP c.658A>G with carnosine content, color b* and L*, drip and cooking losses, pH24h and glycolytic potential values (P≤0.05). The same associations were observed for SNP c.818G>A, but they were not significant after FDR correction. Results suggest that specific SLC15A4 gene variants might increase muscle carnosine content and improve meat quality. With a minor allele frequency of 0.17 for SNP c.658A>G in Yorkshire pigs, selection in favor of the c.658A allele may be considered as a mean to improve pork quality attributes.

Joint QTL mapping and gene expression analysis identify positional candidate genes influencing pork quality traits

Meat quality traits have an increasing importance in the pig industry because of their strong impact on consumer acceptance. Herewith, we have combined phenotypic and microarray expression data to map loci with potential effects on five meat quality traits recorded in the longissimus dorsi (LD) and gluteus medius (GM) muscles of 350 Duroc pigs, i.e. pH at 24 hours post-mortem (pH₂₄), electric conductivity (CE) and muscle redness (a*), lightness (L*) and yellowness (b*). We have found significant genome-wide associations for CE of LD on SSC4 (~104 Mb), SSC5 (~15 Mb) and SSC13 (~137 Mb), while several additional regions were significantly associated with meat quality traits at the chromosome-wide level. There was a low positional concordance between the associations found for LD and GM traits, a feature that reflects the existence of differences in the genetic determinism of meat quality phenotypes in these two muscles. The performance of an eQTL search for SNPs mapping to the regions associated with meat quality traits demonstrated that the GM a* SSC3 and pH₂₄ SSC17 QTL display positional concordance with cis-eQTL regulating the expression of several genes with a potential role on muscle metabolism.

Characterization and management of long runs of homozygosity in parental nucleus lines and their associated crossbred progeny

BACKGROUND

In nucleus populations, regions of the genome that have a high frequency of runs of homozygosity (ROH) occur and are associated with a reduction in genetic diversity, as well as adverse effects on fitness. It is currently unclear whether, and to what extent, ROH stretches persist in the crossbred genome and how genomic management in the nucleus population might impact low diversity regions and its implications on the crossbred genome.

METHODS

We calculated a ROH statistic based on lengths of 5 (ROH5) or 10 (ROH10) Mb across the genome for genotyped Landrace (LA), Large White (LW) and Duroc (DU) dams. We simulated crossbred dam (LA × LW) and market [DU × (LA × LW)] animal genotypes based on observed parental genotypes and the ROH frequency was tabulated. We conducted a simulation using observed genotypes to determine the impact of minimizing parental relationships on multiple diversity metrics within nucleus herds, i.e. pedigree-(A), SNP-by-SNP relationship matrix or ROH relationship matrix. Genome-wide metrics included, pedigree inbreeding, heterozygosity and proportion of the genome in ROH of at least 5 Mb. Lastly, the genome was split into bins of increasing ROH5 frequency and, within each bin, heterozygosity, ROH5 and length (Mb) of ROH were evaluated.

RESULTS

We detected regions showing high frequencies of either ROH5 and/or ROH10 across both LW and LA on SSC1, SSC4, and SSC14, and across all breeds on SSC9. Long haplotypes were shared across parental breeds and thus, regions of ROH persisted in crossbred animals. Averaged across replicates and breeds, progeny had higher levels of heterozygosity (0.0056 ± 0.002%) and lower proportion of the genome in a ROH of at least 5 Mb (-0.015 ± 0.003%) than their parental genomes when genomic relationships were constrained, while pedigree relationships resulted in negligible differences at the genomic level. Across all breeds, only genomic data was able to target low diversity regions.

CONCLUSIONS

We show that long stretches of ROH present in the parents persist in crossbred animals. Furthermore, compared to using pedigree relationships, using genomic information to constrain parental relationships resulted in maintaining more genetic diversity and more effectively targeted low diversity regions.

Genome-wide Association Study to Identify Quantitative Trait Loci for Meat and Carcass Quality Traits in Berkshire

Meat and carcass quality attributes are of crucial importance influencing consumer preference and profitability in the pork industry. A set of 400 Berkshire pigs were collected from Dasan breeding farm, Namwon, Chonbuk province, Korea that were born between 2012 and 2013. To perform genome wide association studies (GWAS), eleven meat and carcass quality traits were considered, including carcass weight, backfat thickness, pH value after 24 hours (pH24), Commission Internationale de l’Eclairage lightness in meat color (CIE L), redness in meat color (CIE a), yellowness in meat color (CIE b), filtering, drip loss, heat loss, shear force and marbling score. All of the 400 animals were genotyped with the Porcine 62K SNP BeadChips (Illumina Inc., USA). A SAS general linear model procedure (SAS version 9.2) was used to pre-adjust the animal phenotypes before GWAS with sire and sex effects as fixed effects and slaughter age as a covariate. After fitting the fixed and covariate factors in the model, the residuals of the phenotype regressed on additive effects of each single nucleotide polymorphism (SNP) under a linear regression model (PLINK version 1.07). The significant SNPs after permutation testing at a chromosome-wise level were subjected to stepwise regression analysis to determine the best set of SNP markers. A total of 55 significant (p<0.05) SNPs or quantitative trait loci (QTL) were detected on various chromosomes. The QTLs explained from 5.06% to 8.28% of the total phenotypic variation of the traits. Some QTLs with pleiotropic effect were also identified. A pair of significant QTL for pH24 was also found to affect both CIE L and drip loss percentage. The significant QTL after characterization of the functional candidate genes on the QTL or around the QTL region may be effectively and efficiently used in marker assisted selection to achieve enhanced genetic improvement of the trait considered.

Genome-wide association analyses for meat quality traits in Chinese Erhualian pigs and a Western Duroc × (Landrace × Yorkshire) commercial population

BACKGROUND

Understanding the genetic mechanisms that underlie meat quality traits is essential to improve pork quality. To date, most quantitative trait loci (QTL) analyses have been performed on F2 crosses between outbred pig strains and have led to the identification of numerous QTL. However, because linkage disequilibrium is high in such crosses, QTL mapping precision is unsatisfactory and only a few QTL have been found to segregate within outbred strains, which limits their use to improve animal performance. To detect QTL in outbred pig populations of Chinese and Western origins, we performed genome-wide association studies (GWAS) for meat quality traits in Chinese purebred Erhualian pigs and a Western Duroc × (Landrace × Yorkshire) (DLY) commercial population.

METHODS

Three hundred and thirty six Chinese Erhualian and 610 DLY pigs were genotyped using the Illumina PorcineSNP60K Beadchip and evaluated for 20 meat quality traits. After quality control, 35 985 and 56 216 single nucleotide polymorphisms (SNPs) were available for the Chinese Erhualian and DLY datasets, respectively, and were used to perform two separate GWAS. We also performed a meta-analysis that combined P-values and effects of 29 516 SNPs that were common to Erhualian, DLY, F2 and Sutai pig populations.

RESULTS

We detected 28 and nine suggestive SNPs that surpassed the significance level for meat quality in Erhualian and DLY pigs, respectively. Among these SNPs, ss131261254 on pig chromosome 4 (SSC4) was the most significant (P = 7.97E-09) and was associated with drip loss in Erhualian pigs. Our results suggested that at least two QTL on SSC12 and on SSC15 may have pleiotropic effects on several related traits. All the QTL that were detected by GWAS were population-specific, including 12 novel regions. However, the meta-analysis revealed seven novel QTL for meat characteristics, which suggests the existence of common underlying variants that may differ in frequency across populations. These QTL regions contain several relevant candidate genes.

CONCLUSIONS

These findings provide valuable insights into the molecular basis of convergent evolution of meat quality traits in Chinese and Western breeds that show divergent phenotypes. They may contribute to genetic improvement of purebreds for crossbred performance.

A Whole Genome Association Study on Meat Quality Traits Using High Density SNP Chips in a Cross between Korean Native Pig and Landrace

A whole genome association (WGA) study was performed to detect significant polymorphisms for meat quality traits in an F₂ cross population (N = 478) that were generated with Korean native pig sires and Landrace dams in National Livestock Research Institute, Songwhan, Korea. The animals were genotyped using Illumina porcine 60k SNP beadchips, in which a set of 46,865 SNPs were available for the WGA analyses on ten carcass quality traits; live weight, crude protein, crude lipids, crude ash, water holding capacity, drip loss, shear force, CIE L, CIE a and CIE b. Phenotypes were regressed on additive and dominance effects for each SNP using a simple linear regression model, after adjusting for sex, sire and slaughter stage as fixed effects. With the significant SNPs for each trait (p<0.001), a stepwise regression procedure was applied to determine the best set of SNPs with the additive and/or dominance effects. A total of 106 SNPs, or quantitative trait loci (QTL) were detected, and about 32 to 66% of the total phenotypic variation was explained by the significant SNPs for each trait. The QTL were identified in most porcine chromosomes (SSCs), in which majority of the QTL were detected in SSCs 1, 2, 12, 13, 14 and 16. Several QTL clusters were identified on SSCs 12, 16 and 17, and a cluster of QTL influencing crude protein, crude lipid, drip loss, shear force, CIE a and CIE b were located between 20 and 29 Mb of SSC12. A pleiotropic QTL for drip loss, CIE L and CIE b was also detected on SSC16. These QTL need to be validated in commercial pig populations for genetic improvement in meat quality via marker-assisted selection.

Genome-wide association study of meat quality traits in a White Duroc×Erhualian F2 intercross and Chinese Sutai pigs

Thousands of QTLs for meat quality traits have been identified by linkage mapping studies, but most of them lack precise position or replication between populations, which hinder their application in pig breeding programs. To localize QTLs for meat quality traits to precise genomic regions, we performed a genome-wide association (GWA) study using the Illumina PorcineSNP60K Beadchip in two swine populations: 434 Sutai pigs and 933 F2 pigs from a White Duroc×Erhualian intercross. Meat quality traits, including pH, color, drip loss, moisture content, protein content and intramuscular fat content (IMF), marbling and firmness scores in the M. longissimus (LM) and M. semimembranosus (SM) muscles, were recorded on the two populations. In total, 127 chromosome-wide significant SNPs for these traits were identified. Among them, 11 SNPs reached genome-wise significance level, including 1 on SSC3 for pH, 1 on SSC3 and 3 on SSC15 for drip loss, 3 (unmapped) for color a*, and 2 for IMF each on SSC9 and SSCX. Except for 11 unmapped SNPs, 116 significant SNPs fell into 28 genomic regions of approximately 10 Mb or less. Most of these regions corresponded to previously reported QTL regions and spanned smaller intervals than before. The loci on SSC3 and SSC7 appeared to have pleiotropic effects on several related traits. Besides them, a few QTL signals were replicated between the two populations. Further, we identified thirteen new candidate genes for IMF, marbling and firmness, on the basis of their positions, functional annotations and reported expression patterns. The findings will contribute to further identification of the causal mutation underlying these QTLs and future marker-assisted selection in pigs.

Association mapping of quantitative trait loci for carcass and meat quality traits at the central part of chromosome 2 in Italian Large White pigs

Association mapping of the central part of porcine chromosome 2 harboring QTLs for carcass and meat quality traits was performed with 17 gene-tagged SNPs located between 44.0 and 77.5 Mb on a physical map (Sscrofa10.2) in Italian Large White pigs. For the analyzed animals records of estimated breeding values for average daily gain, back fat thickness, lean cuts, ham weight, feed conversion ratio, pH1, pHu, CIE L*, CIE a*, CIE b* and drip loss were available. A significant QTL for fat deposition (adjusted P=0.0081) and pH1 (adjusted P=0.0972) to MYOD1 at position 44.4 Mb and a QTL for growth and meatiness (adjusted P=0.0238-0.0601) to UBL5 at position 68.9 Mb were mapped. These results from association mapping are much more accurate than those from linkage mapping and facilitate further search for position candidate genes and causative mutations needed for application of markers through marker assisted selection.

UBE3B and ZRANB1 polymorphisms and transcript abundance are associated with water holding capacity of porcine M. longissimus dorsi

The degradation of myofibrillar proteins during meat maturation affects the water holding capacity (WHC) of meat. Our study sought to identify polymorphisms in UBE3B and ZRANB1, genes encoding proteins involved in ubiquitination, and to evaluate the relationship between genotype, transcript abundance, and WHC of pork. A single SNP of ZRANB1, c.552A>G (p.Ile153Val), and two silent SNPs of UBE3B, c.1921A>T and c.4292C>T, were associated with muscle pH, conductivity, meat colour, or drip loss in German Landrace (GL, n=266) and Pietrain×(Large White×German Landrace) (PiF1, n=316). Further, carriers of the minor alleles at the SNPs tended to have increased transcript abundance. Consistent with the protein degradation promoting and inhibiting effects of UBE3B and ZRANB1, respectively, and the expected impact on WHC, their expressions were positively and negatively associated with WHC. The results implicate that the SNPs in both genes are in linkage with a causal site that affects transcript abundance and WHC.

The structure of a gene co-expression network reveals biological functions underlying eQTLs

What are the commonalities between genes, whose expression level is partially controlled by eQTL, especially with regard to biological functions? Moreover, how are these genes related to a phenotype of interest? These issues are particularly difficult to address when the genome annotation is incomplete, as is the case for mammalian species. Moreover, the direct link between gene expression and a phenotype of interest may be weak, and thus difficult to handle. In this framework, the use of a co-expression network has proven useful: it is a robust approach for modeling a complex system of genetic regulations, and to infer knowledge for yet unknown genes. In this article, a case study was conducted with a mammalian species. It showed that the use of a co-expression network based on partial correlation, combined with a relevant clustering of nodes, leads to an enrichment of biological functions of around 83%. Moreover, the use of a spatial statistics approach allowed us to superimpose additional information related to a phenotype; this lead to highlighting specific genes or gene clusters that are related to the network structure and the phenotype. Three main results are worth noting: first, key genes were highlighted as a potential focus for forthcoming biological experiments; second, a set of biological functions, which support a list of genes under partial eQTL control, was set up by an overview of the global structure of the gene expression network; third, pH was found correlated with gene clusters, and then with related biological functions, as a result of a spatial analysis of the network topology.

A two-step approach to map quantitative trait loci for meat quality in connected porcine F(2) crosses considering main and epistatic effects

The aim of this study was to map QTL for meat quality traits in three connected porcine F(2) crosses comprising around 1000 individuals. The three crosses were derived from the founder breeds Chinese Meishan, European Wild Boar and Pietrain. The animals were genotyped genomewide for approximately 250 genetic markers, mostly microsatellites. They were phenotyped for seven meat quality traits (pH at 45 min and 24 h after slaughter, conductivity at 45 min and 24 h after slaughter, meat colour, drip loss and rigour). QTL mapping was conducted using a two-step procedure. In the first step, the QTL were mapped using a multi-QTL multi-allele model that was tailored to analyse multiple connected F(2) crosses. It considered additive, dominance and imprinting effects. The major gene RYR1:g.1843C>T affecting the meat quality on SSC6 was included as a cofactor in the model. The mapped QTL were tested for pairwise epistatic effects in the second step. All possible epistatic effects between additive, dominant and imprinting effects were considered, leading to nine orthogonal forms of epistasis. Numerous QTL were found. The most interesting chromosome was SSC6. Not all genetic variance of meat quality was explained by RYR1:g.1843C>T. A small confidence interval was obtained, which facilitated the identification of candidate genes underlying the QTL. Epistasis was significant for the pairwise QTL on SSC12 and SSC14 for pH24 and for the QTL on SSC2 and SSC5 for rigour. Some evidence for additional pairwise epistatic effects was found, although not significant. Imprinting was involved in epistasis.

Effects of polymorphisms in the 3' untranslated region of the porcine PPARGC1A gene on muscle fiber characteristics and meat quality traits

Peroxisome proliferator-activated receptor γ coactivator 1 α (PPARGC1A) is a transcriptional coactivator that is involved in a variety of biological processes including muscle fiber type composition. Here, we identified two single nucleotide polymorphisms (SNPs; *2690T>C and *2864T>C) and one insertion/deletion in the 3' untranslated region of porcine PPARGC1A. These SNPs were genotyped by direct sequencing in a total of 439 pigs representing three different pig breeds (Berkshire, n = 156; Yorkshire, n = 163; Landrace, n = 120). We evaluated the effects of diplotypes of individual PPARGC1A 3'UTR SNPs on muscle fiber characteristics and meat quality traits. The *2690T>C polymorphism was significantly associated with the percentage of type I and IIb fibers for both muscle fiber number and area composition (P < 0.05), and also showed a significant association with muscle pH, a parameter of meat quality (P = 0.0188). The *2864T>C polymorphism was also associated with meat quality traits including muscle pH (P = 0.0071), drip loss (P = 0.0006), and lightness (P = 0.0702), but showed no significant association with muscle fiber characteristics. Interestingly, each SNP affected PPARGC1A expression significantly at the protein level but not at the mRNA level, thereby accounting for phenotypic variability among genotypes. Taken together, our data suggest that the *2690T>C and *2864T>C polymorphisms can be used as genetic markers for selection toward improved meat quality.