J. Dairy Sci. 2009. 92:773-776. doi:10.3168/jds.2008-1266
© 2009 American Dairy Science Association ®
Technical note: A new and cost-effective method for detection of the bovine acyl-CoA:diacylglycerol acyltransferase 1 K232A polymorphism in cattle
R. S. Steinberg*,
L. Pereira*,
G. A. Lacorte*,
M. G. C. D. Peixoto
,
R. S. Verneque
,
R. L. Teodoro
,
M. A. Machado
,
C. G. Fonseca* and
M. R. S. Carvalho*,1
* Departamento de Biologia Geral, Universidade Federal de Minas Gerais, Belo Horizonte, 31-270-901, Brazil
Embrapa Gado de Leite, Juiz de Fora, 36038-330, Brazil
1 Corresponding author: mraquel{at}icb.ufmg.br
 |
ABSTRACT
|
|---|
A new, quick, and inexpensive method for detecting the bovine acyl-CoA:diacylglycerol acyltransferase1 (DGAT1) polymorphism (K232A) through tetra-primer amplification refractory mutation system by PCR (ARMS-PCR) was developed in the present investigation. The DGAT1 gene was recently identified as underlying variation in milk production traits. To date, PCR techniques such as PCR-RFLP have been used for detecting the DGAT1 K232A polymorphism, despite being expensive and laborious. The method proposed here, a tetra-primer ARMS-PCR, showed 100% sensitivity and specificity when compared with PCR-RFLP results obtained in a sample of 80 animals tested in a double-blind system. Our results suggest that the use of tetra-primer ARMS-PCR for DGAT1 K232A polymorphism genotyping could greatly reduce costs providing information for both research purposes and for dairy cattle breeders who perform DGAT1 genotyping for gene-assisted selection.
Key Words: DGAT1 polymerase chain reaction-restriction fragment length polymorphism tetra-primer amplification refractory mutation system-polymerase chain reaction marker-assisted selection
The acyl-CoA:diacylglycerol acyltransferase 1 (DGAT1)geneencodesthediacylglycerol-O-transferase (EC 2.3.1.20) microsomal enzyme, which catalyzes the last and limiting step in the synthesis of triglyceride (Cases et al., 1998). In cattle, the QTL associated with milk production traits was mapped to the centromeric region of the bovine chromosome 14 (Riquet et al., 1999). One of the genes in the candidate interval is DGAT1. Therefore, DGAT1 constitutes both a positional and a functional candidate gene that accounts for variation in milk yield and composition. Moreover, DGAT1–/– knockout mice have been reported that do not produce milk, providing additional evidence for the function of DGAT1 in lactation (Smith et al., 2000). Studies have reported correlations between the DGAT1 K232A polymorphism and milk production or composition in some Bos taurus breeds (Grisart et al., 2002; Winter et al., 2002; Kühn et al., 2004), such as the lysine-encoding allele (232K) that was associated with decreasing protein and milk yields and increasing fat yield. The alanine-encoding allele (232A) was associated with increasing protein and milk yields and decreasing fat yield. This polymorphism was also correlated with intramuscular fat deposition as reported for German Holstein and Charolais breeds, for which the lysine-encoding allele was associated with increased marbling score in the musculus semitendinosus (Thaller et al., 2003). Of note, DGAT1 232K was recently related to increased content of saturated fat in milk (Schennink et al., 2007, 2008). Increasing amounts of unsaturated fatty acids in milk is an important selection objective because milk is regarded as one of the major sources of saturated fat in the human diet in Europe (Hulshof et al., 1999). Therefore, selective breeding could make a significant contribution to changing the fat composition of cows milk, which affects the human diet (Schennink et al., 2007).
Three methods have already been described for genotyping the DGAT1 K232A polymorphism. The first is based on PCR-RFLP using the CfrI restriction enzyme (Winter et al., 2002), which is an expensive and time-consuming method that requires a relatively large amount of the restriction enzyme to prevent partial digestion; the second is a laborious technique based on an oligonucleotide ligation assay (Grisart et al., 2002); and the third method is based on quantitative PCR (Schennink et al., 2007), and therefore, a real-time thermocycler is needed.
In the present study, a new, rapid, and inexpensive method for the detection of the DGAT1 K232A polymorphism based on the tetra-primer amplification refractory mutation system-PCR (ARMS-PCR; Newton et al., 1989; Ye et al., 2001) is described. Usually, detection of SNP by tetra-primer ARMS-PCR relies on the pairing of the 3' nucleotide of the primer with the SNP. To enhance the annealing specificity, a mismatch is introduced in 1 of the 3 positions previous to that base-paired with the SNP. The DGAT1 K232A polymorphism is particularly suitable for the development of a tetra-primer ARMS-PCR system because it is produced by a double substitution (GC AA). Therefore, the use of an additional mismatch is not necessary. This method is depicted in Figure 1
.

View larger version (15K):
[in this window]
[in a new window]
|
Figure 1. Amplification of the polymorphism DGAT1 K232A by tetra-primer amplification refractory mutation system-PCR. Continuous lines represent the DNA fragment containing the DGAT1 K232A polymorphism. Arrows indicate the primers used in the PCR reactions. Bent arrows show the annealing failure of the allele specific primers due to 3' mismatches. Dashed lines indicate sizes of PCR products. Tetra-primer reactions produce 3 PCR products with 512 bp (external band), 369 bp (232K allele-specific band), and 181 bp (232A allele-specific band).
|
|
Basically, the system is composedof 4 primers, 2 external (DGAT1-Ala-For and DGAT1-Lys-Rev) and 2 internal (DGAT1-Ala-Rev and DGAT1-Lys-For), each one of the latter annealing at its 3' extremity to one of the alleles GC (232A) or AA (232K). Thus, both alleles are simultaneously amplified and subsequently identified by size after gel electrophoresis (Ye et al., 2001).
Genomic DNA was extracted by the proteinase-K:phenol-chloroform method (Sambrook and Russell, 2001) from blood or semen samples from 80 animals (33 of the Guzerat breed and 47 of the Holstein breed). Four primers were designed based on GenBank sequence AJ318490: DGAT1-Ala-For: 5'-gtc aac ctc tgg tgc cga gag-3'; DGAT1-Ala-Rev: 5'-agc tcc ccc gtt ggc cgc-3'; DGAT1-Lys-For: 5'-tcg tag ctt tgg cag gta aga a-3'; DGAT1-Lys-Rev: 5'-cac ctg gag ctg ggt gag gaa-3'. The PCR-amplified fragments were 512 bp (external band), 369 bp (232K allele), and 181 bp (232A allele). The external band functions as the internal PCR-positive control.
Tetra-primer ARMS-PCR amplifications were carried out in a final volume of 25 µ L with IV-B PCR buffer (10 mM Tris, pH 8.4, 50 mM KCl, 0.1% Triton X-100, 1.75 mM MgCl2), 0.2 mM each dNTP (Invitrogen, Carlsbad, CA), 8% dimethyl sulfoxide (J. T. Baker, Xalostoc, Edo. De Mexico, Mexico), 1 U of Taq DNA polymerase (Phoneutria Biotecnologia and Serviços, Belo Horizonte, MG, Brazil), 25 to 100 ng of genomic DNA, 0.2 µ M DGAT1-Ala-For and DGAT1-Lys-Rev primers, and 0.8 µ M DGAT1-Ala-Rev and DGAT1-Lys-For primers (Integrated DNA Technologies Inc., Coralville, IA). Additional primer proportions tested were 1:3 and 1:2 external to internal primers. The PCR cycling conditions were initial denaturation step for 5 min at 94°C, followed by 25 cycles at 94°C for 30 s, 65°C for 30 s, and 72°C for 30 s, and a final extension step for 5 min at 72°C. Reactions were undertaken using a PxE 0.2 Thermal Cycler (Electron Corporation, Waltham, MA). The PCR products were run and visualized on 2% ethidium bromide-stained agarose gels or on 8% silver stained polyacrylamide gels. Allele sizes were estimated by comparison with a FX174 RF DNA/HaeIII Fragments Ladder (Invitrogen). To confirm the identity and the size of amplified products, 3 independent PCR reactions were set up: 1) including the external primers only, 2) with the lysine-specific primers, and 3) including the alanine-specific primers. The sizes of the products obtained coincided with those observed in tetra-primer ARMS-PCR reactions. The products resulting from the 3 PCR reactions were sequenced using the Big Dye Terminator v3.1 Cycle Sequencing kit and a 3130 Genetic Analyzer (Applied Biosystems, Foster City, CA). Sequence identities were confirmed by homology through a BLAST search (http://www.ncbi.nlm.nih.gov/blast/blast.cgi). Sensitivity and specificity of the method were evaluated in a double-blind experiment, for which 80 animals were simultaneously genotyped by both PCR-RFLP using the CfrI restriction enzyme (Fermentas, Burlington, Ontario, Canada; Winter et al., 2002) and tetra-primer ARMS-PCR.
Tetra-primer ARMS-PCR amplifications yielded 3 bands as follows (Figure 2
): a fragment with 512 bp, produced by the external primers, regarded as the internal control for PCR amplification, a 369-bp fragment (lysine-encoding allele), and a 181-bp fragment (alanine-encoding allele). Concordant results were obtained for all of the 160 alleles tested by PCR-RFLP with CfrI restriction enzyme and tetra-primer ARMS-PCR, which showed 100% sensitivity and specificity. In the BLAST search, the 3 fragments sequenced showed homology with the bovine DGAT1 gene, confirming the identity of the alleles.

View larger version (47K):
[in this window]
[in a new window]
|
Figure 2. Ethidium bromide-stained (2%) agarose electrophoresis gel. Lane 1 = FX174 RF DNA/HaeIII Fragments ladder (Invitrogen, Carlsbad, CA). Molecular size markers are indicated with arrows: 603 and 310 bp. Lanes 2 to 20 = tetra-primer amplification refractory mutation system (ARMS)-PCR products. Individuals in lanes 2 and 13 are heterozygous; those in lanes 3, 4, 7, 8, 9, 10, and 12 are homozygous 232K/232K; and those in lanes 5 and 6 are homozygous 232A/232A. *Amplification failure. The sizes of the tetra-primer ARMS-PCR fragments are 512 bp (external band), 369 bp (232K allele-specific band), and 181 bp (232A allele-specific band). The allele K band is frequently weaker than the other 2. To confirm genotypes in these cases, the reactions were repeated with larger amounts of DNA.
|
|
The alanine-encoding allele has been suggested to be associated with a healthier milk, which means that it contains high protein levels and low fat, particularly low saturated fat. Large frequencies of this allele have been reported for many Bos taurus breeds. However, smaller frequencies of this allele were found in Bos indicus: 4% in Gyr and 2.5% in Red Sindhi, and null in Nellore and Guzerat animals (Lacorte et al., 2006). According to our recent results, alanine-encoding allele in DGAT1 K232A polymorphism is found at frequencies of approximately 0.1% in Nellore (55 individuals) and 2% in Guzerat (352 animals, unpublished data). The 232A allele was not observed in a study including a small number of animals of 6 Bos indicus breeds (Sahiwal, Rathi, Deoni, Tharparkar, Red Kandhari, and Punganur; Tantia et al., 2006). Therefore, a large number of individuals should be tested to identify those bearing the advantageous allele for inclusion in marker-assisted selection programs. The costs of this procedure would be significantly reduced by the tetra-primer ARMS-PCR protocol reported here. Additionally, because of its high sensitivity and specificity, ease, and cost-effectiveness, this new method is appropriate for large-scale DGAT1 K232A genotyping.
 |
ACKNOWLEDGEMENTS
|
|---|
The authors thank Vania Maldini Penna and the Colégio Brasileiro de Criadores de Guzerá for the blood samples. This research was financially supported by Fundação de Amparo à Pesquisa de Minas Gerais (FAPEMIG) and Conselho Nacional de Pesquisa (CNPq). Fellowships were granted to R. S. Steinberg and L. Pereira (FAPEMIG). G. Lacorte was supported by a fellowship from the Coordenadoria de Aperfeiçoamento de Pessoal de Nível Siperior (CAPES).
Received for publication April 11, 2008.
Accepted for publication October 1, 2008.
 |
REFERENCES
|
|---|
Cases, S., S. Smith, Y. Zheng, H. Myers, S. Lear, E. Sande, S. Novak, C. Collins, C. Welch, A. Lusis, S. Erickson, and R. J. Farese. 1998. Identification of a gene encoding an acyl CoA:diacylglycerol acyltransferase, a key enzyme in triacylglycerol synthesis. Proc. Natl. Acad. Sci. USA 95:13018–13023.[Abstract/Free Full Text]
Grisart, B., W. Coppieters, F. Farnir, L. Karim, C. Ford, P. Berzi, N. Cambisano, M. Mni, S. Reid, P. Simon, R. Spelman, M. Georges, and R. Snell. 2002. Positional candidate cloning of a QTL 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.[Abstract/Free Full Text]
Hulshof, K., M. van Erp-Baart, M. Anttolainen, W. Becker, S. Church, C. Couet, E. Hermann-Kunz, H. Kesteloot, T. Leth, I. Martins, O. Moreiras, J. Moschandreas, L. Pizzoferrato, A. Rimestad, H. Thorgeirsdottir, J. van Amelsvoort, A. Aro, A. Kafatos, D. Lanzmann-Petithory, and G. van Poppel. 1999. Intake of fatty acids in western Europe with emphasis on trans fatty acids: The TRANSFAIR Study. Eur. J. Clin. Nutr. 53:143–157.[CrossRef][Medline]
Kühn, C., G. Thaller, A. Winter, O. Bininda-Emonds, B. Kaupe, G. Erhardt, J. Bennewitz, M. Schwerin, and R. Fries. 2004. Evidence for multiple alleles at the DGAT1 locus better explains a quantitative trait locus with major effect on milk fat content in cattle. Genetics 167:1873–1881.[Abstract/Free Full Text]
Lacorte, G., M. Machado, M. Martinez, A. Campos, R. Maciel, R. Verneque, R. Teodoro, M. Peixoto, M. Carvalho, and C. Fonseca. 2006. DGAT1 K232A polymorphism in Brazilian cattle breeds. Genet. Mol. Res. 5:475–482.[Medline]
Newton, C., A. Graham, L. Heptinstall, S. Powell, C. Summers, N. Kalsheker, J. Smith, and A. Markham. 1989. Analysis of any point mutation in DNA. The amplification refractory mutation system (ARMS). Nucleic Acids Res. 17:2503–2516.[Abstract/Free Full Text]
Riquet, J., W. Coppieters, N. Cambisano, J. Arranz, P. Berzi, S. Davis, B. Grisart, F. Farnir, L. Karim, M. Mni, P. Simon, J. Taylor, P. Vanmanshoven, D. Wagenaar, J. Womack, and M. Georges. 1999. Fine-mapping of quantitative trait loci by identity by descent in outbred populations: Application to milk production in dairy cattle. Proc. Natl. Acad. Sci. USA 96:9252–9257.[Abstract/Free Full Text]
Sambrook, J., and D. W. Russell. 2001. Molecular Cloning. A Laboratory Manual. Third edition, ed. Cold Spring Harbor Laboratory Press, New York, NY.
Schennink, A., J. Heck, H. Bovenhuis, M. Visker, H. van Valenberg, and J. van Arendonk. 2008. Milk fatty acid unsaturation: Genetic parameters and effects of stearoyl-CoA desaturase (SCD1) and acyl CoA: diacylglycerol acyltransferase 1 (DGAT1). J. Dairy Sci. 91:2135–2143.[Abstract/Free Full Text]
Schennink, A., W. Stoop, M. Visker, J. Heck, H. Bovenhuis, J. van der Poel, H. van Valenberg, and J. van Arendonk. 2007. DGAT1 underlies large genetic variation in milk-fat composition of dairy cows. Anim. Genet. 38:467–473.[CrossRef][Medline]
Smith, S., S. Cases, D. Jensen, H. Chen, E. Sande, B. Tow, D. Sanan, J. Raber, R. Eckel, and R. J. Farese. 2000. Obesity resistance and multiple mechanisms of triglyceride synthesis in mice lacking Dgat. Nat. Genet. 25:87–90.[CrossRef][Medline]
Tantia, M., R. Vijh, B. Mishra, B. Mishra, S. Kumar, and M. Sodhi. 2006. DGAT1 and ABCG2 polymorphism in Indian cattle (Bos indicus) and buffalo (Bubalus bubalis) breeds. BMC Vet. Res. 2:32.[CrossRef][Medline]
Thaller, G., C. Kühn, A. Winter, G. Ewald, O. Bellmann, J. Wegner, H. Zühlke, and R. Fries. 2003. DGAT1, a new positional and functional candidate gene for intramuscular fat deposition in cattle. Anim. Genet. 34:354–357.[CrossRef][Medline]
Winter, A., W. Krämer, F. Werner, S. Kollers, S. Kata, G. Durstewitz, J. Buitkamp, J. Womack, G. Thaller, and R. Fries. 2002. Association of a lysine-232/alanine polymorphism in a bovine gene encoding acyl-CoA:diacylglycerol acyltransferase (DGAT1) with variation at a quantitative trait locus for milk fat content. Proc. Natl. Acad. Sci. USA 99:9300–9305.[Abstract/Free Full Text]
Ye, S., S. Dhillon, X. Ke, A. Collins, and I. Day. 2001. An efficient procedure for genotyping single nucleotide polymorphisms. Nucleic Acids Res. 29:E88.[CrossRef][Medline]