The objective of this study was to estimate genetic parameters for major milk fatty acids and milk production traits. One morning milk sample was collected from 1,918 Holstein-Friesian heifers located in 398 commercial herds in The Netherlands. Each sample was analyzed for total percentages of fat and protein, and for detailed fatty acid percentages (computed as fatty acid weight as a proportion of total fat weight). Intraherd heritabilities were high for C4:0 to C16:0, ranging from 0.42 for C4:0 to 0.71 for C10:0. Saturated and unsaturated C18 fatty acids had intraherd heritability estimates of approximately 0.25, except for C18:2 cis-9, trans-11, which was 0.42. Standard errors of the heritabilities were between 0.07 and 0.12. Genetic correlations were high and positive among C4:0 to C14:0, as well as among unsaturated C18, but correlations of C4:0 to C14:0 with unsaturated C18 were generally weak. The genetic correlation of C16:0 with fat percentage was positive (0.65), implying that selection for fat percentage should result in a correlated increase of C16:0, whereas unsaturated C18 fatty acids decreased with increasing fat percentage (-0.74). Milk fat composition can be changed by means of selective breeding, which offers opportunities to meet consumer demands regarding health and technological aspects.

With regard to human health aspects of milk fat, increasing the amount of unsaturated fatty acids in milk is an important selection objective. The cow's diet has an influence on the degree of unsaturation, but literature suggests that genetics also plays a role. To estimate genetic variation in milk fatty acid unsaturation indices, milk fatty acid composition of 1,933 Dutch Holstein Friesian heifers was measured and unsaturation indices were calculated. An unsaturation index represents the concentration of the unsaturated product proportional to the sum of the unsaturated product and the saturated substrate. Intraherd heritabilities were moderate, ranging from 0.23 +/- 0.07 for conjugated linoleic acid (CLA) index to 0.46 +/- 0.09 for C16 index. We genotyped the cows for the SCD1 A293V and DGAT1 K232A polymorphisms, which are known to alter milk fatty acid composition. Both genes explain part of the genetic variation in unsaturation indices. The SCD1 V allele is associated with lower C10, C12, and C14 indices, and with higher C16, C18, and CLA indices in comparison to the SCD1 A allele, with no differences in total unsaturation index. In comparison to the DGAT1 K allele, the DGAT1 A allele is associated with lower C10, C12, C14, and C16 indices and with higher C18, CLA, and total indices. We conclude that selective breeding can contribute to higher unsaturation indices, and that selective breeding can capitalize on genotypic information of both the SCD1 A293V and the DGAT1 K232A polymorphism.

Dietary fat may play a role in the aetiology of many chronic diseases. Milk and milk-derived foods contribute substantially to dietary fat, but have a fat composition that is not optimal for human health. We measured the fat composition of milk samples in 1918 Dutch Holstein Friesian cows in their first lactation and estimated genetic parameters for fatty acids. Substantial genetic variation in milk-fat composition was found: heritabilities were high for short- and medium-chain fatty acids (C4:0-C16:0) and moderate for long-chain fatty acids (saturated and unsaturated C18). We genotyped 1762 cows for the DGAT1 K232A polymorphism, which is known to affect milk-fat percentage, to study the effect of the polymorphism on milk-fat composition. We found that the DGAT1 K232A polymorphism has a clear influence on milk-fat composition. The DGAT1 allele that encodes lysine (K) at position 232 (232K) is associated with more saturated fat; a larger fraction of C16:0; and smaller fractions of C14:0, unsaturated C18 and conjugated linoleic acid (P < 0.001). We conclude that selective breeding can make a significant contribution to change the fat composition of cow's milk.

The effects of beta-lactoglobulin (beta-LG), beta-casein (beta-CN), and kappa-CN variants and beta-kappa-CN haplotypes on the relative concentrations of the major milk proteins alpha-lactalbumin (alpha-LA), beta-LG, alpha(S1)-CN, alpha(S2)-CN, beta-CN, and kappa-CN and milk production traits were estimated in the milk of 1,912 Dutch Holstein-Friesian cows. We show that in the Dutch Holstein-Friesian population, the allele frequencies have changed in the past 16 years. In addition, genetic variants and casein haplotypes have a major impact on the protein composition of milk and explain a considerable part of the genetic variation in milk protein composition. The beta-LG genotype was associated with the relative concentrations of beta-LG (A >> B) and of alpha-LA, alpha(S1)-CN, alpha(S2)-CN, beta-CN, and kappa-CN (B > A) but not with any milk production trait. The beta-CN genotype was associated with the relative concentrations of beta-CN and alpha(S2)-CN (A(2) > A(1)) and of alpha(S1)-CN and kappa-CN (A(1) > A(2)) and with protein yield (A(2) > A(1)). The kappa-CN genotype was associated with the relative concentrations of kappa-CN (B > E > A), alpha(S2)-CN (B > A), alpha-LA, and alpha(S1)-CN (A > B) and with protein percentage (B > A). Comparing the effects of casein haplotypes with the effects of single casein variants can provide better insight into what really underlies the effect of a variant on protein composition. We conclude that selection for both the beta-LG genotype B and the beta-kappa-CN haplotype A(2)B will result in cows that produce milk that is more suitable for cheese production.

The effects of lactation stage, negative energy balance (NEB), and milk fat depression (MFD) were estimated on detailed milk fat composition in primiparous Holstein-Friesian cows. One morning milk sample was collected from each of 1,933 cows from 398 commercial Dutch herds in winter 2005. Milk fat composition was measured using gas chromatography, and fat and protein percentage were measured using infrared spectrometry. Each fatty acid changed 0.5 to 1 phenotypic standard deviation over lactation, except odd-chain C5:0 to C15:0, branched-chain fatty acids, and trans-10, cis-12 conjugated linoleic acid (CLA). The greatest change was an increase from 31.2 to 33.3% (wt/wt) for C16:0 from d 80 to 150 of lactation. Energy status was estimated for each cow as the deviation from each average lactation fat-to-protein ratio (FPdev). A high FPdev (>0.12) indicated NEB. Negative energy balance was associated with an increase in C16:0 (0.696 +/- 0.178) and C18:0 (0.467 +/- 0.093), which suggested mobilization of body fat reserves. Furthermore, NEB was associated with a decrease in odd-chain C5:0 to C15:0 (-0.084 +/- 0.020), which might reflect a reduced allocation of C3 components to milk fat synthesis. A low FPdev indicated MFD (<-0.12) and was associated with a decrease in C16:0 (-0.681 +/- 0.255) and C18:0 (-0.128 +/- 0.135) and an increase in total unsaturated fatty acids (0.523 +/- 0.227). The study showed that both lactation stage and energy balance significantly contribute to variation in milk fat composition and alter the activity of different fatty acid pathways.

BackgroundIdentifying genomic regions, and preferably individual genes, responsible for genetic variation in milk fat composition of bovine milk will enhance the understanding of biological pathways involved in fatty acid synthesis and may point to opportunities for changing milk fat composition via selective breeding. An association study of 50,000 single nucleotide polymorphisms (SNPs) was performed for even-chain saturated fatty acids (C4:0-C18:0), even-chain monounsaturated fatty acids (C10:1-C18:1), and the polyunsaturated C18:2cis9,trans11 (CLA) to identify genomic regions associated with individual fatty acids in bovine milk.ResultsThe two-step single SNP association analysis found a total of 54 regions on 29 chromosomes that were significantly associated with one or more fatty acids. Bos taurus autosomes (BTA) 14, 19, and 26 showed highly significant associations with seven to ten traits, explaining a relatively large percentage of the total additive genetic variation. Many additional regions were significantly associated with the fatty acids. Some of the regions harbor genes that are known to be involved in fat synthesis or were previously identified as underlying quantitative trait loci for fat yield or content, such as ABCG2 and PPARGC1A on BTA 6; ACSS2 on BTA 13; DGAT1 on BTA 14; ACLY, SREBF1, STAT5A, GH, and FASN on BTA 19; SCD1 on BTA26; and AGPAT6 on BTA 27.ConclusionsMedium chain and unsaturated fatty acids are strongly influenced by polymorphisms in DGAT1 and SCD1. Other regions also showed significant associations with the fatty acids studied. These additional regions explain a relatively small percentage of the total additive genetic variance, but they are relevant to the total genetic merit of an individual and in unraveling the genetic background of milk fat composition. Regions identified in this study can be fine mapped to find causal mutations. The results also create opportunities for changing milk fat composition through breeding by selecting individuals based on their genetic merit for milk fat composition.

The aim of this study was to estimate genetic parameters for test-day milk urea nitrogen (MUN) and its relationships with milk production traits. Three test-day morning milk samples were collected from 1,953 Holstein-Friesian heifers located on 398 commercial herds in The Netherlands. Each sample was analyzed for somatic cell count, net energy concentration, MUN, and the percentage of fat, protein, and lactose. Genetic parameters were estimated using an animal model with covariates for days in milk and age at first calving, fixed effects for season of calving and effect of test or proven bull, and random effects for herd-test day, animal, permanent environment, and error. Coefficient of variation for MUN was 33%. Estimated heritability for MUN was 0.14. Phenotypic correlation of MUN with each of the milk production traits was low. The genetic correlation was close to zero for MUN and lactose percentage (-0.09); was moderately positive for MUN and net energy concentration of milk (0.19), fat yield (0.41), protein yield (0.38), lactose yield (0.22), and milk yield (0.24), and percentage of fat (0.18), and percentage of protein (0.27); and was high for MUN and somatic cell score (0.85). Herd-test day explained 58% of the variation in MUN, which suggests that management adjustments at herd-level can reduce MUN. This study shows that it is possible to influence MUN by herd practice and by genetic selection.

The objective of this study was to estimate genetic parameters for major milk proteins. One morning milk sample was collected from 1,940 first-parity Holstein-Friesian cows in February or March 2005. Each sample was analyzed with capillary zone electrophoresis to determine the relative concentrations of the 6 major milk proteins. The results show that there is considerable genetic variation in milk protein composition. The intraherd heritabilities for the relative protein concentrations were high and ranged from 0.25 for beta-casein to 0.80 for beta-lactoglobulin. The intraherd heritability for the summed whey fractions (0.71) was higher than that for the summed casein fractions (0.41). Further, there was relatively more variation in the summed whey fraction (coefficient of variation was 11% and standard deviation was 1.23) compared with the summed casein fraction (coefficient of variation was 2% and standard deviation was 1.72). For the caseins and alpha-lactalbumin, the proportion of phenotypic variation explained by herd was approximately 14%. For beta-lactoglobulin, the proportion of phenotypic variation explained by herd was considerably lower (5%). Eighty percent of the genetic correlations among the relative contributions of the major milk proteins were between -0.38 and +0.45. The genetic correlations suggest that it is possible to change the relative proportion of caseins in milk. Strong negative genetic correlations were found for beta-lactoglobulin with the summed casein fractions (-0.76), and for beta-lactoglobulin with casein index (-0.98). This study suggests that there are opportunities to change the milk protein composition in the cow's milk using selective breeding.