Association of a single nucleotide polymorphism in titin gene with marbling in Japanese Black beef cattle
- Takahisa Yamada†1,
- Seiki Sasaki†2,
- Shin Sukegawa†3,
- Sachiyo Yoshioka4,
- Youichi Takahagi3,
- Mitsuo Morita2,
- Hiroshi Murakami3,
- Fumiki Morimatsu3,
- Tatsuo Fujita5,
- Takeshi Miyake1 and
- Yoshiyuki Sasaki4Email author
© Sasaki et al; licensee BioMed Central Ltd. 2009
Received: 27 January 2009
Accepted: 07 May 2009
Published: 07 May 2009
Marbling defined by the amount and distribution of intramuscular fat is an economically important trait of beef cattle in Japan. We have recently reported that single nucleotide polymorphisms (SNPs) in the endothelial differentiation, sphingolipid G-protein-coupled receptor, 1 (EDG1) gene were associated with marbling in Japanese Black beef cattle. As well as EDG1, the titin (TTN) gene, involved in myofibrillogenesis, has been previously shown to possess expression difference in musculus longissimus muscle between low-marbled and high-marbled steer groups, and to be located within genomic region of a quantitative trait locus for marbling. Thus TTN was considered as a positional functional candidate for the gene responsible for marbling. In this study, we explored SNP in TTN and analyzed association of the SNP with marbling.
A SNP in the promoter region of TTN, referred to as g.231054C>T, was the only difference detected between high- and low-marbled steer groups. The SNP was associated with marbling in 3 experiments using 101 sires (P = 0.004), 848 paternal half-sib progeny steers from 5 sires heterozygous for the g.231054C>T (P = 0.046), and 820 paternal half-sib progeny steers from 3 sires homozygous for C allele at the g.231054C>T (P = 0.051), in Japanese Black beef cattle. The effect of genotypes of the SNP on subcutaneous fat thickness was not statistically significant (P > 0.05).
These findings suggest that in addition to the EDG1 SNPs, the TTN SNP polymorphism is associated with marbling and may be useful for effective marker-assisted selection to increase the levels of marbling in Japanese Black beef cattle. Further replicate studies will be needed to confirm the allelic association observed here, and to expand the results to evaluate all possible genotypic combinations of alleles.
Marbling has an important effect on the economics of beef production [1–4]. Thus, it is greatly interesting to obtain better knowledge on the molecular architecture of marbling and to generate new opportunities for more effective marker-assisted selection.
The endothelial differentiation, sphingolipid G-protein-coupled receptor, 1 (EDG1) gene, which is involved in blood vessel formation , exhibited higher expression levels in high-marbled steer group than in low-marbled steer group in musculus longissimus muscle across all ages of the test period [6, 7]. We have recently showed that single nucleotide polymorphisms (SNPs) in the EDG1 were associated with marbling in Japanese Black beef cattle .
From the result of our previous differential-display PCR (ddPCR) analysis in musculus longissimus muscle , the titin (TTN) gene exhibited higher expression levels in low-marbled steer group than in high-marbled steer group in the late stage of the test period. Since TTN is known to be involved in myofibrillogenesis , the differential expression suggested that the decrease in the TTN expression might promote proliferation, differentiation, or maturation of adipocyte-lineage cells by weakening the structural integrity of the sarcomere, thereby resulting in high levels of marbling.
We have also located the TTN gene within genomic region of a quantitative trait locus for marbling [10, 11], which is mapped in a half-sib family of Japanese Black cattle to bovine chromosome 2 region . Thus, the TTN gene was regarded as a positional functional candidate for the gene responsible for marbling.
TTN gene structure
We screened the NCBI bovine genome sequence database (National Center for Biotechnology Information, Bethesda, MD) with the publicly available bovine TTN cDNA sequence [GenBank:XR_027685] including 101,593 bp, and obtained 5,046,004-bp bovine genomic sequence [GenBank:NW_001494580] containing the TTN gene. By comparing the genomic and cDNA sequences, we showed that the bovine TTN gene spans 276,648 bp and comprises 211 bp of the first exon, 1,315 bp of the last exon, 100,067 bp of 313 internal exons, and 175,055 bp of 314 introns.
Two Holstein steers and two somatic nuclear-derived cloned steers  from a Japanese Black Itofuku sire with a very high estimated breeding value for marbling , which were assigned for low-marbled and high-marbled steer groups, respectively, in our previous ddPCR analysis , were used for SNP detection in this study. The two low-marbled Holstein steers and two high-marbled cloned steers were previously shown to have different TTN gene expression patterns . The details of DNA samples from these steers were described previously . We used two high-marbled cloned steers to confirm the authenticity of newly discovered SNP in the TTN gene.
Primers designed for SNP detection in the bovine TTN gene
-7,203 to -7,183
-5,262 to -5,283
-5,428 to -5,406
-3,487 to -3,506
-3,618 to -3,599
-1,674 to -1,693
-1,806 to -1,782
329 to 308
275,190 to 275,209
276,756 to 276,736
We sequenced the 6,615-bp 5' flanking region, the 211-bp first exon region, and the 1,315-bp last exon region of the TTN gene for the 4 DNA samples. This sequence analysis revealed only one SNP in the TTN gene between the 2 steer groups: a C to T transition located 652 bp upstream of the putative transcription initiation site in the promoter region (g.231054C>T). The g.231054C>T has been submitted to the dbSNP with ss accession number 120037472. The two Holstein steers were homozygous for C allele at the g.231054C>T site, whereas the two cloned steers heterozygous for C and T alleles at the g.231054C>T site.
Samples and data
We performed 3 experiments for association of the g.231054C>T with marbling and subcutaneous fat thickness. We used 101 Japanese Black sires in experiment 1. The sires used to be or are used in the Oita Prefectural Institute of Animal Industry (Oita, Japan). There was no strong bias for a specific father or a specific maternal grandfather of the sires, and the sire panel likely represents a variety of the sire lines. In experiment 2, 848 paternal half-sib Japanese Black progeny steers (55 to 445 steers per sire) produced from 5 sires heterozygous for the g.231054C>T, with dams considered to represent a random sample of the female population, were used. We used 820 paternal half-sib Japanese Black progeny steers (88 to 535 steers per sire) produced from 3 sires homozygous for C allele at the g.231054C>T, with dams considered to be a random mating population, in experiment 3. These progeny steers (in experiments 2 and 3) were fattened, and shipped to the carcass market in the Oita prefecture. Semen or blood of the sires and adipose tissues of the progeny steers were collected for SNP genotyping. DNA samples were prepared from the materials according to standard protocols.
The predicted breeding values of the sires and the progeny steers for marbling score and subcutaneous fat thickness were obtained from the Oita recording system for beef cattle reported by Sasaki et al. , and used as phenotypic values in this study. In the recording system, the breeding values were predicted from carcass records of Japanese Black steers and heifers, fattened in the Oita prefecture. The fattened animals were shipped to various carcass markets from 1988 to 2003, where they were slaughtered and their carcasses evaluated. The data were edited to connect across subclasses such that each market-year subclass had 50 or more animals and each farm had 10 or more animals. The final number of animals was 48,045, and there were 89 market-year subclasses, 332 farms, and 228 sires.
Marbling score and subcutaneous fat thickness were measured on carcasses dissected at the sixth and seventh rib section, according to the Japanese meat grading system by certified graders from the Japan Meat Grading Association (Tokyo, Japan) . Marbling score was scored from 1 to 12 (beef marbling standard number) with a standard model panel, in which higher scores correspond to more marbling. Score 1 likely corresponds to "Practically Devoid", "Traces", "Slight" and "Small", score 2 or 3 "Modest", "Moderate" and "Slightly Abundant", and score 4 or 5 "Moderately Abundant". Subcutaneous fat thickness was measured at the rib interface perpendicular to the outside surface at a point 3-fourths the length of the longissimus muscle from its chine bone end.
Data were analyzed by the REML method using the MTDFREML programs , and genetic and environmental variances were estimated. The BLUP option in the programs using the estimated variance components was chosen to predict the breeding values of animals with a single trait model. Sex, market-year, and farm were considered fixed effects. Fattening period and slaughter age were also considered as up to quadratic covariates. The fattening period denotes the period from the start of fattening to shipping to market for each animal. These fixed effects were all significant (P < 0.001). Random effects included additive genetic effect of the individuals; that is, the animal model was adopted to predict the breeding values.
This study conformed to the guidelines for animal experimentation of the Graduate School of Agriculture, Kyoto University (Kyoto, Japan).
The g.231054C>T was genotyped using PCR-restriction fragment length polymorphism (RFLP) method. PCR primers used for PCR-RFLP were 5'-TCATCTCCTAACTACTTCCCA-3' and 5'-ACAAAATCTGAACCTGGCTT-3' (Nucleotide positions relative to the putative transcription initiation site of the TTN gene were -842 to -822 and -616 to -635, respectively). PCR amplifications were carried out as described in SNP detection section, using a final volume of 20 μl and the annealing temperature of 55°C. An aliquot of PCR-amplified products was digested at 37°C for 1 h with restriction enzyme Hpy CH4III, and electrophoresed on a 3.0% agarose gel. Agarose gels were stained with ethidium bromide and photographed under an ultraviolet light. Using this method, 227-bp PCR fragments containing the SNP site were amplified, and Hpy CH4III-digested into 36- and 191-bp fragments at the C allele, but not at the T allele: the CC homozygotes, the TT homozygotes, and the CT heterozygotes resulted in 2 bands (36 and 191 bp), 1 band (227 bp), and 3 bands (36, 191, and 227 bp).
Effect of the SNP genotypes on marbling and subcutaneous fat thickness in experiment 1
No. of animals
Subcutaneous fat thickness (mm)
2.81a ± 0.25
-2.85 ± 0.74
2.13b ± 0.23
-2.84 ± 0.80
Effect of the SNP genotypes on marbling and subcutaneous fat thickness in experiment 2
No. of animals
Subcutaneous fat thickness (mm)
3.42a ± 0.08
-0.33 ± 0.36
3.22b ± 0.06
-0.22 ± 0.25
3.19b ± 0.06
-0.05 ± 0.27
Effect of the SNP genotypes on marbling and subcutaneous fat thickness in experiment 3
No. of animals
Subcutaneous fat thickness (mm)
2.88 ± 0.06
-0.04 ± 0.24
2.76 ± 0.04
0.23 ± 0.18
Based on 3 experiments, we showed that the g.231054C>T is associated with marbling in Japanese Black beef cattle, with the T allele resulting in high levels of marbling. Especially, in experiments 2 and 3, because the dams can be considered to represent a random sample of the Japanese Black population, the association is likely to be true. We used the predicted breeding values as phenotypic values, and then sire contributions to the predicted breeding values were the same for all offspring, even if the heterozygous sires were used in experiment 2. Thus, it is likely to be desirable to use larger numbers of dams in the analyses using paternal half-sib progeny steers such as experiments 2 and 3. We should note that the experiments in this study were the case. Further, the association of the g.231054C>T with marbling was corroborated by an independent study using 104 paternal half-sib Japanese Black families with a total of 821 progeny steers (S. Sukegawa, unpublished data).
Based on the association of the g.231054C>T with marbling, together with TTN expression difference between low-marbled (with g.231054C>T C allele) and high-marbled steer groups (with g.231054C>T C and T alleles) and co-localization of the marbling quantitative trait locus with the TTN, we can hypothesize that the SNP in the promoter region might have an impact on TTN expression and also marbling by affecting TTN promoter activity. However, the SNP is not located within any of as yet identified canonical sequences involved in gene transcription. Thus, a more likely event is that the TTN SNP is in linkage disequilibrium with an unidentified and truly relevant mutation, rather than functional and a causal mutation for marbling.
The effect of genotypes of the SNP was not statistically significant (P > 0.05) for subcutaneous fat thickness. Further, the marbling quantitative trait locus corresponding to chromosomal position of the TTN did not have statistically significant effect on subcutaneous fat thickness [10, 11]. Thus, it is likely that the TTN SNP is not associated with subcutaneous fat thickness in Japanese Black beef cattle. This might be supported by the fact that Japanese Black breed exhibits low genetic correlation between marbling and subcutaneous fat thickness .
Several previous studies have reported polymorphisms associated with beef marbling using beef cattle breeds other than Japanese Black [18–27]. We have recently reported that SNPs in the EDG1 gene were associated with marbling in Japanese Black breed . Thus, our present study seems to be the second report to show polymorphisms associated with marbling using Japanese Black breed.
The information on the TTN SNP obtained in this study, as well as the EDG1 SNPs, may be applied to effective marker-assisted selection to increase the levels of marbling in Japanese Black beef cattle. Further replicate studies will be needed to confirm the allelic association observed here, and to expand the results to evaluate all possible genotypic combinations of alleles.
This work was supported by a Grant-in-Aid for Scientific Research (B) (no. 14360166) from the Ministry of Education, Science, Sports, and Culture of Japan, and by the research funds of Japanese Livestock Technology Association.
- Boylston TD, Morgan SA, Johnson KA, Busboom JR, Wright RW, Reeves JJ: Lipid content and composition of Wagyu and domestic breeds of beef. J Agric Food Chem. 1995, 43: 1202-1207. 10.1021/jf00053a015.View ArticleGoogle Scholar
- Busboom JR, Jeremiah LE, Gibson LL, Johnson KA, Gaskins CT, Reeves JJ, Wright RW: Effects of biological source on cooking and palatability attributes of beef produced for the Japanese market. Meat Sci. 1993, 35: 241-258. 10.1016/0309-1740(93)90054-L.View ArticlePubMedGoogle Scholar
- Matsuishi M, Fujimori M, Okitani A: Wagyu beef aroma in Wagyu (Japanese Black cattle) beef preferred by the Japanese over imported beef. Anim Sci J. 2001, 72: 498-504.Google Scholar
- Japan Meat Grading Association: New Beef Carcass Grading Standards. 1988, Tokyo, Japan: Japan Meat Grading AssociationGoogle Scholar
- Liu Y, Wada R, Yamashita T, Mi Y, Deng CX, Hobson JP, Rosenfeldt HM, Nava VE, Chae SS, Lee MJ, Liu CH, Hla T, Spiegel S, Proia RL: Edg-1, the G protein-coupled receptor for sphingosine-1-phosphate, is essential for vascular maturation. J Clin Invest. 2000, 106: 951-961. 10.1172/JCI10905.PubMed CentralView ArticlePubMedGoogle Scholar
- Sasaki Y, Nagai K, Nagata Y, Doronbekov K, Nishimura S, Yoshioka S, Fujita T, Shiga K, Miyake T, Taniguchi Y, Yamada T: Exploration of genes showing intramuscular fat deposition-associated changes in musculus longissimus muscle. Anim Genet. 2006, 37: 40-46. 10.1111/j.1365-2052.2005.01380.x.View ArticlePubMedGoogle Scholar
- Sasaki Y, Yamada T, Miyake T: Quantitative and molecular genetic approaches for the improvement of carcass traits in the Wagyu cattle. Proceedings of the 8th World Congress on Genetics Applied to Livestock Production: 13–18 August 2006; Belo Horizonte, MG, Brazil. Edited by: Madalena FE. 2006, Belo Horizonte MG, Brazil: Federal University of Minas Gerais, 13-01.Google Scholar
- Yamada T, Itoh M, Nishimura S, Taniguchi Y, Miyake T, Sasaki S, Yoshioka S, Fujita T, Shiga K, Morita M, Sasaki Y: Association of single nucleotide polymorphisms in endothelial differentiation, sphingolipid G-protein-coupled receptor, 1 gene with marbling in Japanese Black beef cattle. Anim Genet. 2009, 40: 209-216. 10.1111/j.1365-2052.2008.01822.x.View ArticlePubMedGoogle Scholar
- Itoh-Satoh M, Hayashi T, Nishi H, Koga Y, Arimura T, Koyanagi T, Takahashi M, Hohda S, Ueda K, Nouchi T, Hiroe M, Marumo F, Imaizumi T, Yasunami M, Kimura A: Titin mutations as the molecular basis for dilated cardiomyopathy. Biochem Biophys Res Commun. 2002, 291: 385-393. 10.1006/bbrc.2002.6448.View ArticlePubMedGoogle Scholar
- Takasuga A, Watanabe T, Mizoguchi Y, Hirano T, Ihara N, Takano A, Yokouchi K, Fujinaka A, Chiba K, Kobayashi N, Tatsuda K, Oe T, Furukawa-Kuroiwa M, Nishimura-Abe A, Fujita T, Inoue K, Mizoshita K, Ogino A, Sugimoto Y: Identification of bovine QTL for growth and carcass traits in Japanese Black cattle by replication and identical-by-descent mapping. Mamm Genome. 2007, 18: 125-136. 10.1007/s00335-006-0096-5.View ArticlePubMedGoogle Scholar
- Watanabe T, Takasuga A, Mizoguchi Y, Hirano T, Ihara N, Takano A, Yokouchi K, Fujinaka A, Chiba K, Kobayashi N, Tatsuda K, Oe T, Furukawa M, Nishimura-Abe A, Imai K, Fujita T, Nagai S, Inoue K, Mizoshita K, Ogino A, Sugimoto Y: A comprehensive QTL map of Japanese Black cattle (Wagyu). Proceedings of the 29th International Conference on Animal Genetics: 11–16 September 2004; Tokyo, Japan. Edited by: Sasaki Y. 2004, Kyoto, Japan: Kyoto University, 136-Google Scholar
- Yamada T, Taniguchi Y, Nishimura S, Yoshioka S, Takasuga A, Sugimoto Y, Sasaki Y: Radiation hybrid mapping of genes showing intramuscular fat deposition-associated expression changes in bovine musculus longissimus muscle. Anim Genet. 2006, 37: 184-185. 10.1111/j.1365-2052.2006.01426.x.View ArticlePubMedGoogle Scholar
- Shiga K, Fujita T, Hirose K, Sasae Y, Nagai T: Production of calves by transfer of nuclei from cultured somatic cells obtained from Japanese Black bulls. Theriogenology. 1999, 52: 527-535. 10.1016/S0093-691X(99)00149-1.View ArticlePubMedGoogle Scholar
- Oita Prefectural Institute of Animal Industry: Sire Summary. 1999, Oita, Japan: Oita Prefectural Institute of Animal IndustryGoogle Scholar
- Nomenclature for the description of sequence variations. [http://www.hgvs.org/mutnomen/]
- Sasaki Y, Miyake T, Gaillard C, Oguni T, Matsumoto M, Ito M, Kurahara T, Sasae Y, Fujinaka K, Ohtagaki S, Dougo T: Comparison of genetic gains per year for carcass traits among breeding programs in the Japanese Brown and the Japanese Black cattle. J Anim Sci. 2006, 84: 317-323.PubMedGoogle Scholar
- Boldman KG, Kriese LA, Van Vleck LD, Van Tassell CP, Kachman SD: A Manual for Use of MTDFREML, A Set of Programs to Obtain Estimates of Variances and Covariances. 1995, Washington, DC: Agricultural Research Service, United States Department of AgricultureGoogle Scholar
- Barendse W: Assessing lipid metabolism. Patent, International Publication Number: WO 99/23248. 1999, World International Property OrganizationGoogle Scholar
- Barendse W, Bunch RJ, Harrison BE, Thomas MB: The growth hormone 1 GH1:c.457C>G mutation is associated with intramuscular and rump fat distribution in a large sample of Australian feedlot cattle. Anim Genet. 2006, 37: 211-214. 10.1111/j.1365-2052.2006.01432.x.View ArticlePubMedGoogle Scholar
- Barendse W, Bunch RJ, Kijas JW, Thomas MB: The effect of genetic variation of the retinoic acid receptor-related orphan receptor C gene on fatness in cattle. Genetics. 2007, 175: 843-853. 10.1534/genetics.106.064535.PubMed CentralView ArticlePubMedGoogle Scholar
- Buchanan FC, Fitzsimmons CJ, Van Kessel AG, Thue TD, Winkelman-Sim DC, Schmutz SM: Association of a missense mutation in the bovine leptin gene with carcass fat content and leptin mRNA levels. Genet Sel Evol. 2002, 34: 105-116. 10.1051/gse:2001006.PubMed CentralView ArticlePubMedGoogle Scholar
- Casas E, Bennett GL, Smith TP, Cundiff LV: Association of myostatin on early calf mortality, growth, and carcass composition traits in crossbred cattle. J Anim Sci. 2004, 82: 2913-2918.PubMedGoogle Scholar
- Hale CS, Herring WO, Shibuya H, Lucy MC, Lubahn DB, Keisler DH, Johnson GS: Decreased growth in Angus steers with a short TG-microsatellite allele in the P1 promoter of the growth hormone receptor gene. J Anim Sci. 2000, 78: 2099-2104.PubMedGoogle Scholar
- Jiang Z, Kunej T, Michal JJ, Gaskins CT, Reeves JJ, Busboom JR, Dovc P, Wright RW: Significant associations of the mitochondrial transcription factor A promoter polymorphisms with marbling and subcutaneous fat depth in Wagyu × Limousin F2 crosses. Biochem Biophys Res Commun. 2005, 334: 516-523. 10.1016/j.bbrc.2005.06.120.View ArticlePubMedGoogle Scholar
- Michal JJ, Zhang ZW, Gaskins CT, Jiang Z: The bovine fatty acid binding protein 4 gene is significantly associated with marbling and subcutaneous fat depth in Wagyu × Limousin F2 crosses. Anim Genet. 2006, 37: 400-402. 10.1111/j.1365-2052.2006.01464.x.View ArticlePubMedGoogle Scholar
- Nkrumah JD, Li C, Yu J, Hansen C, Keisler DH, Moore SS: Polymorphisms in the bovine leptin promoter associated with serum leptin concentration, growth, feed intake, feeding behavior, and measures of carcass merit. J Anim Sci. 2005, 83: 20-28.PubMedGoogle Scholar
- Thaller G, Kuhn C, Winter A, Ewald G, Bellmann O, Wegner J, Zuhlke H, Fries R: DGAT1, a new positional and functional candidate gene for intramuscular fat deposition in cattle. Anim Genet. 2003, 34: 354-357. 10.1046/j.1365-2052.2003.01011.x.View ArticlePubMedGoogle Scholar
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.