Study of DGAT1 Gene Polymorphisms with Carcass Traits in Iranian Zel and Lori-Bakhtiari Sheep Breeds

Document Type: Research Articles

Authors

1 Department of Animal Science, College of Agriculture and Natural Resources, University of Tehran, Karaj, Iran

2 Department of Animal Science, Faculty of Agriculture, Urmia University, Urmia, Iran

3 Department of Animal Science, Faculty of Agriculture, Payame Noor University, Tehran, Iran

Abstract

Diacylglycerol acyltransferase1 (DGAT1) plays an important role in the metabolism of triglycerides which catalyze the final step of triglyceride synthesis in animals. The objective of this study was to investigate the single-nucleotide polymorphisms (SNPs) in 5'UTR, exon-1, and exon-2 of DGAT1 in two Iranian indigenous sheep breeds. A total of 309 animals including fat-tailed Lori-Bakhtiari (n=152) and thin-tailed Zel (n=157) were used in this study. The genotypic patterns were detected by polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP). Five SSCP patterns were detected for 5'UTR and exon-1 fragment by PCR-SSCP and subsequently confirmed by sequencing PCR products. The sequencing results revealed that there are three novel polymorphisms in 5'UTR and exon-1fragment of DGAT1 at the studied breeds. Out of the detected polymorphisms only A277G substitution in exon-1 of DGAT1 leads to the changes in amino acids (p.Arg26Gly). There was significant correlation (p < 0.05) between fat-tail weight (FTW) and back-fat thickness (BFT) and the observed genotypes in Lori-Bakhtiari breed; therefore, animals with G5 pattern had higher FTW and BFT compared to G1 pattern. The G1 and G5 genotypic patterns or haplotypes were different at their position 101 of 5'UTR region. No significant relationship (p < 0.05) was found between the detected genotypes of 5'UTR and exon-1 fragment of DGAT1 in Zel breed and carcass traits. These results revealed that detected DGAT1 novel SNPs had significant effects on carcass traits and they can be used as a marker for these traits.

Keywords


INTRODUCTION

The main goal of sheep breeding programs is to improve meat production and carcass traits (Kosgey et al. 2006). During the last forty years, the genetic improvement of livestock has been achieved by selection based on phenotypic information, but getting fast genetic improvement is difficult by this method. However marker-assisted selection (MAS) can increase the genetic improvement of polygenic traits which results from the presence of major genes and their effect on the performance of these traits. Therefore, this method has provided opportunities to increase response to selection, especially for traits with small heritability (Dekkers, 2004). Lori-Bakhtiari breed, with more than 1.7 million heads, is one of the most common indigenous sheep breeds in Iran. This breed adapted to the hilly and mountainous regions in the west of Iran and has been raised mostly in villages under semi-intensive systems. The Lori-Bakhtiari breed has high body weight and the largest fat tail compared to other Iranian Indigenous breeds (Kianzad et al. 2003; Vatankhah and Talebi, 2008). Zel breed, with about 1.5 million heads, is the only thin-tailed Iranian breed. Zelis raised in the north part of Iran located between the northern slopes of the Alburz mountain and the Caspian Sea (Kamalzadeh et al. 2008). Diacylglycerol acyltransferase1 (DGAT1) gene is involved in the metabolism of triglycerides and catalyzes the last step of the triacylglycerol synthesis (Hatzopoulos et al. 2011; Li et al. 2013). DGAT1 is extensively expressed in almost all tissues including the mammary glands, small intestine, adipose tissue, and skeletal muscle, with the highest expression level in adipose tissue and small intestine (Xu et al. 2008). The DGAT1 gene spans approximately 8.5 kbp and comprises 17 exons 16 introns, located on chromosome 14 of bovine (Grisart et al. 2002). The homologous sequence was found on chromosome 9 of sheep (Ovis aries). This gene contains about 1470 nucleotides that code a polypeptide with 489 amino acids (Xu et al. 2009). Specific polymorphisms in DGAT1 were detected and their association with some production traits were confirmed in farm animals, especially cattle. Then, DGAT1 was introduced as a candidate gene for milk production traits in dairy (Banos et al. 2008; Li et al. 2009; Signorelli et al. 2009; Cerit et al. 2014) and carcass traits (Anton et al. 2010; Souza et al. 2010; Pannier et al. 2010; Li et al. 2013; Tait et al. 2014; Ardicli et al. 2017) in cattle.Moreover, DGAT1 due to its key role in fat metabolism (Scata et al. 2009) is a candidate gene for milk fat content (Scata et al. 2009; Xu et al. 2009; Yang et al. 2011; An et al. 2013; Tabaran et al. 2014; Martin et al. 2017) and carcass traits (Xu et al. 2009; Mohammadi et al. 2013) in dairy sheep and goat. However, there are only a few studies on the association between DGAT1 polymorphisms and carcass traits in sheep (Xu et al. 2009; Mohammadi et al. 2013; Noshahr and Rafat, 2014; Nanekarani et al. 2016). Xu et al. (2009) sequenced the whole coding region of sheep DGAT1 and surveyed the associations between polymorphism of DGAT1 and 11 meat production traits. They reported thepositive effect of the T allele in exon-17 of DGAT1 on meat quality and quantity in three Chinese sheep breeds. Mohammadi et al. (2013) reported that there is a significant correlation betweenexon-17 of DGAT1 polymorphisms (T487C) and carcass weight and dressing percentage, and concluded that the CC genotypes had higher carcass weight and dressing percentage. Few studies have been carried out on the association of DGAT1 polymorphisms with carcass traits in sheep breeds. All of these studies suggested that further investigations are necessary to make a definitive statement about the association of this gene polymorphisms with carcass traits in sheep (Xu et al. 2009; Noshahr and Rafat, 2014; Mohammadi et al. 2013). The aim of this study is to investigate 5'UTR, exon-1 and exon-2 regions of DGAT1 polymorphisms and their association with carcass traits in Lori-Bakhtiari and Zel Iranian indigenous breeds.

 

MATERIALS AND METHODS

Animals and data collection

A total of 309 blood samples were collected from Lori-Bakhtiari (n=152) and Zel (n=157) sheep breeds. Lori-Bakhtiari and Zel samples were taken from two industrial slaughterhouses located in the Qom and Mazandaran provinces of Iran, respectively. Before slaughter, body weight (SBW) of lambs was measured. After slaughter, carcass weight (CW), fat-tail weight (FTW), carcass weight without fat-tail (CWFT), and back fat thickness (BFT) immediately were measured and fat-tail percentage (FTP), dressing percentage (DP) and dressing percentage without fat-tail (DPWFT) were calculated. Genomic DNA was extracted by the Salting-out procedure (Miller et al. 1998). The quality and concentration of extracted DNA were evaluated by visualizing on 0.8% agarose gel under UV light and spectrophotometer (Thermo, NanoDrop 1000).

 

DNA amplification with PCR

Two primer pairs DGAT1 (5'UTR and exon-1) F:5'-GGAACTACGCTTCCCAGGAC-3'; R:5'-ACGTCTCCGTCCTTGTCTGT-3'; and (exon-2) F:5'-GTCTTGCATCACCAGCTCCT-3'; R:5'-CAGGCATCTACTGGGATTCAG-3' were designed based on the published sheep sequences (GenBank accession no.EU178818) to amplify the determined fragment by Primer-BLAST (Ye et al. 2012). The length of the amplified fragments were 360 bp (5'UTR and partial exon-1) and 445 bp (exon-2) in both breeds. PCR was carried out on 25 µL volume of the following ingredients: 1 × PCR buffer (16 mM (NH4)2SO4, 67 mM Tris-HCl pH 8.8, 0.1% Tween-20), 2 mM MgCl2, 200 µM dNTP, 6 pmol of each primer, 1.5 unit Taq polymerase (GenetBio, Korea) and 50-100 ng of genomic DNA. PCR reactions were performed with a Genecycler (Bio-Rad, My Cycler, USA) with the amplification program consisted of an initial denaturation at 96 ˚C for 2 mins, then 32 cycles of 96 ˚C for 1 min, 57 ˚C for 45 seconds and 72 ˚C for 1 min, and a final extension of 72 ˚C was maintained for 10 mins for the exon-1 of DGAT1 gene. For exon-2 fragment, the thermal cycling conditions consisted of initial denaturation at 96 ˚C for 5 mins, followed by 30 cycles of denaturation at 96 ˚C for 1 min, annealing at 60 ˚C for 45 seconds, extension at 72 ˚C for 1 min, with a final extension at 72 ˚C for 10 mins. The PCR products were identified by 2% agarose gel electrophorese and visualized using UV light.

 

Single-strand conformation polymorphism

SSCP method was applied to detect SNPs in 5'UTR, partial exon-1 and exon-2 of the DGAT1 gene. For SSCP, 5µl PCR products were diluted in 5µl denaturing solution; the mixture was centrifuged and denatured at 98 ˚C for 10 mins, then rapidly chilled on ice for 5 mins and detected on 12% polyacrylamide gel. The electrophoresis was performed in a vertical unit (Bio-rad), in 1× TBE buffer at 4 ˚C. The gel running conditions for the 360 bp fragment of DGAT1 was 350 V for 18hours and 445 bp fragment of DGAT1 was 380 V for 21 hours, respectively. Then, DNA fragments were stained with 0.1% silver nitrate and visualized in 3% NaOH solution (Zhang et al. 2007). Two samples from each SSCP patterns were sequenced in both directions (Bioneer, Seoul, South Korea) and were determined by Bio Edit (Version 7.2.0) and Vector NT1 software. Different SSCP pattern sequences were aligned according to Genbank (EU178818.1) by BLAST algorithm to find homologous sequences in NCBI databases (http:// www.ncbi.nlm.hih.gov/).

 

Statistical analysis

Genotypes patterns or haplotypes frequencies were calculated using GENALEX6.4 software (Peakall and Smouse, 2006). Two fixed models were used to make an association between 360 bp fragment of DGAT1 variants and carcass traits using a generalized linear model (GLM) procedure of SAS software (SAS, 2007). Exon-2 of DGAT1 was eliminated from association analysis due to monomorphic state. Tukey–Kramer’s multiple range tests were used for least-squares means (LSM) comparison. The statistical model used for determining CW, DP, BFT and DPWFT traits in two breeds was as follows:

Another model used for FTW, FTP and DP traits in Lori-Bakhtiari sheep breed was as follows:

Where:

yijklm: phenotypic value of carcass traits.

µ: overall population mean for each trait. ­

Ai: fixed effect of the ith age.

Sj: fixed effect of the jth sex.

Bk: fixed effect of the kth breed.

Gl: fixed effect of the lth genotype.

covariate fixed effect of weight at slaughter.

Sj(Bk): effect of the jth sex in kth breed.

(AG­(il: interaction between the ith age and the lth genotype. (SG(jl­: interaction between the jth sex and the lth genotype. (BG)kl: interaction between the kth breed and the lth genotype.

eijklm: residual random effect.

 

RESULTS AND DISCUSSION

Single strand confirmation polymorphism analysis

The results showed that both PCR products including 360 bp and 445 bp fragments, successfully amplified by the used primers. Five SSCP patterns were detected for exon-1 and its 5′UTR (Figure 1).

 

Figure 1 Five PCR-SSCP genotypes of exon-1 and its 5'UTR fragment (360 bp) of DGAT1 gene for Lori-Bakhtiari and Zel breeds

 

 

Furthermore, one SSCP pattern was observed for exon-2 and this segment was monomorphic in both breeds (Figure 2).

 

Figure 2 PCR-SSCP genotype of exon-2 fragment (445 bp) of DGAT1 gene for Lori-Bakhtiari and Zel breeds

 

 

The frequency of the observed genotypic patterns of 360 bp fragment for each breed is shown in Table 1. The G5 and G2 had higher frequency in Lori-Bakhtiari and Zel breeds, respectively. The G3 had the lowest frequency in both breeds. Also, the frequency of G2, G4, and G5 was different between Lori-Bakhtiari and Zelbreeds.

 

DNA sequence analysis

In this study, sequences of PCR amplicons from five different SSCP genotypic patterns of 360 bp fragment variants were analyzed and compared with the NCBI reference sequences (Table 2). For 360 bp fragment (including 5′UTR and exon-1), three novel single nucleotide polymorphisms (SNPs) were detected in the studied breeds. One of the discovered SNPs is located in the coding region (gA277G) and two of them located in the 5′UTR region (g.T101- and C129A) of DGAT1 (Table 2). Moreover, the protein sequences of exon-1 region were aligned to NCBI reference sequence (Figure 3). The only substitution in exon-1 of G3 pattern (A 277G) leads to the changes in amino acids (p.Arg54Gly) (Figure 3).

 

Correlation between DGAT1 polymorphisms and carcass traits in sheep

In the present study, exon-2 of DGAT1 was monomorphic and omitted from association studies. Association of the 360bp fragment of DGAT1 polymorphisms with carcass traits was analyzed (Table 3). Significant relationship (P<0.05) was found between FTW and BFT traits in the Lori-Bakhtiari breed, and animals with G5 had higher FTW and BFT compared to G1 pattern. Although no significant differences (P>0.05) were observed for DP and FTP among five identified genotypic patterns in Lori-Bakhtiari, animals with G3 pattern for DP and G5 pattern for FTP had the highest performance compared to other genotypes. No significant association (P<0.05) was found between observed genotypes of exon-1 of DGAT1 in Zel breed and the studied carcass traits. DGAT1 was selected as a candidate gene due to its effects on carcass traits that had been confirmed in previous studies on cattle (Pannier et al. 2010; Souza et al. 2010; Li et al. 2013; Tait et al. 2014) and sheep (Xu et al. 2009; Mohammadi et al. 2013; Noshahr and Rafat, 2014). The detection of SNPs in 360 bp fragment was consistent with results of Scata et al. (2009) who amplified and sequenced whole of the DGAT1 gene using 18 primer pairs and reported novel SNPs at 5'UTR (g. C127A) in Sarda, Altamurana and Gentile di Puglia sheep breeds. But their result was not consistent with present study for exon-1. The detection of SNPs in 360 bp fragment was consistent with the results of Scata et al. (2009) who identified novel SNP in Sarda, Altamurana and Gentile di Puglia sheep breeds. Furthermore, results indicated that DGAT1 had significant genetic effects on fat content in Lori-Bakhtiari breed are consistent with the findings of Xu et al. (2009), Mohammadi et al. (2013) and Noshahr and Rafat (2014). Based on the sequence analysis, two differences, including one deletion and one substitution, were found in 5'UTR region between detected genotypes (Table 2). Regulation of DGAT1 expression in adipocytes mostly happens at the transcriptional and post-transcriptional steps (Yu et al. 2002), therefore mutations in 5'UTR region may influence DGAT1 expression and consequently fat content.

 

Table 1 Genotypic patterns frequencies of single strand conformation polymorphism (SSCP) variants of 5′ UTR and exon-1 of DGAT1 gene in Zel and Lori-Bakhtiari breeds

 

 

Table 2 Nucleotide substitutions and types of single-nucleotide polymorphisms (SNPs) for exon-1 of DGAT1 gene in Zel and Lori-Bakhtiari breeds

 

 

Table 3 Association of DGAT1 genotypes with carcass traits in Zel and Lori-Bakhtiari breeds

 

The means within the same row with at least one common letter, do not have significant difference (P>0.05).

SE: standard errors.

 

Figure 3 The amino acid sequence of exon-1 of DGAT1 gene for Lori-Bakhtiari and Zel breeds

 

 

Therefore, identified mutations in the 5'UTR region of DGAT1 could affect the function of DGAT1 enzyme in lipid metabolism. For instance, DGAT1 has considerable effect on fat deposition and fatty acid composition in sheep meat and milk (Scata et al. 2009). In the beginning cod on (ATG) has started from position 202, mutations have occurred before it was not able to change the protein sequence and only the mutation occurred in the exon-1 region of the DGAT1 gene which was able to change the amino acid sequence of translated protein derived from this gene. Therefore, mutation only occurred in position 277 of genotypic pattern G3 and led to the transformation of amino acid arginine to glycine (AGG to GGG). However, there was no significant difference (P>0.05) between this genotypic pattern and other genotypes. This result indicated that despite the change in the sequence of the enzyme, performance did not change significantly. Therefore, this SNP could not be introduced as a marker for the studied traits. By the way, future investigations can be useful to detect possible effects of this SNP on the production traits of sheep. Only CW, DP and BFT traits were analyzed in Zel breed because this breed does not have any fat-tail for measuring other traits. In Zel breed, no significant difference (P> 0.05) was found among the observed five genotypes in carcass traits, but G3 pattern showed higher DP and lower BFT compared to others. The lack of significant differences between most of the identified genotypic patterns may be due to the high standard errors found for some genotypic patterns, resulting from the low number of animals for these genotypes. Therefore, increasing the number of animals will probably change the results. The effects of different genotypes on common carcass traits in Zel and Lori-Bakhtiari breeds were similar, when analyzed separately or with each other. Possible contradictory results can be due to the influences of different genomic background and different genotypic frequencies.

 

CONCLUSION

In this study, five different PCR-SSCP genotypic patterns were detected for 360 bp fragment (including 5′UTR and exon-1) and three novel SNPs were detected for this region. There was a significant correlation (P<0.05) between FTW and BFT traits and observed genotypes in Lori-Bakhtiari breed, but there was no significant correlation (P<0.05) between detected genotypes of exon-1 of DGAT1 in Zel breed and carcass traits. These results indicated that DGAT1 is a potential candidate gene for carcass traits in MAS and breeding programs to improve FTW and BFT traits by selection in Lori-Bakhtiari breed. Although significant correlation (P<0.05) was found between carcass fat content and DGAT1 gene, but it can be due to the effect of other mutations in regions that were not investigated in present study (Xu et al. 2008) and also another QTLs close to the DGAT1 (Coppieters et al. 1998). Further investigation is necessary to make a conclusive decision about using the detected polymorphisms of this gene in breeding programs for carcass traits in sheep.

 

ACKNOWLEDGEMENT

The authors thank all the teams who worked on the experiments and provided technical assistance in the laboratory during this study.

An X.P., Hou J.X., Zhao H.B., Bai L., Peng J.Y., Zhu C.M., Yan Q.M., Song Y.X., Wang J.G. and Cao B.Y. (2013). Polymorphism identification in goat DGAT1 and STAT5A genes and association with milk production traits. Czech J. Anim. Sci. 58, 321-327.

Anton I., Kovacs K., Hollo G., Farkas V., Lehel L., Hajda Z. and Zsolnai A. (2010). Effect of leptin, DGAT1 and TG gene polymorphisms on the intramuscular fat of Angus cattle in Hungary. Livest. Sci. 11, 243-249.

Ardicli S., Dincel D., Samli H. and Balci F. (2017). Effects of polymorphisms at LEP, CAST, CAPN1, GHR, FABP4 and DGAT1 genes on fattening performance and carcass traits in Simmental bulls. Arch Tierz. 60(2), 61-70.

Banos G., Woolliams J., Woodward B., Forbes A. and Coffey M. (2008). Impact of single nucleotide polymorphisms in leptin, leptin receptor, growth hormone receptor, and diacylglycerol acyltransferase (DGAT1) gene loci on milk production, feed, and body energy traits of UK dairy cows. J. Dairy Sci. 91, 3190-3200.

Cerit H., Dumen E. and Sezgin F.H. (2014). Comparison of DGAT1 K232A polymorphism and its effects on some milk quality parameters in Holstein and Native Black Race cattles. Kafkas Univ. Vet. Fak. Derg. 20, 301-305.

Coppieters W., Riquet J., Arranz J.J., Berzi P., Cambisano N., Grisart B. Karim L., Marcq F., Moreau L. and Nezer C. (1998). A QTL with major effect on milk yield and composition maps to bovine chromosome 14. Mamm. Genome. 9, 540-544.

Dekkers J.C. (2004). Commercial application of marker and gene-assisted selection in livestock: strategies and lessons. J. Anim. Sci. 82, 313-328.

Grisart B., Coppieters W., Farnir F., Karim L., Ford C., Berzi P., Cambisano N., Mni M., Reid S., Simon P., Spelman R., Georges M. and Snell R. (2002). Positional candidate cloning of a QLT in dairy cattle: Identification of a missense mutation in the bovine DGAT1 gene with major effect on milk yield and composition. Genome Res. 12, 222-231.

Hatzopoulos P., Banilas G., Karampelias M., Makariti I. and Kourti A. (2011). The olive DGAT2 gene is developmentally regu-lated and shares overlapping but distinct expression patterns with DGAT1. J. Exp. Bot. 62, 521-532.

Kamalzadeh A., Rajabaigy M. and Kiasat A. (2008). Livestock production systems and trends in livestock industry in Iran. J. Agric. Soc. Sci. 4, 183-188.

Kianzad M.R., Panandam J.M., Emamjomeh Kashan N., Jelan Z.A. and Dahlan I. (2003). Reproductive performance of three Iranian sheep breeds. Asian-Australian J. Anim. Sci. 16, 11-14.

Kosgey I.S., Baker R.L., Udod H.M.J. and Van-Arendonk J.A.M. (2006). Successes and failures of small ruminant breeding programmes in the tropics: a review. Small Rumin. Res. 61, 13-28.

Li J., Xu X., Zhang Q., Wang X., Deng G., Fang X., Gao X., Ren H. and Xua S. (2009). Association between single nucleotide polymorphisms in the Dgat2 gene and beef carcass and quality traits in commercial feedlot steers. Asian-Australasian J. Anim. Sci. 22, 943-954.

Li X., Ekerljung M., Lundstrom K. and Lunden A. (2013). Association of polymorphisms at DGAT1, leptin, SCD1, CAPN1 and CAST genes with color, marbling and water holding capacity in meat from beef cattle populations in Sweden. Meat Sci. 94, 153-158.

Martin P., Palhière I., Maroteau C., Bardou P., Canale-Tabet K., Sarry J., Woloszyn F., Bertrand-Michel J., Racke I. and Besir H. (2017). A genome scan for milk production traits in dairy goats reveals two new mutations in Dgat1 reducing milk fat content. Sci. Rep. 7(1), 1-13.

Miller S., Dykes D. and Polesky H. (1988). A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res. 16(3), 1215-1224.

Mohammadi H., Moradi-Shahrebabak M. and Sadeghi M. (2013). Association between single nucleotide polymorphism in the ovine DGAT1 gene and carcass traits in two Iranian sheep breeds. Anim. Biotechnol. 24, 159-167.

Nanekarani S., Kolivand M. and Goodarzi M. (2016). Polymorphism of a mutation of DGAT1 gene in Lori sheep breed. J. Adv. Agric. Technol. 3(1), 38-41.

Noshahr F.A. and Rafat A. (2014). Polymorphism of DGAT1 gene and its relationship with carcass weight and dressing percentage in Moghani sheep breed. Iranian J. Appl. Anim. Sci. 4(2), 331-334.

Pannier L., Mullen A.M., Hamill R.M., Stapleton P.C. and Sweeney T. (2010). Association analysis of single nucleotide polymorphisms in DGAT1, TG and FABP4 genes and intramuscular fat in crossbred Bos taurus cattle. Meat Sci. 85, 515-518.

Peakall R. and Smouse P.E. (2006). GENALEX 6: genetic analysis in excel. Population genetic software for teaching and research. Mol. Ecol. Notes. 6, 288-295.

SAS Institute. (2007). SAS®/STAT Software, Release 9.1. SAS Institute, Inc., Cary, NC. USA.

Scata M.C., Napolitano F., Casu S., Carta A., De-Matteis G., Signorelli F., Annicchiarico G., Catillo G. and Moioli B. (2009). Ovine acyl CoA: Diacylglycerol acyltransferase 1 – molecular characterization, polymorphisms and association with milk traits. Anim. Genet. 40, 737-742.

Signorelli F., Orru L., Napllitano F., Matteis G.D., Scata M.C., Catillo G., Marchitelli C. and Moioli, B. (2009). Exploring polymorphisms and effects on milk traits of the DGAT1, SCD1 and GHR genes in four cattle breeds. Livest. Sci. 125, 74-79.

Souza F.R.P., Mercadante M.E.Z., Fonseca L.F.S., Ferreira L.M.S., Regatieri I.C., Ayres D.R., Tonhati H., Silva S.L., Razook A.G. and Albuquerque L.G. (2010). Assessment of DGAT1 and LEP gene polymorphisms in three Nelore (Bos indicus) lines selected for growth and their relationship with growth and carcass traits. J. Anim. Sci. 88, 435-441.

Tăbăran A., Mihaiu M., Dan S.D., Reget O., Pivariu B., Cordis I., Cordea D. and Muresan C. (2014). Identification of polymorphism in goat and sheep DGAT1 gene associated with milk production traits. Bull. UASVM Vet. Med. 71(2), 283-286.

Tait R.G., Shackelford S.D., Wheeler T.L., King D.A., Keele J.W., Casas E., Smith T.P. and Bennett G.L. (2014). CAPN1, CAST, and DGAT1 genetic effects on preweaning performance, carcass quality traits, and residual variance of tenderness in a beef cattle population selected for haplotype and allele equalization. J. Anim. Sci. 92, 5382-5393.

Vatankhah M. and Talebi M.A. (2008). Heritability estimates and correlations between production and reproductive traits in Lori-Bakhtiari sheep in Iran. South African J. Anim. Sci. 38(2), 110-118.

Xu Q.L., Che Y.L., Ma R.X. and Xue P. (2009). Polymorphism of DGAT1 associated with intramuscular fat-mediated tenderness in sheep. J. Sci. Food Agric. 89, 232-237.

Yang J., Manolio T.A., Pasquale L.R., Boerwinkle E., Caporaso N. and Cunningham J.M. (2011). Genome partitioning of genetic variation for complex traits using common SNPs. Nat. Genet. 43, 519-525.

Ye J., Coulouris G., Zaretskaya I., Cutcutache I., Rozen S. and Madden T.L. (2012). Primer-BLAST: A tool to design target-specific primers for polymerase chain reaction. BMC Bioinformatics. 13, 134-142.

Yu Y.H., Zhang Y., Oelkers P., Sturley S.L., Rader D.J. and Ginsberg H.N. (2002). Posttranscriptional control of the expression and function of diacylglycerol acyltransferase-1 in mouse adipocytes. J. Biol. Chem. 277, 50876-50884.

Zhang C.L., Wang Y., Chen H., Lan X. and Lei C. (2007). Enhance the efficiency of single-strand conformation polymorphism analysis by short polyacrylamide gel and modified silver staining. Anal. Biochem. 365(2), 286-287.