|
|
||||||||
ANIMAL GENETICS |

* Department of Animal and Poultry Science, University of Guelph, Ontario, Canada N1G-2W1; and
and
Department of Animal Science, Washington State University, Pullman 99164-6351
| Abstract |
|---|
|
|
|---|
0.05). At 7 d postmortem, the percentage of unacceptably tough steaks decreased by 24 and 35%, respectively, for animals carrying 1 and 2 copies of the C allele relative to animals with no C alleles. However, genotype CC had a greater fat yield (+1.44 ± 0.56%; P = 0.037) than genotype GG, with a corresponding allele substitution effect of 0.67 ± 0.27% (P = 0.015). Therefore, the CAST SNP allele C was associated with increased LM tenderness across days of postmortem aging and, importantly for the beef industry, had a significant reduction in the percentage of steaks rated unacceptably tough by consumers based on an assumed threshold level.
Key Words: beef breed calpastatin gene carcass trait meat quality single nucleotide polymorphism tenderness
| INTRODUCTION |
|---|
|
|
|---|
The physiological change in muscle structure during the postmortem period is complex (Koohmaraie, 1994
). The calpain/calpastatin system is an endogenous, calcium-dependent proteinase system, theorized to mediate the proteolysis of key myofibrillar proteins during postmortem storage of carcass and cuts of meat at refrigerated temperatures (Koohmaraie et al., 1995b
). Calpain is responsible for the breakdown of myofibrillar proteins, which are closely related to meat tenderness (Wheeler and Koohmaraie, 1994
). Calpastatin (CAST) inhibits µ- and m-calpain activity and, therefore, regulates postmortem proteolysis. Increased postmortem CAST activity has been correlated with reduced meat tenderness (e.g., Koohmaraie et al., 1995a
; Pringle et al., 1997
).
The CAST gene, mapped to BTA 7 (Bishop et al., 1993
), is considered a candidate gene for beef tenderness. Initial studies did not detect significant association of CAST polymorphisms with tenderness (e.g., Lonergan et al., 1995
; Chung et al., 1999
). However, in both studies the sample sizes were small, which may have contributed to the failure to identify significant associations between the CAST gene polymorphisms and meat tenderness.
More recently, genetic tests for meat tenderness in beef, which utilize genetic polymorphisms in the CAST and/or Calpain genes, have been made available by private companies. Examples are the IGENITY TenderGENE test (Merial Ltd., Atlanta, GA) and the Gene-STAR Tenderness 2 test (Genetic Solutions Pty. Ltd., Albion, Australia).
The objective of this study was to assess the association of a newly discovered SNP in the CAST gene with meat tenderness and other important meat and carcass quality traits in a large sample of crossbred beef cattle.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Animals were crossbred with breed composition formed by several breeds. The major contributing breeds were Angus, Limousin, Charolais, and Simmental. The average contribution of these 4 breeds to the breed composition of animals having any non-zero fraction of the mentioned breeds was 0.46, 0.50, 0.50, and 0.50 for Angus, Limousin, Charolais, and Simmental, respectively, for Commercial cattle and was 0.51, 0.59, 0.53, and 0.41 for Rockwood cattle. Given the primary use of purebred sires, Commercial cattle were more representative of first generation crossbreds than Rockwood cattle.
Rockwood and Commercial cattle represented AI breeding as well as some herd bulls. Animals were finished in commercial feedlots and marketed based on the feedlot operators definition of finish, which was a visual appraisal of fat cover and weight. Resulting HCW and grade fat levels are listed in Table 1
.
|
Primer Design and Sequencing of Amplified DNA Fragments
One pair of primersforward: 5'-cct cga ctg cgt acc aat tcc gaa gta aag cca aag gaa ca-3' and reverse: - 5'-att tct ctg atg gtg gct gct cac t-3'was designed based on a partial genomic DNA sequence of the bovine CAST gene (Accession number AY008267), using the software Primer 3 (http://www-genome.wi.mit.edu/cgi-bin/primer/primer3.cgi). To simplify the sequencing procedures when performing many reactions for different genes, a common tail was also added to the forward primer: 5'-cct cga ctg cgt acc aat tcc-3'.
The PCR master mix was 6.71 µL of distilled H2O, 1.0 µL of 10x PCR buffer II, 0.73 µL of 25 mM MgCl2 solution, 0.6 µL of 10 mM dNTPs Solution, 0.13 µL of forward primer (20 pmol/µL), 0.13 µL of reverse primer (20 pmol/µL), 0.1 µL of AmpliTaq Gold (5 U/µL; all from Applied Biosystems, Foster City, CA), and 0.6 µL of DNA template (50 ng/µL) in a 10-µL reaction solution. The PCR cycling condition was 95°C (10.0 min) for 1 cycle; then 94°C (30 s), 69 to 62°C (30 s) and 72°C (30 s) for 8 cycles; then 94°C (30 s), 61°C (30 s), and 72°C (30 s) for 27 cycles; and finally 72°C (5.0 min) for 1 cycle, maintaining 4°C thereafter. The PCR reactions were performed on the GenAmp PCR System 9700 (Applied Biosystems). The corresponding PCR products were 523 bp. Sequencing of the PCR products revealed a G/C transversion in the products, resulting in the presence/absence of a RsaI restriction site at the position 257 bp with recognizing sequence gt/ac. Sequencing was performed on the ABI PRISM 377 DNA sequencer (Applied Biosystems) by the Laboratory Service Division, University of Guelph.
PCR-RFLP Genotyping
Individual PCR-RFLP genotypes were distinguished for polymorphisms of CAST RsaI. The PCR products were digested by the restriction endonuclease. Digestion was conducted at 37°C for 4 h and in a 10-µL reaction solution including 2.7 µL of distilled H2O, 1.0 µL of specific NEBuffer, 0.3 µL (3 units) of restriction endonucleases, and 6 µL of PCR product solution. The DNA fragments from the digested PCR products were separated by 2.0%-agarose (Fisher Scientific Ltd., Nepean, Ontario, Canada) gel containing ethidium bromide (0.4 µg/mL). Electrophoresis was performed in 1x TBE buffer including 108 g of Tris, 55 g of boric acid, and 40 mL of 0.5 M EDTA in 1,000 mL of 10x concentrated stock solution, pH 8.0, containing ethidium bromide (3 µL of 10 mg EB/mL per 100 mL of gel solution) under 120 V for about 45 min. Genotype for each individual was read under ultraviolet light, using Molecular Analyst Software (BioRad Laboratories, Molecular Bioscience Group, Hercules, CA).
Phenotypic Information
Information on tenderness (shear force) of LM at 2 (SFL2), 7 (SFL7), 14 (SFL14), and 21 (SFL21) d postmortem and of semitendinosus muscle at 7 (SFS7) d postmortem, chemical fat (CF), grade fat (GFAT), marbling grade (MG), LM area (LMA), percentage of lean (LEANYL), fat (FATYL), and bone (BONEYL) yield, and HCW were available on most of the 628 genotyped animals.
Warner-Bratzler peak shear force measurements (kg) were used as an objective method of assessing tenderness (Shackelford et al., 1999
). The GFAT was the back-fat thickness measurement taken at the 12th and 13th rib interface. The LMA was the measure of the LM area at the 12th and 13th rib interface using a tracing of the muscle. The CF was the chemical analysis on a core meat sample that determined the percentage of i.m. fat in the LM. The LEANYL, FATYL, and BONEYL were determined by dissection of a 4-bone rib section.
Marbling grade was the grade used for grading in Canada (A, AA, AAA, Prime) with most carcasses falling in 1 of the first 3 categories. The assessment of marbling is based on the average amount, size, and distribution of fat particles or deposits in the LM. Canadian beef carcass grading utilizes only 4 of the 9 recognized levels of marbling from the USDA marbling standards with the following minimums for each category: traces (A), slight (AA), small (AAA), and slightly abundant (Prime) marbling (CBGA, 2005
). Because only a few carcasses were classified as Prime, those carcasses were combined with AAA carcasses for the analyses. A complete description of carcass measures available and methods used are given in Laborde et al. (2001)
.
Cattle were slaughtered at the slaughter plant of Better Beef Ltd. in Guelph. Standardized carcass and meat quality measures were made either at the University of Guelph Meat Science Laboratory (FATYL, LEANYL, BONEYL, CF, and tenderness evaluation of LM and semitendinosus muscle) or at the Better Beef plant (GFAT, HCW, LMA, and MG) by University of Guelph meat scientists. Commercial cattle were slaughtered from 1998 to 2000, and Rockwood cattle were slaughtered in 1998.
The genotypes of all animals were used to determine the allele frequencies. For the study of the association of CAST SNP with carcass and meat quality traits, only animals with required phenotypic information were used. The resulting number of records ranged from 466 for SFL2 to 628 for SFL7. Table 1
gives the number of records, means, and SD of the analyzed traits.
Statistical Analyses
All analyses were performed using the statistical software SAS (SAS Institute, Inc., Cary, NC). Descriptive characteristics of quantitative traits were obtained using PROC MEANS. Allele frequencies were tabulated and compared by
2 analysis using PROC FREQ.
Genotype Analyses.
Associations of the genotypes for the SNP in the CAST gene with CF, GFAT, LMA, LEANYL, FATYL, BONEYL, and HCW were evaluated using PROC MIXED, fitting the following model, which included the known SNP genotypes as fixed effects:
![]() |
in which Yijklm is the trait measured in the animal m of the sex j and the slaughter group k having the CAST SNP genotype i and born from the sire l; u is the overall mean for the trait; Geni is the effect of the genotype i for the SNP in the CAST gene (CC, CG, or GG); Sexj is the fixed effect of the sex j (bull, heifer or steer); Slgk is the fixed effect of the slaughter group k (61 levels); ß1, ß2, ß3, ß4 are the regression coefficients of the linear regressions on breed composition (proportion of Angus, Limousin, Charolais, and Simmental, respectively); Sirel is the random effect of the sire l; and eijklm is the residual random effect associated with animal m.
Slaughter groups were defined as animals from the same source (Commercial or Rockwood) with the same slaughter date. The slaughter date accounted for most of the feedlot-of-origin variation in the commercial cattle, because the majority of the slaughter dates (90%) had cattle from a single feedlot. All other slaughter dates (10%) had cattle from only 2 feedlots. If the latter slaughter dates were divided into 2 groups, the resulting number of animals within the corresponding slaughter group would become too small, so they were maintained as a single group.
Animals originated from 78 sires, and all sires were known and assumed unrelated. The average size of the paternal half-sib families was 8.1. The percentages of sires with less than 5, from 6 to 10, from 11 to 15, and more than 15 offspring were 42.3, 23.1, 25.6, and 9.0%, respectively.
For analyzing the shear force measurements, CF was included as a linear covariate in the model to account for possible effects of i.m. fat on meat tenderness:
![]() |
in which Yijklm is the shear force measurement of the animal m of sex j and slaughter group k having the CAST SNP genotype i, and born from the sire l; B5 is the regression coefficient of the linear regression on CF; and all other terms in model 2 are as previously described for model 1.
Following Bertrand et al. (2001)
, phenotypic and genetic relationships between marbling and tenderness are not especially high but show a favorable direction, indicating that greater marbling is slightly associated with greater tenderness.
The repeated LM shear force measurements were analyzed individually within each postmortem aging period and as the average shear force measure across d 5 of postmortem aging, using model 2. Because these analyses ignore the covariance structure due to the sequential nature of the data on each animal, a repeated measures analysis of tenderness was also carried out, which modeled the covariance structure and allowed comparison of genotypes and trends over aging times (Littell et al., 1998
).
The repeated measures analysis was carried out by PROC MIXED of SAS, using the Random and Repeated statements. The same effects as in model 2 were fitted in addition to the fixed effects of aging time (postmortem days) and the genotype by aging time interaction and the random effect of animal, which was the subject on which repeated measures across aging times were taken. The covariance structure within animals was assumed as block diagonal with one 4 x 4 block (4 measures) per animal. A heterogeneous first order autoregressive covariance structure was used:
![]() |
in which
is the variance associated with the postmortem aging period i and
is the autoregressive parameter. The (co)variance components were estimated by the REML method.
A heterogeneous first order autoregressive covariance structure within animal was chosen based on preliminary analyses, in which different covariance structures were modeled and the best fitting model was selected using the Akaike information criterion (Akaike, 1973
).
The numerator df to test the fixed effects was obtained using the option DDFM=BETWITHIN. This option divides the residual df into between-subject and within-subject portions, and a check is made whether a fixed effect changes within any subject. If so, the within-subject df are assigned to the effect; otherwise the between-subject df are assigned. Genotype least squares means for each postmortem aging period were obtained using the Estimate statement.
The effect of CAST SNP genotypes on MG was analyzed by
2 analysis (PROC FREQ) and as a linear trait, applying model 1. In this case, scores of 1, 2, and 3 were assigned to grades A, AA, and AAA, respectively.
Average allele substitution effects (Falconer and Mackay, 1996
) were estimated using the same models used to estimate the genotype effects but replacing the classification effect of genotypes by a linear regression on the number of C alleles (0, 1, or 2).
The association of the number of C alleles with the percentage of unacceptably tough LM steaks (LM shear force > 5.7 kg) was also determined. Researchers at the USDA (Wheeler et al., 1997
) showed that when peak shear force was above the threshold of 5.7 kg, 100% of consumers would rate a steak as unacceptably tough.
The percentage of animals from each genotype having a LM shear force above the threshold was calculated using both the observed shear force and the shear force adjusted for sex, slaughter group, breed, CF, and sire effects. In the latter case, the adjusted shear force was obtained using the residuals of model 2 used to analyze shear force within postmortem aging periods, as previously described, but without fitting the SNP genotype effect. These residuals were added to the expected value of an observation in a given postmortem aging period to produce the adjusted shear force. Expected values were the least squares means for the time (postmortem aging period) effect given by the repeated measures analysis. The association of the number of C alleles with the percentage of unacceptably tough steaks was evaluated by
2 analysis.
To keep reasonable probability values for Type I errors, 2 levels of tests were performed. For initial assessment of the results, an overall value of P < 0.05 (
) was used. For a more detailed review of the results, a modified Bonferroni correction was used (
/
n; Mantel, 1980
) to account for the number of tests. The value of n was determined using a trait-wise approach, grouping traits according to type (Ye, 2003
). Traits were grouped into 2 groups as follows: carcass yield traits (LEANYL, FATYL, BONEYL, GFAT, LMA, and HCW) and meat quality traits (CF, MG, SFL2, SFL7, SFL14, SFL21, and SFS7). So n was equal to 6 and 7 for carcass yield and meat quality traits, respectively, with the corresponding modified Bonferroni significance levels of 0.020 and 0.019.
| RESULTS |
|---|
|
|
|---|
|
2 test, and P = 0.70 by F-test, treating MG as a linear trait).
Results of the analyses of LM tenderness within each postmortem aging period and using the average LM tenderness measure over the 4 measurements are presented in Table 3
. Analyses of the SNP genotypes within each postmortem aging period showed a highly significant effect (P = 0.007) on SFL21 and a trend to affect SFL2 and SFL7 (P
0.07). The CAST genotypes significantly influenced (P = 0.012) the average LM tenderness (SFLavg) over the 21-d postmortem period; genotype CC had the most tender LM compared with genotypes CG and GG. The difference between the homozygous genotypes CC and GG was 0.38 ± 0.13 kg; the heterozygous genotype showed an intermediate tenderness.
|
0.03) for all measures of LM tenderness, except for SFL14, for which there was a trend (P = 0.10). The allele substitution effect ranged from 0.26 ± 0.11 kg to 0.13 ± 0.06 kg for tenderness evaluated at 2 and 21 d of aging, respectively. Genotype effect on SFL21 and SFLavg and average allele substitution effect on SFL2 and SFLavg were significant after modified Bonferroni correction for trait-wise multiple tests (P < 0.05), providing strong evidence for the effect of the SNP on LM tenderness.
In line with the effect on LM tenderness, CAST SNP genotypes tended to affect SFS7 (P = 0.14), with increased tenderness in beef from genotype CC vs. CG and GG (0.17 ± 0.09 kg and 0.10 ± 0.12 kg, respectively). However, the average allele substitution effect on SFS7 was small and did not differ from zero (P = 0.22; Table 3
).
Results of the repeated measures analysis are presented in Table 4
. A highly significant overall effect of genotypes across aging days (P = 0.005) was found; genotype CC had more tender LM (0.32 ± 0.13 kg) than GG, and genotype CG had an intermediate tenderness. The genotype least squares means for each postmortem aging period are also shown in Table 4
. Means were significantly different (P < 0.05) for LM tenderness at 2, 7, and 21 d postmortem; genotype CC showed increased tenderness. However, the difference between the homozygous genotypes CC and GG decreased from 0.53 ± 0.20 kg to 0.21 ± 0.12 kg from 2 to 21 d of aging, respectively. These differences correspond to 25 and 20% of the phenotypic SD of the shear force measurements at 2 and 21 d postmortem, respectively.
|
Table 4
also presents the average allele substitution effect estimated by repeated measures analysis, which was highly significant on the overall tenderness across the postmortem periods (P = 0.002) and significant on the measures of LM tenderness at different postmortem aging periods (P < 0.03), except for SFL14. The allele substitution effect for each postmortem aging period was estimated by fitting the regression on the number of C alleles nested within aging time. Similar to the previous analyses of individual postmortem aging periods, the allele substitution effect given by the repeated measures analysis ranged from 0.27 ± 0.10 kg to 0.13 ± 0.06 kg for tenderness evaluated at 2 and 21 d postmortem, respectively.
The genotype effect on the overall tenderness across postmortem aging periods and average allele substitution effect on SFL2 and SFL7, and on the overall tenderness across postmortem aging periods were significant after modified Bonferroni correction for trait-wise multiple tests (P < 0.05), providing again strong evidence for the effect of the CAST SNP on tenderness.
There was no significant effect of either breed (P > 0.07) or sex (P > 0.20) in all analyses of tenderness (both within postmortem aging periods and by repeated measures analysis). On the contrary, slaughter group showed a highly significant effect (P < 0.001) on tenderness in all analyses (data not shown).
The covariate CF had a highly significant linear effect on SFL7, SFL14, and SFL21 (P < 0.007), with an estimated change in shear force of 0.11 kg, 0.09 kg, and 0.07 kg per 1% increase in CF for SFL7, SFL14, and SFL21, respectively (data not shown).
Table 5
presents the association of the number of C alleles with the percentage of unacceptably tough LM steaks (LM shear force > 5.7 kg). The percentage of unacceptably tough steaks based on the observed shear force values was associated with the number of C alleles at 2 (P = 0.001) and 7 (P = 0.05) d postmortem. At 2 d postmortem, 31.5%, 32.8%, and 18.9% of the steaks from cattle carrying 0, 1, and 2 copies of C allele were unacceptably tough, respectively. The corresponding features for 7 d postmortem were 35.2, 26.0, and 23.3%, respectively.
|
Table 6
shows that CAST SNP tended to influence LMA (P = 0.12) and LEANYL (P = 0.12) and significantly influenced FATYL (P = 0.037) and BONEYL (P = 0.037). The CC genotype tended to have smaller LMA (1.77 ± 1.50 cm2) and lower LEANYL (0.97 ± 0.51%) with greater fat yield (+1.44 ± 0.56%) and lower bone yield (0.47 ± 0.22%) than the GG genotype. The average allele substitution effects on FATYL and BONEYL were significant (P < 0.02), even after modified Bonferroni correction for trait-wise multiple tests (P < 0.05).
|
| DISCUSSION |
|---|
|
|
|---|
5/8 Simmental (n = 33) was genotyped.
The allele C was associated with significant increase in LM tenderness across d of postmortem aging. The magnitude of the effect of the C allele in this study was similar to that reported for the GeneSTAR tender allele (Genetic Solutions Pty. Ltd., Albion, Australia, http://www.geneticsolutions.com.au/files/GeneSTAR/pdf/GeneNOTE_4.pdf), which was 0.37 kg on LM shear force at 7 d postmortem between animals carrying 2 or no tender alleles. The polymorphism in the CAST gene used in the GeneSTAR test (Barendse, 2002
) is, however, different from the one reported in the current study.
The difference in LM shear force at 7 d postmortem between homozygous genotypes CC and GG given by the repeated measures analysis (0.37 ± 0.16 kg) was equivalent to 62% of the estimated difference between the shear force at 2 and 7 d postmortem (time effect, Table 4
), which was a measure of the tenderization of the meat as it aged. The magnitude of the benefit of the favorable genotype can then be put in perspective, relative to the benefits of aging, which are known and are a costly means of increasing beef tenderness. Selection for the beneficial genotype could then either increase beef tenderness, at the same level of aging, or provide a cost savings to attain the same level of tenderness with less aging.
Assuming a heritability of 30% for tenderness (Koch et al., 1982
) evaluated in any of the postmortem aging periods and using the phenotypic SD in Table 1
, the average allele substitution effects given by the repeated measures analysis (which were similar to those from the analyses within postmortem aging periods) would correspond to 29, 22, 14, and 23% of the respective genetic SD for shear force at 2, 7, 14, and 21 d postmortem, respectively, which are large genetic effects for a single SNP.
The number of C alleles was significantly associated with the percentage of unacceptably tough LM steaks at 2 and 7 d postmortem. The percentage of unacceptably tough steaks at 7 d postmortem decreased by 24 and 35% for cattle carrying 1 and 2 copies, respectively, of C allele compared relatively with steaks from animals with 0 copies, regardless if the percentage was calculated based on the observed or adjusted shear force. For 2 d postmortem, the presence of 2 copies of the C allele was associated with lower percentage of unacceptably tough steaks based in both observed (41% lower) and adjusted shear force (25% lower) compared relatively with steaks from animals with 1 or 0 copies.
The need for 2 C alleles to reduce the percentage of unacceptably tough steaks at 2 d of aging was expected because the least squares mean for shear force at 2 d postmortem for the heterozygous genotype CG (Tables 3
and 4
) was above the assumed toughness threshold (shear force >5.7 kg), making the effect of only 1 C allele masked in terms of reduced percentage of unacceptably tough steaks. At 14 and 21 d postmortem, there was not a clear association between the number of C alleles and the percentage of unacceptable tough steaks, likely because, at these 2 d of postmortem aging, the percentage of animals with steaks above the assumed toughness threshold was much smaller than at 2 and 7 d postmortem due to the tenderization of meat as it aged.
The magnitude of the association of number of C alleles with the percentage of unacceptably tough LM steaks at 7 d postmortem was smaller than that reported for the GeneSTAR tender allele (http://www.geneticsolutions.com.au/files/GeneSTAR/pdf/GeneNOTE_4.pdf), which was 33% and 62% for cattle carrying one and two copies of the tender allele relatively to animals with no copies, respectively.
Therefore, the CAST SNP C allele was associated with increased LM tenderness but tended to reduce LM area and lean yield and significantly increase fat yield. This evidence that the CAST SNP may also affect other important carcass traits, such as fat yield, warrants further investigation.
| IMPLICATIONS |
|---|
|
|
|---|
| Footnotes |
|---|
2 Corresponding author: Schenkel{at}uoguelph.ca
Received for publication July 29, 2005. Accepted for publication October 17, 2005.
| LITERATURE CITED |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
W. Barendse, B. E. Harrison, R. J. Hawken, D. M. Ferguson, J. M. Thompson, M. B. Thomas, and R. J. Bunch Epistasis Between Calpain 1 and Its Inhibitor Calpastatin Within Breeds of Cattle Genetics, August 1, 2007; 176(4): 2601 - 2610. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. L. Van Eenennaam, J. Li, R. M. Thallman, R. L. Quaas, M. E. Dikeman, C. A. Gill, D. E. Franke, and M. G. Thomas Validation of commercial DNA tests for quantitative beef quality traits J Anim Sci, April 1, 2007; 85(4): 891 - 900. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |