BackgroundIn the Western world bovine milk products are an important protein source in human diet. The major proteins in bovine milk are the four caseins (CN), αS1-, αS2-, β-, and k-CN and the two whey proteins, β-LG and α-LA. It has been shown that both the amount of specific CN and their isoforms including post-translational modifications (PTM) influence technological properties of milk. Therefore, the aim of this study was to 1) estimate genetic parameters for individual proteins in Danish Holstein (DH) (n = 371) and Danish Jersey (DJ) (n = 321) milk, and 2) detect genomic regions associated with specific milk protein and their different PTM forms using a genome-wide association study (GWAS) approach.ResultsFor DH, high heritability estimates were found for protein percentage (0.47), casein percentage (0.43), k-CN (0.77), β-LG (0.58), and α-LA (0.40). For DJ, high heritability estimates were found for protein percentage (0.70), casein percentage (0.52), and α-LA (0.44). The heritability for G-k-CN, U-k-CN and GD was higher in the DH compared to the DJ, whereas the heritability for the PD of αS1-CN was lower in DH compared to DJ, whereas the PD for αS2-CN was higher in DH compared to DJ. The GWAS results for the main milk proteins were in line what has been earlier published. However, we showed that there were SNPs specifically regulating G-k-CN in DH. Some of these SNPs were assigned to casein protein kinase genes (CSNK1G3, PRKCQ).ConclusionThe genetic analysis of the major milk proteins and their PTM forms revealed that these were heritable in both DH and DJ. In DH, genomic regions specific for glycosylation of k-CN were detected. Furthermore, genomic regions for the major milk proteins confirmed the regions on BTA6 (casein cluster), BTA11 (PEAP), and BTA14 (DGAT1) as important regions influencing protein composition in milk. The results from this study provide confidence that it is possible to breed for specific milk protein including the different PTM forms.Electronic supplementary materialThe online version of this article (doi:10.1186/s12863-016-0421-2) contains supplementary material, which is available to authorized users.

Genetic parameters were estimated for the major milk proteins using bivariate and multi-trait models based on genomic relationships between animals. The analyses included, apart from total protein percentage, αS1-casein (CN), αS2-CN, β-CN, κ-CN, α-lactalbumin, and β-lactoglobulin, as well as the posttranslational sub-forms of glycosylated κ-CN and αS1-CN-8P (phosphorylated). Standard errors of the estimates were used to compare the models. In total, 650 Danish Holstein cows across 4 parities and days in milk ranging from 9 to 481d were selected from 21 herds. The multi-trait model generally resulted in lower standard errors of heritability estimates, suggesting that genetic parameters can be estimated with high accuracy using multi-trait analyses with genomic relationships for scarcely recorded traits. The heritability estimates from the multi-trait model ranged from low (0.05 for β-CN) to high (0.78 for κ-CN). Genetic correlations between the milk proteins and the total milk protein percentage were generally low, suggesting the possibility to alter protein composition through selective breeding with little effect on total milk protein percentage.

BackgroundAccurate genomic prediction requires a large reference population, which is problematic for traits that are expensive to measure. Traits related to milk protein composition are not routinely recorded due to costly procedures and are considered to be controlled by a few quantitative trait loci of large effect. The amount of variation explained may vary between regions leading to heterogeneous (co)variance patterns across the genome. Genomic prediction models that can efficiently take such heterogeneity of (co)variances into account can result in improved prediction reliability. In this study, we developed and implemented novel univariate and bivariate Bayesian prediction models, based on estimates of heterogeneous (co)variances for genome segments (BayesAS). Available data consisted of milk protein composition traits measured on cows and de-regressed proofs of total protein yield derived for bulls. Single-nucleotide polymorphisms (SNPs), from 50K SNP arrays, were grouped into non-overlapping genome segments. A segment was defined as one SNP, or a group of 50, 100, or 200 adjacent SNPs, or one chromosome, or the whole genome. Traditional univariate and bivariate genomic best linear unbiased prediction (GBLUP) models were also run for comparison. Reliabilities were calculated through a resampling strategy and using deterministic formula.ResultsBayesAS models improved prediction reliability for most of the traits compared to GBLUP models and this gain depended on segment size and genetic architecture of the traits. The gain in prediction reliability was especially marked for the protein composition traits β-CN, κ-CN and β-LG, for which prediction reliabilities were improved by 49 percentage points on average using the MT-BayesAS model with a 100-SNP segment size compared to the bivariate GBLUP. Prediction reliabilities were highest with the BayesAS model that uses a 100-SNP segment size. The bivariate versions of our BayesAS models resulted in extra gains of up to 6% in prediction reliability compared to the univariate versions.ConclusionsSubstantial improvement in prediction reliability was possible for most of the traits related to milk protein composition using our novel BayesAS models. Grouping adjacent SNPs into segments provided enhanced information to estimate parameters and allowing the segments to have different (co)variances helped disentangle heterogeneous (co)variances across the genome.