Evaluation of Genome Based Estimated Breeding Values for Meat Quality in a Berkshire Population Using High Density Single Nucleotide Polymorphism Chips

Article information

Asian-Australas J Anim Sci. 2014;27(11):1540-1547
1Dasan Pig Breeding Co., Namwon 590-831, Korea.
2Gyeongbuk Livestock Research Institute, Yeungju 750-780, Korea.
3Department of Animal Science and Technology, Sunchon National University, Suncheon 540-950, Korea.
4Animal Genomics and Bioinformatics Division, National Institute of Animal Science, Suwon 441-706, Korea.
5Genomic Informatics Center, Hankyong National University, Anseong 456-649, Korea.
*Corresponding Author: Jong-Joo Kim. Tel: +82-53- 810-3027, Fax: +82-53-801-3027, E-mail: kimjj@ynu.ac.kr
School of Biotechnology, Yeungnam University, Gyeongsan 712-749, Korea
aThe two authors equally contributed to this work.
Received 2014 May 19; Revised 2014 July 03; Accepted 2014 July 14.

Abstract

The accuracy of genomic estimated breeding values (GEBV) was evaluated for sixteen meat quality traits in a Berkshire population (n = 1,191) that was collected from Dasan breeding farm, Namwon, Korea. The animals were genotyped with the Illumina porcine 62 K single nucleotide polymorphism (SNP) bead chips, in which a set of 36,605 SNPs were available after quality control tests. Two methods were applied to evaluate GEBV accuracies, i.e. genome based linear unbiased prediction method (GBLUP) and Bayes B, using ASREML 3.0 and Gensel 4.0 software, respectively. The traits composed different sets of training (both genotypes and phenotypes) and testing (genotypes only) data. Under the GBLUP model, the GEBV accuracies for the training data ranged from 0.42±0.08 for collagen to 0.75±0.02 for water holding capacity with an average of 0.65±0.04 across all the traits. Under the Bayes B model, the GEBV accuracy ranged from 0.10±0.14 for National Pork Producers Council (NPCC) marbling score to 0.76±0.04 for drip loss, with an average of 0.49±0.10. For the testing samples, the GEBV accuracy had an average of 0.46±0.10 under the GBLUP model, ranging from 0.20±0.18 for protein to 0.65±0.06 for drip loss. Under the Bayes B model, the GEBV accuracy ranged from 0.04±0.09 for NPCC marbling score to 0.72±0.05 for drip loss with an average of 0.38±0.13. The GEBV accuracy increased with the size of the training data and heritability. In general, the GEBV accuracies under the Bayes B model were lower than under the GBLUP model, especially when the training sample size was small. Our results suggest that a much greater training sample size is needed to get better GEBV accuracies for the testing samples.

INTRODUCTION

Pork meat constitutes more than 40% of worldwide meat production, forming a prominent source of human food (Rothschild et al., 2011). Customer demand for pork depends upon the meat quality and its physical and biochemical components (Bonneau et al., 2010). Thus, pork quality is an economically important factor and one of the major selection benchmarks for breeding process in the swine industry (Luo et al., 2012). With the advancements in meat processing technologies, it is now possible to predict ham, loin primal, and sub-primal cut weights (van Wijk et al., 2005). These developments paved the way for the pork industry to adopt more precise value-based classifying systems to satisfy the demands of pork market (Brorsen et al., 1998).

Pig breeding programs have been implemented mainly towards the reduction of production cost in the last several decades. Selection was, therefore, focused on economically important traits such as litter size, weight gain, back fat, and feed conversion. Recently, breeding objectives have also put weight on retail carcass yield and meat quality (van Wijk et al., 2005). Carcass quality traits are highly heritable, causing an efficient selection response (Newcom et al., 2002).

The Berkshire breed (Sus scrofa domesticus) was found 300 years ago in the Berkshire county of United Kingdom. The king of England preferred Berkshire pork for his own personal meat supply, because of excellent meat quality (American Berkshire Association, 2013). They have a dark skin, which protects them from sunburn. Some parts of body such as legs, face and tail, are white pointed. Berkshire pigs are characterized with pink skin color and a strong body type with short neck. Legs of the breed are short and blocky and feet are strong. An adult pig of the breed has an average weight of 272 kg (600 pounds). Individuals of the breed are usually friendly and curious and exhibit excellent disposition (Kawaida, 1993).

Recently, the demand for Berkshire pork has increased in Asia (McLaughlin, 2004), because of excellent meat quality such as richness, texture, marbling, juiciness, tenderness and flavor (Goodwin and Burroughs, 1995; Brewer et al., 2002). The thin muscle fibers and excellent water holding capacity (WHC) of the meat increases its popularity (Goodwin and Burroughs, 1995). The Berkshire pork is exported mainly to Japan and US (Lammers et al., 2011).

The evaluation of breeding values based on pedigree information has many limitations (Dekkers et al., 2010). Some phenotypic traits are difficult and costly to measure, resulting in low accuracy of estimated breeding values (Badke et al., 2014). These shortcomings could be overcome by using a genome based best linear unbiased prediction method (GBLUP), in which genomic estimated breeding values (GEBV) were predicted with a high density marker map covering the porcine whole genome (Meuwissen et al., 2001; Dekkers et al., 2010).

Genomic selection is a selection decision made at an early stage based on GEBV (Hayes et al., 2009). The main benefits of genomic selection are genetic enhancement by minimizing the generation interval and yielding a higher accuracy of estimated breeding values (Hayes et al., 2009; van Raden et al., 2009; Christensen et al., 2010; Wiggans et al., 2011; Boddhireddy et al., 2014).

The GEBVs can be predicted with dense, genome-wide maps of single nucleotide polymorphisms (SNPs), which can lead to a more precise prediction of pig breeding value at a young age. Cleverland et al. (2010) reported that GEBV accuracies in pigs were as good as those in dairy cattle, if the training population size was large enough. Important factors on the accuracy of genomic predictions include number of phenotypes, (i.e. training data, used to form the prediction equation), heritability, effective population size, genome size, marker density, and genetic architecture of the trait, in particular number of loci affecting the trait and distribution of their effects (Daetwyler et al., 2008; Goddard, 2009; Meuwissen, 2009). Recently, Bayesian methods have gained popularity in evaluating genomic selection, due to the fact that different variance is fitted to each SNP (Fernando et al., 2007; Moser et al., 2009).

In this study, we evaluated accuracy of GEBV for port quality traits in a Berkshire population under the GBLUP and Bayes B models.

MATERIALS AND METHODS

Animals and phenotypes

A set of Berkshire samples (n = 1,205) were collected in Dasan breeding farm, Namwon, Cheonbuk province, Korea, between 2008 and 2013. The piglets were weaned at 3 to 4 weeks of age and moved into piglet pens, in each of which about 100 piglets were raised for 60 days. Then, the pigs were placed in growth/fattening pens 20 pigs in size for 90 to 120 days. The pigs were fed with the commercial feeds according to the regiments of Purina Ltd. The samples were slaughtered approximately 211 (±23) days of age in an abattoir in Namwon and cooled at 0°C for 24 h in a chilling room. Among the carcasses, 801 and 404 samples were analyzed for meat quality and composition in the laboratories of National Livestock Research Institute in Suwon and Sunchon National University in Sunchon, Korea, respectively. A total of 16 carcass and meat quality traits were considered: 1) back fat thickness (BF), 2) Commission Internationale de l’Eclairage (CIE) a, 3) CIE b, 4) CIE l 5) collagen, 6) carcass weight (CWT), 7) drip loss, 8) fat, 9) heat (cooking) loss, 10) moisture, 11) National Pork Producers Council (NPCC) color score, 12) NPCC marbling score, 13) pH24, 14) protein, 15) shear force, and 16) WHC. For each individual, slaughter age (sage), gender, and year-season of birth were recorded.

The CWT and BF of each carcass were measured. The parts of loins (longissimus dorsi, LD) on the left side of the cold carcasses were used to determine meat quality parameters. As soon as all samples were placed in vacuum bags, the samples were transported to the laboratory and then frozen at −50°C until they were analyzed. The middle portions of each loin were used for experiment. For analysis of moisture content, fat content, drip loss, and heat loss, the only subcutaneous fat of meat samples was removed. For the others, all visible fat was trimmed off.

The proximate composition of each LD muscle was obtained with a slightly modified method of AOAC (2000). Briefly, moisture content was measured by drying 3 g of samples place in aluminum dishes at 104°C for 15 h. The crude protein contents were measured by the Kjeldahl method (VAPO45, Gerhardt Ltd., Idar-Oberstein, Germany). The crude fat contents were extracted according to the method described by Folch and Sloane-Stanley (1957). The total collagen content was determined by measuring hydroxyproline and using a multiplication factor of 7.14 (Etherington and Sims, 1981). The surface color and marbling scores of each loin was categorized based on NPPC standard (NPPC, 2000).

The surface color value were measured by the CIE L*, a* and b* system using a Minolta colorimeter (Model CR-410, Minolta Co. Ltd., Osaka, Japan). The colorimeter was calibrated against a white reference tile plate (L* = 89.2, a* = 0.921, b* = 0.783), and the diameter size of aperture was 4 cm. The color L* (lightness), a* (redness), and b* (yellowness) values were obtained after 30 min blooming at room temperature. The average value of five random measurements taken from different locations was used for the statistical analysis.

The pH value of each meat sample was determined with a pH meter (Orion 2 Star, Thermo scientific, Beverly, MA, USA). Water holding capacity was determined by centrifugation. Briefly, 5 g of minced meat sample was placed into a centrifuge tube with a filter paper (No. 4, Whatman International Ltd., Maidstone, England), and centrifuged at 3,000×g for 10 min. WHC was calculated as the remaining moisture in the meat sample on the basis of the moisture content of the original meat sample. The drip loss was measured as the percentage weight loss of a standardized (3×3×3 cm) meat sample placed in a sealed petri-dish at 4°C during the storage of 2 d. The heat (cooking) loss was determined as the percentage weight loss of a standardized (3×3×3 cm) meat sample after cooking in an electric grill with double pans (Nova EMG-533, 1,400 W, Evergreen enterprise, Seoul, Korea) for 90 s, until the internal temperature of the meat sample reached 72°C.

The samples were prepared in a cubic form (30×30×20 mm), heated until internal temperature of the samples reached 72°C±2°C, and then cooled for 30 min at room temperature. Each sample was cut perpendicular to the longitudinal orientation of the muscle fiber with a Warner-Bratzler shear attachment on a texture analyzer (TA-XT2, Stable Micro System Ltd., Surrey, UK). The maximum shear force value (kg) was recorded for each sample. Test and pre-test speeds were set at 2.0 mm/s. Post-test speeds were set at 5.0 mm/s.

Molecular data

The 1,205 pigs were genotyped with the Illumina Porcine 62 k SNP chips, in which a total of 62,163 SNPs that covered the entire porcine genome were embedded. To evaluate GEBV, the SNPs on 18 autosomal chromosomes were considered for quality control tests. Any SNP was excluded with <90% call rates, <5% minor allele frequency, or significant departure from Hardy Weinberg equilibrium (p<0.001). Those individuals with less than 90% genotyping call rate were also removed. After the quality control procedures using PLINK v7.0 (Purcell et al., 2007), 36,605 SNPs for 1,191 individuals were used. An imputation procedure was applied to predict missing genotypes of the SNPs with Beagle vs3.3.2 (Browning and Browning, 2007; Nothnagel et al., 2009).

Statistical analysis

To obtain GEBV under the GBLUP model, a linear mixed (Animal) model was fitted with the fixed effects of gender for pH24 and BF, and birth year-season for CWT, pH24, CIE a, CIE b, CIE L, drip loss, heat loss and shear force, and a covariate, slaughter age for fat, protein and BF, respectively. Statistical significance of the fixed factors or covariate was tested using SAS general linear model procedure (vs9.2). For the rest of meat quality traits, the effects of the factors were not significant (p>0.05). The Animal model, then, can be written as

y=Xb+Zg+e

where y is the vector of phenotypic record of the animal, b is the vector of overall mean, fixed and covariate effects, g is the vector of breeding values, X is the design matrix for the fixed and covariate effects, Z is the design matrix allocating records to breeding values and e is the vector of the residual of the phenotype. To construct genome relationship matrix (G), the subroutine in R was used (version 2.15.0), which was then incorporated into the Animal model using ASREML (average sparsity residual maximum likelihood) 3.0 (Gilmour et al., 1995). The mixed model equation was then

[XXXZZXZZ+G-1α][b^g^]=[XyZy]

where α = σe2g2 = (1–h2)/h2, σe2 is the residual variance, σg2 the genetic variance, and h2 heritability.

The breeding values for both phenotyped and non-phenotyped individuals can be predicted by solving the equation:

[g^]=[ZZ+G-1(σe2/σg2)]-1[Zy]

G matrix was calculated based on the observed allele frequencies of the markers. The equation used to calculate G matrix was:

G=(M-P)(M-P)2j=1mpj(1-pj)

The marker matrix, M, had order of n×m, in which n is the number of individuals and m is the number of markers. In the M, alleles were coded as AA (homozygous for the first allele) = −1, AB (heterozygous) = 0, BB (homozygous for the second allele) = 1. The elements of P matrix were calculated using the formula Pj = 2(Pj–0.5), where Pj was the minor allele frequency of the marker locus j. (M-P) is called the incidence matrix (Z) for markers. The P matrix was subtracted from the M matrix to set the mean values of the allele effects to 0, and to give more credit to rare alleles than to common alleles. Genomic inbreeding coefficient would be greater if the individual is homozygous for rare alleles than if homozygous for common alleles.

The meat quality traits composed different sets of training (with phenotype records) and testing (without phenotypes) data. For each individual, GEBV value was predicted, and the expected accuracy of GEBV for the ith individual was calculated using standard errors of GEBV as

1-[(stderror(gi))2/σg2]

The mean and standard deviation of the GEBV accuracies was calculated for each trait.

Evaluation of the GEBV was also carried out under the Bayes B model using Gensel 4.0 (Fernando and Garrick, 2008) software. The Bayes model was:

y=μ+Xb+i=1i=36605ziai+e

where

  • y = the vector of phenotypes

  • μ = overall mean.

  • X = the incidence matrix of the fixed and covariate effects.

  • b = the vector of fixed and covariate effects.

  • zi = a vector of genotypes of a fitted marker i, that is coded as −10, 0, or 10.

  • ai = a random substitution effect of the fitted marker i with its variance, σai2.

  • e = the vector of random residuals that was assumed to be normally distributed.

For the marker effects, a mixture model was applied, i.e. a fraction of markers (π) with zero effect and 1-π of markers with non-zero effects, which was used to predict GEBV. Then, the genetic variances of the markers with non-zero effects would have σai2>0 (Habier et al., 2011).

The π values ranged between 0.996 and 0.999 depending on the traits with different sizes of the training data.

The fixed effects for each trait were fitted as under the GBLUP model. The estimates of genetic and residual variances that were obtained from ASREML analysis were used as prior values for the Bayes B analysis. A total of 41,000 iterations of Markov chain were run for the analyses, with the first 1,000 iterations of burn-in period and each of 100 iterations was selected to calculate posterior mean and variance for the marker effects. The GEBVs were based on a weighted sum of the number of copies of the more frequent allele at each SNP locus, with the weights being the estimated allele substitution effects (β). The sum of all the marker scores for an individual gave the genomic breeding value.

GEBV=i=0i=36605zig^i

RESULTS AND DISCUSSION

Summary statistics for the sixteen meat quality traits were displayed in Table 1. The coefficient of variation were various between traits, e.g. 1.4% for moisture and 51% for drip loss. A set of 36,605 SNPs was chosen from the 62,163 SNPs in Illumina Porcine 60 k Beadchip (Table 2). The number of SNPs (4,426) was the greatest in Sus scrofa chromosome (SSC) 1, while SSC18 had the smallest number of SNPs (886). The SSCs 2, 4, 6, 7, 8, 9, 13, and 14 had more than 2,000 SNPs. The physical map with all of the available SNPs spanned about 2,195 Mb with an average distance of 67.9±106.7 Kb between adjacent SNPs. However, the average distances were various between chromosomes, ranging between 51.9 Kb in SSC14 and 92.4 in SSC15.

Summary statistics for 16 meat quality traits in a Berkshire population

The number of available SNPs and average distances between adjacent SNPs in the 18 Sus scrofa autosomes (SSC) of the Berkshire pig population

The heritabilities that were estimated using genome relationship matrix (G) ranged between 6% and 46% (Table 3). Tomiyama et al. (2011) reported 0.54 and 0.32 for (BF) at finish and CWT, respectively, in a Japanese Berkshire population (n = 4,773). Jung et al. (2011) reported that heritabilities of pH2 4 h, CIE a, b, and L, WHC, NPPC marbling, drip loss, heat loss and shear force ranged between 0.51 and 0.66 in a Berkshire population (n = 808), Korea. Compared with the two reports, the heritability estimates of the meat quality traits in this study were low.

Accuracies of genomic estimated breeding value (GEBV) under the GBLUP and Bayes B models for the training and testing samples in the Berkshire population1

The GEBV accuracies of the testing as well as training data were calculated using the GBLUP and Bayes B methods (Table 3). Under the GBLUP model, the average (±standard deviation) of GEBV accuracy means across the traits was 0.65±0.04 for training data, which ranged from 0.42±0.08 for collagen to 0.75±0.02 for WHC. Under the Bayes B model, the GEBV accuracy mean ranged from 0.10±0.14 for NPCC marbling score to 0.76±0.04 for drip loss with the average of 0.49±0.10 across all the traits. For the testing data sets, the GEBV accuracy was lower, i.e. the overall average of the traits was 0.46±0.10 under the GBLUP model, ranging from 0.20±0.18 for protein to 0.65±0.04 for drip loss, and 0.38±0.13 under the Bayes B model, ranging from 0.04±0.09 for NPCC marbling score to 0.72±0.05 for drip loss (Table 3).

For all the traits, the GEBV accuracies were greater for the training samples than for the testing samples under both the GBLUP and Bayes B models. This makes sense in that, for the training samples, the GEBV prediction was based on both genotypes and phenotypes, while only the genotype information was exploited for the testing samples. However, the GEBV accuracy differences between the training and the testing samples were smaller for the traits with great training sample size. For example, the average GEBV accuracy of CWT (moisture) under the GBLUP model, with 1,051 (693) training individuals, was 0.60 (0.74) for the training and 0.56 (0.54) for the testing samples, respectively (Table 3).

The GEBV accuracy depends on four factors; 1) size of the training population, 2) the heritability of each trait, 3) the extend of linkage disequilibrium between the markers and the QTL, and 4) the distribution of QTL effects (Goddard, 2009; Hayes et al., 2009). Our results supported the first two factors. There was a general tendency of high GEBV accuracy of testing samples with the training sample size (Table 3) under both the GBLUP and Bayes B models (Figure 1). For example, the average GEBV accuracy of NPCC color score (drip loss) for the testing samples, with 358 (1,051) training samples, was 0.30 (0.64) under the GBLUP model. For the traits with the sample training sample size (e.g. CWT and drip loss, n = 1,051), the testing samples had greater GEBV accuracy for drip loss (average was 0.72 under the Bayes B model) than for CWT (0.54), for which heritability of the former (latter) trait was 0.27 (0.13) (Table 3).

Figure 1

Plot of the relationship between the size of training samples and GEBV accuracy of the testing samples under the GBLUP and Bayes B model. GEBV, genomic estimated breeding value; GBLUP, genome based linear unbiased prediction method.

In general, the average GEBV accuracy values were similar between the GBLUP and Bayes B methods. However, for some traits with small sample size, e.g. collagen, moisture and NPCC marbling score, the accuracy of both under the Bayes B model was much smaller than under the GBLUP model (Table 3). This may be partly due to small sample size of the training data, which would cause the estimation of GEBV to be more sensitive to the prior values of the Bayes B method. Cleveland et al. (2010) reported that GEBV accuracy decreased with small training sample size, because of not enough information to accurately estimate SNP effects.

There were a few reports about GEBV accuracy in pig populations. Badke et al. (2014) reported that the GEBV accuracy for BF was 0.45 to 0.47 in a Yorkshire population (965 training samples). Our results showed that, for the trait, 0.59 and 0.65 accuracies were obtained by the GBLUP and Bayes B analyses, respectively, using 1,043 training samples (Table 3).

There are many reports about successful genomic selection procedures in dairy cattle (Hayes et al., 2009; van Raden et al., 2009; Wiggans et al., 2011). The dairy cattle breeding industry is benefitting from genomic selection mainly by reduced generation intervals (Wellmann et al., 2013), while, in the pig industry, the benefit from genome selection was less, partly due to short generation intervals (Wellmann et al., 2013). However, some studies reported that genomic selection was relevant in pig breeding, by improving maternal traits (Simianer, 2009; Lillehammer et al., 2011) or by boosting selection intensities (Tribout et al., 2012; Wellmann et al., 2013).

CONCLUSION

Herein, we presented the first report about GEBV accuracy of a Berkshire population in Korea, and our results were in general agreement with the previous GEBV studies, i.e. GEBV accuracy depends on the size of training data as well as heritability of the tested trait. For some traits such as drip loss, GEBVs were predicted with a good accuracy under the GBLUP or Bayes B model (Table 3). However, more training samples are needed to further improve the GEBV accuracy for pork quality, especially with low heritable traits, to efficiently implement genome selection programs in Berkshire industry in Korea.

ACKNOWLEDGMENTS

This research was supported by a grant (PJ009032) from the Next Generation BioGreen 21 Program, Rural Development Administration, Republic of Korea.

References

American Berkshire Association. 2013. Information on Berkshires Available: http://www.americanberkshire.com. Accessed Nov 1, 2013.
AOAC. 2000. Official Methods of Analysis Gaithersburg, MD: Association of Official Analytical Chemists. 17th edth ed. Arlington, VA, USA:
Badke YM, Bates RO, Ernst CW, Fix J, Steibel JP. 2014;Accuracy of estimation of genomic breeding value in pigs using low density genotypes and imputation. G3 4:623–631.
Boddhireddy P, Kelly MJ, Northcutt S, Prayaga KC, Rumph J, Denise S. 2014;Genomic predictions in Angus cattle: comparisons of sample size, response variables, and clustering methods for cross-validation. J Anim Sci 92:485–497.
Bonneau M, Lebret B. 2010;Production systems and influence on eating quality of pork. Meat Sci 84:293–300.
Brewer MS, Jensen J, Sosnicki AA, Fields B, Wilson E, McKeith FK. 2002;The effect of pig genetics on palatability, color, and physical characteristics of fresh pork loin chops. Meat Sci 61:249–256.
Brorsen BW, Akridge JT, Boland MA, Mauney S, Forrest JC. 1998;Performance and alternative component pricing systems for pork. J Agric Appl Econ 30:313–324.
Browning SR, Browning BL. 2007;Rapid and accurate haplotype phasing and missing data inference for whole-genome association studies by use of localized haplotype clustering. Am J Hum Genet 81:1084–1097.
Christensen OF, Lund MS. 2010;Genomic prediction when some animals are not genotyped. Genet Sel Evol 46:2.
Cleveland MA, Forni S, Garrick DJ, Deeb N. 2010. Prediction of genomic breeding values in a commercial pig population. In : Proceedings of the 9th World Congress on Genetics Applied to Livestock Production. Paper 0266; August 1–6, 2010; Leipzig, Germany.
Daetwyler HD, Villanueva B, Woolliams JA. 2008;Accuracy of predicting the genetic risk of disease using a genome-wide approach. PLoS One 3(10):e3395.
Dekkers JCM, Mathur PK, Knol EF. 2010. Genetic improvement of the pig. The Genetics of the Pig In : Rothschild MF, Ruvinsky A, eds. chapter 162nd edth ed. CABI. Cambridge, MA, USA: p. 390–425.
Erbe M, Hayes BJ, Matukumalli LK, Goswami S, Bowman PJ, Reich CM, Mason BA, Goddard ME. 2012;Improving accuracy of genomic predictions within and between dairy cattle breeds with imputed high-density single nucleotide polymorphism panels. J Dairy Sci 95:4114–4129.
Etherington DE, Sims TJ. 1981;Detection and estimation of collage. J Sci Food Agric 32:539–546.
Fernando RL, Habier D, Stricker C, Dekkers JCM, Totir LR. 2007;Genomic selection. Acta Agric Scand A Anim Sci 57:192–195.
Folch J, Lee M, Sloane-Stanley GH. 1957;A simple method for the isolation and purification of total lipids from animal tissues. J Biol Chem 226:497–507.
Gilmour AR, Thompson R, Cullis BR. 1995;Average information REML; an efficient algorithm for variance parameter estimation in linear mixed models. Biometrics 51:1440–1450.
Goddard ME. 2009;Genomic selection: Prediction of accuracy and maximization of long term response. Genetica 136:245–257.
Goodwin R, Burroughs S. 1995. Genetic Evaluation Terminal Line Program Results National Pork Producers Council. Des Moines, IA, USA:
Habier D, Fernando RL, Kizilkaya K, Garrick DJ. 2011;Extension of the bayesian alphabet for genomic selection. BMC Bioinformatics 12:186.
Hayes BJ, Bowman PJ, Chamberlain AJ, Goddard ME. 2009;Invited review: Genomic selection in dairy cattle: Progresses and challenges. J Dairy Sci 92:433–443.
Jung J-H, Kim C-W, Park B-Y, Choi J-S, Park H-C. 2011;Genetic parameter estimates for meat quality traits in Berkshire pigs. J Anim Sci Technol 53:289–296.
Kawaida H. 1993;Studies on the performance of meat production and meat in pigs. Report of Kagoshima Prefectural Animal Husbandry Experiment Stations 26:1–195.
Lammers PJ, Mabry JW, Honeyman , Mark S, Swantek , Matthew P, Roush WB. 2011;Developing berkshire market pig growth curves. Iowa State Research Farm Progress Reports
Legarra A, Misztal I. 2008;Technical note: Computing strategies in genome-wide selection. J Dairy Sci 91:360–366.
Lillehammer M, Meuwissen THE, Sonesson AK. 2011;Genomic selection for maternal traits in pigs. J Anim Sci 89:3908–3916.
Luo W, Cheng D, Chen S, Wang L, Li Y, Ma X, Song X, Liu X, Li W, Liang J, Yan H, Zhao K, Wang C, Wang L, Zhang L. 2012;Genome-wide association analysis of meat quality traits in a porcine Large White×Minzhu intercross population. Int J Biol Sci 8:580–595.
McLaughlin K. 2004;Now, It’s the Other Red Meat; Atkins Craze Gives a Boost to Fattier, Tastier Pork; Tracing Your Chops’ Pedigree. Wall Street J :D1.
Meuwissen THE, Hayes BJ, Goddard ME. 2001;Prediction of total genetic value using genome-wide dense marker maps. Genetics 157:1819–1829.
Meuwissen THE. 2009;Accuracy of breeding values of ‘unrelated’ individuals predicted by dense SNP genotyping. Genet Sel Evol 41:35.
Moser G, Tier B, Crump RE, Khatkar MS, Raadsma HW. 2009;A comparison of five methods to predict genomic breeding values of dairy bulls from genome-wide SNP markers. Genet Sel Evol 41:56.
Newcom DW, Baas TJ, Mabry JW, Goodwin RN. 2002;Genetic parameters for pork carcass components. J Anim Sci 80:3099–3106.
NPPC. 2000. Pork Composition and Quality Assessment Procedures National Pork Producers Council. Des Moines, IA, USA:
Nothnagel M, Ellinghaus D, Schreiber S, Krawczak M, Franke A. 2009;A comprehensive evaluation of SNP genotype imputation. Hum Genet 125:163–171.
Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC. 2007;PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81:559–575.
Rothschild MF, Ruvinsky A. 2011. The Genetics of the Pig 2nd edth ed. CAB Int. NY, USA:
Simianer H. 2009. The potential of genomic selection to improve litter size in pig breeding programs. In : Proc 60th Annual meeting of the European Association of Animal Production. Barcelona, Spain. p. 210.
Tomiyama M, Kanetani T, Tatsukawa Y, Mori H, Oikawa T. 2011;Genetic relationships and expected responses for genetic improvement of carcass traits of Berkshire pigs. Sci. Agric. (Piracicaba, Braz.) 68:594–597.
Tribout T, Larzul C, Phocas F. 2012;Efficiency of genomic selection in a purebred pig male line. J Anim Sci 90:4164–4176.
Tsuruta S, Misztal I. 2008;Technical note: Computing options for genetic evaluation with a large number of genetic markers. J Anim Sci 86:1514–1518.
van Wijk HJ, Arts DJ, Matthews JO, Webster M, Ducro BJ, Knol EF. 2005;Genetic parameters for carcass composition and pork quality estimated in a commercial production chain. J Anim Sci 83:324–333.
van Raden PM, Van Tassell CP, Wiggans GR, Sonstegard TS, Schnabel RD, Taylor JF, Schenkel FS. 2009;Invited review: Reliability of genomic predictions for North American Holstein bulls. J Dairy Sci 92:16–24.
Wellmann R, Preuß S, Tholen E, Heinkel J, Wimmers K, Bennewitz J. 2013;Genomic selection using low density marker panels with application to a sire line in pigs. Genet Sel Evol 45:28.
Wiggans GR, VanRaden PM, Cooper TA. 2011;The genomic evaluation system in the United States: past, present, future. J Dairy Sci 94:3202–3211.

Article information Continued

Figure 1

Plot of the relationship between the size of training samples and GEBV accuracy of the testing samples under the GBLUP and Bayes B model. GEBV, genomic estimated breeding value; GBLUP, genome based linear unbiased prediction method.

Table 1

Summary statistics for 16 meat quality traits in a Berkshire population

Trait Average SD Min Max CV
Back fat thickness (mm) 25.1 5.2 12 46 20.8
CIE a 9.3 4.6 3.4 19.9 50.1
CIE b 3.7 1.5 0.32 8.3 41.3
CIE L 48.9 2.9 38 58.4 6.1
Collagen 0.89 0.13 0.53 1.39 14.9
Carcass weight (kg) 86.8 5.7 71 111 6.5
Drip loss (%) 40.0 2.0 0.1 14.4 50.9
Fat (%) 2.7 1.2 0.42 10.2 44.6
Heat loss (%) 23.4 6.3 6.8 39 27.1
Moisture (%) 75.4 1.0 69.9 77.6 1.4
NPCC color score 3.5 0.6 1.5 5 17.7
NPCC marbling score 2.2 0.7 1 5 30.6
pH24 5.7 0.2 5.1 6.7 3.6
Protein (%) 23.7 0.9 20.9 26.2 44.6
Shear force (kg) 2.9 0.8 1.3 6.1 26.7
Water holding capacity(%) 58.4 3.4 50.1 67.8 5.8

SD, standard deviation; Min, Minimum; Max, Maximum; CV, coefficient of variation (%), CIE, Commission Internationale de l’Eclairage; NPCC, National Pork Producers Council.

Table 2

The number of available SNPs and average distances between adjacent SNPs in the 18 Sus scrofa autosomes (SSC) of the Berkshire pig population

SSC Number of SNPs1 Average interval size (kb) Standard deviation (kb) Total distance2 (bp)
1 4,426 71.2 135.9 31,518,363
2 2,513 64.6 108.6 162,437,491
3 2,000 72.7 111.8 145,400,228
4 2,420 59.7 104.3 144,372,696
5 1,628 68.4 95.8 111,295,637
6 2,335 73.3 150.1 171,234,551
7 2,347 57.7 76.7 135,453,641
8 2,094 70.9 95 148,611,690
9 2,299 66.7 106.2 153,432,870
10 1,302 60.9 94.6 79,397,854
11 1,119 78.2 120.5 87,533,711
12 1,101 61.9 100.6 68,134,174
13 2,920 74.7 114.9 218,135,680
14 2,958 51.9 58.4 153,403,678
15 1,832 92.4 159.3 169,232,786
16 1,242 69.2 101.2 85,986,771
17 1,183 58.3 86.2 68,944,563
18 886 68.6 100.2 60,833,248
Total: 36,605 Average: 67.85 Average: 106.683 Total: 2,195,359,632

SNP, single nucleotide polymorphism; EBV, estimated breeding value.

1

Among the 62,163 SNPs in the Illumina Porcine 62 k beadchip, those SNPs were selected for genome EBV evaluation after quality control tests; any SNP was excluded with <90% call rates, <5% minor allele frequency, or significant departure from Hardy Weinberg equilibrium (p<0.001).

2

The distances between the first and the last SNPs that were located on their respective chromosomes.

Table 3

Accuracies of genomic estimated breeding value (GEBV) under the GBLUP and Bayes B models for the training and testing samples in the Berkshire population1

Trait Heritability2 Number of samples Number oftraining samples Training data4 Number oftesting samples Testing data3


GBLUP Bayes B GBLUP Bayes B
BF (mm) 0.34 1,191 1,043 0.68±0.03 0.72±0.04 148 0.59±0.05 0.65±0.06
CIE a 0.30 1,191 1,191 0.74±0.03 0.73±0.04 0
CIE b 0.30 1,191 1,191 0.75±0.03 0.74±0.04 0
CIE L 0.11 1,191 1,191 0.59±0.05 0.60±0.07 0
Collagen 0.06 1,191 693 0.42±0.08 0.13±0.15 498 0.25±0.17 0.13±0.15
CWT (kg) 0.13 1,191 1,051 0.60±0.04 0.59±0.07 140 0.56±0.05 0.54±0.08
Drip loss (%) 0.27 1,191 1,051 0.72±0.03 0.76±0.04 140 0.65±0.04 0.72±0.05
Fat (%) 0.46 1,191 686 0.73±0.02 0.74±0.04 505 0.49±0.11 0.57±0.11
Heat loss (%) 0.15 1,191 1,191 0.63±0.04 0.66±0.05 0
Moisture (%) 0.33 1,191 693 0.74±0.02 0.12±0.15 498 0.54±0.10 0.15±0.15
NPCC color score 0.23 1,191 358 0.64±0.03 0.44±0.33 833 0.30±0.17 0.38±0.33
NPCC marbling score 0.30 1,191 358 0.70±0.02 0.10±0.14 833 0.37±0.15 0.04±0.09
pH24 0.13 1,191 1,186 0.60±0.04 0.60±0.07 5 0.56±0.04 0.53±0.05
Protein (%) 0.27 1,191 686 0.51±0.05 0.28±0.17 505 0.20±0.18 0.25±0.17
Shear force (kg) 0.10 1,191 1,191 0.57±0.05 0.53±0.09 0
WHC (%) 0.34 1,191 833 0.75±0.02 0.16±0.16 358 0.53±0.09 0.17±0.16
Average 0.65±0.04 0.49±0.10 0.46±0.10 0.38±0.13

GBLUP, genome based linear unbiased prediction method; BF, backfat thickness; CIE, Commission Internationale de l’Eclairage; CWT, carcass weight; NPCC, National Pork Producers Council; WHC, water holding capacity.

1

The 16 traits have different training (with both genotypes and phenotypes) and testing (genotypes only) samples.

2

The heritabilities were estimated under the GBLUP model.

3

Mean±standard deviation of GEBV.