Role of CSN2, CSN3, and BLG genes and the polygenic background in the cattle milk protein profile

To devise better selection strategies in dairy cattle breeding programs, a deeper knowledge of the role of the major genes encoding for milk protein fractions is required. The aim of the present study was to assess the effect of the CSN2, CSN3, and BLG genotypes on individual protein fractions (αS1-CN, αS2-CN, β-CN, κ-CN, β-LG, α-LA) expressed qualitatively as percentages of total nitrogen content (% N), quantitatively as contents in milk (g/L), and as daily production levels (g/d). Individual milk samples were collected from 1,264 Brown Swiss cows reared in 85 commercial herds in Trento Province (northeast Italy). A total of 989 cows were successfully genotyped using the Illumina Bovine SNP50 v.2 BeadChip (Illumina Inc.), and a genomic relationship matrix was constructed using the 37,519 SNP markers obtained. Milk protein fractions were quantified and the β-CN, κ-CN, and β-LG genetic variants were identified by reversed-phase HPLC (RP-HPLC). All protein fractions were analyzed through a Bayesian multitrait animal model implemented via Gibbs sampling. The effects of days in milk, parity order, and the CSN2, CSN3, and BLG genotypes were assigned flat priors in this model, whereas the effects of herd and animal additive genetic were assigned Gaussian prior distributions, and inverse Wishart distributions were assumed for the respective co-variance matrices. Marginal posterior distributions of the parameters of interest were compared before and after the inclusion of the effects of the 3 major genes in the model. The results showed that a high portion of the genetic variance was controlled by the major genes. This was particularly apparent in the qualitative protein profile, which was found to have a higher heritability than the protein fraction contents in milk and their daily yields. When the genes were included individually in the model, CSN2 was the major gene controlling all the casein fractions except for κ-CN, which was controlled directly by the CSN3 gene. The BLG gene had the most influence on the 2 whey proteins. The genetic correlations showed the major genes had only a small effect on the relationships between the protein fractions, but through comparison of the correlation coefficients of the proteins expressed in different ways they revealed potential mechanisms of regulation and competitive synthesis in the mammary gland. The estimates for the effects of the CSN2 and CSN3 genes on protein profiles showed overexpression of protein synthesis in the presence of the B allele in the genotype. Conversely, the β-LG B variant was associated with a lower concentration of β-LG compared with the β-LG A variant, independently of how the protein fractions were expressed, and it was followed by downregulation (or upregulation in the case of the β-LG B) of all other protein fractions. These results should be borne in mind when seeking to design more efficient selection programs aimed at improving milk quality for the efficiency of dairy industry and the effect of dairy products on human health.