Identification and association of the single nucleotide polymorphisms in calpain3 (CAPN3) gene with carcass traits in chickens
© Zhang et al; licensee BioMed Central Ltd. 2009
Received: 04 January 2009
Accepted: 05 March 2009
Published: 05 March 2009
The aim of this study is to screen single nucleotide polymorphisms (SNP) of chicken Calpain3 (CAPN3) gene and to analyze the potential association between CAPN3 gene polymorphisms and carcass traits in chickens. We screened CAPN3 single nucleotide polymorphisms (SNP) in 307 meat-type quality chicken from 5 commercial pure lines (S01, S02, S03, S05, and D99) and 4 native breeds from Guangdong Province (Huiyang Huxu chicken and Qingyuan Ma chicken) and Sichuan Province (Caoke chicken and Shandi Black-bone chicken), China.
Two SNPs (11818T>A and 12814T>G) were detected by single strand conformation polymorphism (SSCP) method and were verified by DNA sequencing. Association analysis showed that the 12814T>G genotypes were significantly associated with body weight (BW), carcass weight (CW), breast muscle weight (BMW), and leg muscle weight (LMW). Haplotypes constructed on the two SNPs (H1, TG; H2, TT; H3, AG; and H4, AT) were associated with BW, CW (P < 0.05), eviscerated percentage (EP), semi-eviscerated percentage (SEP), breast muscle percentage (BMP), and leg muscle percentage (LMP) (P < 0.01). Diplotype H1H2 was dominant for BW, CW, and LMP, and H2H2 was dominant for EP, SEP, and BMP.
We speculated that the CAPN3 gene was a major gene affecting chicken muscle growth and carcass traits or it was linked with the major gene(s). Diplotypes H1H2 and H2H2 might be advantageous for carcass traits.
The accumulation of muscle mass in the body is determined by the relative rates of muscle protein synthesis and degradation, namely, this process is closely related to the muscle protein turnover. Studies on the rate of muscle growth in layer and broiler chickens have shown that the rate of muscle growth are mainly determined by changes in the rate of muscle protein degradation, with little or no change in the rate of muscle protein synthesis .
Calpains, i.e. intracellular Ca2+-dependent cysteine proteases (EC 220.127.116.11), have long been demonstrated to be involved in muscle growth and development and are proenzymes that are regulated by Ca2+ binding and autoproteolytic modification . The calpains were discovered because the Z-disks in muscle strips incubated in a Ca2+-containing solution disappeared but with no other ultrastructurally detectable change . Four calpain genes (CAPN1, CAPN2, CAPN3, and CAPN1.5) are expressed ubiquitously in chicken . Study have further showed that there was less calpain activity, but more calpastatin, in broiler breast muscle . Among the calpain family members, calpain 3 (Capn3, previously named p94) is particularly interesting and has two insertion sequences, IS1 and IS2, which are involved in the regulation of its function and activity . CAPN3 was identified by positional cloning as the gene responsible for limb-girdle muscular dystrophy 2A (LGMD2A) [7, 8], where it plays an important role in myofibrillar integrity . Capn3 binds specifically to connectin (a striated muscle specific giant filamentous protein) at the N2 line , a region where proteolysis has been linked to post-mortem tenderization . CAPN3 has received considerable attention in recent years because of its potential influence on muscle growth and its role in myofibrillar organization. However, much of the available information about this gene was taken from studies on humans, mice, and other mammalian species, and might not be directly applicable to poultry.
Calpain 3 is ubiquitously expressed in chicken  and is mainly located in skeletal muscle. Compared to the other calpains, Capn3 activity is regulated by different mechanisms, as its activation and function as a muscle cytoskeleton regulator requires very low concentration of calcium and different interacting proteins . We hypothesized that CAPN3 gene may be a major gene affecting chicken carcass traits.
In this study, we screened CAPN3 single nucleotide polymorphisms (SNP) in 307 meat-type quality chickens from 5 commercial pure lines (S01, S02, S03, S05, and D99) and 4 native breeds [Huiyang Huxu chicken (HH), Qingyuan Ma chicken (QY), Caoke chicken (CK) and Mountainous black-bone chicken (MB)] in China, with an intention to evaluate potential association between CAPN3 polymorphisms and the chicken carcass traits.
SNPs of the Chicken CAPN3 Gene
According to UCSC Genome Browser on Chicken May 2006 Assembly http://genome.ucsc.edu/cgi-bin/hgBlat?command=start using Gallus gallus Calpain3 sequence (GenBank accession number XM_421157), chicken CAPN3 gene contains 24 exons. We chose exons 3, 4, 5, 8, 9, 10, 16, 17, and 19 for analyses based on a preliminary analysis of the BGI SNPs, which contain reads from the Broiler, Layer and Silkie chicken strains and are available at the UCSC Genome Browser. A PCR-based SSCP method was successfully developed for screening the nucleotide substitutions in these exons. Two SNPs, a T/A variant at position nt11818 (located in intron 8) and a T/G variant at nt12814 (located in exon 10), were identified in the chicken DNA sequence (GenBank Accession No: NC_006092.2, GI: 118136304). Variant 12814T>G was silent and did not cause the amino acid change.
Genotype and Allele Frequencies in Different Chicken Populations
Additional file 1 showed the allele and genotype frequencies of the two SNPs in nine chicken populations. For SNP 11818T>A, the allele A was the dominant allele and presented the highest allele frequencies among all chicken populations (average = 58.67%). For SNP 12814T>G, the allele frequency of T was higher than that of G in populations S03 and S02, while in S01, D99, S05, QY, HH, and SD, the allele G was identified as the dominant one. The χ2-test for the genotype frequency showed that there were statistically significant differences among the nine chicken populations for 12814T>G (P < 0.01), whereas the difference for 11818T>A was not significant (P > 0.05). Furthermore, variant 11818T>A was significantly deviated from the Hardy-Weinberg equilibrium (HWE) in samples D99, S05, SD, and CK, whereas variant 12814T>G was not violated HWE.
Association of CAPN3 SNPs with Carcass Traits
Least square mean carcass traits, by genotype, of chicken CAPN3 for each locus
1615.80a ± 26.85
1568.22b ± 28.93
1578.44b ± 35.10
72.64a ± 0.27
72.09b ± 0.29
72.82a ± 0.35
12.68a ± 0.13
12.64a ± 0.14
12.56b ± 0.17
18.18b ± 0.15
18.14b ± 0.17
18.31a ± 0.20
1534.53b ± 38.69
1646.77a ± 25.27
1538.59b ± 37.21
1356.31b ± 43.65
1467.70a ± 28.51
1358.04b ± 41.98
86.74a ± 0.41
87.04a ± 0.26
85.53b ± 0.39
153.03b ± 4.74
163.31a ± 3.09
150.03b ± 4.56
220.97b ± 7.39
241.68a ± 4.83
220.35b ± 7.11
For 12814T>G, BW, carcass weight (CW), BMW and LMW of chickens with the TG genotype were significantly higher than those with the GG and TT genotypes, respectively (P < 0.05). The semi-eviscerated percentage (SEP) values of chickens with the GG genotype were significantly lower than those with the TG and TT genotypes, respectively (P < 0.01). No significant difference was detected for other carcass traits.
Construction of Haplotypes and Their Frequencies
Haplotypes inferred based on the two single nucleotide polymorphisms of CAPN3 gene
Diplotypes of chicken CAPN3 gene
Association of the Haplotypes with Carcass Traits
Effects of haplotypes on carcass traits
Associations between the diplotypes and chicken carcass traits
BW (g) *
CW (g) *
SEP (%) **
EP (%) **
BMP (%) **
LMP (%) **
1525.81 ± 67.03
1363.02 ± 74.01
86.10 ± 0.69
72.09 ± 0.68
12.15 ± 0.34
17.72 ± 0.50
1655.48 ± 48.15
1540.09 ± 53.16
86.89 ± 0.50
73.05 ± 0.49
12.56 ± 0.24
18.59 ± 0.29
1432.47 ± 84.64
1254.29 ± 93.46
85.76 ± 0.51
71.63 ± 0.50
12.67 ± 0.25
17.86 ± 0.29
1523.36 ± 49.22
1339.97 ± 54.35
87.65 ± 0.87
73.10 ± 0.86
13.20 ± 0.43
18.17 ± 0.40
1615.98 ± 45.88
1418.61 ± 50.66
87.03 ± 0.47
72.52 ± 0.47
12.86 ± 0.23
18.45 ± 0.27
1567.74 ± 59.06
1390.66 ± 65.22
86.08 ± 0.61
71.91 ± 0.60
12.27 ± 0.30
18.07 ± 0.35
1600.89 ± 48.77
1423.76 ± 53.85
86.49 ± 0.51
72.62 ± 0.50
12.53 ± 0.24
18.12 ± 0.29
1651.25 ± 37.35
1450.81 ± 41.24
87.09 ± 0.38
72.49 ± 0.38
12.57 ± 0.19
18.26 ± 0.22
1535.84 ± 59.08
1349.38 ± 65.24
87.02 ± 0.61
73.08 ± 0.60
13.10 ± 0.30
18.06 ± 0.35
Understanding the genetic basis of protein metabolism in chickens will provide an opportunity for genetic improvement of the muscle growth traits. The study of candidate genes, based on the related results that were well characterized in humans and mice, is one of the primary methods to determine whether certain specific genes are associated to economic traits in farm animals.
Previous studies showed that calpain gene expression was correlated with muscle organization [15–17], functionally related to titin muscle . In poultry, Maeda et al.  measured muscle protein turnover rate and calpain activity in the muscle of two quail lines divergently selected for body size. Zhang et al.  found that three variants in CAPN1 gene are associated with some carcass traits in chickens. Muscle protein turnover rate and calpain activity were significantly higher in line with small body size compared with the large line. In the current study, we analyzed CAPN3 gene polymorphisms and found 2 SNPs (11818T>A and 12814T>G), which showing associations (P < 0.05) with partial carcass traits. The result suggests that CAPN3 may act as a candidate gene of quantitative trait loci (QTL) in regulating muscle growth.
Haplotype or haplotype block analysis provided a practical solution to resolve the innate problems of the single-marker analysis, such as noisy, unsatisfied, and obscured important localization information . Both haplotype diversity and the method of SNP selection based on maximum haplotype diversity were always preferred . In this study, haplotypes were constructed using the two CAPN3 SNPs and were employed to evaluate the potential association between haplotypes and carcass traits. The diplotype H2H2 was found to be associated with higher SEP (%), EP (%), and BMP (%) than other diplotypes, respectively. Chickens with diplotype H1H2 had higher BW, CW, and LMP (%) than those chickens harboring other diplotypes. Apparently, haplotype H2 may be the most advantageous haplotype for carcass traits. Note that our haplotype blocks were based on two SNPs and more tag SNPs should be typed in order to detect the complete haplotype blocks in different chicken breeds.
It should be mentioned that the four native chicken breeds (HH, QY, SD, and CK) and the five commercial pure lines have been subjected to intensive artificial selection on commercial traits, such as carcass weight, eviscerated weight, and breast muscle weight, which might account for the significant deviation from the Hardy – Weinberg equilibrium for variant 11818T>A in some populations. Moreover, it is hard to include a term in the statistical analysis to adjust for relatedness among the individual chickens from the population. Our model also did not include the effects of hatch, location (e.g. cage), and rearing phase (e.g. randomization of breeds to location) that have major effects on growth traits. Therefore, different chicken populations with various domestication background and selection history, further optimized analysis model, more tag SNPs are needed to verify the genetic effects of the CAPN3 gene.
In summary, two SNPs in intron 8 and exon 10 of CAPN3 gene were found in a large cohort of chicken samples from China and there was an association between CAPN3 haplotypes and carcass traits. It is suggestive that CAPN3 may be a potential major gene or in close linkage disequilibrium with a QTL for muscle growth in chicken. Further analyses of the effects of CAPN3 polymorphisms are essential to confirm the association between the alleles or haplotypes and chicken carcass traits.
In total, 307 meat-type quality chicken (with a male/female ratio as 1:1) were used. These samples included 5 commercial pure lines and 4 Chinese domestic chicken breeds. All populations were raised under the same condition and were randomly selected. All the pure lines were officially recognized and were originated from the local breeds in Sichuan (CK and SD) and Guangdong (HH and QY) Provinces. The CK and SD chickens have spotty feathers, black or yellow skin, and a favorable meat quality. The HH and QY chickens are named by their yellow plumage, skin, and shank and have a high quality of meat. Different breeds were hatched in different incubators under the same condition on June 25, 2007. The birds were kept in single breed pens and slaughtered on October 13, 2007. All birds had free access to food and water. Commercial corn-soybean diets that met all NRC requirements were provided in the study. From birth to 3 week (wk) of age, birds received a starter feed (2.90 megacalories (Mcal) of metabolizable energy (ME)/kg and 20.5 g/kg of crude protein (CP)). Birds were fed a grower diet (3.00 Mcal of ME/kg and 18.5 g/kg of CP) from 4 to 6 wk of age and were transferred to the growing pens at the age of 7 wk. Before slaughter, blood was collected and the genomic DNA was isolated by the standard phenol/chloroform method. For each population, 27 to 62 chickens were randomly sampled for collecting blood and slaughtering.
Birds were raised in cage according to the conventional program of commercial broiler. Traits were measured at 90 days age, including BW of the live birds after 12 h with no access to feed, CW, semi-eviscerated weight (SEW), eviscerated weight (EW), BMW, LMW, abdominal fat weight (AFW), subcutaneous fat thickness (SFT). The CW was measured on the chilled carcass removed feathers. Semi-eviscerated weight was measured on the carcass removed trachea, esophagus, gastrointestinal tract, spleen, pancreas, and gonad. Eviscerated weight was measured on the semi-eviscerated weight after removal of head, claws, heart, liver, gizzard, glandular stomach, and abdominal fat. Subcutaneous fat thickness was measured at the caudal spondyle including the skin and fat width with a vernier caliper after dressing. The ratios of each of the above traits to CW were calculated as eviscerated percentage (EP), semi-eviscerated percentage (SEP), breast muscle percentage (BMP), and leg muscle percentage (LMP), respectively.
PCR-SSCP Screening for CAPN3 Mutations
Nine primer pairs were designed to amplify CAPN3 exons according to the sequence of Gallus gallus CAPN3 (GenBank Accession No: NC_006092.2, GI: 118136304) (Additional file 2). Primers synthesis was completed by the TaKaRa Biotechnology Co., Ltd (Dalian, China). The PCR reaction was performed with the following condition: one denaturation cycle at 94°C for 6 min, followed by 35 cycles at 94°C for 45 s, 52°C for 45 s, 72°C for 1 min, and ended with an extension cycle at 72°C for 8 min. The 10 μL reaction volume included 0.8 μL DNA template (50 ng/μL), 5 μL 2 × Taq PCR MasterMix (Beijing TIAN WEI Biology Technique Corporation, Beijing, China), 3.6 μL ddH2O, 0.3 μL of each primer (10 pmol/μL). PCR products were separated on 1% agarose gel and were visualized on gel imaging system (Gel DocTMEQ170-8060) and photographed. 3 μL PCR products was denatured at 98°C for 10 min, then were placed on ice for 5 min. Electrophoresis was run at 10 V/cm on a 12% polyacrylamide gel electrophoresis for 9 to 10 h at 20°C. Gels were stained using the silver staining. Individual single strand conformation polymorphism banding pattern was determined under visible light. Samples showing different bands in the gel were further amplified and purified, and were sequenced by a commercial sequencing company (Shanghai Yingjun Biology Technique Corporation, Shanghai, China).
Data were analyzed using the GLM (General Linear Models Procedures) of SAS 8.0 (SAS Institute Inc., Cary, NC). The genetic effects were analyzed by a general mixed procedure in the SAS package, and the following model was used:
Y = μ + B + S + G + bX + e
where Y = the dependant variable; μ = the population mean; B = fixed effect of the breed; S = fixed effect of the sex; G = genotype value; e = random effect; b = the regression coefficients of the linear regressions on BW (X); BW was included in the model as it was significant for traits (SEW, EW, BMW, LMW, and AFW) at single locus. Multiple comparisons were analyzed with the least square means (LSM), and it was excluded from the model if its effect had a P value larger than 0.05 on a given trait. The significance of the least square means was tested with the Duncan's Multiple Range test (P < 0.01). The data for all carcass traits were verified for normal distribution by using the Shapiro-Wilks test in SAS 8.0. None of them were normally distributed. Thus, the parameters of BW, CW, SEW, EW, BMW, LMW, AFW, and SFT were estimated on the square root scale before they were analyzed.
Haplotypes were constructed based on the 2 SNPs identified in all 307 experimental animals by using the PHASE 2.0 program . The program reconstructed haplotypes from a population data. The genetic status of the subjects was expressed as the combination of two haplotypes (diplotype configuration). Genetic effects of the haplotypes were performed with the mixed model mentioned above.
The authors thank Xiao-Chen Li and Lian Xiao for technical support. This work was financially supported by Chinese national key technologies R & D programme (Grant No: 2008BADB2B01), the National 863 Project of China (No. 2008AA101001) and the provincial Eleventh five-year program on animal breeding of Sichuan Province (2006YZGG-19; 2006YZGG-25). YGY was supported by the "Century Program" (or Hundreds-Talent Program) of the Chinese Academy of Sciences.
- Goll DE, Kleese WC, Szpacenko A: Skeletal muscle proteases and protein turnover from animal growth regulation. 1989, Plenum Publishing Corporation, 141-176.
- Goll DE, Thompson VF, Li H, Wei W, Cong J: The calpain system. Physiol Rev. 2003, 83 (3): 731-801.View ArticlePubMed
- Dayton WR, Schollmeyer JV, Lepley RA, Cortes LR: A calcium-activated protease possibly involved in myofibrillar protein turnover, isolation of a lowcalcium requiring form of the protease. Biochim Biophys Acta. 1981, 659: 48-61.View ArticlePubMed
- Okumura F, Shimogiri T, Shinbo Y, Yoshizawa K, Kawabe K, Mannen H, Okamoto S, Cheng HH, Maeda Y: Linkage mapping of four chicken calpain genes. J Anim Sci. 2005, 76 (2): 121-127. 10.1111/j.1740-0929.2005.00246.x.View Article
- Schreurs F, Heide Van Der JG, Leenstra FR, Dewit W: Endogenous proteolytic enzymes in chicken muscles. Differences among strains with growth rates and protein efficiencies. Poultry Sci. 1995, 74: 523-537.View Article
- Herasse M, Ono Y, Fougerousse F, Kimura E, Stockholm D, Beley C, Montarras D, Pinset C, Sorimachi H, Suzuki K: Expression and functional characteristics of calpain 3 isoforms generated through tissue-specific transcriptional and posttranscriptional events. Mol Cell Biol. 1999, 19 (6): 4047-4055.PubMed CentralView ArticlePubMed
- Richard I, Broux O, Allamand V, Fougerousse F, Chiannilkulchai N, Bourg N, Brenguier L, Devaud C, Pasturaud P, Roudaut C: Mutations in the proteolytic enzyme calpain 3 cause limb-girdle muscular dystrophy type 2A. Cell. 1995, 81 (1): 27-40. 10.1016/0092-8674(95)90368-2.View ArticlePubMed
- Richard I, Roudaut C, Saenz A, Pogue R, Grimbergen JE, Anderson LV, Beley C, Cobo AM, de Diego C, Eymard B: Calpainopathy-a survey of mutations and polymorphisms. Am J Hum Genet. 1999, 64 (6): 1524-1540. 10.1086/302426.PubMed CentralView ArticlePubMed
- Poussard S, Duvert M, Balcerzak D, Ramassamy S, Brustis JJ, Cottin P, Ducastaing A: Evidence for implication of muscle-specific calpain (p94) in myofibrillar integrity. Cell Growth Differ. 1996, 7 (11): 1461-1469.PubMed
- Sorimachi H, Imajoh-Ohmi S, Emori Y, Kawasaki H, Ohno S, Minami Y, Suzuki K: Molecular cloning of a novel mammalian calcium-dependent protease distinct from both m- and μ-types. J Biol Chem. 1989, 264 (33): 20106-20111.PubMed
- Taylor RG, Geesink GH, Thompson VF, Koohmaraie M, Goll DE: Is Z-disk degradation responsible for postmortem tenderization?. J Anim Sci. 1995, 73 (5): 1351-1367.PubMed
- Taveau M, Bourg M, Sillon G, Roudaut C, Bartoli M, Richard I: Calpain 3 Is Activated through Autolysis within the Active Site and Lyses Sarcomeric and Sarcolemmal Components. Mol Cell Biol. 2003, 23 (24): 9127-9135. 10.1128/MCB.23.24.9127-9135.2003.PubMed CentralView ArticlePubMed
- Wang J, He XM, Ruan J, Dai MT, Chen J, Zhang Y, Hu YF, Ye C, Li ST, Cong LJ: ChickVD: a sequence variation database for the chicken genome. Nucleic Acids Res. 2005, 33: D438-D441. 10.1093/nar/gki092.PubMed CentralView ArticlePubMed
- 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.PubMed CentralView ArticlePubMed
- Ono Y, Kakinuma K, Torii F, Irie A, Nakagawa K, Labeit S, Abe K, Suzuki K, Sorimachi H: Possible regulation of the conventional calpain system by skeletal muscle-specific calpain, p94/calpain 3. J Biol Chem. 2004, 279 (4): 2761-2771. 10.1074/jbc.M308789200.View ArticlePubMed
- Stockholm D, Herasse M, Marchand S, Praud C, Roudaut C, Richard I, Sebille A, Beckmann JS: Calpain 3 mRNA expression in mice after denervation and during muscle regeneration. Am J Physiol Cell Physiol. 2001, 280 (6): C1561-1569.PubMed
- Walder K, McMillan J, Lapsys N, Kriketos A, Trevaskis J, Civitarese A, Southon A, Zimmet P, Collier G: Calpain 3 gene expression in skeletal muscle is associated with body fat content and measures of insulin resistance. Int J Obes Relat Metab Disord. 2002, 26 (4): 442-449. 10.1038/sj.ijo.0801946.View ArticlePubMed
- Garvey SM, Rajan C, Lerner AP, Frankel WN, Cox GA: The muscular dystrophy with myositis (mdm) mouse mutation disrupts a skeletal muscle-specific domain of titin. Genomics. 2002, 79 (2): 146-149. 10.1006/geno.2002.6685.View ArticlePubMed
- Maeda Y, Kawabe K, Okamoto S, Hashiguchi T: Genetical studies on muscle protein turnover rate and calcium activated neutral protease activity in the skeletal muscle of the Japanese quail, Coturnix coturnix japonica. Anim Sci and Tech. 1991, 62 (9): 813-821.
- Zhang ZR, Liu YP, Jiang XS, Du HR, Zhu Q: Study on association of single nucleotide polymorphism of CAPN1 gene with muscle fibre and carcass traits in quality chicken populations. J Anim Breed Genet. 2008, 125: 258-264. 10.1111/j.1439-0388.2008.00723.x.View ArticlePubMed
- Daly MJ, Rioux JD, Schaffner SF, Hudson TJ, Lander ES: High-resolution haplotype structure in the human genome. Nat Genet. 2001, 29 (2): 229-232. 10.1038/ng1001-229.View ArticlePubMed
- Huang Q, Fu YX, Boerwinkle E: Comparison of strategies for selecting single nucleotide polymorphisms for case/control association studies. Hum Genet. 2003, 113 (3): 253-257. 10.1007/s00439-003-0965-x.View ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.