|
|
||||||||
Department of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, Box 7023, SE-750 07 Uppsala, Sweden
1 Corresponding author: Jessica.naslund{at}hgen.slu.se
| ABSTRACT |
|---|
|
|
|---|
Key Words: DGAT1 genetic polymorphism cattle milk fat
| INTRODUCTION |
|---|
|
|
|---|
Milk fat contains approximately 98% triglycerides. Characteristic of milk fat from ruminants is the occurrence of short-chain fatty acids comprising 4 to 6 carbon atoms. When triglycerides are synthesized, fatty acids are first attached to positions sn-1 and sn-2 (according to the stereospecific numbering) on the glycerol molecule under the catalyzing effects of enzymes. When a fatty acid is attached to the third position, the enzyme acyl-CoA:diacylglycerol acyltransferase 1 (DGAT1) acts as a catalyst. It is interesting to note that it is only in this position that we find the short-chain fatty acids with 4 to 6 carbon atoms (Marshall and Knudsen, 1977).
A bovine QTL affecting milk yield and composition was reported to be located on bovine chromosome 14 near the centromeric end (Coppieters et al., 1998; Heyen et al., 1999; Looft et al., 2001). Through fine mapping of the QTL region (Riquet et al., 1999; Farnir et al., 2002), the DGAT1 gene emerged as a strong candidate for the QTL effect. The gene encoding the DGAT enzyme had previously been identified in studies of mice (Cases et al., 1998). The A-A to G-C dinucleotide substitution in exon VIII of DGAT1, which gives rise to an amino acid substitution from lysine (K) to alanine (A), was postulated to be the causative mutation (Grisart et al., 2002; Winter et al., 2002). The allele frequency and the effect of the K232A polymorphism have been characterized in dairy cattle populations in New Zealand (Grisart et al., 2002; Spelman et al., 2002), Israel (Weller et al., 2003), the Netherlands (Grisart et al., 2002), Germany (Thaller et al., 2003; Sanders et al., 2006), Poland (Pareek et al., 2005; Strzalkowska et al., 2005), and France (Gautier et al., 2007). The effect of the amino acid substitution is an increase in milk and protein yield, and a decrease in fat yield and fat and protein percentage. In a gene expression study it was shown that the lysine variant had a greater enzyme activity level (Vmax) compared with the alanine variant (Grisart et al., 2004).
The existence of multiple alleles at the DGAT1 locus or a second linked QTL was suggested by Bennewitz et al. (2004), and soon after a variable number of tandem repeats (VNTR) polymorphism was identified in the promoter region of the DGAT1 gene that was associated with variation in milk fat percentage (Kühn et al., 2004). The effect was not as marked as the DGAT1 K232A effect, and in a study in French dairy cattle, the VNTR explained only a small fraction of the QTL variance (Gautier et al., 2007). Moreover, a polymorphism in the gene coding for the enzyme CYP11B1 showed association with milk production and functional traits in German Holstein cattle (Kaupe et al., 2007). However, when accounting for the effects of the DGAT1 K232A and CYP11B1 V30A polymorphisms simultaneously in the statistical analysis, the latter contributed only marginally to the genetic variation in the analyzed milk production traits, compared with the DGAT1 K232A polymorphism alone.
The dairy herd at the Swedish University of Agricultural Sciences includes cows of the SRB and the SLB breeds. The SRB cows belong to 1 of 2 selection lines producing milk with a high or low percentage of milk fat at equal levels of milk energy content, such that cows from both lines have similar nutritional requirements. This material offers an interesting opportunity to explore the biology and genetics behind variation in milk composition and especially fat content. The SLB breed of today is a Swedish counterpart to the various Holstein populations in Europe and North America, and as such, functions as a reference when comparing results between studies. The aims of the present study were to investigate how different breeding strategies affect the allele frequencies at the K232A polymorphism and to estimate the effect of the DGAT1 K232A in the 2 major Swedish dairy breeds.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Phenotypic Data
Data comprised 16,866 test-day records from 462 lactations (1 to 3) of the 239 cows. Milk samples were collected weekly from April 1990 through December 2003. Milk samples were treated with bronopol (Boots Microcheck, Nottingham, UK) immediately after milking, and samples of fresh milk were analyzed for SCC by flow cytometry (Fossomatic 5200, A/S Foss Electric, Hillerød, Denmark). Infrared technique was used to determine the concentrations of fat, protein, and lactose (MilkoScan 93; A/S Foss Electric). Milk data included milk yield, milk yield expressed as kilograms of ECM, concentration and yield of fat and protein, and SCC transformed to the natural logarithm scale (lnSCC). Observations with an SCC >300,000 cells/mL were not included in the statistical analysis to avoid milk with characteristics of mastitis. Lactation number, lactation week, and season and year of milk sampling were registered for each individual milk sample. Only registrations from the first 305 d in lactation were used. The mean values and standard deviations for the analyzed milk traits for each parity are in Table 1
.
|
Polymerase chain reaction amplification was carried out on a PTC-200 DNA Engine (MJ Research, Waltham, MA) in a 25-µL reaction volume using the AmpliTaq Gold (Applied Biosystems, Foster City, CA) PCR kit. The reactions were performed with the addition of betaine (Sigma, St. Louis, MO) to improve the PCR reaction. To distinguish between the 2 variants a single PCR was performed in which part of exon 8 was amplified.
The samples were analyzed by using the pyrosequencing method (Ronaghi et al., 1998). The biotinylated PCR product (20 µL) was immobilized onto streptavidin-coated paramagnetic beads (Dynal AS, Oslo, Norway) using binding buffer (5 mM Tris-HCl, 1 M NaCl, 0.5 mM EDTA, 0.05% Tween 20, pH 7.6) in a total volume of 80 µL during a 10-min incubation on a vortex mixer (Vortex Genie 2, Scientific Industries, Bohemia, NY). Biotinylated single-stranded DNA was obtained by washing the immobilized PCR product in 0.2 M NaOH and washing the beads in washing buffer (10 mM Tris-acetate, pH 7.6). A total of 15 pmol of detection primer, designed with its 3' end immediately upstream of the polymorphic site 5'-GCT CGT AGC TTT GGC AGG TA-3', was allowed to hybridize onto single-stranded DNA in 40 µL of annealing buffer (20 mM Tris-acetate, 2 mM MgAc2, pH 7.6) at 80°C for 2 min with subsequent cooling to room temperature. Pyrosequencing was carried out on the PSQ96 Pyrosequencer instrument using the PSQ96 SNP Reagent kit (Biotage AB, Uppsala, Sweden) containing dATP
S, dCTP, dGTP, dTTP, enzyme mixture (DNA polymerase, ATP sulfurylase, luciferase, and apyrase), and substrate mixture (adenosine 5'-phospho-sulfate and luciferin).
Statistical Analysis
Data were analyzed with the following model:
![]() | [1] |
where yijklmno = test-day record; µ = the overall mean; ysk = the fixed effect of year-season of testing (k = 1, 2...50); groupj = the fixed effect of group, (j = SLB, SRB/HF, SRB/HF); parityl = the fixed effect of parity number, (l = 1, 2, 3); wim = week in milk; bm1 ... bm4 = the regression coefficients associated with the fixed lactation function (m = 1, 2+); Pjn = random environmental effect of cow n; Ajn = random additive genetic effect of polygenic background for cow n, and eijklmno = the random residual effect.
Separate lactation curves (Ali and Schaeffer, 1987) were fitted for all parities. The shape of the lactation curve was the same for lactations 2 and 3; that is, only the level (intercept) of the lactation curve was different for lactations 2 and 3. Using a parametric curve has several advantages compared with the modeling of subclasses; there are fewer parameters to be estimated, one overcomes the arbitrary definition of the classes, and the model takes into account that the residual variance changes continuously over time (Jaffrezic et al., 2000).
The model included a relationship matrix whereby the additive polygenetic background is accounted for to avoid bias in the estimates of DGAT1 K232A gene effects. Parity and week in milk (wim) were used as repeated measures variables, where parity identifies the multivariate factor and wim the levels of time according to the unstructured and first-order autoregressive structure, denoted UN@AR(1), for multivariate repeated measures. The parameters for the random effects were estimated from the data.
Three different alternatives for f(gtp) were considered in model 1:
![]() |
where x1 = the number of copies (0, 1, or 2) of the K variant of the DGAT1 gene for the cow n; a = the additive effect; x2 = indicator variable equal to 1 for the AK genotype and 0 otherwise; and d = the dominance effect;
![]() |
where genotypei = the fixed effect of DGAT1 K232A genotype (i = AA, AK, or KK); and
![]() |
where
= the average allele substitution effect (Falconer and Mackay, 1996).
Data were analyzed with the MIXED procedure of the SAS system (SAS Institute, 2006), and traits were analyzed separately. The significance of fixed effects was tested using the F-value computed by the MIXED procedure and the containment method for computing the denominator degrees of freedom.
Differences in allele frequencies between the groups were tested using Fishers exact test.
| RESULTS |
|---|
|
|
|---|
|
|
Effect of the DGAT1 K232A Polymorphism
The additive effect of DGAT1 was highly significant for fat concentration in both the SLB and SRB/HF groups. According to the genotypic values in Table 4
, the increase in fat percentage was 0.52 and 0.51, respectively, for each copy of the K variant. The genotypic value of the KK genotype was also positive for fat yield and protein concentration in the SLB group, whereas it was not significantly associated with protein yield in any of the groups. There was no significant dominance effect observed for any of the traits. With regard to yield and concentration of lactose and SCC, there were no effects of the DGAT1 polymorphism (data not shown). Imprecise estimates of the additive and dominance effects of DGAT1 were obtained within the SRB/LF line because of the lack of SRB/LF cows with the KK genotype and the low number of AK heterozygotes.
|
The allele substitution effect of changing an A allele for a K allele was in accordance with the estimated allele effects based on genotypic values. However, the standard errors of the estimates were lower, resulting in significant allele substitution effects for fat yield and protein yield in the SRB/HF group of cows.
| DISCUSSION |
|---|
|
|
|---|
In Sweden the milk price to the farmer has, for decades, favored milk volume rather than milk components, giving the DGAT1 A variant a selective advantage. This might be a reason for the high frequency of this variant in the present material, in which the frequency of the A variant in the SLB breed (0.86) was considerably higher than in other Holstein populations (Grisart et al., 2002; Spelman et al., 2002; Thaller et al., 2003; Weller et al., 2003; Kaupe et al., 2004) in which selection targeted fat and protein concentrations. Although the A allele also is clearly predominant among the SRB cows in this study, the frequency of the K variant in the 2 selection lines was somewhat higher than observed in most of the previous studies on Red dairy breeds (Winter et al., 2002; Thaller et al., 2003; Kaupe et al., 2004). This could be a result of the reinforced selection on composition traits. In New Zealand, both Ayrshire and Holstein-Friesian cows carry a relatively high frequency of the K variant (0.60 and 0.22, respectively), possibly because of the strong selection for high DM content of milk (Spelman et al., 2002). Because the SRB breed has a reputation of producing milk with relatively high fat and protein percentages and the SLB breed is known to have a high milk yield, the 2 lines of SRB cows were expected to have a greater frequency of the K variant than the SLB cows. However, there was no difference in allele frequencies between the breeds in this limited material.
Effect of the DGAT1 K232A Polymorphism
In this study we used a unique data set consisting of a large set of phenotypic records of genotyped daughters. In contrast to using granddaughter designs, this approach enabled estimation of additive and dominance effects. The accuracy of the estimates regarding the effect of the K allele suffered from the small number of KK homozygous individuals in the material. Apart from this, the estimated effects of the DGAT1 K232A polymorphism in this study were similar to previously published results (Grisart et al., 2002; Spelman et al., 2002; Thaller et al., 2003; Sanders et al., 2006). The allele substitution effects for fat and protein concentration were, however, larger in this study than those reported by Thaller et al. (2003) but were similar to the allele substitution effects in the daughter design by Grisart et al. (2002). The lack of dominance effects were in accordance with previous results (Grisart et al., 2002).
The DGAT1 K232A effect on fat concentration has been reported to differ in magnitude between breeds (Spelman et al., 2002; Thaller et al., 2003; Gautier et al., 2007). However, the estimated allele substitution effect was not significantly different for the SRB/HF compared with the SLB cows, and the genotypic values for fat concentration in this material were similar for the SRB/HF and SLB cows. This contradicts the existence of breed differences regarding the DGAT1 effects in this material. Differences in effects on fat concentration may be because of the interaction with background genes, which might differ between breeds. Although DGAT1 plays a highly significant role in the formation of triglycerides, other genes can be of importance in the various pathways involved in the milk fat synthesis. However, accounting for a random polygenic component in the model and using a known pedigree relationship matrix is expected to reduce the effect of the background genome. Because this study is part of a larger project aiming to identify genetic factors affecting off-flavors in milk we chose not to genotype for the VNTR reported by Kühn et al. (2004), because its effect would be expected to be similar but less pronounced compared with the effect of the DGAT1 K232A polymorphism. However, if the effect of the DGAT1 K232A polymorphism on occurrence of off-flavors is found to be marginal, genotyping for additional polymorphisms in the chromosomal region around the DGAT1 gene may be performed.
By creating divergent lines for milk fat concentration, the 2 groups of cows are more likely to segregate for genes with major effects on the trait, which was indeed shown to be the case with the DGAT1 polymorphism in which the applied selection was reflected in the frequency of DGAT1 K232A variants. A previous study with a limited number of cows from the 2 selection lines indicated that the applied selection has altered not only the fat concentration of the milk, but also fat composition (Agenas et al., 2003). Knowing the function of the DGAT1 enzyme and the special features of bovine milk fatty acid composition, it could be expected that the DGAT1 K232A polymorphism would also have an effect on fat composition. Our initial results from fatty acid analysis of a subsample of the milk samples included in the present study indicate that the DGAT1 genotype AA is associated with an increased proportion of C18:2 (Näslund et al., 2006). Because milk fat composition is known to influence the quality of dairy products, further studies on the effect of DGAT1 are warranted.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication May 2, 2007. Accepted for publication January 7, 2008.
| REFERENCES |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |