Original Research ARTICLE
Identification of carcass and meat quality QTL in an F2 Duroc × Pietrain pig resource population using different least-squares analysis models
- Department of Animal Science, Michigan State University, East Lansing, MI, USA
A three-generation resource population was constructed by crossing pigs from the Duroc and Pietrain breeds. In this study, 954 F2 animals were used to identify quantitative trait loci (QTL) affecting carcass and meat quality traits. Based on results of the first scan analyzed with a line-cross (LC) model using 124 microsatellite markers and 510 F2 animals, 9 chromosomes were selected for genotyping of additional markers. Twenty additional markers were genotyped for 954 F2 animals and 20 markers used in the first scan were genotyped for 444 additional F2 animals. Three different Mendelian models using least-squares for QTL analysis were applied for the second scan: a LC model, a half-sib (HS) model, and a combined LC and HS model. Significance thresholds were determined by false discovery rate (FDR). In total, 50 QTL using the LC model, 38 QTL using the HS model, and 3 additional QTL using the combined LC and HS model were identified (q < 0.05). The LC and HS models revealed strong evidence for QTL regions on SSC6 for carcass traits (e.g., 10th-rib backfat; q < 0.0001) and on SSC15 for meat quality traits (e.g., tenderness, color, pH; q < 0.01), respectively. QTL for pH (SSC3), dressing percent (SSC7), marbling score and moisture percent (SSC12), CIE a* (SSC16), and carcass length and spareribs weight (SSC18) were also significant (q < 0.01). Additional marker and animal genotypes increased the statistical power for QTL detection, and applying different analysis models allowed confirmation of QTL and detection of new QTL.
Quantitative trait loci (QTL) mapping has been conducted using numerous pig populations to identify genomic regions controlling phenotypic variation for hundreds of traits (http://www.animalgenome.org/cgi-bin/QTLdb/SS/index). Nevertheless, the implementation of QTL into breeding programs which is a major goal of QTL mapping has been limited not only due to insufficient numbers of identified causative mutations, but because of unknown linkage disequilibrium (LD) phase between markers and QTL resulting from cross breeding systems (Spelman and Bovenhuis, 1998; Hayes et al., 2009). We have developed a F2 Duroc × Pietrain resource population at Michigan State University (Edwards et al., 2008b) and reported QTL for carcass merit and meat quality traits (Edwards et al., 2008a). The Duroc and Pietrain breeds are used in breeding programs as sire breeds worldwide, and these breeds exhibit variation in carcass merit and meat quality phenotypes. Pietrain pigs have been shown to have less backfat (Affentranger et al., 1996; Edwards et al., 2003) and larger longissimus muscle area (LMA; Edwards et al., 2003). Duroc and Duroc-sired pigs generally have more favorable meat quality (Langlois and Minvielle, 1989; Affentranger et al., 1996; Jeremiah et al., 1999; Edwards et al., 2003), whereas Pietrain and Pietrain-sired pigs are leaner with average meat quality (Edwards et al., 2003).
A line-cross (LC) model, which assumes the founder lines to be fixed for alternative QTL alleles, has been most commonly used to identify QTL for F2 population designs (Haley et al., 1994). The first genome scan for our Duroc × Pietrain population was performed using a LC analysis (Edwards et al., 2008a, b). However, for crosses between outbred lines such as domestic animals, not all QTL alleles are completely fixed so effects under the LC model can be biased downwards (Pérez-Enciso and Varona, 2000). To identify QTL segregating within parental breeds, a half-sib (HS) model that does not assume fixation of QTL alleles in the founder lines was introduced by Knott et al. (1996), and Kim et al. (2005) developed a combined line-cross and half-sib (CB) model that accounts for both line and HS effects. We have recently utilized LC, HS, and CB models to identify QTL for growth traits in our Duroc × Pietrain population (Choi et al., 2010). The objective of this study was to confirm previously identified carcass merit and meat quality QTL regions with addition of new marker genotypes and additional F2 animals, and to detect new QTL for carcass merit and meat quality traits using three different least-squares models under different assumptions; (1) founders fixed for alternative QTL alleles (LC model), (2) segregation of QTL alleles at similar frequencies in founders (HS model), and (3) segregation of QTL alleles at different frequencies in founders (CB model).
Materials and Methods
Animals and Phenotypic Data
A three-generation resource population developed at Michigan State University was used for this study. A detailed description of the animals and phenotypic data was previously reported (Edwards et al., 2008a, b). All grandparents were confirmed to be homozygous normal for the polymorphism at position 1843 in the RYR1 gene (Edwards et al., 2008b). Animal protocols were approved by the Michigan State University All University Committee on Animal Use and Care (AUF# 09/03-114-00). A total of 954 F2 pigs were used which included the 510 animals evaluated in the first genome scan. These pigs were produced from 6 F1 boars and 50 F1 sows which were retained from 4 F0 Duroc sires and 15 F0 Pietrain dams. The F2 pigs were analyzed for 38 carcass and meat quality traits.
Details of carcass and meat quality phenotype collection were reported in Edwards et al., 2008a. Briefly, animals were slaughtered at the Michigan State University Meat Laboratory (East Lansing, MI, USA) or a federally inspected commercial plant (DeVries Meats, Coopersville, MI, USA). Slaughter age was 165.8 ± 9.2 days and the minimum off-farm body weight (BW) for slaughter was 82.54 kg. Hot carcass weight (HCW), and pH and temperature of the longissimus muscle (LM) at 45-min and 24-h postmortem were obtained. After overnight chilling, backfat thickness, number of ribs and carcass length were measured, and the weights of primal cuts were recorded. A single trained evaluator scored color, marbling, and firmness using two 2.54-cm thick chops cut from the LM, and objective color scores of CIE L*, a*, and b* were obtained using a Minolta colorimeter. The remaining section of the LM was used to determine drip loss, cook yield, Warner-Bratzler shear force, proximate analysis measures, and sensory attributes. A trained sensory panel evaluated juiciness, tenderness, overall tenderness, connective tissue, and off-flavor using an 8-point hedonic scale. Descriptive statistics for phenotypes used in this study are presented in Table 1.
Nine chromosomes (SSC3–7, 12, 15, 16, and 18) were selected based on results of the first genome scan (Edwards et al., 2008a, b) which had been completed using 510 F2 animals and 124 microsatellite markers. For the second scan 20 additional microsatellite markers were selected on these chromosomes (1–4 markers per chromosome; Choi et al., 2010) in order to increase the power of QTL detection and to narrow the QTL locations. All F0, F1, and the 954 F2 pigs were genotyped for the 20 new markers, and the 444 additional F2 pigs were also genotyped for 20 markers flanking the QTL regions on the 9 selected chromosomes. Sex-averaged genetic linkage maps were estimated for all autosomes using CRI-MAP version 2.4 (Green et al., 1990) and converted to the Haldane map function (Choi et al., 2010).
Three different models using least-squares (LC, HS, and CB models) were adopted for QTL analysis (Kim et al., 2005) and analyses were performed using the methods described in Choi et al. (2010). Significance thresholds were determined by false discovery rate (FDR; Weller et al., 1998).
The LC analysis assumes the QTL to be fixed for alternative alleles in the founder lines. Probabilities of each F2 individual being homozygous for two Duroc alleles (P11), homozygous for two Pietrain alleles (P22), or heterozygous (P12 or P21) were estimated at fixed 1-cM intervals across the genome using the QTL Express software (Seaton et al., 2002). By denoting the mean of homozygous animals for the Duroc allele as positive additive (a), the mean of heterozygous animals as dominance (d) and the mean of homozygous animals for the Pietrain allele as negative additive (−a), the following linear model was fitted at every cM across the genome.
yj = Xjb + aPaj + dPdj + ej
Where yj is the phenotype of F2 progeny j, Xj, and b are the design matrix and solution vector for the fixed effects, respectively, a and d are the estimated additive and dominance effects of a putative QTL at the given location, respectively, Paj = P11 − P22 is the conditional expectation of the number of Duroc alleles carried by animal j, Pdj = P12 + P21 is the conditional probability of animal j to be heterozygous, and ej is the residual error.
The HS analysis assumes the QTL to be segregating in the parental breeds, and the 6 F1 sires were regarded as common parents. QTL Express (Seaton et al., 2002) was used to calculate the probabilities of individuals inheriting allele (A1) or allele (A2) from the common F1 sire (A1 or A2) at fixed 1-cM intervals (Knott et al., 1996). In these analyses contrasts were made between the two haplotypes of every F1 sire.
yij = Xijb + si + αHSiPSij + eij
Where yij is the phenotype of F2 progeny j of F1 sire i, Xij, and b are the design matrix and the solution vector for fixed effects, respectively, si is the effect of the ith F1 sire, αHSi is the substitution effect for the two putative QTL alleles (A1 or A2) carried by the ith F1 sire, PSij is the probability that the F2 individuals inherited the arbitrary allele (Ai1) from F1 sire i, and eij is the residual error.
The CB model assumes the QTL to be segregating in the parental breeds.
yij = Xijb + sij + aPaij + dPdij + αCBiPSij + eij
Where yij is the phenotype of F2 progeny j of F1 sire i, Xij, and b are the design matrix and the solution vector for fixed effects, respectively, si is the effect of the ith F1 sire, a and d are the additive and dominance effects of breed-origin alleles, respectively, Paij and Pdij are the corresponding breed-origin coefficients as described above, αCBi is the substitution effect for the two putative QTL alleles carried by the ith F1 sire, PSij is the probability that the F2 individuals inherited the arbitrary allele (Ai1) from F1 sire i, and eijis the residual error. In this model, a and d account for the average effects of breed-origin alleles through both the F1 sire and the F1 dam and αCBi represents the difference between the two QTL alleles that a given F1 sire received from the two parental breeds as a deviation from their average additive effect (Kim et al., 2005). To avoid increasing Type I error rate due to multiple testing, a significance threshold of q < 0.05 was used, where q is the FDR corrected p-value. QTL detected using the LC, HS, or CB models were declared using the following criteria:
(1) LC QTL declared if qLC = min(qLC, qHS) < 0.05
(2) HS QTL declared if qHS = min(qLC, qHS) < 0.05
(3) CB QTL declared if qCB < 0.05 and qLC > 0.05 and qHS > 0.05
A QTL was declared under the CB model only if it had not been previously detected using the LC or HS models.
Three different models for QTL analysis revealed a total of 91 QTL for carcass and meat quality traits on all autosomes except SSC11 and 17. The LC analysis revealed 50 QTL (Table 2) including 14 new QTL on 6 chromosomes (SSC3, 6, 7, 12, 16, and 18) which had not been identified in the first genome scan of this population (Edwards et al., 2008a). The HS analysis revealed 38 QTL, and 3 additional QTL were detected using the CB model (Table 2). The thresholds used in this study were −log10(P) = 3.78 and −log10(P) = 2.88 at the 1 and 5% FDR levels, respectively. As an example, the genome scan for ham weight is shown in Figure 1. At the 1% FDR level, two QTL were identified using the LC model on SSC6 and 7, and one QTL was identified using the HS model on SSC7. At the 5% FDR level, additional QTL were revealed on SSC3 with the LC model and on SSC5, 8, and 9 with the HS model.
Figure 1. Genome scan results for ham weight determined using different analysis models. A whole genome scan to identify QTL for the trait ham weight was performed using three different analysis models (line-cross, red line; half-sib, blue line; combined line-cross and half-sib, green line). The X-axis indicates positions of chromosomes 1–18. Horizontal lines indicate significance thresholds (lower line, 5% FDR; upper line, 1% FDR).
A total of 50 significant QTL were identified on SSC1, 3–10, 12, 14, 16, and 18 using the LC model (Table 2). Of these, 29 QTL were below the 1% FDR threshold on SSC1, 3, 6, 7, 12, 14, and 18. On SSC1, QTL affecting LMA and spareribs weight detected at 12 and 236 cM supported our previous results, but a QTL for dressing percent which was significant at the 1% chromosome-wise level in the first scan of this population (Edwards et al., 2008a) did not reach significance in the second scan. On SSC3, a QTL for 45-min pH was significant at the 1% FDR level, also confirming results from our first scan (Edwards et al., 2008a), and a QTL for ham weight was newly identified at the 5% FDR level.
On SSC6, QTL for moisture and firmness were located in the S0087–S0220 interval, QTL influencing meat quality traits were mapped to the SW2173–SW1647 interval, and QTL affecting fat deposition and carcass traits were identified in the SW1647–SW1881–SW322 interval (Figure 2). The QTL detected in these marker intervals showed additive pleiotropic effects indicating that the Duroc allele contributed to increased fat deposition and reduced muscularity. In contrast to SSC6, QTL affecting muscle mass located in the SW2019–SW859 interval on SSC7 showed negative additive effects, and the Pietrain allele was associated with higher muscularity. The incorporation of the new SSC7 marker SW2019 in the SW1369–SW850 marker interval allowed refining the QTL position detected in the first scan, as well as increasing the statistical power and narrowing the QTL interval. A QTL for LMA detected in the SW859–S0115 interval in the first scan was repositioned at 86 cM in the SW2019–SW859 interval in the second scan.
Figure 2. Line-cross analysis indicated strong evidence for QTL influencing fat deposition traits on SSC6. Highly significant QTL for traits related to fat deposition were identified on SSC6. Confidence intervals for fat percentage and 10th-rib backfat (BF10) were estimated using 10,000 bootstrap permutations as 136–146 cM (blue bar) and 159–165 cM (red bar), respectively. Marker positions are shown as triangles on the X-axis (black, markers used for both QTL scans and genotyped only in 510 animals; green, markers used for both QTL scans and genotyped in all animals; red, markers used for second scan only and genotyped in all animals). Horizontal lines indicate significance thresholds (lower line, 5% FDR; upper line, 1% FDR).
On SSC12, QTL for fat related traits including marbling score, belly weight, and intramuscular fat percent detected in the SW874–S0090 interval in the first scan were identified in the second scan in the SW957–SW874 interval at the 1% FDR level. In addition, at the 5% FDR level, QTL for a* and b* not identified in the first scan were mapped to 93 and 110 cM of SSC12, respectively. A QTL for LMA mapped to 42 cM and QTL for intramuscular fat and moisture percent located at 143 cM were newly discovered on SSC16 in the second scan. In the SW2540–SW1023 interval of SSC18, not only was a QTL for spareribs weight confirmed from the first scan, but QTL for carcass length and last-lumber backfat were also newly identified in the second scan.
Half-sib analysis revealed a total of 38 QTL on SSC2, 3, 5, 7–10, 13–16, and 18 (Table 2). Of these, 20 QTL identified on SSC3, 5, 7, 8, 15, 16, and 18 were significant at the 1% FDR level including 13 QTL detected on SSC15.
A QTL affecting 45-min carcass temperature (q < 0.01) was detected at 47 cM on SSC3. On SSC5, a QTL for first-rib backfat was declared as a HS QTL (q < 0.01) in the second scan, whereas a first-rib backfat QTL had previously been identified in this location with the LC analysis in the first scan (Edwards et al., 2008a). On SSC8, a QTL affecting ham weight was identified at 39 cM, and QTL for off-flavor and cook yield were mapped to the distal region of SSC8 near S0178. In addition, a QTL for ham weight on SSC9 and a QTL for protein percent on SSC10 were identified (q < 0.01). On SSC16, HS QTL were identified for L* (q < 0.02), a* (q < 0.01), and moisture percent (q < 0.02). A highly significant HS QTL influencing spareribs weight (q < 0.0007) was detected on SSC18 with an estimated location at 70 cM. The location of the LC QTL for spareribs weight on SSC18 was estimated at 24 cM so these QTL were considered to be unique QTL.
On SSC7, QTL affecting ham weight and number of ribs were identified in the SW859–SW2040–S0115 interval (q < 0.01), and QTL for marbling score and loin weight significant at the 5% FDR level were located in the same interval. For ham weight, the QTL identified with the HS analysis was mapped to 139 cM (q < 0.0002), whereas the ham weight QTL revealed with the LC analysis was mapped to 104 cM (q < 0.006). Since these QTL detected by the different models mapped to distinct locations, they were considered to be separate unique QTL.
The HS analysis revealed evidence for QTL influencing meat quality traits in the SW1683–SW906–SW1983 interval on SSC15 (Figure 3). In the SW1683–SW906 interval, a QTL for protein percent had the highest test statistic [–log10(P) = 27.55; q < 0.0001] among the QTL detected on SSC15. In addition, a QTL for 24 h pH, a trait that is associated with many other meat quality traits, was highly significant (q < 0.0001; Figure 3). The LC analysis also revealed significant QTL for these traits in the same interval, but the HS model showed much higher statistical evidence.
Figure 3. Half-sib analysis indicated strong evidence for QTL influencing meat quality traits on SSC15. Highly significant QTL for meat quality traits were identified on SSC15. Confidence interval for protein percentage was estimated by 10,000 bootstrap permutations as 77–85 cM (gray bar). Marker positions are shown as triangles on the X-axis (black, markers used for both QTL scans and genotyped only in 510 animals; green, markers used for both QTL scans and genotyped in all animals; red, markers used for second scan only and genotyped in all animals). Horizontal lines indicate significance thresholds (lower line, 5% FDR; upper line, 1% FDR).
In addition to QTL identified with the LC and HS analyses, three additional QTL exceeded the 5% FDR significance threshold using the CB analysis. A QTL for pH decline from 45-min to 24-h was mapped to 117 cM on SSC3, a QTL for spareribs weight was detected in the SW2019–SW859 interval on SSC7 and a QTL for 24-h carcass temperature was found in the SW2540–SW1023 interval on SSC18. Although the statistical power was sufficient to detect QTL, the CB model revealed a small number of additional QTL because most QTL had been declared using either the LC or HS models due to higher test statistics with these analyses.
Effect of Additional Markers and Animals on QTL Detection
The QTL analyses under three different models revealed QTL for pH associated traits at different positions on SSC3 (Figure 4). The first scan of this chromosome using the LC model had revealed QTL for 45-min pH and pH decline from 45-min to 24-h postmortem (Edwards et al., 2008a). The second scan included two additional markers and genotypes for additional F2 pigs. A QTL for 45-min pH using the LC model (q ≤ 0.0076) was mapped at 135 cM near marker SW2047 (134.8 cM) and a QTL for pH decline from 45-min to 24-h was detected using the CB model (q ≤ 0.0469) at 117 cM near marker ACTG2 (116.5 cM) which did not reach the significance threshold in the LC analysis (q ≤ 0.065). In addition, a QTL for 45-min carcass temperature was detected using the HS model (q ≤ 0.002) located at 47 cM in the SW2021–S0206 marker interval. These results confirm results of the first scan for 45-min pH and pH decline, and add new results for 45-min carcass temperature.
Figure 4. Quantitative trait loci results determined by different models for pH related traits on SSC3. Line-cross model detected a QTL (q < 0.01) for 45 min pH (red line) and half-sib model detected a QTL (q < 0.01) for 45 min carcass temperature (blue line). Combined model identified a QTL (q < 0.01) affecting pH decline from 45 min to 24 h (green line). Marker positions are shown as triangles on the X-axis (black, markers used for both QTL scans and genotyped only in 510 animals; green, markers used for both QTL scans and genotyped in all animals; red, markers used for second scan only and genotyped in all animals). Horizontal lines indicate significance thresholds (lower line, 5% FDR; upper line, 1% FDR).
We have recently used LC, HS, and CB models to identify QTL for growth traits in our Duroc × Pietrain resource population, and we reported that additional markers and animals contributed to reduce the confidence intervals and increase the test statistics for QTL detection (Choi et al., 2010). For the present study, QTL affecting the a* and b* objective color measures were newly detected on SSC12 (q < 0.04). In order to determine how the QTL peaks for these traits were changed, analyses were performed under 4 different scenarios; 5 and 7 markers with 510 and 948 animals (Figure 5). The results indicated that increasing the number of animals or increasing the number of markers was effective in increasing the power to detect QTL on this chromosome, and that increasing the numbers of both animals and markers allowed detection of the a* and b* QTL.
Figure 5. Effect of additional markers and animals for detecting meat color QTL on SSC12. Effects of additional marker genotypes and animals for detecting QTL for a* and b* objective meat color measures on SSC12 were compared under 4 different scenarios (5 and 7 markers with 510 F2 animals, 5 and 7 markers with 948 F2 animals). Blue lines indicate 5 markers (black and green triangles on the X-axis) and red lines indicate 7 markers (black, green and red triangles on the X-axis). Solid lines indicate 948 animals and dotted lines indicate 510 animals. Marker positions are shown as triangles on the X-axis (black, markers used for both QTL scans and genotyped only in 510 animals; green, markers used for both QTL scans and genotyped in all animals; red, markers used for second scan only and genotyped in all animals). Horizontal lines indicate significance thresholds (lower line, 5% FDR; upper line, 1% FDR).
This study identified 91 QTL for pig carcass and meat quality traits located on all autosomes except SSC11 and 17 using three least-squares Mendelian analysis models. The LC analysis, which detected QTL segregating between breeds, revealed 50 QTL including 13 new QTL on 6 chromosomes (SSC3, 6, 7, 12, 16, and 18) that had not been identified in the first genome scan of this population (Edwards et al., 2008a). The HS analysis, which detected QTL segregating within breeds, revealed 38 QTL including 18 on SSC15. Three additional QTL were detected using the CB model (Kim et al., 2005).
Application of the three different models for SSC3 identified not only QTL influencing muscularity under the LC model, but also QTL affecting pH and carcass temperature using all three models. The LC QTL for 45-min pH detected at 135 cM near SW2047 confirmed the 45-min pH QTL observed in the first scan (Edwards et al., 2008a). Beeckmann et al. (2003) reported a QTL for 45-min pH at the same interval in a Wild boar × Meishan F2 population. Several studies (Óvilo et al., 2002a; de Koning et al., 2003; Evans et al., 2003; Wimmers et al., 2006) reported QTL affecting muscle pH in the SW2021–S0206 marker interval, a region where we identified a QTL for 45-min carcass temperature under the HS model. In addition, Duan et al. (2009) reported a QTL for pH decline from 45-min to 3-h in the SW2021–S0206 interval in a White Duroc × Chinese Erhualian population. We also detected a QTL for pH decline from 45-min to 24-h under the CB model, however, our QTL was located at 117 cM near ACTG2.
Significant QTL affecting backfat thickness were located on SSC6 within the SW1647–SW1881–SW322 marker interval at 160–174 cM. A 10,000 bootstrap permutation analysis showed the average QTL positions for each backfat trait to be located in the 160.12–167.96 cM region. The 95% confidence interval for 10th-rib backfat was estimated to be 159–165 cM (6 Haldane cM), which was considerably narrowed from the 38.5 Haldane cM interval observed for the first scan (Edwards et al., 2008a). Not only were QTL affecting fat deposition traits observed in this region, but QTL influencing muscularity were also identified at the same marker interval since Duroc alleles contributed to both fat accumulation and reduced muscle content.
Our results for backfat thickness traits were in agreement with other studies (Malek et al., 2001b; Óvilo et al., 2002b; Varona et al., 2002) that identified QTL for fatness traits in this region of SSC6. This region includes the leptin receptor (LEPR) gene which is considered as a potential candidate gene for fatness (Ernst et al., 1997; Óvilo et al., 2005; Mohrmann et al., 2006), and studies to identify a causal mutation in LEPR have been conducted (Mackowski et al., 2005; Muñoz et al., 2009). We also observed a QTL for intramuscular fat percent using the LC model in a position more proximal to this region of SSC6 at 141 cM, which coincided with a backfat thickness QTL detected with the HS model in a Duroc × Pietrain population by Liu et al. (2008). The SSC6 region affecting intramuscular fat percent also included QTL for marbling score and a*, which were all mapped to 141–152 cM in the SW2173–SW1647 interval. The confidence interval for these QTL did not overlap with the confidence interval for 10th-rib backfat. This result was consistent with previous studies (Szyda et al., 2003; Óvilo et al., 2005) which reported that QTL effects for backfat and intramuscular fat content resulted from different closely linked loci on SSC6. QTL affecting intramuscular fat content have been reported (de Koning et al., 2000; Grindflek et al., 2001) in the same region where we detected a QTL for marbling score, although no other reports of subjective marbling score QTL in this SSC6 region have been reported. Also, Harmegnies et al. (2006) identified QTL for a* as well as fat thickness in this same region.
The different models revealed distinct QTL regions on SSC7 with LC and HS QTL identified at 84–104 cM and at 130–178 cM, respectively. A highly significant QTL influencing muscle mass identified in the SW2019–SW859 interval had an additive effect for which Duroc alleles increased carcass length and decreased dressing percent, LMA and ham weight. In this region, Yue et al. (2003) found a 1% genome-wide level significant QTL influencing carcass composition traits such as carcass length in a Wild boar × Meishan population. Liu et al. (2008) reported QTL for carcass length and dressing percent with similar allelic substitution effect in their Duroc × Pietrain population as we observed in our study. However, Nezer et al. (2002) identified a QTL for carcass length at the more distal position from our QTL in a Pietrain × Large White population. In addition, Sato et al. (2003) detected a QTL for dressing percent in a Duroc × Meishan population in the same region as our study. A QTL for number of ribs was detected using the HS analysis. A QTL for number of ribs had been detected in this position at the 1% genome-wise significance level using the LC analysis in the first scan (Edwards et al., 2008a), however, evidence from the second scan suggests the HS model better describes the QTL allele frequency in the parental breeds. Also on SSC7, analyses using both the LC and HS models identified QTL for ham weight at different locations, which were in the SW2019–SW859 interval with the LC model and in the SW859–SW2040 interval with the HS model. Similarly, Milan et al. (2002) also reported suggestive QTL for ham weight at different positions using LC and HS models, and their LC QTL detected in the SLA–S0102 marker interval was in a similar region to our LC QTL.
We have recently used LC, HS, and CB models to identify QTL for growth traits in our Duroc × Pietrain resource population, and we reported that additional markers and animals contributed to reduce the confidence intervals and increase the test statistics for QTL detection (Choi et al., 2010). In the present study, genotyping of additional markers and animals increased the statistical power and facilitated discovery of new QTL which had not been observed in the first scan (Edwards et al., 2008a). For example, QTL for the objective color measures of a* and b* were identified on SSC12 with the addition of more F2 pigs and more marker genotypes using the LC analysis. The LC analysis also identified QTL on SSC12 related to intramuscular fat percent and moisture at 47–50 cM and at 69 cM, respectively. A QTL for marbling was located in the SW957–SW874 marker interval, whereas the position of this QTL had been more distal for the first scan (Edwards et al., 2008a). The additive effects of these QTL indicated that Duroc alleles increased marbling and intramuscular fat percent, and decreased moisture percent. Harmegnies et al. (2006) also reported a QTL for a* although at a more distal position than our current result, and Malek et al. (2001a) detected a QTL for subjective color score in the same region as our result.
The HS analysis revealed strong evidence for QTL affecting meat quality traits on SSC15 at 74–90 cM in the SW1683–SW1983 marker interval where 13 and 2 QTL were significant at the 1 and 5% FDR levels, respectively, including a highly significant QTL for protein percent. Significant QTL had been identified in this region using the LC analysis in the first scan (Edwards et al., 2008a), and a negative additive effect had been seen for protein percent, color, and tenderness traits suggesting contributions from segregation of Pietrain alleles. The pleiotropic effects of Pietrain alleles contributing to leanness resulted in effects on other meat quality traits resulting in more muscularity, paler muscle color, and less tenderness. QTL for 24-h pH, L*, and tenderness significant at the 1% genome-wise level were identified in this region of SSC15 in a Berkshire × Yorkshire population (Malek et al., 2001a; Thomsen et al., 2004; Kim et al., 2005). Very few studies have measured the trait of protein percent and no QTL for protein percent have been reported on SSC15.
Several candidate genes such as myostatin (MSTN), Titin (TTN), and protein kinase AMP-activated gamma 3 (PRKAG3) are located in the SW1683–SW1983 interval (Sonstegard et al., 1998; Milan et al., 2000; Davoli et al., 2003). Stinckens et al. (2008) reported that pigs of the Pietrain breed had higher muscularity as a result of association between a polymorphism in the MSTN MEF3 binding site and muscle mass. Wimmers et al. (2007) reported that a polymorphism in TTN was associated with leanness in a Pietrain population. Also, Milan et al. (2000) mapped the PRKAG3 gene in the SW1683–SW1983 marker interval and Ciobanu et al. (2001) identified three polymorphic sites in the PRKAG3 gene that affect meat quality traits including 24-h pH.
This study used a Duroc × Pietrain F2 resource population and identified a total of 91 QTL for carcass merit and meat quality phenotypes. Three different least-squares models were applied under different assumptions; (1) founders fixed for alternative QTL alleles (LC model), (2) segregation of QTL alleles at similar frequencies in founders (HS model), and (3) segregation of QTL alleles at different frequencies in founders (combined model). The addition of new marker and animal genotypes contributed to increasing the statistical power for QTL detection, and the application of alternative models allowed confirmation of QTL and detection of new QTL segregating either between or within breeds.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This project was supported by National Research Initiative Competitive Grant no. 2004-35604-14580 from the USDA National Institute of Food and Agriculture. We thank D. Edwards for technical contributions in development of the pig resource population.
Affentranger, P., Gerwig, C., Seewer, G. J. F., Schworer, D., and Kunzi, N. (1996). Growth and carcass characteristics as well as meat and fat quality of three types of pigs under different feeding regimens. Livest. Prod. Sci. 45, 187–196.
Choi, I., Steibel, J. P., Bates, R. O., Raney, N. E., Rumph, J. M., and Ernst, C. W. (2010). Application of alternative models to identify QTL for growth traits in an F2 Duroc x Pietrain pig resource population. BMC Genet. 11, 97. doi: 10.1186/1471-2156-11-97
Ciobanu, D., Bastiaansen, J., Malek, M., Helm, J., Woollard, J., Plastow, G., and Rothschild, M. (2001). Evidence for new alleles in the protein kinase adenosine monophosphate-activated (3-subunit gene associated with low glycogen content in pig skeletal muscle and improved meat quality. Genetics 159, 1151–1162.
Davoli, R., Braglia, S., Lama, B., Fontanesi, L., Buttazzoni, L., Baiocco, C., and Russo, V. (2003). Mapping, identification of polymorphisms and analysis of allele frequencies in the porcine skeletal muscle myopalladin and titin genes. Cytogenet. Genome Res. 102, 152–156.
de Koning, D. J., Pong-Wong, R., Varona, L., Evans, G. J., Giuffra, E., Sanchez, A., Plastow, G., Noguera, J. L., Andersson, L., and Haley, C. S. (2003). Full pedigree quantitative trait locus analysis in commercial pigs using variance components. J. Anim. Sci. 81, 2155–2163.
de Koning, D. J., Rattink, A. P., Harlizius, B., van Arendonk, J. A., Brascamp, E. W., and Groenen, M. A. (2000). Genome-wide scan for body composition in pigs reveals important role of imprinting. Proc. Natl. Acad. Sci. U.S.A. 97, 7947–7950.
Duan, Y. Y., Ma, J. W., Yuan, F., Huang, L. B., Yang, K. X., Xie, J. P., Wu, G. Z., and Huang, L. S. (2009). Genome-wide identification of quantitative trait loci for pork temperature, pH decline, and glycolytic potential in a large-scale White Duroc x Chinese Erhualian resource population. J. Anim. Sci. 87, 9–16.
Edwards, D. B., Ernst, C. W., Raney, N. E., Doumit, M. E., Hoge, M. D., and Bates, R. O. (2008a). Quantitative trait locus mapping in an F2 Duroc x Pietrain resource population: II. Carcass and meat quality traits. J. Anim. Sci. 86, 254–266.
Edwards, D. B., Ernst, C. W., Tempelman, R. J., Rosa, G. J. M., Raney, N. E., Hoge, M. D., and Bates, R. O. (2008b). Quantitative trait loci mapping in an F2 Duroc x Pietrain resource population: I. Growth traits. J. Anim. Sci. 86, 241–253.
Evans, G. J., Giuffra, E., Sanchez, A., Kerje, S., Davalos, G., Vidal, O., Illán, S., Noguera, J. L., Varona, L., Velander, I., Southwood, O. I., de Koning, D. J., Haley, C. S., Plastow, G. S., and Andersson, L. (2003). Identification of quantitative trait loci for production traits in commercial pig populations. Genetics 164, 621–627.
Harmegnies, N., Davin, F., De Smet, S., Buys, N., Georges, M., and Coppieters, W. (2006). Results of a whole-genome quantitative trait locus scan for growth, carcass composition and meat quality in a porcine four-way cross. Anim. Genet. 37, 543–553.
Jeremiah, L. E., Gibson, J. P., Gibson, L. L., Ball, R. O., Aker, C., and Fortin, A. (1999). The influence of breed, gender, and PSS (Halothane) genotype on meat quality, cooking loss, and palatability of pork. Food Res. Int. 32, 59–71.
Liu, G., Kim, J. J., Jonas, E., Wimmers, K., Ponsuksili, S., Murani, E., Phatsara, C., Tholen, E., Juengst, H., Tesfaye, D., Chen, J. L., and Schellander, K. (2008). Combined line-cross and half-sib QTL analysis in Duroc-Pietrain population. Mamm. Genome 19, 429–438.
Mackowski, M., Szymoniak, K., Szydlowski, M., Kamyczek, M., Eckert, R., Rozycki, M., and Switonski, M. (2005). Missense mutations in exon 4 of the porcine LEPR gene encoding extracellular domain and their association with fatness traits. Anim. Genet. 36, 135–137.
Malek, M., Dekkers, J. C., Lee, H. K., Baas, T. J., Prusa, K., Huff-Lonergan, E., and Rothschild, M. F. (2001a). A molecular genome scan analysis to identify chromosomal regions influencing economic traits in the pig. II. Meat and muscle composition. Mamm. Genome 12, 637–645.
Malek, M., Dekkers, J. C., Lee, H. K., Baas, T. J., and Rothschild, M. F. (2001b). A molecular genome scan analysis to identify chromosomal regions influencing economic traits in the pig. I. Growth and body composition. Mamm. Genome 12, 630–636.
Milan, D., Bidanel, J. P., Iannuccelli, N., Riquet, J., Amigues, Y., Gruand, J., Le Roy, P., Renard, C., and Chevalet, C. (2002). Detection of quantitative trait loci for carcass composition traits in pigs. Genet. Sel. Evol. 34, 705–728.
Milan, D., Jeon, J. T., 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, P. L., Chardon, P., and Andersson, L. (2000). A mutation in PRKAG3 associated with excess glycogen content in pig skeletal muscle. Science 288, 1248–1251.
Mohrmann, M., Roehe, R., Knap, P. W., Looft, H., Plastow, G. S., and Kalm, E. (2006). Quantitative trait loci associated with AutoFOM grading characteristics, carcass cuts and chemical body composition during growth of Sus scrofa. Anim. Genet. 37, 435–443.
Muñoz, G., Óvilo, C., Silió, L., Tomás, A., Noguera, J. L., and Rodríguez, M. C. (2009). Single- and joint-population analyses of two experimental pig crosses to confirm quantitative trait loci on Sus scrofa chromosome 6 and leptin receptor effects on fatness and growth traits. J. Anim. Sci. 87, 459–468.
Nezer, C., Moreau, L., Wagenaar, D., and Georges, M. (2002). Results of a whole genome scan targeting QTL for growth and carcass traits in a Pietrain x Large White intercross. Genet. Sel. Evol. 34, 371–387.
Óvilo, C., Clop, A., Noguera, J. L., Oliver, M. A., Barragán, C., Rodríguez, C., Silió, L., Toro, M. A., Coll, A., Folch, J. M., Sánchez, A., Babot, D., Varona, L., and Pérez-Enciso, M. (2002a). Quantitative trait locus mapping for meat quality traits in an Iberian x Landrace F2 pig population. J. Anim. Sci. 80, 2801–2808.
Óvilo, C., Oliver, A., Noguera, J. L., Clop, A., Barragán, C., Varona, L., Rodríguez, C., Toro, M., Sánchez, A., Pérez-Enciso, M., and Silió, L. (2002b). Test for positional candidate genes for body composition on pig chromosome 6. Genet. Sel. Evol. 34, 465–479.
Óvilo, C., Fernández, A., Noguera, J. L., Barragán, C., Letón, R., Rodríguez, C., Mercadé, A., Alves, E., Folch, J. M., Varona, L., and Toro, M. (2005). Fine mapping of porcine chromosome 6 QTL and LEPR effects on body composition in multiple generations of an Iberian by Landrace intercross. Genet. Res. 85, 57–67.
Sato, S., Oyamada, Y., Atsuji, K., Nade, T., Kobayashi, E., Mitsuhashi, T., Nirasawa, K., Komatsuda, A., Saito, Y., Terai, S., Hayashi, T., and Sugimoto, Y. (2003). Quantitative trait loci analysis for growth and carcass traits in a Meishan x Duroc F2 resource population. J. Anim. Sci. 81, 2938–2949.
Stinckens, A., Luyten, T., Bijttebier, J., Van den Maagdenberg, K., Dieltiens, D., Janssens, S., De Smet, S., Georges, M., and Buys, N. (2008). Characterization of the complete porcine MSTN gene and expression levels in pig breeds differing in muscularity. Anim. Genet. 39, 586–596.
Thomsen, H., Lee, H. K., Rothschild, M. F., Malek, M., and Dekkers, J. C. (2004). Characterization of quantitative trait loci for growth and meat quality in a cross between commercial breeds of swine. J. Anim. Sci. 82, 2213–2228.
Varona, L., Ovilo, C., Clop, A., Noguera, J. L., Perez-Enciso, M., Coll, A., Folch, J. M., Barragan, C., Toro, M. A., Babot, D., and Sanchez, A. (2002). QTL mapping for growth and carcass traits in an Iberian by Landrace pig intercross: additive, dominant and epistatic effects. Genet. Res. 80, 145–154.
Wimmers, K., Fiedler, I., Hardge, T., Murani, E., Schellander, K., and Ponsuksili, S. (2006). QTL for microstructural and biophysical muscle properties and body composition in pigs. BMC Genet. 7, 15. doi: 10.1186/1471-2156-7-15
Wimmers, K., Murani, E., Te Pas, M. F. W., Chang, K. C., Davoli, R., Merks, J. W. M., Henne, H., Muraniova, M., Da Costa, N., Harlizius, B., Schellander, K., Böll, I., Braglia, S., De Wit, A. A. C., Cagnazzo, M., Fontanesi, L., Prins, D., and Ponsuksili, S. (2007). Associations of functional candidate genes derived from gene-expression profiles of prenatal porcine muscle tissue with meat quality and muscle deposition. Anim. Genet. 38, 474–484.
Yue, G., Stratil, A., Cepica, S., Schrffel, J., Schrffelova, D., Fontanesi, L., Cagnazzo, M., Moser, G., Bartenschlager, H., Reiner, G., and Geldermann, H. (2003). Linkage and QTL mapping for Sus scrofa chromosome 7. J. Anim. Breed. Genet. 120, 56–65.
Keywords: pig, QTL, carcass merit, meat quality
Citation: Choi I, Steibel JP, Bates RO, Raney NE, Rumph JM and Ernst CW (2011) Identification of carcass and meat quality QTL in an F2 Duroc × Pietrain pig resource population using different least-squares analysis models. Front. Gene. 2:18. doi: 10.3389/fgene.2011.00018
Received: 08 February 2011; Accepted: 13 April 2011;
Published online: 02 May 2011.
Edited by:Jean Michel Claude Elsen, Institut National de la Recherche Agronomique, France
Reviewed by:Tosso Leeb, University of Bern, Switzerland
Juliette Riquet, Institut National de la Recherche Agronomique, France
Copyright: © 2011 Choi, Steibel, Bates, Raney, Rumph and Ernst. This is an open-access article subject to a non-exclusive license between the authors and Frontiers Media SA, which permits use, distribution and reproduction in other forums, provided the original authors and source are credited and other Frontiers conditions are complied with.
*Correspondence: Catherine W. Ernst, Department of Animal Science, Michigan State University, 3385 Anthony Hall, East Lansing, MI 48824-1225, USA. e-mail: email@example.com