Overlap in genomic variation associated with milk fat composition in Holstein Friesian and Dutch native dual-purpose breeds.

The aim of this study was to identify if genomic variations associated with fatty acid (FA) composition are similar between the Holstein-Friesian (HF) and native dual-purpose breeds used in the Dutch dairy industry. Phenotypic and genotypic information were available for the breeds Meuse-Rhine-Yssel (MRY), Dutch Friesian (DF), Groningen White Headed (GWH), and HF. First, the reliability of genomic breeding values of the native Dutch dual-purpose cattle breeds MRY, DF, and GWH was evaluated using single nucleotide polymorphism (SNP) effects estimated in HF, including all SNP or subsets with stronger associations in HF. Second, the genomic variation of the regions associated with FA composition in HF (regions on Bos taurus autosome 5, 14, and 26), were studied in the different breeds. Finally, similarities in genotype and allele frequencies between MRY, DF, GWH, and HF breeds were assessed for specific regions associated with FA composition. On average across the traits, the highest reliabilities of genomic prediction were estimated for GWH (0.158) and DF (0.116) when the 8 to 22 SNP with the strongest association in HF were included. With the same set of SNP, GEBV for MRY were the least reliable (0.022). This indicates that on average only 2 (MRY) to 16% (GWH) of the genomic variation in HF is shared with the native Dutch dual-purpose breeds. The comparison of predicted variances of different regions associated with milk and milk fat composition showed that breeds clearly differed in genomic variation within these regions. Finally, the correlations of allele frequencies between breeds across the 8 to 22 SNP with the strongest association in HF were around 0.8 between the Dutch native dual-purpose breeds, whereas the correlations between the native breeds and HF were clearly lower and around 0.5. There was no consistent relationship between the reliabilities of genomic prediction for a specific breed and the correlation between the allele frequencies of this breed and HF. In conclusion, most of the genomic variation associated with FA composition in the Dutch dual-purpose breeds appears to be breed-specific. Furthermore, the minor allele frequencies of genes having an effect on the milk FA composition in HF were shown to be much smaller in the breeds MRY, DF, and GWH, especially for the MRY breed.

Selection of SNP from 50K and 777K arrays to predict breed of origin in cattle.

Reliable breed assignment can be performed with SNP. Currently, high density SNP chips are available with large numbers of SNP from which the most informative SNP can be selected for breed assignment. Several methods have been published to select the most informative SNP to distinguish among breeds. In this study, we evaluated Delta, Wright's FST, and Weir and Cockerham's FST, and extended these methods by adding a rule to avoid selection of sets of SNP in high linkage disequilibrium (LD) providing the same information. The SNP that had a r2 value>0.3 with any of the SNP already selected were discarded. The different selection methods were evaluated for both the 50K SNP and 777K Bovine BeadChip. Animals from 4 cattle breeds (989 Holstein Friesian, 97 Groningen White headed, 137 Meuse-Rhine-Yssel, and 64 Dutch Friesian) were genotyped. After editing 30,447 and 452,525 SNP were available for the 50K and 777K SNP chip, respectively. All selection methods showed that only a small set of SNP is needed to differentiate among the 4 Dutch cattle breeds, whereas comparison of the selection methods showed only small differences. In general, the 777K performed marginally better than the 50K BeadChip, especially at higher confidence thresholds. The rule to avoid selection of SNP in high LD reduced the required number of SNP to achieve correct breed assignment. The Global Weir and Cockerham's FST performed marginally better than other selection methods. There was little overlap in the SNP selected from the 2 BeadChips, whereas the number of SNP selected was about the same.

Differences in milk fat composition predicted by mid-infrared spectrometry among dairy cattle breeds in the Netherlands.

The aim of this study was to estimate breed differences in milk fatty acid (FA) profile among 5 dairy cattle breeds present in the Netherlands: Holstein-Friesian (HF), Meuse-Rhine-Yssel (MRY), Dutch Friesian (DF), Groningen White Headed (GWH), and Jersey (JER). For this purpose, total fat percentage and detailed FA contents in milk (14 individual FA and 14 groups of FA) predicted from mid-infrared spectra were used. Mid-infrared spectrometry profiles were collected during regular milk recording from a range of herds with different combinations of breeds, including both purebred and crossbred cows. The data set used for the analyses contained 41,404 records from a total of 24,445 cows. In total 7,626 cows were crossbreds belonging to the breeds HF, MRY, DF, GWH, and JER; 1,769 purebreds (≥87.5%) belonging to the breeds MRY, DF, GWH, and JER; and the other 15,050 cows were HF. Breed effects were estimated using a single-trait animal model. The content in milk of short-chain FA C4:0, C6:0, C8:0, C10:0, C12:0, C14:0, and C16:0 was higher for JER and the content in milk of C16:0 was lower for GWH compared with the other breeds; when adjusting for breed differences in fat percentage, however, not all breed differences were significant. Breed differences were also found for cis-9 C14:1, cis-9 C16:1, C18:0, and a number of C18 unsaturated FA. In general, differences in fat composition in milk between HF, MRY, and DF were not significant. Jerseys tended to produce more saturated FA, whereas GWH tended to produce relatively less saturated FA. After adjusting for differences in fat percentage, breed differences in detailed fat composition disappeared or became smaller for several short- and medium-chain FA, whereas for several long-chain unsaturated FA, more significant breed differences were found. This indicates that short- and medium-chain FA are for all breeds more related to total fat percentage than long-chain FA. In conclusion, between breed differences were found in detailed FA composition and content of individual FA. Especially, for FA produced through de novo synthesis (short-chain FA, C12:0, C14:0, and partly C16:0) differences were found for JER and GWH, compared with the breeds HF, MRY, and DF.

Validation of fatty acid predictions in milk using mid-infrared spectrometry across cattle breeds.

The aim of this study was to investigate the accuracy to predict detailed fatty acid (FA) composition of bovine milk by mid-infrared spectrometry, for a cattle population that partly differed in terms of country, breed and methodology used to measure actual FA composition compared with the calibration data set. Calibration equations for predicting FA composition using mid-infrared spectrometry were developed in the European project RobustMilk and based on 1236 milk samples from multiple cattle breeds from Ireland, Scotland and the Walloon Region of Belgium. The validation data set contained 190 milk samples from cows in the Netherlands across four breeds: Dutch Friesian, Meuse-Rhine-Yssel, Groningen White Headed (GWH) and Jersey (JER). The FA measurements were performed using gas-liquid partition chromatography (GC) as the gold standard. Some FAs and groups of FAs were not considered because of differences in definition, as the capillary column of the GC was not the same as used to develop the calibration equations. Differences in performance of the calibration equations between breeds were mainly found by evaluating the standard error of validation and the average prediction error. In general, for the GWH breed the smallest differences were found between predicted and reference GC values and least variation in prediction errors, whereas for JER the largest differences were found between predicted and reference GC values and most variation in prediction errors. For the individual FAs 4:0, 6:0, 8:0, 10:0, 12:0, 14:0 and 16:0 and the groups' saturated FAs, short-chain FAs and medium-chain FAs, predictions assessed for all breeds together were highly accurate (validation R 2 > 0.80) with limited bias. For the individual FAs cis-14:1, cis-16:1 and 18:0, the calibration equations were moderately accurate (R 2 in the range of 0.60 to 0.80) and for the individual FA 17:0 predictions were less accurate (R 2 < 0.60) with considerable bias. FA concentrations in the validation data set of our study were generally higher than those in the calibration data. This difference in the range of FA concentrations, mainly due to breed differences in our study, can cause lower accuracy. In conclusion, the RobustMilk calibration equations can be used to predict most FAs in milk from the four breeds in the Netherlands with only a minor loss of accuracy.

Short communication: milk fat composition of 4 cattle breeds in the Netherlands.

Milk fatty acid (FA) composition was compared among 4 cattle breeds in the Netherlands: Dutch Friesian (DF; 47 animals/3 farms), Meuse-Rhine-Yssel (MRY; 52/3), Groningen White Headed (GWH; 45/3), and Jersey (JER; 46/3). Each cow was sampled once between December 2008 and March 2009 during the indoor housing season, and samples were analyzed using gas chromatography. Significant breed differences were found for all traits including fat and protein contents, 13 major individual FA, 9 groups of FA, and 5 indices. The saturated fatty acid proportion, which is supposed to be unfavorable for human health, was smaller for GWH (68.9%) compared with DF (74.1%), MRY (72.3%), and JER (74.3%) breeds. The proportion of conjugated linoleic acid and the unsaturation index, which are associated positively with human health, were both highest for GWH. Differences in milk fat composition can be used in strategies to breed for milk with a FA profile more favorable for human health. Our results support the relevance of safeguarding the local Dutch breeds.