- Research article
Genome-wide identification of quantitative trait loci in a cross between Hampshire and Landrace II: Meat quality traits
BMC Geneticsvolume 9, Article number: 22 (2008)
Meat quality traits are important in pig breeding programs, but they are difficult to include in a traditional selection program. Marker assisted selection (MAS) of meat quality traits is therefore of interest in breeding programs and a Quantitative Trait Locus (QTL) analysis is the key to identifying markers that can be used in MAS. In this study, Landrace and Hampshire intercross and backcross families were used to investigate meat quality traits. Hampshire pigs are commonly used as the sire line in commercial pig breeding. This is the first time a pedigree including Hampshire pigs has been used for a QTL analysis of meat quality traits.
In total, we analyzed 39 meat quality traits and identified eight genome-wide significant QTL peaks in four regions: one on chromosome 3, two on chromosome 6 and one on chromosome 16. At least two of the QTLs do not appear to have been detected in previous studies. On chromosome 6 we identified QTLs for water content in M. longissimus dorsi (LD), drip loss in LD and post mortem pH decline in LD. On chromosomes 3 and 16 we identified previously undetected QTLs for protein content in LD and for freezing and cooking loss respectively.
We identified at least two new meat quality trait QTLs at the genome-wide significance level. We detected two QTLs on chromosome 6 that possibly coincide with QTLs detected in other studies. We were also able to exclude the C1843T mutation in the ryanodine receptor (RYR1) as a causative mutation for one of the chromosome 6 QTLs in this cross.
Since the first Quantitative Trait Locus (QTL) analysis in pigs was published in 1994 , QTL analyses have been widely used to identify chromosomal regions harbouring genes for various complex traits in the pig such as growth, carcass composition and meat quality . Meat quality traits have been studied before using crosses between Wild Boar and Large White [3, 4], Meishan and Yorkshire , Meishan and Large White/Landrace , Duroc and Landrace/Yorkshire , Berkshire and Yorkshire , Iberian and Landrace , Pietrain and Meishan and Wild Boar , and between Duroc and Berlin Miniature pig .
In this study we used a cross between Finnish Landrace and Swedish Hampshire set up by the Swedish breeding company Quality Genetics, as a combined intercross/backcross design. Landrace has been used in several QTL crosses before, but so far the Hampshire breed has not been used in any QTL intercross which provided an opportunity to detect specific QTL alleles that have been selected in this breed. Landrace and Hampshire pigs differ in a number of traits including coat colour, body composition, fertility and meat quality. Landrace has a long body compared to the shorter more compact Hampshire pig and Hampshire is more muscular than Landrace .
A mutation in PRKAG3, the RN-mutation (RN-), has a large impact on the technological yield and meat quality and has been widespread amongst Hampshire pigs . Its high frequency was most likely the result of its ability to increase the lean meat content of pigs. The effect of the RN-mutation in this cross on traits such as technological yield, meat quality and colour characteristics of pork has been published elsewhere [14–16].
A genome scan detecting QTLs for carcass traits in this cross was published previously , and in this study we report the results for meat quality traits. We identified four QTL regions on three different chromosomes that reached genome-wide significance. At least two of these have not been detected in previous studies.
The total length of the linkage map including all autosomes was estimated to be 23.8 Morgans (M). The average distance between markers was 23.1 cM, with five telomeric regions on SSC 3, 5, 7, 8 and 16 exceeding 50 cM between markers. The linkage map is presented in Table 1.
The genome-wide significance thresholds were F = 8.3 and F = 10.2 at the 5% and 1% level respectively, for traits analysed using a combined F2 and backcross analysis. The sensory traits were analysed using only F2 animals and for these the significance thresholds were F = 9.9 and F = 12.9 at the 5% and 1% genome-wide significance level, respectively.
Three PRKAG3 alleles, RN- (R225Q), rn+ (wild type) and rn* (V224I), were segregating in this family material. Allele frequencies in the F2 and backcross generations are presented in Table 2. As an internal positive control, we performed a QTL analysis of muscle glycogen content without including the PRKAG3-genotype as a fixed effect. This gave a significant QTL effect with an F-value for muscle glycogen content of over 100, which completely disappears when the PRKAG3-genotype is included as a fixed effect (data not shown). This confirms the high quality of the phenotypic data and an excellent matching of genotype and phenotype data.
We performed genome scans for 39 meat quality traits (Table 3) and observed eight genome-wide significant QTL tests (Table 4; Figure 1). This is more than the expected number of Type I errors, given the fact that we have carried out genome scans for 39 traits and used the 5% significance level. All QTL tests that reached chromosome-wise significance are compiled in Table 5. We observed 46 suggestive QTLs which are only slightly more than the 39 expected Type I errors (~1 per genome scan/trait). Thus, a large proportion of the suggestive QTLs is expected to be false positives and further studies are needed to sort out which ones are true positives. The suggestive QTLs are therefore not further discussed here except those that were co-localized with QTLs showing genome-wide significance. No QTL showing genetic imprinting was detected in this study (data not shown).
On SSC3, a QTL for protein content in M. Longissimus dorsi (LD) was detected with a peak at 156 cM. The QTL showed an additive effect and the Hampshire allele was associated with a higher protein content. In the same region of chromosome 3 we detected QTLs for glycogen content in LD, water content in LD and colour a* (redness), all of which reached chromosome-wise significance (Table 5). The Hampshire allele was associated with reduced glycogen and water content and higher degree of redness (colour a*).
QTLs affecting water content in LD, drip loss in LD during four days and pH decline in LD between 45 min and 3 hours post mortem were detected between positions 51 and 69 cM on SSC6 (Table 4; Figure 1). It is likely that these significant effects reflect the action of a single QTL. The QTL showed an additive effect and the Hampshire allele was associated with higher water content, drip loss and pH decline after slaughter. In the same interval, QTL tests with chromosome-wise significance for freezing and cooking loss, drip loss during 24 hours and intramuscular fat in LD were obtained (Table 5). We excluded the previously published porcine C1843T mutation in the ryanodine receptor gene (RYR1)  as a causative mutation for this QTL since it did not segregate in the pedigree discussed herein.
Another QTL on SSC6, with its peak at position 119 cM, was identified for three highly correlated traits, chewing resistance, tenderness and chewing time, scored by a trained sensory panel. The significance and estimated effects of this QTL must be interpreted with caution since only 53 animals were scored for these traits. This QTL showed overdominance, which means that the heterozygous class had the most extreme phenotypic value, and was associated with higher tenderness (Table 4). A suggestive QTL (1% chromosome-wise significance) for the total impression of the meat was found in the same region of chromosome 6 (Table 5).
We identified a QTL for freezing and cooking loss at 41 cM on SSC16 showing overdominance; the heterozygotes showed reduced freezing and cooking losses (Table 4). In the same region, suggestive QTLs for drip loss during 4 days and 24 hours were also identified (Table 5).
Meat quality is obviously of great importance in commercial pig breeding and it is a trait that is difficult and expensive to measure accurately on a large number of pigs in a progeny testing scheme. It is therefore of considerable interest to identify QTLs in experimental populations and exploit such loci by marker assisted selection (MAS) in breeding programs. Furthermore, the molecular characterization of genes controlling meat quality and meat content can provide new insights into muscle metabolism. This is illustrated by the identification of missense mutations in RYR1  and PRKAG3  that have major effects on lean meat content and meat quality, as well as by the point mutation in intron 3 of IGF2  underlying a major QTL for muscle growth and lean meat content. Thus, further genetic studies of the QTLs reported here may lead to new basic knowledge as well as practical applications.
We identified two QTL regions on SSC6. The first QTL region, located at position 51–69 cM, affects water content in LD, drip loss in LD over four days and pH decline in LD between 45 min and 3 hours post mortem. Several other studies have also identified QTLs for meat quality traits in this region. QTLs for meat quality, stress resistance and carcass composition were mapped to SSC6 in crosses including the Piétrain breed. These QTLs are most likely explained by a mutation in RYR1 occurring at a high frequency in the Piétrain breed [10, 20]. The pigs in our cross do not carry this mutation. Another study using non-carriers of the RYR1 mutation has also identified QTLs for meat quality traits in this region. Malek et al. identified a suggestive QTL for pH 24 hours post mortem in loin using a Berkshire x Yorkshire cross . The location of this QTL is in the same region as our QTLs and they both showed an additive effect. However, Malek et al. did not detect QTLs for drip loss and cooking loss in this region even though these traits were included in their study.
The second QTL region on SSC6, with a peak at 119 cM, influenced chewing resistance, chewing time and tenderness. These traits are highly correlated and we assume that it is a single QTL that influences all three traits. A panel of individuals subjectively scored these traits and only 53 pigs were included. The small sample size reduces the power to detect QTLs for these traits and the results should be interpreted with caution. The QTL showed overdominance and was estimated to explain an astonishing ~35% of the residual variance, which could be an overestimation due to the few number of pigs analyzed. To put these results in perspective, we performed a QTL analysis for muscle glycogen content on chromosome 15 using the same 53 pigs to test if we could detect the segregation at the RN locus with this small number of pigs (PRKAG3-genotype was not included as a fixed effect and PRKAG3 was excluded as a marker in the linkage map). We obtained a statistically significant F-value of 11.0 at approximately the correct position (data not shown). This demonstrates that we can detect loci with major phenotypic effects using only 53 animals. Interestingly, Szyda et al. have reported a QTL for tenderness with an overlapping location to our QTL using a Norwegian commercial slaughter pig cross including Duroc, Norweigan Landrace and Yorkshire . Further studies are required to find out whether our observation reflects a Type I error or a new major locus with an important effect on meat quality.
A QTL for protein content in LD was found at 156 cM on SSC3 and, to the best of our knowledge, no QTL with similar effect has previously been reported in this region. Similarly, we are not aware of any previously reported QTL with a strikingly similar effect to the one for freezing and cooking loss in LD that we mapped to position 41 cM on SSC16. Pierzchala et al. identified QTLs for conductivity, pH measurements and stress response on SSC16 in crosses between Meishan, Wild Boar and Piétrain but they did not detect a QTL for cooking loss even though this trait was scored . Paszek et al. also detected a QTL for pH in muscle on SSC16 but did not see any QTL for muscle moisture in that region . Our QTL had no significant effect on pH values.
In this study 39 meat quality traits were analyzed and we identified eight QTLs at the genome-wide significance level. The QTLs were located in four regions, one on chromosome 3, two on chromosome 6 and one on chromosome 16. This was the first time the Hampshire breed was used in a QTL study of meat quality traits and it enabled us to detect two previously undetected QTLs on chromosome 3 and 16. We also identified two QTLs on chromosome 6 that coincide with QTLs detected in previous studies. One of the chromosome 6 QTLs is located in the same region as QTLs explained by the C1843T mutation in the ryanodine receptor (RYR1), however we have been able to exclude this as a causative mutation for our QTL. Several interesting QTL regions have been identified in this study and, although they require further investigation, they may be interesting for Marker Assisted Selection (MAS) in the future.
Animals and genotyping
The breeding company Quality Genetics established a three generation-cross between Finnish Landrace and Swedish Hampshire for commercial reasons. A combined intercross and backcross design was used. Eight Landrace boars (L) were crossed with 41 Hampshire sows (H) generating 52 animals in the F1-generation (LH). F1-animals were then intercrossed to produce 136 F2-animals. LH animals from the F1-generation were also reciprocally backcrossed to 42 purebred Hampshire pigs producing 112 (LH × H) and 72 (H × LH) offspring. Including the parental generation of the purebred Hampshire pigs the pedigree comprised a total of 527 animals. The offspring represented 86 full-sib families.
Husbandry and slaughtering as well as the phenotypic measurements have previously been described in detail [14–16, 23]. The pigs were raised at three different breeding farms referred to as herd, but they were slaughtered at the same commercial slaughterhouse. During the experiment, the stunning procedure at the slaughterhouse changed, from individual stunning with CO2 to stunning in groups of five pigs. The traits analyzed in the current study are listed in Table 1.
A total of 120 microsatellite markers covering the autosomes were PCR amplified in 450 animals (excluding the 77 purebred Hampshire sows in the parental generation) and genotyped using either an ABI PRISM® 3100 Genetic Analyzer and ABI GeneMapper™ Genotyping Software in Copenhagen or a MegaBACE™ 1000 DNA Analysis System and Genetic Profiler (Amersham Biosciences) in Uppsala.
The three alleles, denoted rn+ (wild type), RN- (R225Q) and rn* (V224I), at the PRKAG3/RN- locus were scored according to a previously described method using pyrosequencing . The single point mutation (C → T) in the pig ryanodine receptor (ryr1) gene changing an arginine to a cysteine at amino acid 615  was genotyped with pyrosequencing using the following primers: forward primer with an M13-tag sequence CACGACGTTGTAAAACGACAGTGCCCTCACACCTTGAC, reverse primer CCAGGGAGCAAGTTCTCAGT, M13-biotinylated primer CACGACGTTGTAAAACGAC and sequencing primer AGTAATGAGATCTTGGTTGGAG. A 20 μl PCR reaction with 1× PCR Buffer II, 2.5 mM MgCl2, 0.3 mM dNTP, 0.03 μM forward primer, 0.3 μM of each reverse and M13-biotinylated primer, 0.75 U of AmpliTaq GOLD polymerase (Applied Biosystems) and 50–100 ng DNA was run using a standard touch-down PCR protocol. Starting with 95°C for 15 min, then 14 touch down cycles 95°C 30 s, 65–52°C 30 s, 72°C 30 s, followed by 30 cycles 95°C 30 s, 52°C 30 s, 72°C 30 s and ending with 72°C for 10 min. A standard pyrosequencing protocol was employed (Biotage).
Linkage maps were built using the CRI-MAP program version 2.1 or 2.4 . The sex average maps were used in the QTL analyses. Phenotypes were checked for normal distribution using the Ryan-Joiner normality test in MiniTab and transformed when necessary. The QTL analyses were performed using QTL express  and the combined F2 and backcross option as described in detail in our previous study . The model including additive and dominance effects was compared with a model also including a parent-of-origin effect for all traits. For meat quality traits, except male hormones, the fixed effects herd, sex, stunning procedure and PRKAG3/RN-genotype were used; six different PRKAG3/RN genotypes were observed. For sex hormones, only herd, stunning procedure and PRKAG3/RN-genotype were included as fixed effects. The QTL analysis of sensory traits only included F2 progeny and these were therefore analyzed using the F2 design in QTL express. The model included herd, stunning procedure and PRKAG3/RN-genotype. Genome-wide significant thresholds were determined by permutation tests . One thousand permutations were performed for all traits and an average calculated. Two different thresholds were permutated for traits analysed using the different options in QTL express. Chromosome-wise significant thresholds were also determined by permutation tests using thousand permutations. Confidence intervals (CI) were estimated for each genome-wide significant QTL using the bootstrap method  and 10,000 iterations were performed. For the genome-wide significant QTLs, the residual variance explained by the QTL was computed as ((Residual sums of squares reduced model – Residual sums of squares full model)/Residual sums of squares reduced model) ×100.
Andersson L, Haley CS, Ellegren H, Knott SA, Johansson M, Andersson K, Andersson-Eklund L, Edfors-Lilja I, Fredholm M, Hansson I, et al: Genetic mapping of quantitative trait loci for growth and fatness in pigs. Science. 1994, 263 (5154): 1771-1774. 10.1126/science.8134840.
Hu ZL, Dracheva S, Jang W, Maglott D, Bastiaansen J, Rothschild MF, Reecy JM: A QTL resource and comparison tool for pigs: PigQTLDB. Mamm Genome. 2005, 16 (10): 792-800. 10.1007/s00335-005-0060-9.
Andersson-Eklund L, Marklund L, Lundstrom K, Haley CS, Andersson K, Hansson I, Moller M, Andersson L: Mapping quantitative trait loci for carcass and meat quality traits in a wild boar x Large White intercross. J Anim Sci. 1998, 76 (3): 694-700.
Nii M, Hayashi T, Mikawa S, Tani F, Niki A, Mori N, Uchida Y, Fujishima-Kanaya N, Komatsu M, Awata T: Quantitative trait loci mapping for meat quality and muscle fiber traits in a Japanese wild boar x Large White intercross. J Anim Sci. 2005, 83 (2): 308-315.
Paszek AA, Wilkie PJ, Flickinger GH, Miller LM, Louis CF, Rohrer GA, Alexander LJ, Beattie CW, Schook LB: Interval mapping of carcass and meat quality traits in a divergent swine cross. Anim Biotechnol. 2001, 12 (2): 155-165. 10.1081/ABIO-100108342.
de Koning DJ, Harlizius B, Rattink AP, Groenen MA, Brascamp EW, van Arendonk JA: Detection and characterization of quantitative trait loci for meat quality traits in pigs. J Anim Sci. 2001, 79 (11): 2812-2819.
Grindflek E, Szyda J, Liu Z, Lien S: Detection of quantitative trait loci for meat quality in a commercial slaughter pig cross. Mamm Genome. 2001, 12 (4): 299-304. 10.1007/s003350010278.
Malek M, Dekkers JC, Lee HK, Baas TJ, Prusa K, Huff-Lonergan E, Rothschild MF: A molecular genome scan analysis to identify chromosomal regions influencing economic traits in the pig. II. Meat and muscle composition. Mamm Genome. 2001, 12 (8): 637-645. 10.1007/s003350020019.
Ovilo C, Clop A, Noguera JL, Oliver MA, Barragan C, Rodriguez C, Silio L, Toro MA, Coll A, Folch JM, Sanchez A, Babot D, Varona L, Perez-Enciso M: Quantitative trait locus mapping for meat quality traits in an Iberian x Landrace F2 pig population. J Anim Sci. 2002, 80 (11): 2801-2808.
Geldermann H, Müller E, Moser G, Reiner G, Bartenschlager H, Cepica S, Stratil A, Kuryl J, Moran C, Davoli R, Brunsch C: Genome-wide linkage and QTL mapping in porcine F2 families generated from Pietrain, Meishan and Wild Boar crosses. J Anim Breed Genet. 2003, 120: 363-393. 10.1046/j.0931-2668.2003.00408.x.
Wimmers K, Fiedler I, Hardge T, Murani E, Schellander K, Ponsuksili S: QTL for microstructural and biophysical muscle properties and body composition in pigs. BMC Genet. 2006, 7: 15-10.1186/1471-2156-7-15.
Whittemore C: The Science and Practice of Pig Production. 1998, Oxford , Blackwell Science Ltd, 2
Milan D, Jeon JT, Looft C, Amarger V, Robic A, Thelander M, Rogel-Gaillard C, Paul S, Iannuccelli N, Rask L, Ronne H, Lundstrom K, Reinsch N, Gellin J, Kalm E, Roy PL, Chardon P, Andersson L: A mutation in PRKAG3 associated with excess glycogen content in pig skeletal muscle. Science. 2000, 288 (5469): 1248-1251. 10.1126/science.288.5469.1248.
Lindahl G, Enfalt AC, von Seth G, Josell A, Hedebro-Velander I, Andersen HJ, Braunschweig M, Andersson L, Lundstrom K: A second mutant allele (V199I) at the PRKAG3 (RN) locus--I. Effect on technological meat quality of pork loin. Meat Science. 2004, 66 (3): 609-619. 10.1016/S0309-1740(03)00179-7.
Lindahl G, Enfalt AC, von Seth G, Josell A, Hedebro-Velander I, Andersen HJ, Braunschweig M, Andersson L, Lundstrom K: A second mutant allele (V199I) at the PRKAG3 (RN) locus--II. Effect on colour characteristics of pork loin. Meat Science. 2004, 66 (3): 621-627. 10.1016/S0309-1740(03)00180-3.
Zamaratskaia G, Madej A, Babol J, Squires E, Lundstrom K: Free oestrone in adipose tissue and its relation to androstenone and skatole in entire male pigs. Reprod Domest Anim. 2005, 40 (2): 156-160. 10.1111/j.1439-0531.2005.00573.x.
Karlskov-Mortensen P, Bruun CS, Braunschweig MH, Sawera M, Markljung E, Enfalt AC, Hedebro-Velander I, Josell A, Lindahl G, Lundstrom K, von Seth G, Jorgensen CB, Andersson L, Fredholm M: Genome-wide identification of quantitative trait loci in a cross between Hampshire and Landrace I: carcass traits. Anim Genet. 2006, 37 (2): 156-162. 10.1111/j.1365-2052.2005.01405.x.
Fujii J, Otsu K, Zorzato F, de Leon S, Khanna VK, Weiler JE, O'Brien PJ, MacLennan DH: Identification of a mutation in porcine ryanodine receptor associated with malignant hyperthermia. Science. 1991, 253 (5018): 448-451. 10.1126/science.1862346.
Van Laere AS, Nguyen M, Braunschweig M, Nezer C, Collette C, Moreau L, Archibald AL, Haley CS, Buys N, Tally M, Andersson G, Georges M, Andersson L: A regulatory mutation in IGF2 causes a major QTL effect on muscle growth in the pig. Nature. 2003, 425 (6960): 832-836. 10.1038/nature02064.
Yue G, Stratil A, Kopecny M, Schröffelova D, Schröffel Jr J, Hojny J, Cepica S, Davoli R, Zambonelli P, Brunsch C, Sternstein I, Moser G, Bartenschlager H, Reiner G, Geldermann H: Linkage and QTL mapping for Sus scrofa chromosome 6. J Anim Breed Genet. 2003, 120 (Suppl. 1): 45-55. 10.1046/j.0931-2668.2003.00423.x.
Szyda J, Grindflek E, Liu Z, Lien S: Multivariate mixed inheritance models for QTL detection on porcine chromosome 6. Genet Res. 2003, 81 (1): 65-73. 10.1017/S0016672302006043.
Pierzchala M, Kuryl J, Reiner G, Bartenschlager H, Moser G, Geldermann H: Linkage and QTL mapping for Sus scrofa chromosome 16. J Anim Breed Genet. 2003, 120 (Suppl. 1): 126-131. 10.1046/j.0931-2668.2003.00433.x.
Josell A, von Seth G, Tornberg E: Sensory and meat quality traits of pork in relation to post-slaughter treatment and RN genotype. Meat Science. 2004, 66 (1): 113-124. 10.1016/S0309-1740(03)00067-6.
Green P, Falls K, Crook S: Documentation for CRIMAP, ver. 2.4. Washington University School of Medicine, St Louis. 1990
Seaton G, Haley CS, Knott SA, Kearsey M, Visscher PM: QTL Express: mapping quantitative trait loci in simple and complex pedigrees. Bioinformatics. 2002, 18 (2): 339-340. 10.1093/bioinformatics/18.2.339.
Churchill GA, Doerge RW: Empirical threshold values for quantitative trait mapping. Genetics. 1994, 138 (3): 963-971.
Visscher PM, Thompson R, Haley CS: Confidence intervals in QTL mapping by bootstrapping. Genetics. 1996, 143 (2): 1013-1020.
The project was funded by the Nordic Joint Committee for Agricultural Research, the Swedish Research Council for Environment, Agricultural Sciences and Spatial Planning, the Danish Agricultural and Veterinary Research Council and PhD stipends from the Royal Veterinary and Agricultural University, Denmark to Peter Karlskov-Mortensen, Camilla S. Bruun and Milena Sawera. Animal material was provided by the Swedish pig breeding company Quality Genetics. We thank Ulla Gustafsson, Majken Madvig Jansen, and Tina Neergård Mahler for excellent technical assistance and Greger Larsen for valuable comments on the manuscript.
EM compiled the Uppsala genotyping data, contributed to the QTL analyses, carried out the genotyping of the RYR1 mutation, summarized the data and drafted the manuscript. MHB carried out the DNA extractions and the genotyping in Uppsala. PKM contributed to the genotyping in Copenhagen and to the QTL analyses. CSB contributed to the genotyping in Copenhagen and to the QTL analyses. MS contributed to the genotyping in Copenhagen and to the QTL analyses. ICC contributed to retyping of some markers in Uppsala. IHV was responsible for the animal material and collection of samples for DNA extraction. ÅJ contributed to collection of the phenotypic data. KL contributed to the design of the study and to the collection of the phenotypic data. GvS contributed to collection of the phenotypic data. CBJ contributed to the genotyping in Copenhagen and to the QTL analyses. MF contributed to the design of the study and supervised this study in Copenhagen. LA designed the study and supervised the molecular analysis in Uppsala, edited and made final improvements to the manuscript. All authors read and approved the final manuscript.
Ellen Markljung, Martin H Braunschweig contributed equally to this work.