J. Anim Sci.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Page, B. T.
Right arrow Articles by Smith, T. P. L.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Page, B. T.
Right arrow Articles by Smith, T. P. L.
J. Anim. Sci. 2002. 80:3077-3085
© 2002 American Society of Animal Science

Evaluation of single-nucleotide polymorphisms in CAPN1 for association with meat tenderness in cattle1,2

B. T. Page*, E. Casas*, M. P. Heaton*, N. G. Cullen{dagger}, D. L. Hyndman{ddagger}, C. A. Morris{dagger}, A. M. Crawford{ddagger}, T. L. Wheeler*, M. Koohmaraie*, J. W. Keele* and T. P. L. Smith*,3

* USDA, ARS, U.S. Meat Animal Research Center, Clay Center, NE; and {dagger} AgResearch, Ruakura Research Center, PB3123, Hamilton, New Zealand; and and {ddagger} AgResearch, Molecular Biology Unit, Dunedin, New Zealand

3 Correspondence:
Spur 18D, P.O. Box 166 (phone: 402-762-4366; fax: 402-762-4390; E-mail:
smith{at}email.marc.usda.gov).


    Abstract
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Implications
 References
 
Micromolar calcium activated neutral protease (CAPN1) was evaluated as a candidate gene for a quantitative trait locus (QTL) on BTA29 affecting meat tenderness by characterization of nucleotide sequence variation in the gene. Single-nucleotide polymorphisms (SNP) were identified by sequencing all 22 exons and 19 of the 21 introns in two sires (Piedmontese x Angus located at the U.S. Meat Animal Research Center in Clay Center, NE; Jersey x Limousin located at AgResearch in New Zealand) of independent resource populations previously shown to be segregating meat tenderness QTL on BTA29. The majority of the 38 SNP were found in introns or were synonymous substitutions in the coding regions, with two exceptions. Exons 14 and 9 contained SNP that were predicted to alter the protein sequence by the substitution of isoleucine for valine in Domain III of the protein, and alanine for glycine in Domain II of the protein. The resource populations were genotyped for these two SNP in addition to six intronic polymorphisms and two silent substitutions. Analysis of genotypes and shear force values in both populations revealed a difference between paternal CAPN1 alleles in which the allele encoding isoleucine at position 530 and glycine at position 316 associated with decreased meat tenderness (increased shear force values) relative to the allele encoding valine at position 530 and alanine at position 316 (P < 0.05). The association of maternal alleles with meat tenderness phenotypes is consistent with the hypothesis of CAPN1 as the gene underlying the QTL effect in two independent resource populations and presents the possibility of using these markers for selective breeding to reduce the numbers of animals with unfavorable meat tenderness traits.

Key Words: Calpain • Genetic Polymorphism • Meat • Quantitative Trait Loci


    Introduction
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Implications
 References
 
Variation in meat tenderness has significant impact on consumer satisfaction with beef, however, classical selection techniques have not been effective in eliminating the small fraction of animals yielding undesirable palatability traits. Establishing the genetic basis for variation in meat tenderness would likely aid in the development of selection criteria for improving meat tenderness in cattle.

The micromolar calcium-activated neutral protease (CAPN1) gene encodes a cysteine protease, µ-calpain, that degrades myofibrillar proteins under postmortem conditions and appears to be the primary enzyme in the postmortem tenderization process (Koohmaraie, 1992; 1994; 1996). Regulation of µ-calpain activity has been correlated with variation in meat tenderness (Geesink and Koohmaraie, 1999). Bovine CAPN1 has been mapped to the telomeric end of BTA29 (Smith et al., 2000a), and recently, a QTL for tenderness was found to be segregating in this region of BTA29 in two resource populations (Casas et al., 2000; Morris et al., 2001). The moderate effect (0.4 standard deviation) of this QTL on Warner-Bratzler shear force, a mechanical measurement of tenderness (AMSA, 1995) in the U.S. Meat Animal Research Center (MARC) population supported the evaluation of CAPN1 as a candidate gene for tenderness within this population. Our objective was to determine if genetic variation within the candidate gene CAPN1 could be associated with the meat tenderness values of two resource populations exhibiting meat tenderness QTL on BTA29.


    Materials and Methods
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Implications
 References
 
Bacterial Artificial Chromosome Sequencing

A random shotgun library was prepared from a Bacterial Artificial Chromosome (BAC) clone (Smith et al., 2000a; Warren et al., 2000) containing the CAPN1 gene by partial digestion of ~5 µg of BAC DNA, with the restriction enzyme CviJI (0.3 units per µL) essentially as described (Gingrich et al., 1996). Briefly, the DNA was incubated with enzyme in a 50-µL volume for 20 min, the resulting DNA smear was separated on a 1% agarose gel, and fragments in the 1 to1.5 kb range were isolated using a commercial kit (Novagen, Madison, WI). The eluted fragments were ligated to a pBluescript vector (Stratagene, La Jolla, CA) that had been linearized with restriction enzyme EcoRV and treated with alkaline phosphatase. Ligated fragments were transformed into Top 10 competent cells (Invitrogen, Valencia, CA). Sequencing was performed on PCR-derived templates as described (Smith et al., 2000b), using T7 and T3 as universal sequencing primers on an ABI 3700 and ABI 377 sequencers.

The sequences from BAC subclones were analyzed with the sequence assembly algorithm Phrap (Ewing and Green 1998; Ewing et al., 1998) to detect overlaps and form contigs. The previously published (Smith et al., 2000a) partial sequence of the bovine CAPN1 gene was included in the data set to seed formation of contigs surrounding the gene. Sequence comparisons of the contig sequences with human genomic sequence surrounding the human CAPN1 gene (GenBank accession number NT_030106.2) and the E. coli genome were performed using BLASTN (Altshul et al., 1990; 1997).

Animal Sequencing and Single Nucleotide Polymorphisms Identification

Primers (Table 1Go) for PCR amplification were designed from the CAPN1 gene sequence (GenBank accession AF252504 and AF248054) to obtain overlapping DNA fragments of 1 to 2 kb in predicted length. Production of amplicons for direct sequencing was performed in a 20-µL volume with Amplitaq Gold or HotStar Taq polymerase enzyme following the protocols supplied by the manufacturers (Applied Biosystems, Inc.; Qiagen, Inc). Sequencing was performed using one-fourth reaction volumes with BigDye Terminator as recommended by the manufacturer (Applied Biosytems, Inc.) and analyzed on an ABI 3700 sequencer. Amplicons were sequenced with the amplification primers.


View this table:
[in this window]
[in a new window]
 
Table 1. Primer pairs used for amplification of genomic CAPN1
 
Sequence data were entered into the MARC database, base calls were made using Phred, and the sequences were aligned using Phrap (Ewing and Green, 1998; Ewing et al., 1998). Consed viewer was used to identify and tag single nucleotide polymorphisms (SNP) (Gordon et al., 1998) as described (Heaton et al., 2001). Genotyping was performed using a mass spectrometry-based MassArray system as suggested by the manufacturer (Sequenom, Inc., San Diego, CA), employing a PCR primer tailed with universal primer sequence in combination with a universal primer carrying a biotin tag as described (Stone et al., 2002). Polymerase chain reaction amplification for MassArray assays was carried out as outlined above except that reactions were carried out in 50-µL volumes.

Resource Populations

The MARC resource family used in this study has been previously described (Casas et al., 1998). Briefly, a half-sib family was developed using a Piedmontese x Angus sire. A total of 209 offspring were produced by matings primarily to MARC III dams (1/4 Angus, 1/4 Hereford, 1/4 Pinzgauer, and 1/4 Red Poll). Dams were artificially inseminated, and the offspring were born during the spring of 1995. Calves were weaned at an average of 200 d and raised from weaning to slaughter on a corn-corn silage diet. Steers were slaughtered at a commercial beef processing facility after 194 to 312 d on feed, and heifers were slaughtered after 207 to 287 d on feed. Meat tenderness was measured in all animals on the longissimus thoracis as Warner-Bratzler shear force (kg) at 3 and 14 d postmortem after belt cooking as described by Wheeler et al. (1998). This family was used to detect QTL based on differences between the two paternal alleles for growth, carcass composition, and meat quality traits (Casas et al., 2000; 2001).

The New Zealand AgResearch resource family used in this study was the progeny of a Limousin x Australian Jersey sire. The bull’s semen was imported into New Zealand for insemination of New Zealand Limousin and New Zealand Jersey cows (mainly grade animals), thus providing backcross animals of two types for a joint New Zealand-Australian DNA marker study. Importation into New Zealand of semen from this and two other bulls from the Australian "Mapping Project" (Pitchford et al., 1998), provided common genetic links between New Zealand and Australian grandparents for the joint study. The cows in New Zealand were inseminated in 1995 and 1996, and calves were born in 1996 and 1997 in industry and research herds. Heifer and steer calves were collected after weaning and grown out on a research property (Tokanui Station) near AgResearch’s headquarters. Further details of trial design and animal management (live calves and at slaughter) are given by Morris et al. (2001). Briefly, animals were preallocated to slaughter groups at 22 to 28 mo of age, over 18 kill days in 1998 and 10 kill days in 1999, with preallocation based on breed of calf, sire, and balanced as far as possible for live weight before the first slaughter day. The right striploin (longissimus lumbarum) was removed for tenderness measurement before stimulation of the remainder of the carcass. Each striploin was maintained at 15°C after removal. Five steak portions from the unstimulated striploin were cut for the cooking and shear force measurements with the first steak processed at rigor mortis and the remainder at 1.3, 2.0, 2.3, and 4.0 d postmortem. A different technique for measuring shear force was used in the AgResearch population as shear force was measured by MIRINZ tenderometer (Fraserhurst and MacFarlane, 1983).

Haplotype Analysis

Haplotypes inherited from the sires were established based on a selection of 10 SNP representing genetic variation within the Piedmontese x Angus sire. Two SNP representing predicted amino acid changes, two SNP representing silent substitutions within the coding region, and six SNP representing intron variation were selected out of 38 total SNP heterozygous in the Piedmontese x Angus sire. Six of the SNP selected reside on the half of the gene 5' to the approximately 100-kb intron 10, while the remaining four are located in the half 3' to intron 10 (Figure 1Go). The haplotypes inherited by the sire from the Piedmontese grandsire and Angus grandam were identified by inferring the haplotypes based on offspring that were homozygous for all of the SNP tested and comparing these genotype patterns with markers used in the QTL analysis. All offspring from the MARC family were genotyped for the 10 SNP, offspring from the New Zealand family were genotyped for eight of these SNP informative in the Jersey x Limousin sire.



View larger version (12K):
[in this window]
[in a new window]
 
Figure 1. Schematic of the calcium-activated neutral protease gene. The boxes represent exon sequence and the connecting lines represent the intron sequence. The approximately 100-kb intron 10 is shown as a break in the sequence. The open boxes represent coding region and the hatched boxes represent untranslated sequence. The black boxes located below the schematic of the gene indicate sequence generated from this study. Contigs 1 and 2 from this study are indicated and contig 3 is located within the ~100-kb intron. The domains of the protein are marked by roman numerals I to IV and the set of 10 single-nucleotide polymorphisms (SNP) for haplotype construction are shown as stars. The first star on the left is SNP 1 followed by the rest of the SNP in numerical order. The first two SNP were included only in the U.S. Meat Animal Research Center quantitative trait locus population whereas 3 to 10 were included in both resource populations.

 
Dam Genotypes

Dam alleles were predicted based on the sire’s haplotype and the SNP genotypes of the progeny. The effect of the allele inherited from the dam was estimated simultaneously with the effect of the paternally inherited haplotype in the MARC population using multiple regression of SAS (SAS Inst. Inc., Cary, NC). Variation associated with the maternally inherited alleles was partitioned separately from that associated with paternal alleles in order to isolate novel variation from previously reported differences associated with paternal alleles (Casas et al., 2000). The two SNP predicting valine or isoleucine (V530/I530) and glycine or alanine (G316/A316) were evaluated. The analysis provided least squares means by using a model that included the effects of sex (steers or heifers) and days on feed as a covariate.


    Results
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Implications
 References
 
Bacterial Artificial Chromosome Sequencing.

The first step in the evaluation of CAPN1 as a candidate was to obtain sufficient sequence to enable identification of SNP within the entire coding region of the bovine gene. Partial sequence (11,055 bp) of the bovine CAPN1 gene (Smith et al., 2000a) was available; however, the reported CAPN1 sequence excluded sequences between exon 9 and exon 11. To obtain this sequence, a subclone library was made from a BAC containing the CAPN1 gene. Clones were sequenced from both ends to generate 5,947 sequence reads, a predicted 10-fold coverage of the BAC. Analysis of the sequences assembled 4,133 of the reads into 362 sequence contigs, with the remainder being low quality sequence reads or failing to show overlap with other sequences.

The three largest bovine contigs (containing 1,776 reads), were examined by BLAST analysis. The first contig contained exons 1 to 10 of the bovine CAPN1 gene, the second contig contained exons 11 to 22 of the bovine CAPN1 gene, and the third contig shared homology to the human CAPN1 intron 10 (estimated to be approximately 100 kb based on the draft human sequence). The three contigs containing CAPN1 genomic sequence represented a total length of 43,848 bp. No matches or overlaps were found from the remaining contigs that added sequence information.

The two contigs containing coding portions of the bovine CAPN1 gene added significant new sequence, relative to the previous report by Smith et al. (2000a), within and around the exons as indicated by filled rectangles in Figure 1Go. Furthermore, two significant aspects of CAPN1 gene architecture were revealed that had not been identified in the previous study. The additional sequence upstream of the 100-kb intron revealed a 130-bp intron in exon 9, that had not been detected by the exon primer-based PCR approach used previously (Smith et al., 2000a). Thus, the portion of the bovine CAPN1 gene upstream from the 100-kb intron includes 10 exons, rather than the nine as previously reported. The second aspect involved the first intron downstream of the large intron which had been reported to be approximately 4 kb, but the new sequence data demonstrated that this intron is divided into a 2,791-bp intron, 12-bp exon, and an 891-bp intron (Figure 1Go). Therefore, the half of the gene downstream of the large intron includes 12 exons, rather than 10 as previously reported. In summary, the data indicate that the bovine CAPN1 gene consists of 22 exons and 21 introns. This is the same arrangement observed for the human gene by comparison of the human cDNA (accession NM_005186) and the draft genome (Lander et al., 2001) sequences (data not shown). The updated sequence information was sufficient to support primer design to scan all exons of the gene. No effort was made to completely sequence the ~100-kb intron 10. In all, more than 28 kb of sequence containing the entire coding region of the bovine CAPN1 gene was determined.

Single-Nucleotide Polymorphism Iidentification in MARC Quantitative Trait Locus Population.

The hypothesis that CAPN1 variation underlies the QTL effect observed in the MARC resource population predicts the presence of sequence variation between the two alleles of the Piedmontese x Angus sire of the population (sire PA). To identify this variation, a set of PCR primers was designed to amplify and sequence the coding region of the gene (Table 1Go). A total of 12.5 kb of the gene was encompassed by 35 amplicons. All 22 exons were completely sequenced to identify coding region variation, and all but the two largest introns were included to identify sufficient numbers of SNP to support construction of haplotypes. This procedure identified 38 SNP in sire PA, of which five were in exons. Two of the five exon SNP predicted amino acid variation in the µ-calpain protein, including a GCC alanine codon in exon 9 vs a GGC glycine codon (amino acid number 316), and a GTC valine codon in exon 14 vs an ATC isoleucine codon (amino acid number 530).

A set of microsatellite markers had been used previously to define the sire allele with effect on shear force measurement of meat tenderness. However, mapping the CAPN1 gene to the QTL interval (Smith et al., 2000a) did not define which allele of the gene would associate with decreased meat tenderness (increased shear force), nor did the previous analysis establish haplotypes to permit this analysis. Therefore, the haplotypes of the sire were determined by genotyping of the progeny to identify homozygotes of each haplotype and permit unambiguous assignment of SNP alleles to each haplotype. This analysis demonstrated that one sire allele contained alanine at position 316 (allele A316) and valine at position 530 (allele V530), whereas the other sire allele contained glycine (allele G316) and isoleucine (allele I530) at those same two positions. A set of 10 SNP heterozygous in sire PA were chosen to genotype the entire resource population, to determine the haplotype associated with increased shear force. The SNP chosen and their alternative alleles in this bull are shown in Table 2Go, with the positions in the gene shown graphically in Figure 1Go. We arbitrarily define allele 1 of this sire as the GTTTCGATCG haplotype containing G316 and I530, and allele 2 as the ACCCTCGCTA haplotype containing A316 and V530 with each base designation representing the allele of each SNP in the order shown in Table 2Go. Alleles 1 and 2 were inherited from the Piedmontese and Angus parent of the sire, respectively. Consistent with the original microsatellite-based QTL analysis, the haplotype analysis revealed that the Piedmontese allele had an effect on meat tenderness in this population (Table 3Go). This result suggests that the G316 I530 haplotype is a marker for a functional allele of CAPN1 associated with increased shear force in this population.


View this table:
[in this window]
[in a new window]
 
Table 2. Ten SNP used to generate the haplotypes for the two resource population sires, in the order of location starting from the 5' end of the calcium-activated neutral protease gene1
 

View this table:
[in this window]
[in a new window]
 
Table 3. Least squares mean, standard error, and probability values for meat tenderness measured as Warner-Bratzler shear force (kg) at 3 and 14 d postmortem
 
An alternative hypothesis for the observed effect is that CAPN1 markers are in linkage disequilibrium with a functional allele of another gene mapping to this part of BTA29. In this case, recombination will eventually break the association of the functional allele with CAPN1 haplotypes. To test whether CAPN1 alleles are the source of the meat tenderness variation, the contribution of alleles from the dams of the population were examined. If one or both of the predicted amino acid variations are markers for functional alleles with respect to shear force, then maternal inheritance should have a similar impact on phenotype. On the other hand, the association will be unlikely if the functional allele is in a different gene, since the dams are distantly related to the Piedmontese breed and an ample number of generations has occurred for recombination of the CAPN1 alleles and putative functional alleles of other genes. The dam haplotypes were inferred by subtraction of the known sire haplotype from the genotypes of the calves, and analyzed simultaneously with the sire allele as a fixed effect. Table 3Go summarizes the analyses of the sire alleles and the dam alleles with regard to shear force mean differences and their significance as measured by the P value. The magnitude of effect is measured by both shear force mean differences in the two alleles tested and the calculated effect in standard deviations. The dam alleles containing I530 or V530 were analyzed separately from the dam alleles containing G316 or A316 as numbers of dams containing inferred genotypes for both markers were insufficient for haplotype analysis. The G316 allele inherited from the dams was present at 86% frequency and did not show an effect on meat tenderness. The I530 allele was present in the dams at 18% frequency and showed an effect on meat tenderness (increase in shear force d 14, P < 0.04), providing strong support for the hypothesis that this SNP is in linkage disequilibrium with a functional allele.

Single Nucleotide Polymorphism Identification in AgResearch Quantitative Trait Locus Population.

A second resource population at AgResearch in New Zealand had been developed that showed evidence of a QTL on BTA29 (C. Morris, unpublished data). A reciprocal backcross population using a Jersey x Limousin sire (sire JL) revealed a significant QTL approximately centered at the map position of CAPN1, affecting MIRINZ tenderometer shear force. This resource population had a different breed composition and was raised in a different environment (New Zealand vs Nebraska), by a different production protocol, suggesting it represented a rigorous test for the association of particular CAPN1 haplotypes with meat tenderness across populations. The CAPN1 gene of the Jersey x Limousin sire was sequenced. No additional variation within the coding region was identified that would be predicted to alter the amino acid sequence of the protein. Genotypes from the 10 SNP used in the MARC population demonstrated that sire JL was heterozygous for eight of the 10, including the two predicting amino acid variation. Genotypes for the eight informative SNP were collected on the progeny of the resource population to define haplotypes of the bull and to determine which haplotype had the effect of increased shear force. The two haplotypes are shown in Table 2Go, with allele 1 arbitrarily defined as the allele containing G316 and I530. Use of these SNP markers in a QTL analysis of the New Zealand population demonstrated that allele 1 of sire JL is associated with decreased meat tenderness (i.e., increased shear force, P < 0.0001) in this Jersey x Limousin cross (Table 4Go). Allele 1 was derived from the Limousin grandsire and is highly similar to allele 1 of the MARC sire, providing support for the hypothesis that these SNP represent markers of functional alleles of CAPN1.


View this table:
[in this window]
[in a new window]
 
Table 4. Least squares mean, standard error, and probability values for meat tenderness measured as MIRINZ Tenderometer shear force (kg) at 1.3, 2.0, and 2.3 d postmortem
 
The dam contributions for the two amino acid variants were inferred using the known sire contributions to the progeny. The shear force values as measured by the MIRINZ tenderometer test and results of the association analyses for both the dam and sire contributions to the progeny are summarized in Table 4Go. The analysis of the dam contributions revealed an effect for the G316 allele similar to that observed in the sire alleles of both resource populations. No effect was observed for the I530 allele when the dam contributions were analyzed as independent variables.

The identification of two µ-calpain variants that show effects on meat tenderness in unrelated QTL populations led to examination of the gene in a wider source of germplasm to ascertain the depth of CAPN1 diversity in cattle. The same primers previously used to amplify overlapping portions of the gene were applied to sequence a sampling of 16 beef breeds and the Holstein dairy breed represented in the MARC bovine diversity panel (described in Heaton et al., 2001). An additional 134 SNP were detected in this panel, but none of the eleven present in exons were predicted to result in amino acid variation in the protein, nor were there any other obviously significant changes to the gene. The I530 allele was observed in all of the breeds in the diversity panel except Brangus and was present at an overall frequency of 30%. The G316 allele was observed in all of the breeds represented by the diversity panel and was present at an overall frequency of 86%. The A316 allele was not observed in Brahman, Simmental, Gelbvieh, Salers, Maine-Anjou, or Chianina breeds.


    Discussion
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Implications
 References
 
We have characterized variation in the bovine CAPN1 gene to identify markers for alleles affecting meat tenderness. In the Piedmontese x Angus sire of the MARC resource population, a haplotype was defined by 10 SNP that were representative of the variation observed by virtue of exonic and intronic SNP lying in both halves of the gene separated by an intron of over 100 kb. This group of SNP successfully determined the sire contribution to progeny, even in animals heterozygous at multiple positions and in the absence of dam genotypes. The initial study indicated that the paternally inherited haplotype containing G316 and I530 was associated with increased shear force in the MARC resource population. This conclusion is limited, since these were the same animals as the original quantitative study and, therefore, the two sire alleles of any gene in the QTL interval would be expected to show contrast in the population. However, the dam alleles of the resource population represented an ideal resource to address the possibility that the observed CAPN1 sequence variation represents useful markers for functional alleles affecting meat tenderness. The dams were from a different genetic background than the sire (Piedmontese x Angus sire vs Hereford, Angus, Pinzgauer, and Red Poll composite dams), making it unlikely that an association of phenotype with CAPN1 markers in the dam alleles would occur as a result of population stratification or breed effect. Subtraction of the sire haplotype from the progeny genotypes was an efficient means to determine the dam allele contribution, even in the absence of DNA samples from the dams.

No association with decreased meat tenderness was detected for the G316 allele contributed from the dams in the MARC resource population, in contrast to the effect observed in the sire alleles. The low frequency of A316 in the dams resulted in many of the progeny inheriting a V530 and G316 allele from the dam which is out of phase with the inheritance pattern from the sire. The inheritance of these two opposing alleles from the dam limits the conclusion that can be made when an effect is not observed. The I530 allele, however, showed a similar effect on shear force when the dam alleles were analyzed as had been observed in the sire alleles, providing strong support for the utility of this marker. Further strong support is provided by the discovery that the I530 allele also is correlated with increased shear force in the AgResearch population, which is comprised of different breeds, raised on a separate continent, and with different management variables. The association of I530 with decreased tenderness in three different genetic backgrounds as demonstrated by the MARC sire (Piedmontese x Angus), the MARC dams (Red Poll, Pinzgauer, Hereford, Angus) and the AgResearch sire (Jersey x Limousin) provides strong evidence for I530 as a functional marker. The association of the G316 allele with increased shear force in the AgResearch population was demonstrated when inherited from either the sire or the dams providing support for the utility of this marker in addition to the I530 marker. The association of G316 with increased shear force is demonstrated in three different genetic backgrounds as demonstrated by the MARC sire (Piedmontese x Angus), the AgResearch sire (Jersey x Limousin) from Australia, and the AgResearch dams (Jersey and Limousin) from New Zealand.

Although differences were not detected in the AgResearch dams between the I530 and V530 alleles, the low frequency of dams contributing I530 in the AgResearch population may have limited the ability to detect a difference in shear force values. In addition, tenderness measurements were taken only up to d 4 postmortem in this population. It is possible that differences are only detectable at d 14 postmortem as demonstrated by the MARC resource population. Similarly the low frequency of A316 in the MARC resource population may have hindered the detection of a difference between A316 and G316 alleles in the MARC dams.

The overall results provide strong evidence that I530 and G316 are informative markers for meat tenderness variation. Both I530 and V530 alleles were detected in all breeds present in the diversity panel except for Brangus (note the panel has only four Brangus bull samples, and the Brahman and Angus bull samples had both alleles), with no breed appearing fixed for the I530 allele. The observed frequency of I530 across the entire panel was 30%, sufficiently high to suggest that selection could impact variation in meat tenderness. The observed frequency of G316 across the entire panel was 86% with Brahman, Simmental, Gelbvieh, Salers, Maine-Anjou, and Chianina appearing fixed for the G316 allele. However, the low number of samples from each breed does not support the precise calculation of frequency within breeds.

In addition to utility as markers, it is possible that one or both of the amino acid variations reflect a mutation that causes a functional change in the µ-calpain protease. The µ-calpain isoform containing I530, G316, or both may be a functionally different protein that led to variation in myofibrillar proteolysis and resulted in a difference in shear force in the two QTL populations. The I530/V530 variation represents a relatively conservative substitution of nonpolar amino acids in Domain III of the protein, which currently has no known function in terms of the activity of the enzyme. However, it could potentially alter protein folding or stability, impacting autocatalysis or proteolytic activity. Moreover, a major change in the protein is unlikely to be necessary to effect the relatively small change in postmortem proteolysis underlying the variation in shear force measurement. A peptide search using NCBI BLAST of bovine Domain III revealed that human, mice, rat, rabbit, chicken, and monkey versions of µ-calpain all have isoleucine at position 530. In contrast, sequences for swine, sheep, zebrafish, and frog all have valine at this position, which is the most common allele in the beef diversity panel. The region of comparison was well conserved among species (> 90%) except for chicken (66%) and zebrafish (56%). More sequencing of these other species is necessary before any conclusions could be made regarding these potential differences. Similar to the I530 allele, the variation represented by the G316 and A316 alleles represent a relatively conservative substitution of non-polar amino acids. However, this variation is located in Domain II of the protein, which has been identified as the proteolysis domain. An alteration in the proteolysis domain could alter the activity of the protein, and thus effect meat tenderization. A peptide search using the region around this variation revealed that monkey, human, mice, rat, chicken, and frog versions of µ-calpain do not contain glycine or alanine. The region in general is less conserved (< 90%) for the species mentioned, and the presence of glycine, alanine, serine, asparagine, lysine, and glutamic acid demonstrate that this position is not highly conserved among species.

Our results demonstrate the association of the I530 and G316 alleles with meat tenderness in two separate populations as well as in the dam alleles of the MARC population and the AgResearch population, and are consistent with the hypothesis that variation in CAPN1 is responsible for the detection of the BTA29 QTL in the two resource populations.


    Implications
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Implications
 References
 
This report presents evidence that single nucleotide polymorphisms determining amino acid variation of valine or isoleucine at position 530 and glycine or alanine at position 316 in the micromolar calcium-activated neutral protease gene, are markers for meat tenderness variation. If the alternative genotypes can successfully predict meat tenderness phenotypes in larger, unrelated populations, it would represent the first example of a deoxyribonucleic acid test for genetic merit directly related to this trait and would represent a valuable tool for cattle breeders to improve tenderness in cattle.


    Footnotes
 
1 The authors thank K. Tennill, R. Godtel, L. Flathman, and S. Simcox for technical assistance; S. Kluver for manuscript preparation; and D. Nonneman and R. Stone for manuscript review. Back

2 Mention of trade names, proprietary products, or specific equipment does not constitute a guarantee or warranty of the product by the USDA and does not imply its approval to the exclusion of other products that also may be suitable. Back

Received for publication May 30, 2002. Accepted for publication July 16, 2002.


    References
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Implications
 References
 


Altschul, S. F., W. Gish, W. Miller, E. W. Myers, and D. J. Lipman. 1990. Basic local alignment search tool. J. Mol. Biol. 215:403–410.[Medline]

Altschul, S. F., T. L. Madden, A. A. Schaffer, J. Zhang, Z. Zhang, W. Miller, and D. J. Lipman. 1997. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25:3389–3402.[Abstract/Free Full Text]

AMSA. 1995. Research guidelines for cookery, sensory evaluation and instrumental tenderness measurements of fresh meat. Am. Meat Sci. Assoc., Chicago, IL.

Casas, E., R. T. Stone, J. W. Keele, S. D. Shackelford, S. M. Kappes, and M. Koohmaraie. 2001. A comprehensive search for quantitative loci affecting growth and carcass composition of cattle segregating alternative forms of the myostatin gene. J. Anim. Sci. 79:854–860.[Abstract/Free Full Text]

Casas, E., S. D. Shackelford, J. W. Keele, R. T. Stone, S. M. Kappes, and M. Koohmaraie. 2000. Quantitative trait loci affecting growth and carcass composition of cattle segregating alternate forms of myostatin. J. Anim. Sci. 78:560–569.[Abstract/Free Full Text]

Casas, E., J. W. Keele, S. D. Shackelford, M. Koohmaraie, T. S. Sonstegard, T. P. L. Smith, S. M. Kappes, and R. T. Stone. 1998. Association of the muscle hypertrophy locus with carcass traits in beef cattle. J. Anim. Sci. 76:468–473.[Abstract/Free Full Text]

Ewing, B., and P. Green. 1998. Base-calling of automated sequencer traces using Phred. II. Error probabilities. Genome Res. 8:186–194.[Abstract/Free Full Text]

Ewing, B., L. Hillier, M. C. Wendl, and P. Green. 1998. Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 8:175–185.[Abstract/Free Full Text]

Fraserhurst, L., and P. MacFarlane. 1983. A device for measuring the tenderness of meat. N. Z. Patent 190945, Auckland, NZ.

Geesink, G. H., and M. Koohmaraie. 1999. Effect of calpastatin on degradation of myofibrillar proteins by µ-calpain under postmortem conditions. J. Anim. Sci. 77:2685–2692.[Abstract/Free Full Text]

Gingrich, J. C., D. M. Boehrer, and S. B. Basu. 1996. Partial CviJI digestion as an alternative approach to generate cosmid sublibraries for large-scale sequencing projects. Biotechniques 21:99–104.[Medline]

Gordon, D., C. Abajian, and P. Green. 1998. Consed: a graphical tool for sequence finishing. Genome Res. 8:195–202.[Abstract/Free Full Text]

Heaton, M. P., C. G. Chitko-McKnown, W. M. Grosse, J. W. Keele, J. E. Keen, and W. W. Laegreid. 2001. Interleukin-8 haplotype structure from nucleotide sequence variation in commercial populations of U.S. beef cattle. Mamm. Genome 12:219–226.[Medline]

Koohmaraie, M. 1992. The role of Ca2+-dependent proteases (calpains) in postmortem proteolysis and meat tenderness. Biochemie 74: 239–245.[Medline]

Koohmaraie, M. 1994. Muscle proteinases and meat aging. Meat Sci. 36:93–104.

Koohmaraie, M. 1996. Biochemical factors regulating the toughening and tenderization process of meat. Meat Sci. 43:S193–S201.

Lander, E. S., J. Rogers, R. H. Waterston, T. Hawkins, R. A. Gibbs, et al. 2001. Initial sequencing and analysis of the human genome. Nat. Genet. 409:860–922.

Morris, C. A., N. G. Cullen, S. M. Hickey, A. M. Crawford, D. L. Hyndman, C. D. K. Bottema, and W. S. Pitchford. 2001. Progress in DNA marker studies of beef carcass composition and meat quality in New Zealand and Australia. In: Proc. Assoc. Advancement Anim. Breed. Genet. Queenstown, NZ. 14:17–22.

Pitchford, W. S., A. L. Ewers, R. W. Ponzoni, and M. P. B. Deland. 1998. Breed and sire effects on saleable beef yield. Proc. 6th World Cong. Genet. Appl. Livest. Prod., Armidale, NSW, Australia. 23:117–120.

Smith, T. P. L., E. Casas, C. E. Rexroad III, S. M. Kappes, and J. W. Keele. 2000a. Bovine CAPN1 maps to a region of BTA29 containing a quantitative trait locus for meat tenderness. J. Anim. Sci. 78:2589–2594.[Abstract/Free Full Text]

Smith, T. P. L., R. A. Godtel, and R. T. Lee. 2000b. PCR-Based setup for high-throughput cDNA library sequencing on the ABI 3700 automated DNA sequencer. Biotechniques 29:698–700.[Medline]

Stone, R. T., W. M. Grosse, E. Casas, Smith T. P. L., J. M. Keele, and G. L. Bennett. 2002. Use of bovine EST data and human genetic sequence to map 100 gene-specific bovine markers. Mamm. Genome 13:211–215.[Medline]

Warren, W., T. P. L. Smith, C. E. Rexroad III, S. C. Fahrenkrug, T. Allison, C. L. Shu, J. Catanese, and P. J. de Jong. 2000. Construction and characterization of a new bovine bacterial artificial chromosome library with 10 genome-equivalent coverage. Mamm. Genome 11:662–663.[Medline]

Wheeler, T. L, S. D. Shackelford, and M. Koohmaraie. 1998. Cooking and palatability traits of beef longissimus steaks cooked with belt grill on an open hearth electric broiler. J. Anim. Sci. 76:2805–2810.[Abstract/Free Full Text]


This article has been cited by other articles:


Home page
J ANIM SCIHome page
X. Jia, E. Veiseth-Kent, H. Grove, P. Kuziora, L. Aass, K. I. Hildrum, and K. Hollung
Peroxiredoxin-6--A potential protein marker for meat tenderness in bovine longissimus thoracis muscle
J Anim Sci, July 1, 2009; 87(7): 2391 - 2399.
[Abstract] [Full Text] [PDF]


Home page
Endocr. Rev.Home page
B. D. Rodgers and D. K. Garikipati
Clinical, Agricultural, and Evolutionary Biology of Myostatin: A Comparative Review
Endocr. Rev., August 1, 2008; 29(5): 513 - 534.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
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]


Home page
GeneticsHome page
W. Barendse, A. Reverter, R. J. Bunch, B. E. Harrison, W. Barris, and M. B. Thomas
A Validated Whole-Genome Association Study of Efficient Food Conversion in Cattle
Genetics, July 1, 2007; 176(3): 1893 - 1905.
[Abstract] [Full Text] [PDF]


Home page
J ANIM SCIHome page
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 page
J ANIM SCIHome page
E. Casas, S. N. White, T. L. Wheeler, S. D. Shackelford, M. Koohmaraie, D. G. Riley, C. C. Chase Jr., D. D. Johnson, and T. P. L. Smith
Effects of calpastatin and {micro}-calpain markers in beef cattle on tenderness traits
J Anim Sci, March 1, 2006; 84(3): 520 - 525.
[Abstract] [Full Text] [PDF]


Home page
J ANIM SCIHome page
R. T. Stone, E. Casas, T. P. L. Smith, J. W. Keele, G. Harhay, G. L. Bennett, M. Koohmaraie, T. L. Wheeler, S. D. Shackelford, and W. M. Snelling
Identification of genetic markers for fat deposition and meat tenderness on bovine chromosome 5: Development of a low-density single nucleotide polymorphism map
J Anim Sci, October 1, 2005; 83(10): 2280 - 2288.
[Abstract] [Full Text] [PDF]


Home page
J ANIM SCIHome page
S. N. White, E. Casas, T. L. Wheeler, S. D. Shackelford, M. Koohmaraie, D. G. Riley, C. C. Chase Jr., D. D. Johnson, J. W. Keele, and T. P. L. Smith
A new single nucleotide polymorphism in CAPN1 extends the current tenderness marker test to include cattle of Bos indicus, Bos taurus, and crossbred descent
J Anim Sci, September 1, 2005; 83(9): 2001 - 2008.
[Abstract] [Full Text] [PDF]


Home page
J ANIM SCIHome page
E. Casas, S. N. White, D. G. Riley, T. P. L. Smith, R. A. Brenneman, T. A. Olson, D. D. Johnson, S. W. Coleman, G. L. Bennett, and C. C. Chase Jr.
Assessment of single nucleotide polymorphisms in genes residing on chromosomes 14 and 29 for association with carcass composition traits in Bos indicus cattle
J Anim Sci, January 1, 2005; 83(1): 13 - 19.
[Abstract] [Full Text] [PDF]


Home page
J ANIM SCIHome page
B. T. Page, E. Casas, R. L. Quaas, R. M. Thallman, T. L. Wheeler, S. D. Shackelford, M. Koohmaraie, S. N. White, G. L. Bennett, J. W. Keele, et al.
Association of markers in the bovine CAPN1 gene with meat tenderness in large crossbred populations that sample influential industry sires
J Anim Sci, December 1, 2004; 82(12): 3474 - 3481.
[Abstract] [Full Text] [PDF]


Home page
J ANIM SCIHome page
E. Casas, S. D. Shackelford, J. W. Keele, M. Koohmaraie, T. P. L. Smith, and R. T. Stone
Detection of quantitative trait loci for growth and carcass composition in cattle
J Anim Sci, December 1, 2003; 81(12): 2976 - 2983.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Page, B. T.
Right arrow Articles by Smith, T. P. L.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Page, B. T.
Right arrow Articles by Smith, T. P. L.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS