TY - JOUR AU - Bittante, G. AB - ABSTRACT The aims of this study were 1) to investigate the potential application of near-infrared spectroscopy (NIRS) to predict beef quality (BQ) traits, 2) to assess genetic variations of BQ measures and their predictions obtained by NIRS, and 3) to infer the genetic relationship between measures of BQ and their predictions. Young Piedmontese bulls (n = 1,230) were raised and fattened on 124 farms and slaughtered at the same commercial abattoir. The BQ traits evaluated were shear force (SF, kg), cooking loss (CL, %), drip loss (DL, %), lightness (L*), redness (a*), yellowness (b*), saturation index (SI), and hue angle. Near-infrared spectra were collected using a Foss NIRSystems 5000 instrument over a spectral range of 1,100 to 2,498 nm every 2 nm, in reflectance mode. After editing, prediction models were developed on a calibration subset (n = 268) using partial least squares regressions, followed by application of these models to the validation subset (n = 940). Estimations of (co)variance for measures of BQ and NIRS-based predictions were obtained through a set of bivariate Bayesian analyses on the validation subset. Near-infrared predictions were satisfactory for measurements of L* (R2 = 0.64), a* (R2 = 0.68), hue angle (R2 = 0.81), and saturation index (R2 = 0.59), but not for b*, DL, CL, and SF. The loss of additive genetic variance of predicted vs. measured L*, a*, DL, CL, and SF was generally high and was similar to the loss of residual variance, being a function of the calibration parameter R2. As a consequence, estimated heritabilities of measures and predictions of BQ were similar for traits with high calibration R2 values. Genetic correlations between BQ measures and predictions were high for all color traits and DL, and were greater than the corresponding phenotypic correlations, whereas both the phenotypic and genetic correlations for SF and CL were nil. Results suggest that NIRS-based predictions for color features and DL may be used as indicator traits to improve meat quality of the Piedmontese breed. INTRODUCTION Beef quality (BQ) traits are important aspects of beef production from the consumer viewpoint and with respect to related effects on herd profitability (Dikeman, 1990). Because exploitable additive genetic variation for BQ exists (Burrow et al., 2001; Johnston et al., 2003; Warner et al., 2010), it is expected that such traits may be enhanced through selection strategies. However, large-scale recording of individual BQ phenotypes remains a critical issue because the available techniques are time-consuming and, as yet, no high-throughput automated measuring device has been developed. Because opportunities for breeding based on direct measures of BQ phenotypes are limited, optimal selection to enhance BQ remains under investigation. Alternative methodologies for measuring BQ rely on the use of infrared spectroscopy (Prevolnik et al., 2004; Prieto et al., 2009). Because the technique is rapid and accurate, infrared spectroscopy has been successfully used in the beef industry for chemical analysis (Mitsumoto et al., 1991; Hildrum et al., 1995). Conversely, near-infrared spectroscopy (NIRS) is less reliable in predicting technological variables of meat and meat products (Abeni and Bergoglio, 2001; Geesink et al., 2003; Hoving-Bolink et al., 2005). From a genetic point of view, the relevance of NIRS in programs focusing on selection for improved BQ relies on genetic variation of NIRS-based predictions of BQ and on the magnitude of genetic correlation between NIRS-based predictions and BQ measured by defined reference methods. No estimates of such genetic parameters are currently available for BQ. We therefore investigated the ability of NIRS to predict BQ, the genetic variations in BQ measures and their predictions obtained by NIRS, and the genetic relationship between measures of BQ and predictions. MATERIALS AND METHODS Animal Care and Use Committee approval was not obtained for this study because the data were obtained from an existing database; the analyzed records were collected after slaughtering of animals in a commercial abattoir (Carrù, Piemonte, Italy) from March 2005 to February 2007. The authors did not have direct control over the care of the animals included in this study. Animals, Samples, and Data The present study was part of a wider project conducted on 1,230 young Piedmontese bulls marketed as “Vitellone Piemontese della Coscia” under a Protected Geographical Indication, as defined by the European Union. Animals were fattened on 124 farms and slaughtered at the same commercial abattoir from March 2005 to February 2007. The average age at slaughter (SD) was 523 d (73 d). Young bulls were progeny of 109 AI purebred sires and 1,170 dams, all registered in the Italian Piedmontese herdbook. Sire-offspring relationships were ascertained using DNA testing based on 19 microsatellites (Budowle et al., 2005). Twenty-four hours after slaughter, individual samples (1 per animal) of the longissimus thoracis muscle were collected between the fifth and sixth thoracic vertebrae, weighed, individually vacuum-packed, and transferred to the meat laboratory of the Department of Animal Science of the University of Padova (Legnaro, Italy) in a portable cooler (4°C). On arrival, samples were stored at 4°C in a chilling room until measurement of BQ traits. Physical Analysis and NIRS After aging for 8 d, the BQ of meat samples was assessed by measurement of color, drip loss (DL, %), cooking loss (CL, %), and shear force (SF, kg). Drip loss was computed as the difference between the weight of the packaged sample and the weight of the sample dried using blotting paper plus the weight of the heat-dried bag, and was expressed as a percentage of sample weight. Two slices of meat, each 20 mm thick, were cut from each sample. The first slice was exposed to air at 4°C for 1 h (ASPA, 1996), and the color of the exposed surface was determined using a Minolta CM-508 spectrophotometer (Konica Minolta, Milan, Italy) equipped with a D65 light source. The reflectance coordinates [L* (lightness), a* (redness), b* (yellowness)] of CIE (1976) were measured at 5 random positions. Saturation index (SI) was calculated using the formula whereas hue angle (HA) was calculated from the formula HA = tan−1(b*/a*) (AMSA, 1991). The variables L*, a*, b*, SI, and HA were measured 5 times each, and the measures were averaged before statistical analysis. The second slice was weighed, sealed in a polyethylene bag, heated in a water bath at 75°C for 55 min, and weighed again for measurement of CL (ASPA, 1996). Cooking loss was computed as the amount of weight lost after cooking as a percentage of the weight of the uncooked slice. Shear force was measured on 5 1-cm2 cross-sectional round cores, taken at approximately the same location from each cooked slice and running parallel to the longitudinal axis of the muscle fibers, by using a TA-HDi Texture Analyzer (Stable Micro System Ltd., Godalming, UK) equipped with a Warner-Bratzler shear device (100-kg load cell and crosshead speed of 2 mm/s). Both CL and SF were measured 5 times each, and measures were averaged before statistical analysis. Near-infrared spectra were collected on fresh minced samples using a Foss NIRSystems 5000 instrument (FOSS NIRSystems Inc., Laurel, MD) over a spectral range of 1,100 to 2,498 nm, in reflectance mode, every 2 nm. Duplicate spectra were captured for each sample and averaged before data analysis. Statistical Analysis Validation Procedure. The original data set was edited to discard records with errors (e.g., individual identification spectra not matching reference samples). After editing, a total of 1,208 records, including measures of BQ and NIRS spectra, were available for statistical analysis. To evaluate the predictive ability of NIRS measurements for BQ traits and the magnitude of genetic correlations between BQ measurements and predictions from the calibration equations based on NIRS spectra, a holdout validation procedure was carried out. This method involved the partition of data into 2 subsets, a calibration and a test subset. The former was used to develop a calibration equation that could predict individual BQ phenotypes by using NIRS spectra, whereas the latter was employed to validate the calibration equation and to estimate heritabilities and genetic correlations for the measured BQ traits and predictions obtained from NIRS spectra and the calibration equations. Observations to be included in the calibration subset were randomly sampled from the set of available data, with the restriction that at least 2 observations per herd and week of laboratory analysis were present in the subset (268 records). Records not included in the calibration subset were included in the test subset (940 records). Multivariate Data Analysis and Predictive Ability of NIRS. Partial least squares regressions (Unscrambler software, Camo A/S, Oslo, Norway) were used to establish calibration models (Hubert and Van den Branden, 2003). Partial least squares regressions have been used to estimate correlations between reference data and values predicted from spectral information in meat (Prieto et al., 2008) and milk (Soyeurt et al., 2006; Dal Zotto et al., 2008; De Marchi et al., 2009). Models were developed using untreated spectral data (absorbance spectra) for CL, DL, and color traits, and untreated spectral data plus first-derivative spectra for SF. To assess the adequacy of the calibration models, we calculated the root mean square error of calibration (RMSEC) and the coefficient of determination (R2). To evaluate the practical utility of the models, the range error ratio (RER) was calculated as the ratio between the range and the RMSEC of the trait (Williams, 1987). Models with RER values less than 3 have little practical utility; RER values between 3 and 10 indicate the models are limited to good practical utility; and RER values above 10 indicate a high utility value of the model (Williams, 1987). Estimates of (Co)Variance Components. For test subset (n = 940 records), (co)variance components and related parameters were estimated using a Bayesian approach and Markov-chain Monte Carlo methods (Sorensen and Gianola, 2002). A Bayesian technique was used because this offers some advantages over classical statistical methods (Blasco, 2005); in particular, Bayesian inference is based on probabilities, providing great flexibility in the construction of all types of confidence intervals. (Co)variance components for measures of BQ and predictions by NIRS were estimated using 8 separate bivariate Bayesian analyses on the test subset. For all traits, the model included the effects of fattening farms (124 levels), week of BQ laboratory analysis (92 wk), and carcass weight (class 1: <387 kg; class 2: 387 to 410 kg; class 3: 411 to 430 kg; class 4: 431 to 450 kg; class 5: 451 to 474 kg; class 6: >474 kg). All traits were continuous variables, and the values were assumed to be sampled from the following normal multivariate distribution:  where yBQ is a vector of phenotypic observations on measures of BQ; ypBQ is a vector of phenotypic observations on NIRS-based predictions of BQ; X, W1, W2, and Z are incidence matrices relating systematic fattening herd week of laboratory analysis and additive genetic effects to the vector of observations; I is an identity matrix of appropriate order; and R is a 2 × 2 matrix of residual (co)variances. The systematic effect considered in b is the carcass weight. Flat priors were used for systematic effect and dispersion parameters. Prior distributions for the additive genetic effects in u, herd effects in c, and weeks of laboratory analysis in q were normal densities. In a Bayesian setting, we assumed  where G is a 2 × 2 additive genetic (co)variance matrix and A is the numerator relationship matrix between individuals. Likewise, herd effects and week of laboratory analysis were assumed to follow normal bivariate distributions:    where P1 is a 2 × 2 herd (co)variance matrix and P2 is a 2 × 2 week of the laboratory analysis (co)variance matrix. Gibbs Sampler. Marginal posterior distributions of unknown parameters were estimated by performing numerical integration using the Gibbs sampler (Gelfand and Smith, 1990), as implemented in the TM program (http://snp.toulouse.inra.fr/~alegarra/). This was employed to obtain autocorrelated samples from the joint posterior distributions and subsequently from the marginal posterior distributions of all unknowns in the model. The lengths of the chain and of the burn-in period were assessed by visual inspection of trace plots, as well as by the diagnostic tests of Geweke (1992) and Gelman and Rubin (1992). After a preliminary run, it was decided to construct a single chain consisting of 850,000 iterations and to discard the first 50,000 iterations as a very conservative burn-in. Subsequently, 1 of every 200 successive samples was retained, to store draws that were more loosely correlated. Thus, 4,000 samples were used to determine posterior distributions of unknown parameters. The lower and upper bounds of the highest 95% probability density regions for h2, and additive genetic and residual variances were obtained from the estimated marginal densities as well as from a posterior probability for h2 > 0.10. The posterior median was used as the point for estimating (co)variance components and related parameters. Autocorrelations between samples and estimates of Monte Carlo SE (Geyer, 1992) were calculated. To facilitate comparisons with previously reported results, we calculated intraherd heritability as  where is the additive genetic variance and is the residual variance. Genetic correlations were computed as   where is the additive genetic covariance between traits 1 and 2, and and are the additive genetic SD of traits 1 and 2, respectively. Phenotypic correlations were computed as   where is the phenotypic covariance between traits 1 and 2, and and are the phenotypic SD of traits 1 and 2, respectively. RESULTS AND DISCUSSION The Piedmontese (Piemontese), the most important Italian beef cattle breed, is characterized by muscular hypertrophy (Kambadur et al., 1997) and is actively selected for traits such as daily BW gain, live fleshiness (Albera et al., 2001; Mantovani et al., 2010), and direct and maternal calving performance (Carnier et al., 2000; Kizilkaya et al., 2002, 2003). Previous investigations have assessed phenotypic variation in carcasses (Biagini and Lazzaroni, 2005; Dal Maso et al., 2009; Schiavon et al., 2010) and BQ traits (Destefanis et al., 2003; Schiavon et al., 2011) of purebred Piedmontese cattle. Conversely, genetic aspects of BQ traits have been studied only in crossbreeding experiments (Wheeler et al., 2001b), thus not at a population level. To address this issue, a research project has been established in cooperation with the National Association of Piedmontese Breeders (ANABoRaPi, Carrù, Italy) to study the genetics of BQ traits in this breed, both quantitatively (Boukha et al., 2007) and at a molecular level (Ribeca et al., 2009). The possibility of estimating BQ traits using NIRS forms the topic of the present report. The data obtained from BQ analyses (Table 1) were similar to those obtained on other hypertrophic breeds (Vincenti et al., 2007), especially Belgian Blue cattle (Fiems et al., 2003). Table 1. Descriptive statistics of beef quality (BQ) traits for the calibration set (n = 268) and test set (n = 940)1 Item2  Calibration set3  Test set  Mean  CV, %  Minimum  Maximum  R2  RMSEC  RER  Mean  CV, %  Minimum  Maximum  L*  34.28  10.33  27.36  44.54  0.64  2.12  8.10  34.67  9.76  21.08  37.46  a*  16.46  15.38  9.45  22.69  0.68  1.41  9.39  16.39  15.77  12.96  27.71  b*  14.85  14.89  9.39  20.92  0.44  1.65  6.99  15.09  14.60  8.85  20.68  HA  42.21  9.26  32.61  53.19  0.81  1.69  12.18  42.63  8.60  29.79  54.19  SI  22.28  13.64  15.04  31.03  0.59  1.94  8.24  22.36  13.94  14.51  30.94  DL, %  3.98  33.23  1.46  8.05  0.17  1.19  5.53  4.28  32.58  1.37  9.04  CL, %  24.33  14.94  13.72  32.01  0.04  3.55  5.15  24.22  14.38  12.65  33.22  SF, kg  2.65  20.55  1.57  4.41  0.21  0.48  5.92  2.68  20.84  1.51  4.48  Item2  Calibration set3  Test set  Mean  CV, %  Minimum  Maximum  R2  RMSEC  RER  Mean  CV, %  Minimum  Maximum  L*  34.28  10.33  27.36  44.54  0.64  2.12  8.10  34.67  9.76  21.08  37.46  a*  16.46  15.38  9.45  22.69  0.68  1.41  9.39  16.39  15.77  12.96  27.71  b*  14.85  14.89  9.39  20.92  0.44  1.65  6.99  15.09  14.60  8.85  20.68  HA  42.21  9.26  32.61  53.19  0.81  1.69  12.18  42.63  8.60  29.79  54.19  SI  22.28  13.64  15.04  31.03  0.59  1.94  8.24  22.36  13.94  14.51  30.94  DL, %  3.98  33.23  1.46  8.05  0.17  1.19  5.53  4.28  32.58  1.37  9.04  CL, %  24.33  14.94  13.72  32.01  0.04  3.55  5.15  24.22  14.38  12.65  33.22  SF, kg  2.65  20.55  1.57  4.41  0.21  0.48  5.92  2.68  20.84  1.51  4.48  1Calibration set = samples used to develop a calibration equation to predict BQ phenotypes using near-infrared spectra; test set = independent sample used to validate the calibration equation and to estimate heritabilities and genetic correlations for measured BQ and their predictions obtained from near-infrared spectra and the calibration equation. 2L* = lightness; a* = redness; b* = yellowness; HA = hue angle [HA = tan−1(b*/a*)]; SI = saturation index DL = drip loss; CL = cooking loss; SF = shear force. 3R2 = coefficient of determination of calibration; RMSEC = root mean square error of calibration; RER = ratio performance deviation. View Large Table 1. Descriptive statistics of beef quality (BQ) traits for the calibration set (n = 268) and test set (n = 940)1 Item2  Calibration set3  Test set  Mean  CV, %  Minimum  Maximum  R2  RMSEC  RER  Mean  CV, %  Minimum  Maximum  L*  34.28  10.33  27.36  44.54  0.64  2.12  8.10  34.67  9.76  21.08  37.46  a*  16.46  15.38  9.45  22.69  0.68  1.41  9.39  16.39  15.77  12.96  27.71  b*  14.85  14.89  9.39  20.92  0.44  1.65  6.99  15.09  14.60  8.85  20.68  HA  42.21  9.26  32.61  53.19  0.81  1.69  12.18  42.63  8.60  29.79  54.19  SI  22.28  13.64  15.04  31.03  0.59  1.94  8.24  22.36  13.94  14.51  30.94  DL, %  3.98  33.23  1.46  8.05  0.17  1.19  5.53  4.28  32.58  1.37  9.04  CL, %  24.33  14.94  13.72  32.01  0.04  3.55  5.15  24.22  14.38  12.65  33.22  SF, kg  2.65  20.55  1.57  4.41  0.21  0.48  5.92  2.68  20.84  1.51  4.48  Item2  Calibration set3  Test set  Mean  CV, %  Minimum  Maximum  R2  RMSEC  RER  Mean  CV, %  Minimum  Maximum  L*  34.28  10.33  27.36  44.54  0.64  2.12  8.10  34.67  9.76  21.08  37.46  a*  16.46  15.38  9.45  22.69  0.68  1.41  9.39  16.39  15.77  12.96  27.71  b*  14.85  14.89  9.39  20.92  0.44  1.65  6.99  15.09  14.60  8.85  20.68  HA  42.21  9.26  32.61  53.19  0.81  1.69  12.18  42.63  8.60  29.79  54.19  SI  22.28  13.64  15.04  31.03  0.59  1.94  8.24  22.36  13.94  14.51  30.94  DL, %  3.98  33.23  1.46  8.05  0.17  1.19  5.53  4.28  32.58  1.37  9.04  CL, %  24.33  14.94  13.72  32.01  0.04  3.55  5.15  24.22  14.38  12.65  33.22  SF, kg  2.65  20.55  1.57  4.41  0.21  0.48  5.92  2.68  20.84  1.51  4.48  1Calibration set = samples used to develop a calibration equation to predict BQ phenotypes using near-infrared spectra; test set = independent sample used to validate the calibration equation and to estimate heritabilities and genetic correlations for measured BQ and their predictions obtained from near-infrared spectra and the calibration equation. 2L* = lightness; a* = redness; b* = yellowness; HA = hue angle [HA = tan−1(b*/a*)]; SI = saturation index DL = drip loss; CL = cooking loss; SF = shear force. 3R2 = coefficient of determination of calibration; RMSEC = root mean square error of calibration; RER = ratio performance deviation. View Large Predictions of Beef Quality Traits Descriptive statistics for the calibration and test subsets are presented in Table 1. The means and CV of BQ traits in the 2 subsets were similar. The CV were quite large and ranged from 9.26 (HA) to 33.23% (DL), facilitating the development of robust calibration equations. The coefficients of determination of calibration (R2) varied from 0.44 to 0.81 for color indices, but were very small (R2 < 0.25) for DL, CL, and SF (Table 1). The RMSEC varied from 1.41 to 2.12 for color indexes and were 1.19%, 3.55%, and 0.48 kg for DL, CL, and SF, respectively. The RER values varied from 5.15 (CL) to 12.18 (HA), showing limited to good practical utility of the prediction models. Except for b*, the R2 values of all color indices were quite large (0.59 to 0.81), indicating that NIRS could predict these traits satisfactorily (Williams, 2003). The NIRS-based predictions were much less accurate for DL, CL, and SF (R2 < 0.25). In general, predictions from the present work were less accurate than those reported by Leroy et al. (2003), Andrés et al. (2008), and Prieto et al. (2008, 2009), but were more accurate than the values of Hoving-Bolink et al. (2005). The low accuracy of NIRS-based predictions for DL and CL is in agreement with data reported previously (Leroy et al., 2003; Meulemans et al., 2003; Andrés et al., 2008). The poor accuracy of such reference methods may be attributable to the heterogeneity of meat samples, the variability in spectrum wavelengths, the techniques used for the acquisition of spectra, or all 3. Such variation can adversely affect the results achievable when these traits are studied (Prieto et al., 2009). Good predictions for DL and CL have been reported in only 1 study on pork (R2 = 0.71; Forrest et al., 2000) and in only 1 study on beef (R2 = 0.86; Ripoll et al., 2008). We also found that NIRS was limited in its ability to predict SF (R2 = 0.21), confirming previous findings in beef (Leroy et al., 2003; Andrés et al., 2008), pork (Chan et al., 2002), and poultry (Liu et al., 2004). This may be attributable to the use of homogenized samples, a suboptimal spectral wavelength, limited variation in the available reference data for SF, or all 3. Variance Components and Heritability Point estimates (medians of the marginal posterior density of each parameter) for the additive genetic and residual SD of BQ traits measured by reference methods and predicted by NIRS are shown in Table 2. The additive genetic and residual SD were always less for the predicted than the measured BQ traits, even though the difference varied greatly according to the specific trait. Table 2. Posterior median (SD) for additive genetic (σA) and residual (σE) SD of beef quality (BQ) traits measured and predicted by near-infrared spectroscopy (NIRS) Item1  σA  σE  Measures  Predictions  ΔσA,2 %  P3  Measures  Predictions  ΔσE,4 %  P3  L*  1.66 (0.28)  1.25 (0.24)  −25  93.7  2.48 (0.17)  2.09 (0.14)  −16  98.6  a*  1.17 (0.20)  0.94 (0.14)  −20  94.0  1.70 (0.13)  1.24 (0.09)  −27  99.7  b*  0.64 (0.18)  0.60 (0.12)  −6  60.5  1.63 (0.07)  0.94 (0.07)  −42  100  HA  2.37 (0.21)  2.16 (0.21)  −9  93.7  1.81 (0.24)  1.69 (0.24)  −7  73.2  SI  0.95 (0.24)  0.92 (0.20)  −3  55.9  2.26 (0.11)  1.68 (0.10)  −26  99.9  DL, %  0.61 (0.15)  0.16 (0.04)  −74  99.6  1.07 (0.07)  0.40 (0.02)  −62  100  CL, %  0.65 (0.30)  0.18 (0.05)  −72  95.3  2.70 (0.09)  0.40 (0.02)  −85  100  SF, kg  0.15 (0.05)  0.06 (0.02)  −60  92.1  0.46 (0.02)  0.18 (0.01)  −61  100  Item1  σA  σE  Measures  Predictions  ΔσA,2 %  P3  Measures  Predictions  ΔσE,4 %  P3  L*  1.66 (0.28)  1.25 (0.24)  −25  93.7  2.48 (0.17)  2.09 (0.14)  −16  98.6  a*  1.17 (0.20)  0.94 (0.14)  −20  94.0  1.70 (0.13)  1.24 (0.09)  −27  99.7  b*  0.64 (0.18)  0.60 (0.12)  −6  60.5  1.63 (0.07)  0.94 (0.07)  −42  100  HA  2.37 (0.21)  2.16 (0.21)  −9  93.7  1.81 (0.24)  1.69 (0.24)  −7  73.2  SI  0.95 (0.24)  0.92 (0.20)  −3  55.9  2.26 (0.11)  1.68 (0.10)  −26  99.9  DL, %  0.61 (0.15)  0.16 (0.04)  −74  99.6  1.07 (0.07)  0.40 (0.02)  −62  100  CL, %  0.65 (0.30)  0.18 (0.05)  −72  95.3  2.70 (0.09)  0.40 (0.02)  −85  100  SF, kg  0.15 (0.05)  0.06 (0.02)  −60  92.1  0.46 (0.02)  0.18 (0.01)  −61  100  1L* = lightness; a* = redness; b* = yellowness; HA = hue angle [HA = tan−1(b*/a*)]; SI = saturation index DL = drip loss; CL = cooking loss; SF = shear force. 2Median of the marginal posterior density of the difference between the genetic SD of BQ measured and predicted by NIRS. 3P = posterior probability for values of the difference less than zero. 4Median of the marginal posterior density of the difference between the residual SD of BQ measured and predicted by NIRS. View Large Table 2. Posterior median (SD) for additive genetic (σA) and residual (σE) SD of beef quality (BQ) traits measured and predicted by near-infrared spectroscopy (NIRS) Item1  σA  σE  Measures  Predictions  ΔσA,2 %  P3  Measures  Predictions  ΔσE,4 %  P3  L*  1.66 (0.28)  1.25 (0.24)  −25  93.7  2.48 (0.17)  2.09 (0.14)  −16  98.6  a*  1.17 (0.20)  0.94 (0.14)  −20  94.0  1.70 (0.13)  1.24 (0.09)  −27  99.7  b*  0.64 (0.18)  0.60 (0.12)  −6  60.5  1.63 (0.07)  0.94 (0.07)  −42  100  HA  2.37 (0.21)  2.16 (0.21)  −9  93.7  1.81 (0.24)  1.69 (0.24)  −7  73.2  SI  0.95 (0.24)  0.92 (0.20)  −3  55.9  2.26 (0.11)  1.68 (0.10)  −26  99.9  DL, %  0.61 (0.15)  0.16 (0.04)  −74  99.6  1.07 (0.07)  0.40 (0.02)  −62  100  CL, %  0.65 (0.30)  0.18 (0.05)  −72  95.3  2.70 (0.09)  0.40 (0.02)  −85  100  SF, kg  0.15 (0.05)  0.06 (0.02)  −60  92.1  0.46 (0.02)  0.18 (0.01)  −61  100  Item1  σA  σE  Measures  Predictions  ΔσA,2 %  P3  Measures  Predictions  ΔσE,4 %  P3  L*  1.66 (0.28)  1.25 (0.24)  −25  93.7  2.48 (0.17)  2.09 (0.14)  −16  98.6  a*  1.17 (0.20)  0.94 (0.14)  −20  94.0  1.70 (0.13)  1.24 (0.09)  −27  99.7  b*  0.64 (0.18)  0.60 (0.12)  −6  60.5  1.63 (0.07)  0.94 (0.07)  −42  100  HA  2.37 (0.21)  2.16 (0.21)  −9  93.7  1.81 (0.24)  1.69 (0.24)  −7  73.2  SI  0.95 (0.24)  0.92 (0.20)  −3  55.9  2.26 (0.11)  1.68 (0.10)  −26  99.9  DL, %  0.61 (0.15)  0.16 (0.04)  −74  99.6  1.07 (0.07)  0.40 (0.02)  −62  100  CL, %  0.65 (0.30)  0.18 (0.05)  −72  95.3  2.70 (0.09)  0.40 (0.02)  −85  100  SF, kg  0.15 (0.05)  0.06 (0.02)  −60  92.1  0.46 (0.02)  0.18 (0.01)  −61  100  1L* = lightness; a* = redness; b* = yellowness; HA = hue angle [HA = tan−1(b*/a*)]; SI = saturation index DL = drip loss; CL = cooking loss; SF = shear force. 2Median of the marginal posterior density of the difference between the genetic SD of BQ measured and predicted by NIRS. 3P = posterior probability for values of the difference less than zero. 4Median of the marginal posterior density of the difference between the residual SD of BQ measured and predicted by NIRS. View Large The median of the marginal posterior density of the difference (loss of variability) between the genetic SD for L* measured and predicted by NIRS (ΔσA) was −25%, with a posterior probability of being less than zero (P) close to 100%. A similar result was observed for a*, with a ΔσA of −20%, whereas reductions in genetic variability were much less pronounced for b*, HA, and SI (ΔσA < 10%). The loss of genetic variability caused by NIRS estimation was ≥60% for SF, DL, and CL. The loss of residual variation (ΔσE) was generally similar to that of ΔσA, with the exception of b* and SI, which showed greater phenotypic losses. A study on milk coagulation properties measured using a reference method and predicted by infrared spectroscopy showed that both additive genetic and residual SD were less when estimated from infrared-based predictions, with average ΔσA and ΔσE values being −14 and −27%, respectively, for rennet coagulation time, and −6 and −37%, respectively, for curd firmness (Cecchinato et al., 2009). These findings are similar to those observed for color traits in the present work. Considering all traits together, the losses of phenotypic, additive genetic, and residual SD caused by NIRS prediction are strictly associated with the R2 values of the calibration equations (Figure 1). Figure 1. View largeDownload slide Relationships between losses of a) phenotypic, b) additive genetic, and c) residual SD of beef quality (BQ) traits (shear force, cooking loss, drip loss, lightness, redness, yellowness, hue angle [HA = tan−1(b*/a*)], saturation index) predicted by near-infrared spectroscopy and the coefficient of determination (R2) of each calibration equation. Figure 1. View largeDownload slide Relationships between losses of a) phenotypic, b) additive genetic, and c) residual SD of beef quality (BQ) traits (shear force, cooking loss, drip loss, lightness, redness, yellowness, hue angle [HA = tan−1(b*/a*)], saturation index) predicted by near-infrared spectroscopy and the coefficient of determination (R2) of each calibration equation. Features of the marginal posterior distributions of heritability of BQ measured and predicted by NIRS are shown in Table 3. Point estimates of heritability for BQ measures differed principally when color-related traits were examined, ranging from 0.13 to 0.63. The heritabilities of L* and a* were moderate and were in agreement with previous estimates (Aass, 1996; Johnston et al., 2003). Saturation index and b* showed the smallest estimates of heritability (0.15 and 0.13, respectively), whereas HA was the most heritable trait (0.63). Point estimates of heritability were low for CL (0.05) and SF (0.10) and were moderate to low for DL (0.24). Only a few studies have investigated genetic aspects of water-holding capacity, and the results have been contradictory, with heritability estimates for CL of 0.25 (Riley et al., 2003), 0.15 (Johnston et al., 2003), and 0.01 (Cecchi et al., 2004). We found that the heritability estimate for SF was 0.10, similar to that reported by Riley et al. (2003) but less than the values obtained by Splan et al. (1998), Wheeler et al. (2001a), and Dikeman et al. (2005). Table 3. Features of the marginal posterior distributions of heritability of beef quality (BQ) traits measured and predicted by near-infrared spectroscopy (NIRS)1 Item2  h2 measures  h2 predictions  Median  HPD95  P (h2 > 0.10)  Median  HPD95  P (h2 > 0.10)  L*  0.31  0.06; 0.47  99.9  0.26  0.05; 0.34  99.6  a*  0.32  0.14; 0.57  99.9  0.36  0.17; 0.43  100  b*  0.13  0.03; 0.30  67.7  0.29  0.10; 0.52  97.8  HA  0.63  0.45; 0.85  100  0.62  0.43; 0.86  100  SI  0.15  0.04; 0.33  81.8  0.23  0.08; 0.45  95.7  DL, %  0.24  0.05; 0.47  93.8  0.14  0.05; 0.32  71.3  CL, %  0.05  0.01; 0.20  24.2  0.17  0.05; 0.40  83.7  SF, kg  0.10  0.01; 0.28  49.0  0.10  0.02; 0.24  47.2  Item2  h2 measures  h2 predictions  Median  HPD95  P (h2 > 0.10)  Median  HPD95  P (h2 > 0.10)  L*  0.31  0.06; 0.47  99.9  0.26  0.05; 0.34  99.6  a*  0.32  0.14; 0.57  99.9  0.36  0.17; 0.43  100  b*  0.13  0.03; 0.30  67.7  0.29  0.10; 0.52  97.8  HA  0.63  0.45; 0.85  100  0.62  0.43; 0.86  100  SI  0.15  0.04; 0.33  81.8  0.23  0.08; 0.45  95.7  DL, %  0.24  0.05; 0.47  93.8  0.14  0.05; 0.32  71.3  CL, %  0.05  0.01; 0.20  24.2  0.17  0.05; 0.40  83.7  SF, kg  0.10  0.01; 0.28  49.0  0.10  0.02; 0.24  47.2  1Median = median of the marginal posterior density of the parameter; HPD95 = highest posterior density region at 95%; P = posterior probability for values of h2 greater than 0.10. 2L* = lightness; a* = redness; b* = yellowness; HA = hue angle [HA = tan−1(b*/a*)]; SI = saturation index DL = drip loss; CL = cooking loss; SF = shear force. View Large Table 3. Features of the marginal posterior distributions of heritability of beef quality (BQ) traits measured and predicted by near-infrared spectroscopy (NIRS)1 Item2  h2 measures  h2 predictions  Median  HPD95  P (h2 > 0.10)  Median  HPD95  P (h2 > 0.10)  L*  0.31  0.06; 0.47  99.9  0.26  0.05; 0.34  99.6  a*  0.32  0.14; 0.57  99.9  0.36  0.17; 0.43  100  b*  0.13  0.03; 0.30  67.7  0.29  0.10; 0.52  97.8  HA  0.63  0.45; 0.85  100  0.62  0.43; 0.86  100  SI  0.15  0.04; 0.33  81.8  0.23  0.08; 0.45  95.7  DL, %  0.24  0.05; 0.47  93.8  0.14  0.05; 0.32  71.3  CL, %  0.05  0.01; 0.20  24.2  0.17  0.05; 0.40  83.7  SF, kg  0.10  0.01; 0.28  49.0  0.10  0.02; 0.24  47.2  Item2  h2 measures  h2 predictions  Median  HPD95  P (h2 > 0.10)  Median  HPD95  P (h2 > 0.10)  L*  0.31  0.06; 0.47  99.9  0.26  0.05; 0.34  99.6  a*  0.32  0.14; 0.57  99.9  0.36  0.17; 0.43  100  b*  0.13  0.03; 0.30  67.7  0.29  0.10; 0.52  97.8  HA  0.63  0.45; 0.85  100  0.62  0.43; 0.86  100  SI  0.15  0.04; 0.33  81.8  0.23  0.08; 0.45  95.7  DL, %  0.24  0.05; 0.47  93.8  0.14  0.05; 0.32  71.3  CL, %  0.05  0.01; 0.20  24.2  0.17  0.05; 0.40  83.7  SF, kg  0.10  0.01; 0.28  49.0  0.10  0.02; 0.24  47.2  1Median = median of the marginal posterior density of the parameter; HPD95 = highest posterior density region at 95%; P = posterior probability for values of h2 greater than 0.10. 2L* = lightness; a* = redness; b* = yellowness; HA = hue angle [HA = tan−1(b*/a*)]; SI = saturation index DL = drip loss; CL = cooking loss; SF = shear force. View Large Heritability estimates for BQ predicted by NIRS were similar to or larger than estimates for measured BQ, except for L* and DL. Considering all BQ traits, heritabilities of measured and predicted traits showed a correlation of 0.78. When we compared the statistics in Table 1 with the estimates in Table 3, we found that neither the absolute value of heritability of the predicted BQ traits nor the difference between heritabilities of the predicted and measured values was directly linked to the R2 value of the NIRS calibration. The largest variations between values of heritability estimated using measured and predicted traits were observed for b* (0.13 to 0.29), CL (0.05 to 0.17), DL (0.24 to 0.14), and SI (0.15 to 0.23). These traits also exhibited the largest differences between ΔσA and ΔσE (Table 2). Correlations Between Measured and Predicted BQ Traits Point estimates (posterior medians) and lower and upper bounds of the 95% highest posterior density intervals for additive genetic and phenotypic correlations between measured and predicted BQ traits are shown in Table 4. The estimates were much greater for genetic than for phenotypic correlations, with the exception of CL. The genetic correlation for DL was fairly high (0.72), although the phenotypic correlation (0.24) and R2 of the calibration equation (0.17) were low. For all color traits, the estimates between predicted and measured traits were very high and ranged from 0.85 (L*) to 0.99 (HA). The estimated posterior densities of the genetic correlations between measures and predictions of BQ were skewed (data not shown), and were similar among all color traits. The estimated symmetric 95% posterior density regions for color aspects showed that, in the least favorable situation (L*), the genetic correlation between measured and predicted values had a 97.5% posterior probability of being >0.53 (data not shown). The lower boundary of the highest 95% posterior density interval was much less for DL than for color values, but the probability of the genetic correlation being greater than zero was 98% (data not shown). Table 4. Posterior median and bounds of the 95% high posterior density region (HPD95) for the genetic (rA) and phenotypic (rp) correlations between beef quality (BQ) traits measured and predicted by near-infrared spectroscopy (NIRS)1 Correlation2  rA  rp  Median  HPD95  Median  HPD95  L* with pL*  0.85  0.53; 0.98  0.60  0.55; 0.65  a* with pa*  0.98  0.83; 0.99  0.67  0.62; 0.71  b* with pb*  0.93  0.44; 0.99  0.47  0.38; 0.54  HA with pHA  0.99  0.96; 0.99  0.76  0.72; 0.79  SI with pSI  0.95  0.61; 0.99  0.61  0.55; 0.66  DL with pDL  0.72  −0.07; 0.98  0.24  0.17; 0.32  CL with pCL  −0.04  −0.83; 0.84  0.04  −0.03; 0.11  SF with pSF  −0.10  −0.85; 0.86  0.02  −0.05; 0.09  Correlation2  rA  rp  Median  HPD95  Median  HPD95  L* with pL*  0.85  0.53; 0.98  0.60  0.55; 0.65  a* with pa*  0.98  0.83; 0.99  0.67  0.62; 0.71  b* with pb*  0.93  0.44; 0.99  0.47  0.38; 0.54  HA with pHA  0.99  0.96; 0.99  0.76  0.72; 0.79  SI with pSI  0.95  0.61; 0.99  0.61  0.55; 0.66  DL with pDL  0.72  −0.07; 0.98  0.24  0.17; 0.32  CL with pCL  −0.04  −0.83; 0.84  0.04  −0.03; 0.11  SF with pSF  −0.10  −0.85; 0.86  0.02  −0.05; 0.09  1Median = median of the marginal posterior density of the parameter; HPD95 = highest posterior density region at 95%. 2L* = lightness; pL* = lightness predicted by NIRS; a* = redness; pa* = redness predicted by NIRS; b* = yellowness; pb* = yellowness predicted by NIRS; HA = hue angle [HA = tan−1(b*/a*)]; pHA = HA predicted by NIRS; SI = saturation index pSI = SI predicted by NIRS; DL = drip loss; pDL = drip loss predicted by NIRS; CL = cooking loss; pCL = cooking loss predicted by NIRS; SF = shear force; pSF = shear force predicted by NIRS. View Large Table 4. Posterior median and bounds of the 95% high posterior density region (HPD95) for the genetic (rA) and phenotypic (rp) correlations between beef quality (BQ) traits measured and predicted by near-infrared spectroscopy (NIRS)1 Correlation2  rA  rp  Median  HPD95  Median  HPD95  L* with pL*  0.85  0.53; 0.98  0.60  0.55; 0.65  a* with pa*  0.98  0.83; 0.99  0.67  0.62; 0.71  b* with pb*  0.93  0.44; 0.99  0.47  0.38; 0.54  HA with pHA  0.99  0.96; 0.99  0.76  0.72; 0.79  SI with pSI  0.95  0.61; 0.99  0.61  0.55; 0.66  DL with pDL  0.72  −0.07; 0.98  0.24  0.17; 0.32  CL with pCL  −0.04  −0.83; 0.84  0.04  −0.03; 0.11  SF with pSF  −0.10  −0.85; 0.86  0.02  −0.05; 0.09  Correlation2  rA  rp  Median  HPD95  Median  HPD95  L* with pL*  0.85  0.53; 0.98  0.60  0.55; 0.65  a* with pa*  0.98  0.83; 0.99  0.67  0.62; 0.71  b* with pb*  0.93  0.44; 0.99  0.47  0.38; 0.54  HA with pHA  0.99  0.96; 0.99  0.76  0.72; 0.79  SI with pSI  0.95  0.61; 0.99  0.61  0.55; 0.66  DL with pDL  0.72  −0.07; 0.98  0.24  0.17; 0.32  CL with pCL  −0.04  −0.83; 0.84  0.04  −0.03; 0.11  SF with pSF  −0.10  −0.85; 0.86  0.02  −0.05; 0.09  1Median = median of the marginal posterior density of the parameter; HPD95 = highest posterior density region at 95%. 2L* = lightness; pL* = lightness predicted by NIRS; a* = redness; pa* = redness predicted by NIRS; b* = yellowness; pb* = yellowness predicted by NIRS; HA = hue angle [HA = tan−1(b*/a*)]; pHA = HA predicted by NIRS; SI = saturation index pSI = SI predicted by NIRS; DL = drip loss; pDL = drip loss predicted by NIRS; CL = cooking loss; pCL = cooking loss predicted by NIRS; SF = shear force; pSF = shear force predicted by NIRS. View Large The additive genetic correlation between predicted and measured BQ traits was the most important criterion for using NIRS-based predictions as indicator traits in breeding programs (Cecchinato et al., 2009; Rutten et al., 2010). Considering all the BQ traits, the additive genetic correlations exhibited a very low relationship with the heritabilities of predicted traits (Figure 2). A comparison of CL and DL is particularly significant, in that these parameters showed similar heritabilities of predicted values (0.17 and 0.14, respectively) but completely different additive genetic correlations between predicted and measured values (−0.04 and 0.72, respectively). Figure 2. View largeDownload slide Relationships of genetic correlations between measured and predicted beef quality (BQ) traits (shear force, cooking loss, drip loss, lightness, redness, yellowness, hue angle [HA = tan−1(b*/a*)], saturation index) with a) heritabilities of the predicted traits and b) the coefficient of determination (R2) of each calibration equation. Figure 2. View largeDownload slide Relationships of genetic correlations between measured and predicted beef quality (BQ) traits (shear force, cooking loss, drip loss, lightness, redness, yellowness, hue angle [HA = tan−1(b*/a*)], saturation index) with a) heritabilities of the predicted traits and b) the coefficient of determination (R2) of each calibration equation. The additive genetic correlations seemed to be positively correlated with the R2 values of the calibration equations (R2 = 0.67), but a comparison between CL and DL clarifies the range of possible variations. The R2 value of the calibration equation was much more directly correlated with the phenotypic correlation between predicted and measured BQ values (R2 = 0.94) than with the genetic correlation. To our knowledge, this study is the first one dealing with the use of NIRS to assess indicator traits in a breeding program for enhanced BQ. Our findings indicate that, although the R2 values of some calibration equations were low, the additive genetic correlations between predicted and measured values of BQ remained quite high. Because the use of NIRS calibration equations resulted in losses of additive genetic and phenotypic variability in an unpredictable manner, the heritability estimates of predicted values may be greater or less than those of the measured values. Thus, the heritability of predictions alone cannot be considered a good indicator of the suitability of a calibration equation. A combination of the R2 values of calibration, the heritability of predicted values, and the loss of phenotypic variability is required to determine the utility of the technique used. This is particularly important if the available database of predicted and measured values is not large enough to directly estimate additive genetic correlations between measured and predicted traits. Our findings support the possibility of using NIRS spectra and calibration equations to genetically improve color traits, whereas DL, CL, and SF need to be investigated using larger data sets. LITERATURE CITED Aass L. 1996. Variation in carcass and meat quality traits and their relations to growth in dual purpose cattle. Livest. Prod. Sci.  46: 1– 12. Google Scholar CrossRef Search ADS   Abeni F. Bergoglio G. 2001. Characterization of different strains of boiler chicken by carcass measurements, chemical and physical parameters and NIRS on breast muscle. Meat Sci.  57: 133– 137. Google Scholar CrossRef Search ADS PubMed  Albera A. Mantovani R. Bittante G. Groen A. F. Carnier P. 2001. Genetic parameters for daily live-weight gain, live fleshiness and bone thinness in station tested Piemontese young bulls. Anim. Sci.  72: 449– 456. Google Scholar CrossRef Search ADS   AMSA (American Meat Science Association) 1991. Guidelines for Meat Color Evaluation.  44th Recip. Meat Conf. Natl. Live Stock and Meat Board, Chicago, IL. Andrés S. Silva A. Soares-Pereira A. L. Martins C. Bruno-Soares A. M. Murray I. 2008. The use of visible and near infrared reflectance spectroscopy to predict beef M. Longissimus thoracis et lumborum quality attributes. Meat Sci.  78: 217– 224. Google Scholar CrossRef Search ADS PubMed  ASPA (Associazione Scientifica di Produzione Animale) 1996. Metodiche per la determinazione delle caratteristiche qualitative della carne.  Università degli studi di Perugia, Perugia, Italy. Biagini D. Lazzaroni C. 2005. Carcass dissection and commercial meat yield in Piemontese and Belgian Blue double-muscled young bulls. Livest. Prod. Sci.  98: 199– 204. Google Scholar CrossRef Search ADS   Blasco A. 2005. The use of Bayesian statistics in meat quality analyses: A review. Meat Sci.  69: 115– 122. Google Scholar CrossRef Search ADS PubMed  Boukha A. De Marchi M. Albera A. Bittante G. Gallo L. Carnier P. 2007. Genetic parameters of beef quality traits for Piemontese cattle. Ital. J. Anim. Sci.  6( Suppl. 1): 53– 55. Budowle B. Garofano P. Hellman A. Ketchum M. Kanthaswamy S. Parson W. van Haeringen W. Fain S. Broad T. 2005. Recommendations for animal DNA forensic and identity testing. Int. J. Legal Med.  119: 295– 302. https://doi.org/15834735 Google Scholar CrossRef Search ADS PubMed  Burrow H. M. Moore S. S. Johnston D. J. Barendse W. Bindon B. M. 2001. Quantitative and molecular genetic influences on properties of beef: A review. Aust. J. Exp. Agric.  41: 893– 919. Google Scholar CrossRef Search ADS   Carnier P. Albera A. Dal Zotto R. Groen A. F. Bona M. Bittante G. 2000. Genetic parameters for direct and maternal calving ability over parity in Piedmontese cattle. J. Anim. Sci.  78: 2532– 2539. https://doi.org/11048917 Google Scholar CrossRef Search ADS PubMed  Cecchi F. Russo C. Prezioso G. Cianci D. 2004. Genetic variability of meat quality traits in Chinina beef cattle. Ital. J. Anim. Sci.  3: 191– 198. Google Scholar CrossRef Search ADS   Cecchinato A. De Marchi M. Gallo L. Bittante G. Carnier P. 2009. Mid-infrared spectroscopy predictions as indicator traits in breeding programs for enhanced coagulation properties of milk. J. Dairy Sci.  92: 5304– 5313. https://doi.org/19762848 Google Scholar CrossRef Search ADS PubMed  Chan D. E. Walker P. N. Mills E. W. 2002. Prediction of pork quality characteristics using visible and near-infrared spectroscopy. J. ASAE  45: 1519– 1527. CIE (Commission Internationale de l'Éclairage) 1976. Official Recommendations on Uniform Color Spaces, Color Differences Equations and Metric Color Terms. Suppl. No. 2, CIE Publ. No. 15. Colorimetry.  Commission Internationale de l'Éclairage, Paris, France. Dal Maso M. Schiavon S. Tagliapietra F. Simonetto A. Bittante G. 2009. Growth performance and N excretion of double muscled Piemontese bulls fed low protein rations with or without the addition of rumen protected conjugated linoleic acid. Ital. J. Anim. Sci.  8( Suppl. 3): 175– 177. Google Scholar CrossRef Search ADS   Dal Zotto R. De Marchi M. Cecchinato A. Penasa M. Cassandro M. Carnier P. Gallo L. Bittante G. 2008. Reproducibility and repeatability of milk coagulation properties and predictive ability of mid-infrared reflectance spectroscopy. J. Dairy Sci.  91: 4103– 4112. https://doi.org/18832238 Google Scholar CrossRef Search ADS PubMed  De Marchi M. Fagan C. C. O'Donnell C. P. Cecchinato A. Dal Zotto R. Cassandro M. Penasa M. Bittante G. 2009. Prediction of coagulation properties, titratable acidity, and pH of bovine milk using mid-infrared spectroscopy. J. Dairy Sci.  92: 423– 432. https://doi.org/19109300 Google Scholar CrossRef Search ADS PubMed  Destefanis G. Brugiapaglia A. Barge M. T. Lazzaroni C. 2003. Effect of castration on meat quality in Piemontese cattle. Meat Sci.  64: 215– 218. Google Scholar CrossRef Search ADS PubMed  Dikeman M. E. 1990. Genetic effects on the quality of meat from cattle. Proc. 4th World Congr. Genet. Appl. Livest. Prod.  15: 521– 530. Dikeman M. E. Pollak E. J. Zhang Z. Moser D. W. Gill C. A. Dressler E. A. 2005. Phenotypic ranges and relationships among carcass and meat palatability traits for fourteen cattle breeds, and heritabilities and expected progeny differences for Warner-Bratzler shear force in three cattle breeds. J. Anim. Sci.  83: 2461– 2467. https://doi.org/16160060 Google Scholar CrossRef Search ADS PubMed  Fiems L. O. De Campeneere S. Van Caelenbergh W. De Boever J. L. Vanacker J. M. 2003. Carcass and meat quality in double-muscled Belgian Blue bulls and cows. Meat Sci.  63: 345– 352. Google Scholar CrossRef Search ADS PubMed  Forrest J. C. Morgan M. T. Borggaard C. Rasmussen A. J. Jespersen B. L. Andersen J. R. 2000. Development of technology for the early post mortem prediction of water holding capacity and drip loss in fresh pork. Meat Sci.  55: 115– 122. Google Scholar CrossRef Search ADS PubMed  Geesink G. H. Schreutelkamp F. H. Frankhuizen R. Vedder H. W. Faber N. M. Kranen R. W. Gerritzen M. A. 2003. Prediction of pork quality attributes from near infrared reflectance spectra. Meat Sci.  65: 661– 668. Google Scholar CrossRef Search ADS PubMed  Gelfand A. Smith A. F. M. 1990. Sampling based approaches to calculating marginal densities. J. Am. Stat. Assoc.  85: 398– 409. Google Scholar CrossRef Search ADS   Gelman A. Rubin D. B. 1992. Inference from iterative simulation using multiple sequences. Stat. Sci.  7: 457– 472. Google Scholar CrossRef Search ADS   Geweke, J. 1992. Evaluating the accuracy of sampling-based approaches to the calculation of posterior moments (with discussion). Pages 164–193 in Bayesian Statistics.  J. O. Berger, J. M. Bernardo, A. P. Dawid, and A. F. M. Smith ed. Oxford University Press, Oxford, UK. Geyer C. J. 1992. Practical Markov chain Monte Carlo. Stat. Sci.  7: 473– 483. Google Scholar CrossRef Search ADS   Hildrum K. I. Isaksson T. Naes T. Nilsen B. N. Rødbotten M. Lea P. 1995. Near infrared reflectance spectroscopy in the prediction of sensory properties of beef. J. Near Infrared Spectrosc.  3: 81– 87. Google Scholar CrossRef Search ADS   Hoving-Bolink A. H. Vedder H. W. Merks J. W. M. de Klein W. J. H. Reimert H. G. M. Frankhuizen R. van den Broek W. H. A. M. Lambooij E. E. 2005. Perspective of NIRS measurements early post mortem for prediction of pork quality. Meat Sci.  69: 417– 423. Google Scholar CrossRef Search ADS PubMed  Hubert M. Van den Branden K. 2003. Robust methods for partial least squares regression. J. Chemometr.  17: 537– 549. Google Scholar CrossRef Search ADS   Johnston D. J. Reverter A. Ferguson D. M. Thompson J. M. Burrow H. M. 2003. Genetic and phenotypic characterisation of animal, carcass, and meat quality traits from temperate and tropically adapted beef breeds. 3. Meat quality traits. Aust. J. Agric. Res.  54: 135– 147. Google Scholar CrossRef Search ADS   Kambadur R. Sharma M. Smith T. P. L. Bass J. J. 1997. Mutations in myostatin (GDF8) in double-muscled Belgian Blue and Piedmontese cattle. Genome Res.  7: 910– 916. https://doi.org/9314496 Google Scholar CrossRef Search ADS PubMed  Kizilkaya K. Banks B. D. Carnier P. Albera A. Bittante G. Templeman R. J. 2002. Bayesian inference strategies for the prediction of genetic merit using threshold models with an application to calving ease scores in Italian Piemontese cattle. J. Anim. Breed. Genet.  119: 209– 220. Google Scholar CrossRef Search ADS   Kizilkaya K. Carnier P. Albera A. Bittante G. Templeman R. J. 2003. Cumulative t-link threshold models for the genetic analysis of calving ease scores. Genet. Sel. Evol.  35: 489– 512. https://doi.org/12939202 Google Scholar CrossRef Search ADS PubMed  Leroy B. Lambotte S. Dotreppe O. Lecocq H. Istasse L. Clinquart A. 2003. Prediction of technological and organoleptic properties of beef Longissimus thoracis from near-infrared reflectance and transmission spectra. Meat Sci.  66: 45– 54. Google Scholar CrossRef Search ADS   Liu Y. Lyon B. G. Windham W. R. Lyon C. E. Savage E. M. 2004. Prediction of physical, color, and sensory characteristics of broiler breasts by visible/near infrared reflectance spectroscopy. Poult. Sci.  83: 1467– 1474. https://doi.org/15339027 Google Scholar CrossRef Search ADS PubMed  Mantovani R. Cassandro M. Contiero B. Albera A. Bittante G. 2010. Genetic evaluation of type traits in hypertrophic Piemontese cows. J. Anim. Sci.  88: 3504– 3512. https://doi.org/20675599 Google Scholar CrossRef Search ADS PubMed  Meulemans A. Dotreppe O. Leroy B. Istasse L. Clinquart A. 2003. Prediction of organoleptic and technological characteristics of pork meat by near infrared spectroscopy. Sci. Aliments  23: 159– 162. Google Scholar CrossRef Search ADS   Mitsumoto M. Maeda S. Mitsuhashi T. Ozawa S. 1991. Near-infrared spectroscopy determination of physical and chemical characteristics in beef cuts. J. Food Sci.  56: 1493– 1496. Google Scholar CrossRef Search ADS   Prevolnik M. Candek-Potokar M. Škorjanc D. 2004. Ability of NIR spectroscopy to predict meat chemical composition and quality: A review. Czech J. Anim. Sci.  49: 500– 510. Prieto N. Andrés S. Giráldez F. J. Mantecón A. R. Lavín P. 2008. Ability of near infrared reflectance spectroscopy (NIRS) to estimate physical parameters of adult steers (oxen) and young cattle meat samples. Meat Sci.  79: 692– 699. Google Scholar CrossRef Search ADS PubMed  Prieto N. Ross D. W. Navajas E. A. Nute G. R. Richardson R. I. Hyslop J. J. Simm G. Roehe R. 2009. On-line application of visible and near infrared reflectance spectroscopy to predict chemical–physical and sensory characteristics of beef quality. Meat Sci.  83: 96– 103. https://doi.org/20416617 Google Scholar CrossRef Search ADS PubMed  Ribeca C. Bittante G. Albera A. Bonfatti V. Maretto F. Gallo L. 2009. Investigation on variability of candidate genes for meat quality traits in Piemontese cattle. Ital. J. Anim. Sci.  8( Suppl. 2): 132– 134. Google Scholar CrossRef Search ADS   Riley D. G. Chase C. C.Jr. Hammond A. C. West R. L. Johnson D. D. Olson T. A. Coleman S. W. 2003. Estimated genetic parameters for palatability traits of steaks from Brahman cattle. J. Anim. Sci.  81: 54– 60. https://doi.org/12597372 Google Scholar CrossRef Search ADS PubMed  Ripoll G. Alberti P. Panea B. Olletaq J. L. Sanudo C. 2008. Near-infrared reflectance spectroscopy for predicting chemical, instrumental and sensory quality of beef. Meat Sci.  80: 697– 702. Google Scholar CrossRef Search ADS PubMed  Rutten M. J. M. Bovenhuis H. van Arendonk J. A. M. 2010. The effect of the number of observations used for Fourier transform infrared model calibration for bovine milk fat composition on the estimated genetic parameters of the predicted data. J. Dairy Sci.  93: 4872– 4882. https://doi.org/20855022 Google Scholar CrossRef Search ADS PubMed  Schiavon S. De Marchi M. Tagliapietra F. Bailoni L. Bittante G. 2011. Effect of high or low protein ration combined or not with rumen protected conjugated linoleic acid (CLA) on meat CLA content and quality traits of double-muscled Piemontese bulls. Meat Sci.  89: 133– 142. https://doi.org/21561723 Google Scholar CrossRef Search ADS PubMed  Schiavon S. Tagliapietra F. Dal Maso M. Bailoni L. Bittante G. 2010. Effect of low-protein diets and rumen-protected conjugated linoleic acid on production and carcass traits of growing double-muscled Piemontese bulls. J. Anim. Sci.  88: 3372– 3383. https://doi.org/20601521 Google Scholar CrossRef Search ADS PubMed  Sorensen, D., and D. Gianola 2002. Likelihood, Bayesian, and MCMC Methods in Quantitative Genetics.  Springer-Verlag, New York, NY. Google Scholar CrossRef Search ADS   Soyeurt H. Dardenne P. Gillon A. Crocquet C. Vanderick S. Mayeres P. Bertozzi C. Gengler N. 2006. Variation in fatty acid contents of milk and fat within and across breeds. J. Dairy Sci.  89: 4858– 4865. https://doi.org/17106116 Google Scholar CrossRef Search ADS PubMed  Splan R. K. Cundiff L. V. Van Vleck L. D. 1998. Genetic parameters for sex-specific traits in beef cattle. J. Anim. Sci.  76: 2272– 2278. https://doi.org/9781482 Google Scholar CrossRef Search ADS PubMed  Vincenti F. Failla S. Gigli S. Lasagna E. Landi V. Mangione A. Berti C. Sarti F. M. 2007. The hypertrophic Marchigiana: Physical and biochemical parameters for meat quality evaluation. Ital. J. Anim. Sci.  6( Suppl. 1): 491– 493. Warner R. D. Greenwood P. L. Pethick D. W. Ferguson D. M. 2010. Genetic and environmental effects on meat quality. Meat Sci.  86: 171– 183. https://doi.org/20561754 Google Scholar CrossRef Search ADS PubMed  Wheeler T. L. Cundiff L. V. Shackelford S. D. Koohmaraie M. 2001 a. Characterization of biological types of cattle (Cycle V): Carcass traits and longissimus palatability. J. Anim. Sci.  79: 1209– 1222. https://doi.org/11374541 Google Scholar CrossRef Search ADS PubMed  Wheeler T. L. Shackelford S. D. Casas E. Cundiff L. V. Koohmaraie M. 2001 b. The effects of Piedmontese inheritance and myostatin genotype on the palatability of longissimus thoracis, gluteus medius, semimembranosus, and biceps femoris. J. Anim. Sci.  79: 3069– 3074. https://doi.org/11811461 Google Scholar CrossRef Search ADS PubMed  Williams, P. 2003. Near-infrared technology getting the best out of light. Page 109 in A Short Course in the Practical Implementation of Near Infrared Spectroscopy for the User.  1.1 ed. PDKProjects Inc., Nanaimo, Canada. Williams, P. C. 1987. Implementation of near-infrared technology. Pages 143–167 in Near-Infrared Technology in the Agricultural and Food Industries. 2nd ed. P. C. Williams and K. Norris, ed. Am. Assoc. Cereal Chem., St. Paul, MN. American Society of Animal Science This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/), which permits non-commercial reuse, distribution, and reproduction in any medium, provided the original work is properly cited. For commercial re-use, please contact journals.permissions@oup.com TI - Near-infrared reflectance spectroscopy predictions as indicator traits in breeding programs for enhanced beef quality JF - Journal of Animal Science DO - 10.2527/jas.2010-3740 DA - 2011-09-01 UR - https://www.deepdyve.com/lp/oxford-university-press/near-infrared-reflectance-spectroscopy-predictions-as-indicator-traits-K7Pk66viqy SP - 2687 EP - 2695 VL - 89 IS - 9 DP - DeepDyve ER -