Validation of a paternally imprinted QTL affecting pH24h distinct from PRKAG3 on SSC15

We previously performed a genome scan in a White Duroc × Erhualian F(2) population and identified a QTL on SSC15 with strong effect on pH24h in M. Longissimus dorsi and M. Semimembranous muscle tissue. At that time, the mode of inheritance of this QTL was not clarified, and it was also unclear whether the observed QTL effect was completely or partially caused by mutations in the PRKAG3 gene, which is the only major gene on SSC15 so far known to influence pH24h. In this study, effects of the PRKAG3 gene on meat quality traits were estimated by association analyses. Two substitutions in PRKAG3, p.Ile199Val (p.I199V) and p.Thr30Asn (p.T30N), were found to be segregating in the F(2) population and to significantly affect pH24h and total glycogen in meat, respectively. However, we excluded PRKAG3 as a causative gene for the detected QTL based on the following reasons: (i) the gene was located outside of the QTL confidence interval; (ii) when the PRKAG3 substitution was included as a fixed effect in the QTL model, the F-ratio for the QTL increased rather than decreased; (iii) favourable alleles for pH24h at the QTL and at the PRKAG3 p.I199V locus originated from Erhualian and White Duroc founders, respectively; (iv) more importantly, this QTL showed exclusive maternal expression, differing from the Mendelian expression of PRKAG3. In conclusion, this study is the first to report a maternally-expressed QTL for pH24h on SSC15, which is distinct from PRKAG3.     

Characterization of the bovine PRKAG3 gene: structure, polymorphism, and alternative transcripts

The bovine PRKAG3 gene encodes the AMPK gamma3 subunit, one isoform of the regulatory gamma subunit of the AMP-activated protein kinase (AMPK). The AMPK plays a major role in the regulation of energy metabolism and mutations affecting the genes encoding the gamma subunits have been shown to influence AMPK activity. The gamma3 subunit is involved in the regulation of AMPK activity in skeletal muscle and strongly influences glycogen metabolism. Glycogen content in muscle is correlated to meat quality in livestock because it influences postmortem maturation process and ultimate pH. Naturally occurring mutations in the porcine PRKAG3 gene highly affect meat quality by influencing glycogen content before slaughter. We present the characterization of the bovine PRKAG3 gene and a polymorphism analysis in three cattle breeds. Thirty-two SNPs were identified among which 13 are in the coding region, one is in the 3' UTR, and 18 are in the introns. Five of them change an amino acid in the PRKAG3 protein sequence. Allelic frequencies were determined in the three breeds considered, and mutant alleles affecting the coding sequence are found at a very low frequency. Alternative splicing sites were identified at two positions of the gene, introducing heterogeneity in the population of proteins translated from the gene.     

Genetic variations of the porcine PRKAG3 gene in Chinese indigenous pig breeds

Four missense substitutions (T30N, G52S, V199I and R200Q) in the porcine PRKAG3 gene were considered as the likely candidate loci affecting meat quality. In this study, the R200Q substitution was investigated in a sample of 62 individuals from Hampshire, Chinese Min and Erhualian pigs, and the genetic variations of T30N, G52S and V199I substitutions were detected in 1505 individuals from 21 Chinese indigenous breeds, 5 Western commercial pig breeds, and the wild pig. Allele 200R was fixed in Chinese Min and Erhualian pigs. Haplotypes II-QQ and IV-QQ were not observed in the Hampshire population, supporting the hypothesis that allele 200Q is tightly linked with allele 199V. Significant differences in allele frequencies of the three substitutions (T30N, G52S and V199I) between Chinese indigenous pigs and Western commercial pigs were observed. Obvious high frequencies of the "favorable" alleles 30T and 52G in terms of meat quality were detected in Chinese indigenous pigs, which are well known for high meat quality. However, the frequency of the "favorable" allele 199I, which was reported to have a greater effect on meat quality in comparison with 30T and 52G, was very low in all of the Chinese indigenous pigs except for the Min pig. The reasons accounting for this discrepancy remain to be addressed. The presence of the three substitutions in purebred Chinese Tibetan pigs indicates that the three substitutions were ancestral mutations. A novel A/G substitution at position 51 in exon 1 was identified. The results suggest that further studies are required to investigate the associations of these substitutions in the PRKAG3 gene with meat quality of Chinese indigenous pigs, and to uncover other polymorphisms in the PRKAG3 gene with potential effects on meat quality in Chinese indigenous pigs.     

Functional implication of an Arg307Gly substitution in corticosteroid-binding globulin, a candidate gene for a quantitative trait locus associated with cortisol variability and obesity in pig

We previously reported that corticosteroid-binding globulin gene (Cbg) may be the causal gene of a quantitative trait locus associated with cortisol levels, fat deposition, and muscle content in a pig intercross. Sequence analysis of parental animals allowed us to identify four amino-acid substitutions. Here we have examined if any of these single amino acid substitutions could be responsible for the difference in CBG binding and affinity for cortisol between the parental breeds, using in vitro assays of Cbg variants after transfection of mammalian cells. Additionally, the Cbg coding region was analyzed in samples from a synthetic pig line to study association between polymorphism and CBG biochemical properties, carcass composition, and meat quality. Both in vitro transfection assays and the association studies suggest a role of the Arg307Gly mutation in increasing CBG capacity (by >70%) and decreasing CBG affinity for cortisol (by 30%). The Ile265Val substitution may also have an effect on decreasing CBG affinity for cortisol by 25%. The mutations Ser15Ile and Thr257Met do not seem to have an effect on CBG parameters. The Arg307Gly substitution was the only mutation associated with a parameter of meat quality and no mutation was linked to carcass composition.     

Determination of PRKAG1 coding sequence and mapping of PRKAG1 and PRKAG2 relatively to porcine back fat thickness QTL

PRKAG1, PRKAG2 and PRKAG3 encode three isoforms of AMP-activated protein kinase gamma chain. A major effect on meat quality and a medium effect on back fat thickness of the RN- mutation in the PRKAG3 gene has previously been reported. We have now mapped PRKAG1 and PRKAG2 at expected locations on SSC5 and SSC18 by analysis of radiation hybrids (IMpRH panel). PRKAG2 has been mapped in a region where no quantitative trait loci (QTL) has been reported. PRKAG1 has been mapped close to (but probably outside) a region containing a QTL influencing fatness traits. We have determined the full coding sequence of PRKAG1. No missense mutation was identified when comparing the coding sequence of one Meishan and one Large White boars. Further work is, however, required to determine if a polymorphism in PRKAG1 could be responsible for a part of the variability observed on fatness traits.     

UDP-glucose pyrophosphorylase is upregulated in carriers of the porcine RN- mutation in the AMP-activated protein kinase

The AMP-activated protein kinase (AMPK) plays a key role in the regulation of energy metabolism in eukaryotic cells acting as a metabolic sensor. In its activated form AMPK inhibits ATP consuming pathways and stimulates ATP generating pathways. A dominant mutation, denoted RN(-), in the porcine PRKAG3 gene, encoding the regulatory gamma3 subunit of AMPK, results in hyperaccumulation of glycogen in glycolytic skeletal muscle cells. To study the effects of this mutation on protein expression patterns in skeletal muscle, comparative proteome analysis of muscle samples from 12 animals (6 rn (+)/rn (+) and 6 RN(-)/rn (+)) was performed. The major finding of the proteome analysis was that the key enzyme in the synthesis of glycogen, UDP-glucose pyrophosphorylase, was significantly up-regulated in RN(-) carriers. This observation was subsequently supported by studies of enzyme activity and Northern blot analysis. Furthermore, the expression patterns of enzymes related to glycolysis and the citric acid cycle were also affected. Our data suggests that hyperaccumulation of glycogen mediated by the RN(-) mutation is due to an increased synthesis of glycogen.     

Evaluation of the porcine melanocortin 4 receptor (MC4R) gene as a positional candidate for a fatness QTL in a cross between Landrace and Hampshire

Melanocortin 4 receptor (MC4R) is expressed in the appetite-regulating areas of the brain where it is central in the regulation of feed intake and energy balance. A mutation in MC4R causing an Asp298Asn substitution has been associated with fatness, high daily gain and feed intake in the pig. In a previously performed genome scan based on a Hampshire x Landrace cross, we detected one quantitative trait loci (QTL) affecting carcass fat/meat ratio and one QTL affecting the biceps femoris muscle, both close to the position of MC4R on porcine chromosome 1. In this study, the two lines were found to be close to fixation for alternative alleles of the Asp298Asn polymorphism. Additional QTL analyses supported our hypothesis of MC4R as a positional candidate gene but only for the fat/meat QTL. The Asp298Asn polymorphism was also evaluated as a selection target for daily gain in a Danish pig breeding population that included four breeds (Hampshire, Duroc, Landrace and Yorkshire). Over a 12-year period (1990-2002), a significant increase in the allele frequency of 298Asn was found in Landrace and Duroc, whereas a non-significant decrease in the 298Asn allele frequency was observed in Yorkshire. The Hampshire breed was fixed for the 298Asn allele in 1990. The high 298Asn allele frequencies in Hampshire, Landrace and Duroc are most likely due to selection for daily gain, whereas selection for daily gain in the Yorkshire breed apparently focuses on other loci.     

Corticosteroid binding globulin: a new target for cortisol-driven obesity

We present data suggesting that corticosteroid-binding globulin (CBG) may be the causal gene of a previously identified quantitative trait locus (QTL) associated with cortisol levels, fat, and muscle content in a pig intercross. Because Cbg in human and mouse maps in the region orthologous to the pig region containing this QTL, we considered Cbg as an interesting positional candidate gene because CBG plays a major role in cortisol bioavailability. Firstly, we cloned pig Cbg from a bacterial artificial chromosome library and showed by fluorescent in situ hybridization and radiation hybrid mapping that it maps on 7q26 at the peak of the QTL interval. Secondly, we detected in a subset of the pig intercross progeny a highly significant genetic linkage between CBG plasma binding capacity values and the chromosome 7 markers flanking the cortisol-associated QTL. In this population, CBG capacity is correlated positively to fat and negatively to muscle content. Thirdly, CBG capacity was three times higher in Meishan compared with Large White parental breeds and a 7-fold difference was found in Cbg mRNA expression between the two breeds. Overall, the data accumulated in this study point to Cbg gene as a key regulator of cortisol levels and obesity susceptibility.     

Study of candidate genes for glycolytic potential of porcine skeletal muscle: identification and analysis of mutations, linkage and physical mapping and association with meat quality traits in pigs

Several genes (PRKAA2, PRKAB1, PRKAB2, PRKAG3, GAA, GYS1, PYGM, ALDOA, GPI, LDHA, PGAM2 and PKM2), chosen according to their role in the regulation of the energy balance and in the glycogen metabolism and glycolysis of the skeletal muscle, were studied. Eleven single nucleotide polymorphisms (SNPs) were identified in six of these genes (PRKAB1, GAA, PYGM, LDHA, PGAM2 and PKM2). Allele frequencies were analyzed in seven different pig breeds for these loci and for a polymorphism already described for GPI and for three polymorphic sites already reported at the PRKAG3 locus (T30N, G52S and I199V). Linkage mapping assigned PYGM and LDHA to porcine chromosome (SSC) 2, PKM2 to SSC7, GAA to SSC12, PRKAB1 to SSC14 and PGAM2 to SSC18. Physical mapping, obtained by somatic cell hybrid panel analysis, confirmed the linkage assignments of PRKAB1 and GAA and localized ALDOA, PRKAB2 and GYS1 to SSC3, SSC4 and SSC6, respectively. Pigs selected for the association study, for which several meat quality traits were measured, were first genotyped at the PRKAG3 R200Q polymorphic site (RN locus), in order to exclude carriers of the 200Q allele, and then were genotyped for all the mutations considered in this work. Significant associations (P < or = 0.001) were observed for the PRKAG3 T30N and G52S polymorphic sites with meat colour (L* at 24 h post mortem). PGAM2 and PKM2 were significantly associated (P = 0.01) with drip loss percentage and glycogen content at one hour post mortem, respectively.     

RNA‐seq analysis reveals new candidate genes for drip loss in a Pietrain × Duroc × Landrace × Yorkshire population

Drip loss, one of the most important meat quality traits, is characterized by low heritability. To date, the genetic factors affecting the drip loss trait have not been clearly elucidated. The objective of this study was to identify critical candidate genes affecting drip loss. First, we generated a Pietrain × Duroc × Landrace × Yorkshire commercial pig population and obtained phenotypic values for the drip loss trait. Furthermore, we constructed two RNA libraries from pooled samples of longissimus dorsi muscles with the highest (H group) and lowest (L group) drip loss and identified the differentially expressed genes (DEGs) between these extreme phenotypes using RNA‐seq technology. In total, 25 883 genes were detected in the H and L group libraries, and none was specifically expressed in only one library. Comparative analysis of gene expression levels found that 150 genes were differentially expressed, of which 127 were upregulated and 23 were downregulated in the H group relative to the L group. In addition, 68 drip loss quantitative trait loci (QTL) overlapping with 63 DEGs were identified, and these QTL were distributed mainly on chromosomes 1, 2, 5 and 6. Interestingly, the triadin (TRDN) gene, which is involved in muscle contraction and fat deposition, and the myostatin (MSTN) gene, which has a role in muscle growth, were localized to more than two drip loss QTL, suggesting that both are critical candidate genes responsible for drip loss.     

Fine mapping of porcine chromosome 6 QTL and LEPR effects on body composition in multiple generations of an Iberian by Landrace intercross

The leptin receptor gene (LEPR) is a candidate for traits related to growth and body composition, and is located on SSC6 in a region where fatness and meat composition quantitative trait loci (QTL) have previously been detected in several F2 experimental designs. The aims of this work were: (i) to fine map these QTL on a larger sample of animals and generations (F3 and backcross) of an Iberian x Landrace intercross and (ii) to examine the effects of LEPR alleles on body composition traits. Eleven single nucleotide polymorphisms (SNPs) were detected by sequencing LEPR coding regions in Iberian and Landrace pig samples. Three missense polymorphisms were genotyped by pyrosequencing in 33 F0, 70 F1, 418 F2, 86 F3 and 128 individuals coming from the backcross of four F2 males with 24 Landrace females. Thirteen microsatellites and one SNP were also genotyped. Traits analysed were: backfat thickness at different locations (BF(T)), intramuscular fat percentage (IMF(P)), eye muscle area (EM(A)), loin depth (LO(D)), weight of shoulder (SH(W)), weight of ribs (RIB(W)) and weight of belly bacon (BB(W)). Different statistical models were applied in order to evaluate the number and effects of QTL on chromosome 6 and the possible causality of the LEPR gene variants with respect to the QTL. The results support the presence of two QTL on SSC6. One, at position 60-100 cM, affects BF(T) and RIB(W). The other and more significant maps in a narrow region (130-132 cM) and affects BF(T), IMF(P), EM(A), LO(D), SH(W), RIB(W) and BB(W). Results also support the association between LEPR alleles and BF(T) traits. The possible functional implications of the analysed polymorphisms are considered.     

Association of the polymorphism in GYS1 and ACOX1 genes with meat quality traits in pigs

Phenotypic information about several pig meat quality traits on 334 Large White × Meishan F2 pigs was collected. Effects of the association of the FokI variants in the seventh intron of the skeletal muscle glycogen synthase (GYS1) gene and the PstI variants in the ninth intron of the palmitoyl acyl-CoA oxidase 1 (ACOX1) gene on the meat quality traits were examined on all pigs. The FokI variants of the GYS1 gene showed significant effects on pH of m. semipinalis capitis (P < 0.05). Linkage analysis indicated that the peak of the quantitative trait loci (QTL) curve was located around this marker for pH, but it did not reach significance (P > 0.05). The results may be due to several reasons such as linkage disequilibrium to the causal mutations, the limited number of animals or balance of another QTL or marker with negative effects. Significant effects of PstI variants of ACOX1 gene were also found on meat colour value and meat marbling score of both m. longissimus dorsi and m. biceps femoris (P < 0.05). Dominant effects for the affected traits at those two loci were significant except for meat marbling score of m. biceps femoris (P < 0.05). The results of this study give us some evidence for the potential of those dominant markers used in the marker-assisted selection of crossbreeding of the Large White pig sire lines and Meishan-derived synthetic dam lines.     

An improved method for quantitative trait loci detection and identification of within-line segregation in F2 intercross designs

We present a new flexible, simple, and powerful genome-scan method (flexible intercross analysis, FIA) for detecting quantitative trait loci (QTL) in experimental line crosses. The method is based on a pure random-effects model that simultaneously models between- and within-line QTL variation for single as well as epistatic QTL. It utilizes the score statistic and thereby facilitates computationally efficient significance testing based on empirical significance thresholds obtained by means of permutations. The properties of the method are explored using simulations and analyses of experimental data. The simulations showed that the power of FIA was as good as, or better than, Haley-Knott regression and that FIA was rather insensitive to the level of allelic fixation in the founders, especially for pedigrees with few founders. A chromosome scan was conducted for a meat quality trait in an F(2) intercross in pigs where a mutation in the halothane (Ryanodine receptor, RYR1) gene with a large effect on meat quality was known to segregate in one founder line. FIA obtained significant support for the halothane-associated QTL and identified the base generation allele with the mutated allele. A genome scan was also performed in a previously analyzed chicken F(2) intercross. In the chicken intercross analysis, four previously detected QTL were confirmed at a 5% genomewide significance level, and FIA gave strong evidence (P < 0.01) for two of these QTL to be segregating within the founder lines. FIA was also extended to account for epistasis and using simulations we show that the method provides good estimates of epistatic QTL variance even for segregating QTL. Extensions of FIA and its applications on other intercross populations including backcrosses, advanced intercross lines, and heterogeneous stocks are also discussed.     

QTL analysis of white blood cell, platelet and red blood cell-related traits in an F2 intercross between Landrace and Korean native pigs

Haematological traits play important roles in disease resistance and defence functions. The objective of this study was to locate quantitative trait loci (QTL) and the associated positional candidate genes influencing haematological traits in an F₂ intercross between Landrace and Korean native pigs. Eight blood‐related traits (six erythrocyte traits, one leucocyte trait and one platelet trait) were measured in 816 F₂ progeny. All experimental animals were genotyped with 173 informative microsatellite markers located throughout the pig genome. We report that nine chromosomes harboured QTL for the baseline blood parameters: genomic regions on SSC 1, 4, 5, 6, 8, 9, 11, 13 and 17. Eight of twenty identified QTL reached genome‐wide significance. In addition, we evaluated the KIT locus, an obvious candidate gene locus affecting variation in blood‐related traits. Using dense single nucleotide polymorphism marker data on SSC 8 and the marker‐assisted association test, the strong association of the KIT locus with blood phenotypes was confirmed. In conclusion, our study identified both previously reported and novel QTL affecting baseline haematological parameters in pigs. Additionally, the positional candidate genes identified here could play an important role in elucidating the genetic architecture of haematological phenotype variation in swine and in humans.     

Comparative sequence analysis of the PRKAG3 region between human and pig: evolution of repetitive sequences and potential new exons

The PRKAG3 gene encodes the gamma3 chain of AMP-activated protein kinase (AMPK). A non-conservative missense mutation in the PRKAG3 gene causes a dominant phenotype involving abnormally high glycogen content in pig skeletal muscle. We have determined >126 kb (in 13 contigs) of porcine genomic sequence surrounding the PRKAG3 gene and the corresponding mouse region covering the gene. A comparison of these PRKAG3 sequences and the human sequence was conducted and used to predict evolutionarily conserved regions, including regulatory regions. A comparison of the human genomic sequence and a porcine BAC sequence containing the PRKAG3 gene, revealed a conserved organization and the presence of three additional genes, CYP27A1 (cytochrome P450, family 27, subfamily A, polypeptide 1), STK36 (Serine Threonine Kinase 36), and the homolog of the unidentified human mRNA KIAA0173. Interspersed repetitive elements constituted 51.4 and 38.6% of this genomic region in human and pig, respectively. We were able to reliably align 12.6 kb of orthologous repeats shared between pig and human and these showed an average sequence identity of 72.4%. Our analysis revealed that the human KIAA0173 gene harbors alternative 5' untranslated exons originating from repetitive elements. This provides an obvious example how transposable elements may affect gene evolution.     

Detection of quantitative trait loci for meat quality in a commercial slaughter pig cross

The primary goal of this study was to localize quantitative trait loci (QTL) affecting meat quality traits in swine. In total, 42 traits were scored on 305 F₂ individuals from a commercial slaughter pig cross in Norway. F₁ and F₂ individuals were genotyped for 29 markers on Chromosomes (Chrs) 4, 6, and 7, since previous studies had revealed QTL affecting meat quality traits on these chromosomes. The most evident result was detection of a QTL affecting amount of intramuscular fat on Chr 6. The QTL might also influence tenderness, whereas no effect was observed for back-fat thickness. Additionally, suggestive evidence for QTL affecting other meat quality traits was found on Chr 4 and Chr 7. Received: 16 June 2000 / Accepted: 5 December 2000