A single nucleotide polymorphism in CAPN1 associated with marbling score in Korean cattle
BMC Genetics volume 9, Article number: 33 (2008)
Marbling score (MS) is the major quantitative trait that affects carcass quality in beef cattle. In this study, we examined the association between genetic polymorphisms of the micromolar calcium-activated neutral protease gene (micro-calpain, CAPN1) and carcass traits in Korean cattle (also known as Hanwoo).
By direct DNA sequencing in 24 unrelated Korean cattle, we identified 39 sequence variants within exons and their flanking regions in CAPN1. Among them, 12 common polymorphic sites were selected for genotyping in the beef cattle (n = 421). Statistical analysis revealed that a polymorphism in the 3'UTR (c.2151*479C>T) showed significant association with MS (Pcor.= 0.02).
Our findings suggest that polymorphisms in CAPN1 might be one of the important genetic factors involved in carcass quality in beef cattle, although it could be false positive association.
Genetic improvement has long been considered an important factor in the competitiveness of beef cattle production. Identification of the genes and/or polymorphisms underlying quantitative/qualitative traits, and an understanding of how these genes/polymorphisms interact with the environment or with other genes affecting economic traits might be the keys to successful application of marker-assisted selection in the commercial animal population. As one of these economic traits, marbling is intramuscular fat that gives meat flavor and tenderness. Thus, an increase in the degree of marbling raises the level of meat quality.
Calpain is a ubiquitous cytoplasmic cysteine protease, the activity of which is absolutely dependent on calcium . Two genes of calpain (CAPN1 [macro-calpain] and CAPN2 [mili-calpain]) have been identified . CAPN1 degrades myofibrillar proteins under postmortem conditions and appears to be the primary enzyme in the postmortem tenderization process [3–6]. Regulation of CAPN1 activity has been correlated with variation in meat tenderness, and previous studies also identified a quantitative trait locus influencing meat tenderness on chromosome 29 where CAPN1 lies [7–9].
In the CAPN1 gene, more than 100 single nucleotide polymorphisms (SNPs) have been identified in Bos indicus or Bos taurus [8, 10, 11] (S. N. White and T. Smith, unpublished data). Among them, four polymorphisms, two non-synonymous SNPs (G316A and V530I), and two intronic SNPs (C4685T and C4751T) have been found to have significant effects on meat tenderness [10, 12–16].
Association studies of SNPs in CAPN1 with carcass traits have been performed, but there is no publicly available evaluation of the association of these SNPs in Korean cattle (Bos taurus coreanae, also known as Hanwoo). Thus, the objective of this study was to discover polymorphisms and assess the association, including that of reported SNPs, in Korean cattle.
Results and Discussion
By direct DNA sequencing, 39 polymorphisms were identified within exons and their flanking regions of CAPN1: 6 in coding exons, 10 in 3'UTR, and 23 in introns. The locations and allele frequencies of the polymorphisms are shown in Table 1 and Figure 1A. Among those identified, 12 polymorphisms (c.579G>A [K193K], c.630A>G [T210T], c.760-24T>C, c.843+330A>G, c.1199G>A [R400Q], c.1588G>A [V530I], c.1611+104C>T, c.1869+235C>G, c.2151*479C>T, c.2151*765A>G, c.2151*832G>A, and c.2151*845A>G) were selected for larger-scale genotyping based on location (polymorphisms in exons were preferred), minor allele frequency exceeding 0.05, LD (a polymorphism was chosen if it was in absolute LD [r2 = 1] with one or more other polymorphisms), and haplotype tagging status (Table 2 and Figure 1A). The minor allele frequencies of genotyped polymorphisms in Korean cattle (n = 421) are described in Table 1. Pairwise linkage disequilibrium analysis with the 12 polymorphisms showed that CAPN1 can be parsed into two LD blocks (Block1 and Block2), the 1.5 kb region spanning from exon5 to intron7 (Block1) and the 9 kb region from exon11 to exon22 (Block2). The distance between two adjacent SNPs (c.843+330A>G in Block1 and c.1199G>A in Block2) is approximately 13 kb (Figures 1A and 1C). There were three and five common haplotypes (freq. > 0.05) in Block1 and Block2, respectively (Figures 1B and 1C).
Associations of CAPN1 polymorphisms with cold carcass weight (CW) and marbling score (MS) were analyzed using the mixed-effect model. The obtained P- values were corrected for multiple testing using the Bonferroni correction method. Two common haplotypes (Block2_ht2 and Block2_ht3) were also tested with the same covariates in a similar manner (Table 2). Block1_ht1, Block1_ht2, Block1_ht3, Block2_ht1, Block2_ht4, and Block2_ht5 were not analyzed because they were equivalent with c.630A>G (T210T), c.579G>A (K193K), c.843+330A>G, c.760-24T>C, c.1869+235C>G, c.2151*832G>A, and c.2151*479C>T, respectively (Figure 1B).
No positive associations were detected with CW. However, one polymorphism (c.2151*479C>T) in 3'UTR showed significant associations with MS. The "T" allele of c.2151*479C>T revealed an additive effect on MS, i.e., the lowest MS was found in "T/T" (MS = 0.94), intermediate in "C/T" (MS = 1.56), and the highest in "C/C" (MS = 2.34) (Pcor.= 0.02). In further haplotype association analysis, no positive associations were detected with CW and MS (Table 2).
Four polymorphisms, G316A, V530I, C4685T and C4751T have shown significant effects on meat tenderness [10, 12–16]. However, one potential SNP marker with a reported effect on meat tenderness, specifically G316A in exon 9 of the CAPN1 gene, which produces a substitution of alanine for glycine at position 316 of the CAPN1 protein, was homozygous in our population (data not shown). In addition, V530I, C4685T and C4751T were discovered in Korean cattle (c.1588G>A [V530I], c.1611+104C>T and c.1800+169C>T respectively), and C4685T and C4751T are in absolute LD with each other in this study (Figure 1A). But these SNPs showed no association with MS. Interestingly, the results of our study showed that one polymorphism (c.2151*479C>T) in 3'UTR was associated with MS.
Although the mechanisms involved in the association of polymorphisms in the UTR region with MS are not currently understood, our study clearly indicates that c.2151*479C>T has an influence on MS in Korean cattle in light of the fact that strong association (Pcor.= 0.02), high enough to overcome the correction for multiple testing, was detected. We can hypothesize that polymorphisms in the UTR region might impact gene function by affecting the mRNA stability as well as regulatory motifs within UTRs. In addition, although it is true that the allele itself may be functional and directly affect the expression of the phenotype, a more likely event is that the allele is in linkage disequilibrium with another allele at a nearby locus that is the true causal allele.
Our results may be a false positive test because there is no quantitative trait locus (QTL) evidence for marbling in the respective genome region and CAPN1 is not a convincing functional candidate gene for marbling. However, similar to other association studies using the candidate gene strategies, although not complete, would facilitate further studies and contribute to understanding the genetic background of important traits. Further biological and/or functional evidence is needed to confirm the genetic effects of CAPN1 polymorphisms on MS.
In summary, we have identified 39 polymorphisms in CAPN1 in Korean cattle, and 12 common polymorphic sites were selected for genotyping. Statistical analysis revealed that one polymorphism (c.2151*479C>T) in 3'UTR showed significant associations with the important carcass quality trait, marbling score, in Korean cattle.
Animals and phenotypic data
The Korean native cattle genomic DNA samples were obtained from 421 steers produced from 76 sires used in the progeny testing program of the National Livestock Research Institute (NLRI) of Korea. All steers were fed for 731.39 ± 16.53 days under a tightly controlled feeding program in the Daekwanryeong and Namwon branches. The average birth date was May 10 [range: April 15 – June 3]. After two years, animals were slaughtered in May 3 to 18. They were weaned at a mean age of 3 months and fed until they were 6 months old with 30% concentrates and 70% roughage. After 6 months of age, they were fed with concentrates consisting of 15% crude protein (CP)/71% totally digestible nutrients (TDN) until they were 14 months old; 13% CP/72% TDN until 20 months; and 11% CP/73% TDN until 24 months of age. The roughage was offered ad libitum, and steers had free access to fresh water during the whole period. Live weights were determined before slaughter. The mean of the live weights was 539.93 ± 51.96 kg. Yield grades for carcasses were determined by cold carcass weight (CW). After a 24-h chill, CW was measured, and then the left side of each carcass was cut between the last rib and the first lumbar vertebra to determine marbling score (MS). The mean of the CW was 311.47 ± 33.39 kg. MS was determined by assessing the degree of marbling in the cut surface of the ribeye. The degree of marbling was evaluated according to the Korean Beef Marbling Standard (1 = trace, 7 = very abundant) . The mean of the MS was 2.25 ± 1.36.
Sequencing analysis of CAPN1
We have sequenced exons and their flanking regions to discover variants in 24 unrelated Korean native cattle using the ABI PRISM 3730 DNA analyzer (Applied Biosystems, Foster City, CA). Sixteen primer sets for the amplification and sequencing analysis were designed based on GenBank sequences (Ref. Genome seq. NC_007330 released on 26 Dec. 2006 and Ref. mRNA seq. AF221129 released on 9 Feb. 2000). Primer information is given in supplementary data. Sequence variants were verified by chromatograms [see Additional file 1].
Genotyping by single-base extension (SBE) and electrophoresis
For genotyping of polymorphic sites, amplifying and extension primers were designed for single-base extension (SBE) . Primer extension reactions were performed with the SNaPshot ddNTP Primer Extension Kit (Applied Biosystems). To stop the primer extension reaction, one unit of SAP (shrimp alkaline phosphatase) was added to the reaction mixture, and the mixture was incubated at 37°C for 1 hour, followed by 15 min at 72°C for enzyme inactivation. The DNA samples containing extension products and GeneScan 120 Liz size standard solution were added to Hi-Di formamide (Applied Biosystems) according to the recommendation of the manufacturer. The mixture was incubated at 95°C for 5 min, followed by 5 min on ice, and then electrophoresis was performed using the ABI Prism 3100 Genetic Analyzer. The results were analyzed using the programs ABI Prism GeneScan and Genotyper (Applied Biosystems). Probe information is given in supplementary data [see Additional file 1].
The χ2 tests were used to determine whether the individual variant was in equilibrium at each locus in the population (Hardy-Weinberg equilibrium). We examined a widely used measure of linkage disequilibrium between all pairs of biallelic loci, D' (the correlation coefficient [Delta, |D'|]), LOD (logarithm of odds), and r2. Strength of LD between pairs of SNPs was measured as D' using Haploview. Regions of strongly associated markers (LD blocks) were inferred by Gabriel's method as implemented in Haploview [19, 20]. Gabriel's method defines pairs of SNPs to be in strong LD if the one-sided 95% D' confidence bound is between 0.7 and 0.98. The method defines a block if 95% of pairwise SNP comparisons are in strong LD. r2 was also used to determine whether the pairs of sites were in absolute LD. Haplotypes and their frequencies were inferred using the algorithm developed by Stephens et al. . Phase probabilities of each site were calculated for each individual using this software (PHASE) (input option: ignoring families). Phase probabilities of all polymorphic sites for haplotypes were calculated for each individual using this software. Because 95% of samples had phase probabilities greater than 97%, 97% was chosen as threshold for phase probability. Associations between individual SNPs and carcass traits (CW and MS) were determined by the mixed effect model, treating "sire" as a random effect, and "age" at slaughter was also included in the model as a covariate using the library (nlme) in the R statistical package (see Availability and requirements for URL). Other covariates were not available for this analysis. We used a single SNP model. Single SNP/haplotype effects were tested in the mixed effect model. P-values were corrected by Bonferroni correction for multiple testing (the threshold of raw P-value significance was 0.0018 [14 polymorphisms and haplotypes, and 2 phenotypes analyzed]). A Type I error of 5% was used to obtain the Bonferroni corrected p-value. For the haplotype analyses, we fit the model with the same covariates in a similar manner.
Availability and requirements
R statistical package: http://www.r-project.org
Sorimachi H, Freiburg A, Kolmerer B, Ishiura S, Stier G, Gregorio CC, Labeit D, Linke WA, Suzuki K, Labeit S: Tissue-specific expression and alpha-actinin binding properties of the Z-disc titin: implications for the nature of vertebrate Z-discs. Journal of molecular biology. 1997, 270 (5): 688-695. 10.1006/jmbi.1997.1145.
Suzuki K, Sorimachi H: A novel aspect of calpain activation. FEBS letters. 1998, 433 (1-2): 1-4. 10.1016/S0014-5793(98)00856-4.
Beltran JA, Jaime I, Santolaria P, Sanudo C, Alberti P, Roncales: Effect of stress-induced high post-mortem pH on protease activity and tenderness of beef. Meat Sci. 1997, 45: 201-207. 10.1016/S0309-1740(96)00103-9.
Koohmaraie M: The role of Ca(2+)-dependent proteases (calpains) in post mortem proteolysis and meat tenderness. Biochimie. 1992, 74 (3): 239-245. 10.1016/0300-9084(92)90122-U.
Koohmaraie M: Muscle proteinases and meat aging. Meat Sci. 1994, 36: 93-104. 10.1016/0309-1740(94)90036-1.
Koohmaraie M: Biochemical factors regulating the toughening and tenderization process of meat. Meat Sci. 1996, 43: S193-S201. 10.1016/0309-1740(96)00065-4.
Geesink GH, Koohmaraie M: Effect of calpastatin on degradation of myofibrillar proteins by mu-calpain under postmortem conditions. Journal of animal science. 1999, 77 (10): 2685-2692.
Smith TP, Casas E, Rexroad CE, Kappes SM, Keele JW: Bovine CAPN1 maps to a region of BTA29 containing a quantitative trait locus for meat tenderness. Journal of animal science. 2000, 78 (10): 2589-2594.
Casas E, Shackelford SD, Keele JW, Koohmaraie M, Smith TP, Stone RT: Detection of quantitative trait loci for growth and carcass composition in cattle. Journal of animal science. 2003, 81 (12): 2976-2983.
Page BT, Casas E, Heaton MP, Cullen NG, Hyndman DL, Morris CA, Crawford AM, Wheeler TL, Koohmaraie M, Keele JW, Smith TP: Evaluation of single-nucleotide polymorphisms in CAPN1 for association with meat tenderness in cattle. Journal of animal science. 2002, 80 (12): 3077-3085.
Juszczuk-Kubiak E, Sakowski T, Flisikowski K, Wicinska K, Oprzadek J, Rosochacki SJ: Bovine mu-calpain (CAPN1) gene: new SNP within intron 14. Journal of applied genetics. 2004, 45 (4): 457-460.
Page BT, Casas E, Quaas RL, Thallman RM, Wheeler TL, Shackelford SD, Koohmaraie M, White SN, Bennett GL, Keele JW, Dikeman ME, Smith TP: Association of markers in the bovine CAPN1 gene with meat tenderness in large crossbred populations that sample influential industry sires. Journal of animal science. 2004, 82 (12): 3474-3481.
White SN, Casas E, Wheeler TL, Shackelford SD, Koohmaraie M, Riley DG, Chase CC, Johnson DD, Keele JW, Smith TP: A new single nucleotide polymorphism in CAPN1 extends the current tenderness marker test to include cattle of Bos indicus, Bos taurus, and crossbred descent. Journal of animal science. 2005, 83 (9): 2001-2008.
Morris CA, Cullen NG, Hickey SM, Dobbie PM, Veenvliet BA, Manley TR, Pitchford WS, Kruk ZA, Bottema CD, Wilson T: Genotypic effects of calpain 1 and calpastatin on the tenderness of cooked M. longissimus dorsi steaks from Jersey x Limousin, Angus and Hereford-cross cattle. Animal genetics. 2006, 37 (4): 411-414. 10.1111/j.1365-2052.2006.01483.x.
Rincon G, Medrano JF: Assays for genotyping single nucleotide polymorphisms in the bovine CAPN1 gene. Animal genetics. 2006, 37 (3): 294-295. 10.1111/j.1365-2052.2006.01430.x.
Van Eenennaam AL, Li J, Thallman RM, Quaas RL, Dikeman ME, Gill CA, Franke DE, Thomas MG: Validation of commercial DNA tests for quantitative beef quality traits. Journal of animal science. 2007, 85 (4): 891-900. 10.2527/jas.2006-512.
APGS: Report of business for animal products grading. Animal products grading system, National Livestock Co-operatives Federation, Korea. 1995
Vreeland WN, Meagher RJ, Barron AE: Multiplexed, high-throughput genotyping by single-base extension and end-labeled free-solution electrophoresis. Anal Chem. 2002, 74 (17): 4328-4333. 10.1021/ac0258094.
Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics (Oxford, England). 2005, 21 (2): 263-265. 10.1093/bioinformatics/bth457.
Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, Higgins J, DeFelice M, Lochner A, Faggart M, Liu-Cordero SN, Rotimi C, Adeyemo A, Cooper R, Ward R, Lander ES, Daly MJ, Altshuler D: The structure of haplotype blocks in the human genome. Science. 2002, 296 (5576): 2225-2229. 10.1126/science.1069424.
Stephens M, Smith NJ, Donnelly P: A new statistical method for haplotype reconstruction from population data. Am J Hum Genet. 2001, 68 (4): 978-989. 10.1086/319501.
This study was supported by the Technology Development Program for Agriculture and Forestry, Ministry of Agriculture and Forestry, Republic of Korea.
HSC carried out the statistical analyses and drafted the manuscript. LHK, SN, HWL, CSH, and JOK carried out the sequencing, SNP discovery, and genotyping. D-HY and I-CC participated in the animal and phenotypic data collection. HDS participated in the design and coordination, and helped to draft the manuscript. BLP and JSB carried out further statistical analyses. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Supplementary information to CAPN1. The data provided represent primer information for CAPN1 sequencing, genotyping probe and chromatograms of each SNP. (DOC 2 MB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Cheong, H.S., Yoon, DH., Park, B.L. et al. A single nucleotide polymorphism in CAPN1 associated with marbling score in Korean cattle. BMC Genet 9, 33 (2008). https://doi.org/10.1186/1471-2156-9-33