In vitro gastrointestinal digestion of bovine αS1- and αS2-casein variants gives rise to different IgE-binding epitopes

The occurrence and differences of resistant regions containing IgE-binding epitopes of α_S1-casein (α_S1-CN) variants B and C, as well as α_S2-CN A and B, after in vitro gastrointestinal digestion was investigated using mass spectrometry. The amino acid substitutions characterising the genetic variants affected the peptide pattern arising from the caseins and thus modifications in their allergenic epitopes occurred. Peptides f174–193 in α_S1-CN B and f179–198 in α_S1-CN C correspond to the IgE-binding epitope f173–194, which has been reported as one of the major epitopes in α_S1-CN B. Within α_S2-CN, the two variant-specific peptides, f7–29 from variant A and f1–22 from variant B, contain the previously identified IgE-binding epitope f1–20. These peptides, and in consequence the protein variants, may exhibit different immunoreactions, which could be significant in the production of milk with improved nutritional properties, such as hypoallergenic quality, by selection and breeding of cows with particular milk protein genotypes.     

Genetic variants of bovine β- and κ-casein result in different immunoglobulin E-binding epitopes after in vitro gastrointestinal digestion

Immunoglobulin E-mediated allergy to cow milk is a common allergy in industrialized countries, mainly affecting young children and infants. β-Casein (CN) and κ-CN belong to the major allergens in cow milk. Within these milk proteins, genetic polymorphisms occur, which are characterized by substitutions or deletions of AA, resulting in different variants for each protein. Until now, these variants have not been considered when discussing the allergenic potential of bovine milk. In this study, the focus was placed on the arising peptide pattern after in vitro gastrointestinal digestion of several β- and κ-CN variants to determine resistant fragments containing IgE-binding epitopes and to identify potential differences between these variants. β-Casein A¹, A², and B, as well as κ-CN A, B, and E, were separated and isolated from milk of cows homozygous for these variants and digested with an in vitro gastrointestinal digestion model. The resulting peptides were identified using mass spectrometry and compared with previously determined epitopes. Seven β-CN and 4 κ-CN peptides, common in all β- or κ-CN variants, remained of sufficient size to harbor IgE-binding epitopes. In addition, some peptides and, consequently, epitopes differ from each other due to the AA substitution occurring in the individual variants. The distinct peptides AA 108 to 129 of β-CN A¹ and A², AA 103 to 123 of β-CN B, as well as AA 59 to 72, AA 59 to 80, and AA 58 to 80 of all 3 β-CN variants correspond to the IgE-binding epitopes AA 107 to 120 and AA 55 to 70, respectively. In κ-CN, the 2 variant-specific peptides AA 136 to 149 (κ-CN A, E) and AA 134 to 150 (κ-CN B) are congruent with the IgE-binding epitope AA 137 to 148. The present study shows that genetic polymorphisms affected the arising peptide pattern of the caseins and thus modifications in the IgE-binding epitopes occurred. As a consequence, the casein variants could show differences in their allergenicity. Studies investigating the allergenic potential of these different peptides are currently in progress.     

Comprehensive analysis of proteolysis during 8 months of ripening of high-cooked Old Saare cheese

We applied capillary electrophoresis, liquid chromatography coupled with tandem mass-spectrometry (MS/MS), and ultra-performance liquid chromatography to determine the composition of water-insoluble and water-soluble proteinaceous fractions of the cheese and to study in detail the degradation of caseins during 8 mo of ripening of Estonian high-temperature cooked hard cheese Old Saare. The application of high-resolution and high-accuracy MS/MS enabled identification of more than 3,000 small peptides, representing a fairly full casein peptidome containing peptides of 4 to 25 AA in length: 1,049 from β-casein (CN), 944 from α_S1-CN, 813 from α_S2-CN, and 234 from κ-CN. The majority of β-CN- and α_S1-CN-derived peptides originated from the N-terminal parts of the molecule, f6-93 and f1-124, respectively; peptides from α_S2-CN arose predominantly from the C-terminal end f100-162. At the beginning of ripening, we found a relatively high amount of peptides originating from the glycomacropeptide part of κ-CN, whereas peptides from para-κ-CN prevailed during the later stages of ripening of the cheese. The cleavage patterns of β-CN, α_S2-CN, as well as α_S1-CN, showed that primary proteolysis was started mainly by plasmin, although a low proteolytic activity of chymosin was also evident. Based on the analysis of cleavage sites, we observed a significant participation of proteolytic enzymes, including amino- and carboxypeptidases, of both mesophilic and thermophilic starter bacteria in further hydrolysis of oligopeptides during the ripening. Several new phosphopeptides were detected in the result of MS/MS data analysis. The profiles of the estimated concentrations of phosphopeptides revealed that those originating from β-CN and α_S1-CN accumulated during cheese maturation. In contrast, we did not notice any generation of phosphopeptides from the highly phosphorylated part of α_S2-CN, f25-80, presumably due to the inaccessibility of this region to the action of plasmin and chymosin. The analysis of cleavage sites and the combination of principal component and clustering analyses provided a characterization of the complex dynamics of formation and degradation of peptides during cheese maturation. We made an attempt to obtain a comprehensive picture of proteolysis during Old Saare cheese ripening on the basis of the detailed peptidomic data, including also the less abundant peptides determined by MS/MS, and complemented by the data on intact caseins and free AA and reported the results in the paper.     

Qualitative and Quantitative Determination of Caseins with Reverse-Phase and Anion-Exchange HPLC

Whole native caseinate (WNC) and casein (CN) fractions from preparative DE-52 cellulose urea columns were chromatographed using C-8 reverse-phase (RP) and DEAE-type anion-exchange (AEx) HPLC systems. With RP, α_S2-CN and κ-CN eluted first as several small peaks; α_S1-CN eluted later as two peaks, followed by β-CN peaks. With AEx, κ-CN eluted early as a group of peaks, β-CN eluted next, and α_S1-CN and α_S2-CN coeluted last. Standard curves were prepared for α_S1-CN and β-CN using RP-PHLC and showed correlation coefficients of 0.99 and 0.98, respectively. The caseins in WNC, nonfat dry milk casein, commercial casein(ates) and caseins from milks of individual cows were determined.     

Comparison of milk protein composition and rennet coagulation properties in native Swedish dairy cow breeds and high-yielding Swedish Red cows

Recent studies have reported a very high frequency of noncoagulating milk in Swedish Red cows. The underlying factors are not fully understood. In this study, we explored rennet-induced coagulation properties and relative protein profiles in milk from native Swedish Mountain and Swedish Red Polled cows and compared them with a subset of noncoagulating (NC) and well-coagulating (WC) milk samples from modern Swedish Red cows. The native breeds displayed a very low prevalence of NC milk and superior milk coagulation properties compared with Swedish Red cows. The predominant variants in both native breeds were α_S1-casein (α_S1-CN) B, β-CN A² and β-lactoglobulin (β-LG) B. For κ-CN, the B variant was predominant in the Swedish Mountain cows, whereas the A variant was the most frequent in the Swedish Red Polled. The native breeds displayed similar protein composition, but varied in content of α_S1-CN with 9 phosphorylated serines (9P) form. Within the Swedish Mountain cows, we observed a strong inverse correlation between the relative concentration of κ-CN and micelle size and a positive correlation between ionic calcium and gel firmness. For comparison, we investigated a subset of 29 NC and 28 WC milk samples, representing the extremes with regard to coagulation properties based on an initial screening of 395 Swedish Red cows. In Swedish Red, NC milk properties were found to be related to higher frequencies of β-CN A², κ-CN E and A variants, as well as β-LG B, and the predominant composite genotype of β- and κ-CN in the NC group was A²A²/AA. Generally, the A²A²/AA composite genotype was related to lower relative concentrations of κ-CN isoforms and higher relative concentrations of α_S1-, α_S2-, and β-CN. Compared with the group of WC milk samples, NC milk contained a higher fraction of α_S2-CN and α-lactalbumin (α-LA) but a lower fraction of α_S1-CN 9P. In conclusion, milk from native Swedish breeds has good characteristics for cheese milk, which could be exploited in niche dairy products. In milk from Swedish Mountain cows, levels of ionic calcium seemed to be more important for rennet-induced gel firmness than variation in the relative protein profile. In Swedish Red, lower protein content as well as higher fraction of α_S2-CN and lower fraction of α_S1-CN 9P were related to NC milk. Further, a decrease in the frequency of the composite β-κ-CN genotype A²A²/AA through selective breeding could have a positive effect on milk coagulation properties.     

Effects of breed and casein genetic variants on protein profile in milk from Swedish Red,Danish Holstein,and Danish Jersey cows

In selecting cows for higher milk yields and milk quality, it is important to understand how these traits are affected by the bovine genome. The major milk proteins exhibit genetic polymorphism and these genetic variants can serve as markers for milk composition, milk production traits, and technological properties of milk. The aim of this study was to investigate the relationships between casein (CN) genetic variants and detailed protein composition in Swedish and Danish dairy milk. Milk and DNA samples were collected from approximately 400 individual cows each of 3 Scandinavian dairy breeds: Swedish Red (SR), Danish Holstein (DH), and Danish Jersey (DJ). The protein profile with relative concentrations of α-lactalbumin, β-lactoglobulin, and α_S1-, α_S2-, κ-, and β-CN was determined for each milk sample using capillary zone electrophoresis. The genetic variants of the α_S1- (CSN1S1), β- (CSN2), and κ-CN (CSN3) genes for each cow were determined using TaqMan SNP genotyping assays (Applied Biosystems, Foster City, CA). Univariate statistical models were used to evaluate the effects of composite genetic variants, α_S1-β-κ-CN, on the protein profile. The 3 studied Scandinavian breeds differed from each other regarding CN genotypes, with DH and SR having similar genotype frequencies, whereas the genotype frequencies in DJ differed from the other 2 breeds. The similarities in genotype frequencies of SR and DH and differences compared with DJ were also seen in milk production traits, gross milk composition, and protein profile. Frequencies of the most common composite α_S1-β-κ-CN genotype BB/A²A²/AA were 30% in DH and 15% in SR, and cows that had this genotype gave milk with lower relative concentrations of κ- and β-CN and higher relative concentrations of α_S-CN, than the majority of the other composite genotypes in SR and DH. The effect of composite genotypes on relative concentrations of the milk proteins was not as pronounced in DJ. The present work suggests that a higher frequency of BB/A¹A²/AB, together with a decrease in BB/A²A²/AA, could have positive effects on DH and SR milk regarding, for example, the processing of cheese.     

Identification of immunoglobulin E epitopes on major allergens from dairy products after digestion and transportation in vitro

Dairy processing can alter the digestion stability and bioavailability of cow milk proteins in the gastrointestinal tract. However, analysis of stable linear epitopes on cow milk allergens that could enter into intestinal mucosal is limited. Thus, this study aimed to investigate the digestion and transportation properties and residual allergen epitopes entering into gastrointestinal mucosa of 3 commercial dairy products, including pasteurized milk (PM), ultra-heat-treated milk (UHTM), and dried skim milk (DSM). In this work, the digestive stability of the 3 kinds of dairy products has been performed in a standard multistep static digestion model in vitro and characterized by Tricine-SDS-polyacrylamide gel electrophoresis and reversed-phase HPLC. With respect to gastrointestinal digestion in vitro, the main allergens including β-lactoglobulin (β-LG), α-lactalbumin (α-LA), and caseins were degraded gradually, and the resistance peptides remained in the PM with a molecular weight of range from 3.4 to 5.0 kDa. Simultaneously, the potential allergenicity of the cow milk proteins was diminished gradually and is basically consistent after 60 min of gastrointestinal digestion. After gastrointestinal digestion, the remaining peptides were transported via an Ussing chamber and identified by liquid chromatography-MS/MS. By alignment, 10 epitopes peptides were identified from 16 stable peptides, including 5 peptides (AA 92–100, 125–135, 125–138, and 149–162) in β-LG, 2 peptides in α-LA (AA 80–93 and 63–79), 2 peptides in α_S1-casein (AA 84–90 and 125–132), and 1 peptide (AA 25–32) in α_S2-casein were identified by dot-blotting mainly exist in UHTM and PM. This study demonstrates dairy processing can affect the digestion and transport characteristics of milk proteins and in turn alter epitope peptides release.     

Genome-wide association study for αS1- and αS2-casein phosphorylation in Dutch Holstein Friesian

Phosphorylation of caseins (CN) is a crucial post-translational modification that allows caseins to form colloid particles known as casein micelles. Both α_S1- and α_S2-CN show varying degrees of phosphorylation (isoforms) in cow milk and were suggested to be more relevant for stabilizing internal micellar structure than β- and κ-CN. However, little is known about the genetic background of individual α_S2-CN phosphorylation isoforms and the phosphorylation degrees of α_S1- and α_S2-CN (α_S1-CN PD and α_S2-CN PD), defined as the proportion of isoforms with higher degrees of phosphorylation in total α_S1- and α_S2-CN, respectively. We aimed to identify genomic regions associated with these traits using 50K single nucleotide polymorphisms for 1,857 Dutch Holstein Friesian cows. A total of 10 quantitative trait loci (QTL) regions were identified for all studied traits on 10 Bos taurus autosomes (BTA1, 2, 6, 9, 11, 14, 15, 18, 24, and 28). Regions associated with multiple traits were found on BTA1, 6, 11, and 14. We showed 2 QTL regions on BTA1, one affecting α_S2-CN production and the other harboring the SLC37A1 gene, which encodes a phosphorus antiporter and affects α_S1- and α_S2-CN PD. The QTL on BTA6 harbors the casein gene cluster and affects individual α_S2-CN phosphorylation isoforms. The QTL on BTA11 harbors the PAEP gene that encodes for β-lactoglobulin and affects relative concentrations of α_S2-CN-10P and α_S2-CN-11P as well as α_S1-CN PD and α_S2-CN PD. The QTL on BTA14 harbors the DGAT1 gene and affects relative concentrations of α_S2-CN-10P and α_S2-CN-11P as well as α_S1-CN PD and α_S2-CN PD. Our results suggest that effects of identified genomic regions on phosphorylation of α_S1- and α_S2-CN are related to changes in milk synthesis and phosphorus secretion in milk. The actual roles of SLC37A1, PAEP, and DGAT1 in α_S1- and α_S2-CN phosphorylation in Dutch Holstein Friesian require further investigation.     

Identification of rare genetic variants of the αS-caseins in milk from native Norwegian dairy breeds and comparison of protein composition with milk from high-yielding Norwegian Red cows

Several factors influence the composition of milk. Among these, genetic variation within and between cattle breeds influences milk protein composition, protein heterogeneity, and their posttranslational modifications. Such variations may further influence technological properties, which are of importance for the utilization of milk into dairy products. Furthermore, these potential variations may also facilitate the production of differentiated products (e.g., related to specific breeds or specific genetic variants). The objective of this study was to investigate the genetic variation and relative protein composition of the major proteins in milk from 6 native Norwegian dairy breeds representing heterogeneity in geographical origin, using the modern Norwegian breed, Norwegian Red, as reference. In total, milk samples from 144 individual cows were collected and subjected to liquid chromatography-electrospray ionization/mass spectrometry–based proteomics for identification of genetic and posttranslational modification isoforms of the 4 caseins (α_S1-CN, α_S2-CN, β-CN, κ-CN) and the 2 most abundant whey proteins (α-lactalbumin and β-lactoglobulin). Relative quantification of these proteins and their major isoforms, including phosphorylations of α_S1-CN and glycosylation of κ-CN, were determined based on UV absorbance. The presence and frequency of genetic variants of the breeds were found to be very diverse and it was possible to identify rare variants of the CN, which, to our knowledge, have not been identified in these breeds before. Thus, α_S1-CN variant D was identified in low frequency in 3 of the 6 native Norwegian breeds. In general, α_S1-CN was found to be quite diverse between the native breeds, and the even less frequent A and C variants were furthermore detected in 1 and 5 of the native breeds, respectively. The α_S1-CN variant C was also identified in samples from the Norwegian Red cattle. The variant E of κ-CN was identified in 2 of the native Norwegian breeds. Another interesting finding was the identification of α_S2-CN variant D, which was found in relatively high frequencies in the native breeds. Diversity in more common protein genetic variants were furthermore observed in the protein profiles of the native breeds compared with milk from the high-yielding Norwegian Reds, probably reflecting the more diverse genetic background between the native breeds.     

Relation between αS1-casein,genotype, and quality traits of milk from Swedish dairy goats

Locally produced food is becoming popular among Swedish consumers. One product that has increased in popularity is artisan-manufactured goat cheese, and although the dairy goat industry in Sweden is small-scale, production is gradually increasing. In goats, the CSN1S1 gene regulates expression of the protein α_S1-casein (α_S1-CN), which has been found to be important for cheese yield. Over the years, breeding animals have been imported to Sweden from Norway. Historically, a high frequency of the Norwegian goat population carried a polymorphism at the CSN1S1 gene. This polymorphism, called the Norwegian null allele (D), leads to zero or significantly reduced expression of α_S1-CN. Using milk samples from 75 goats, this study investigated associations between expression of α_S1-CN and genotype at the CSN1S1 gene on milk quality traits from Swedish Landrace goats. Milk samples were grouped according to relative level of α_S1-CN (low: 0–6.9% of total protein; medium-high: 7–25% of total protein) and genotype (DD, DG, DA/AG/AA). While the D allele leads to extremely low expression of α_S1-CN, the G allele is low expressing and the A allele is highly expressing for this protein. Principal component analysis was used to explore the total variation in milk quality traits. To evaluate the effect of different allele groups on milk quality attributes, 1-way ANOVA and Tukey pairwise comparison tests were used. The majority (72%) of all goat milk samples investigated showed relative α_S1-CN content of 0% to 6.82% of total protein. The frequency of individuals homozygous for the Norwegian null allele (DD) was 59% in the population of sampled goats, and only 15% carried at least one A allele. A low relative concentration of α_S1-CN was associated with lower total protein, higher pH, and higher relative concentration of β-casein and levels of free fatty acids. Milk from goats homozygous for the null allele (DD) showed a similar pattern as milk with low relative concentration of α_S1-CN, but total protein was only numerically lower, and somatic cell count and α_S2-CN were higher than for the other genotypes. The associations between levels of α_S1-CN and the investigated genotype at the CSN1S1 gene indicate a need for a national breeding program for Swedish dairy goats.     

Whole-genome association study for milk protein composition in dairy cattle

Our objective was to perform a genome-wide association study for content in bovine milk of α_S1-casein (α_S1-CN), α_S2-casein (α_S2-CN), β-casein (β-CN), κ-casein (κ-CN), α-lactalbumin (α-LA), β-lactoglobulin (β-LG), casein index, protein percentage, and protein yield using a 50K single nucleotide polymorphism (SNP) chip. In total, 1,713 Dutch Holstein-Friesian cows were genotyped for 50,228 SNP and a 2-step association study was performed. The first step involved a general linear model and the second step used a mixed model accounting for all family relationships. Associations with milk protein content and composition were detected on 20 bovine autosomes. The main genomic regions associated with milk protein composition or protein percentage were found on chromosomes 5, 6, 11, and 14. The number of chromosomal regions showing significant (false discovery rate <0.01) effects ranged from 3 for β-CN and 3 for β-LG to 12 for α_S2-CN. A genomic region on Bos taurus autosome (BTA) 6 was significantly associated with all 6 major milk proteins, and a genomic region on BTA 11 was significantly associated with the 4 caseins and β-LG. In addition, regions were detected that only showed a significant effect on one of the milk protein fractions: regions on BTA 13 and 22 with effects on α_S1-CN; regions on BTA 1, 9, 10, 17, 19, and 28 with effects on α_S2-CN; a region on BTA 6 with an effect on β-CN; regions on BTA 13 and 21 with effects on κ-CN; regions on BTA 1, 5, 9, 16, 17, and 26 with effects on α-LA; and a region on BTA 24 with an effect on β-LG. The proportion of genetic variance explained by the SNP showing the strongest association in each of these genomic regions ranged from <1% for α_S1-CN on BTA 22 to almost 100% for casein index on BTA 11. Variation associated with regions on BTA 6, 11, and 14 could in large part but not completely be explained by known protein variants of β-CN (BTA 6), κ-CN (BTA 6), and β-LG (BTA 11) or DGAT1 variants (BTA 14). Our results indicate 3 regions with major effects on milk protein composition, in addition to several regions with smaller effects involved in the regulation of milk protein composition.     

Effects of diet on casein and fatty acid profiles of milk from goats differing in genotype for αS1-casein synthesis

This study investigated the interactions between nutrition and the genotype at α_S1-CN loci (CSN1S1) in goats, evaluating the impact of fresh forage-based diets and an energy supplement on the casein and fatty acid (FA) profiles of milk from Girgentana goats. Twelve goats were selected for having the same genotype at the α_S2-CN, β-CN, and κ-CN loci and differing in the CSN1S1 genotype: homozygous for strong alleles (AA) or heterozygous for strong and weak alleles (AF). Goats of each genotype were divided into three groups and, according to a 3 × 3 Latin square design, fed ad libitum three diets: sulla fresh forage (SFF), SFF plus 800 g/day of barley (SFB), and mixed hay plus 800 g/day of barley (MHB). The SFB diet led to higher-energy intake and milk yield. The energy-supplemented diets (SFB, MHB) reduced milk fat and urea and increased coagulation time. The fresh forage diets (SFF, SFB) increased dry matter (DM) and crude protein (CP) intake and milk β-CN. Diet had a more pronounced effect than CSN1S1 genotype on milk FA profile, which was healthier from goats fed the SFF diet, due to the higher content of rumenic acid, polyunsaturated, and omega-3 FAs. The AA milk had longer coagulation time and higher curd firmness, higher short- and medium-chain FAs (SMFA), and lower oleic acid than AF milk. Significant diet by genotype interactions indicated the higher milk yield of AA goats than AF goats with the higher-energy SFB diet and the lower synthesis of SMFA in AF than in AA goats with the SFF diet.     

Effects of milk protein variants on the protein composition of bovine milk

The effects of β-lactoglobulin (β-LG), β-casein (β-CN), and κ-CN variants and β-κ-CN haplotypes on the relative concentrations of the major milk proteins α-lactalbumin (α-LA), β-LG, α_S1-CN, α_S2-CN, β-CN, and κ-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 β-LG genotype was associated with the relative concentrations of β-LG (A » B) and of α-LA, α_S1-CN, α_S2-CN, β-CN, and κ-CN (B > A) but not with any milk production trait. The β-CN genotype was associated with the relative concentrations of β-CN and α_S2-CN (A² > A¹) and of α_S1-CN and κ-CN (A¹ > A²) and with protein yield (A² > A¹). The κ-CN genotype was associated with the relative concentrations of κ-CN (B > E > A), α_S2-CN (B > A), α-LA, and α_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 β-LG genotype B and the β-κ-CN haplotype A²B will result in cows that produce milk that is more suitable for cheese production.     

Genetic and nongenetic factors contributing to differences in αS-casein phosphorylation isoforms and other major milk proteins

Relative concentrations of α_S-casein (α_S-CN) phosphorylation isoforms vary considerably among milk of individual cows. We aimed to explore to what extent genetic and other factors contribute to the variation in relative concentrations of α_S-CN phosphorylation isoforms and the phosphorylation degree of α_S-CN defined as the proportion of isoforms with higher degrees of phosphorylation. We also investigated the associations of genetic variants of milk proteins and casein haplotypes with relative concentrations of α_S-CN phosphorylation isoforms and with the phosphorylation degree of α_S-CN in French Montbéliarde cattle from the cheese production area of Franche-Comté. Detailed milk protein composition was determined by liquid chromatography coupled with electrospray ionization mass spectrometry from 531 test-day morning milk samples. Parity, lactation stage, and genetic variation of cows contributed to the phenotypic variation in relative concentrations of individual α_S-CN phosphorylation isoforms and in the phosphorylation degree of α_S-CN. As lactation progressed, we observed a significant increase for relative concentrations of α_S-CN isoforms with higher degrees of phosphorylation (α_S1-CN-9P, α_S2-CN-13P, and α_S2-CN-14P) as well as for the phosphorylation degree of both α_S1-CN and α_S2-CN. Furthermore, the β-CN I variant was associated with a greater proportion of isoforms with lower degrees of phosphorylation (α_S1-CN-8P, α_S2-CN-10P, and α_S2-CN-11P); the β-CN B variant was associated with a greater proportion of isoforms with higher degrees of phosphorylation (α_S1-CN-9P, α_S2-CN-12P to α_S2-CN-14P). The heritability estimates were low to moderate for relative concentrations of α_S2-CN phosphorylation isoforms (0.07 to 0.32), high for relative concentrations of α_S1-CN-8P (0.84) and α_S1-CN-9P (0.56), and moderate for phosphorylation degrees of α_S1-CN (0.37) and α_S2-CN (0.23). Future studies investigating relations between the phosphorylation degree of α_S-CN and technological properties of milk will be beneficial for the dairy industry.     

Effect of protein concentrate mixtures and dietary addition of exogenous phytase on major milk minerals and proteins,including casein phosphorylation

Variations in major milk minerals, proteins, and their posttranslational modifications are largely under genetic influence, whereas the effect of nongenetic factors is less studied. Through a controlled feeding experiment (incomplete balanced Latin square design), the effect of concentrate mixtures, based on fava beans, rapeseed meal, or soybean meal as main P and protein sources, on milk composition was examined under typical Danish management conditions. Concentrations of P, Ca, and Mg, together with proteomics for relative quantification of major milk proteins and their isoforms, were analyzed in milk samples from 24 cows sampled in 4 periods. Each cow was fed 1 of the 3 diets in each period with or without addition of exogenous phytase. Cows were blocked by lactation stage into early and mid-lactation (23.3 ± 6.7 and 176 ± 15 d in milk, respectively, at the beginning of the experiment, mean ± standard deviation). Significant effects of feed concentrate mixture were observed for milk protein concentration, milk urea nitrogen, citrate, and the percentage of mixed and preformed fatty acids as well as mineral composition, and their distributions within micellar or serum phases. Furthermore, relative contents of α_S1-casein (CN) 9P form and unglycosylated κ-CN and thereby phosphorylation degree of α_S1-CN (PD) and the glycosylation degree of κ-CN were found to be significantly affected by these diets. To our knowledge, we are the first to document that feed concentrate mixture can affect the relative concentrations of α_S1-CN phosphorylation isoforms in milk, and the results suggested an effect on α_S1-CN 9P and PD, but not on α_S1-CN 8P. Furthermore, although only significant for α_S1-CN 8P, we found a lower relative concentration of α_S1-CN 8P and higher α_S1-CN 9P (and thus higher PD) in milk from cows in mid compared with early lactation. Also, protein concentration and concentration of Mg in skim milk and serum as well as relative concentration of α-lactalbumin were found to be significantly affected by lactation stage. Addition of dietary exogenous phytase only had a minor effect on milk composition or functionality with significant effect detected for α-lactalbumin and micellar Mg concentration.     

Effect of protein composition on the cheese-making properties of milk from individual dairy cows

The objective of this study was to evaluate the effect of variations in milk protein composition on milk clotting properties and cheese yield. Milk was collected from 134 dairy cows of Swedish Red and White, Swedish Holstein, and Danish Holstein-Friesian breed at 3 sampling occasions. Concentrations of α_S1-, β-, and κ-casein (CN), α-lactalbumin, and β-lactoglobulin (LG) A and B were determined by reversed phase liquid chromatography. Cows of Swedish breeds were genotyped for genetic variants of β- and κ-CN. Model cheeses were produced from individual skimmed milk samples and the milk clotting properties were evaluated. More than 30% of the samples were poorly coagulating or noncoagulating, resulting in weak or no coagulum, respectively. Poorly and noncoagulating samples were associated with a low concentration of κ-CN and a low proportion of κ-CN in relation to total CN analyzed. Furthermore, the κ-CN concentration was higher in milk from cows with the AB genotype than the AA genotype of κ-CN. The concentrations of α_S1-, β-, and κ-CN and of β-LG B were found to be significant for the cheese yield, expressed as grams of cheese per one hundred grams of milk. The ratio of CN to total protein analyzed and the β-LG B concentration positively affected cheese yield, expressed as grams of dry cheese solids per one hundred grams of milk protein, whereas β-LG A had a negative effect. Cheese-making properties could be improved by selecting milk with high concentrations of α_S1-, β-, and κ-CN, with high κ-CN in relation to total CN and milk that contains β-LG B.     

Structural equation modeling for unraveling the multivariate genomic architecture of milk proteins in dairy cattle

The aims of this study were to investigate potential functional relationships among milk protein fractions in dairy cattle and to carry out a structural equation model (SEM) GWAS to provide a decomposition of total SNP effects into direct effects and effects mediated by traits that are upstream in a phenotypic network. To achieve these aims, we first fitted a mixed Bayesian multitrait genomic model to infer the genomic correlations among 6 milk nitrogen fractions [4 caseins (CN), namely κ-, β-, α_S1-, and α_S2-CN, and 2 whey proteins, namely β-lactoglobulin (β-LG) and α-lactalbumin (α-LA)], in a population of 989 Italian Brown Swiss cows. Animals were genotyped with the Illumina BovineSNP50 Bead Chip v.2 (Illumina Inc.). A Bayesian network approach using the max-min hill-climbing (MMHC) algorithm was implemented to model the dependencies or independence among traits. Strong and negative genomic correlations were found between β-CN and α_S1-CN (?0.706) and between β-CN and κ-CN (?0.735). The application of the MMHC algorithm revealed that κ-CN and β-CN seemed to directly or indirectly influence all other milk protein fractions. By integrating multitrait model GWAS and SEM-GWAS, we identified a total of 127 significant SNP for κ-CN, 89 SNP for β-CN, 30 SNP for α_S1-CN, and 14 SNP for α_S2-CN (mostly shared among CN and located on Bos taurus autosome 6) and 15 SNP for β-LG (mostly located on Bos taurus autosome 11), whereas no SNP passed the significance threshold for α-LA. For the significant SNP, we assessed and quantified the contribution of direct and indirect paths to total marker effect. Pathway analyses confirmed that common regulatory mechanisms (e.g., energy metabolism and hormonal and neural signals) are involved in the control of milk protein synthesis and metabolism. The information acquired might be leveraged for setting up optimal management and selection strategies aimed at improving milk quality and technological characteristics in dairy cattle.     

Effectiveness of mid-infrared spectroscopy for the prediction of detailed protein composition and contents of protein genetic variants of individual milk of Simmental cows

Mid-infrared (MIR) spectroscopy was used to predict the detailed protein composition of 1,517 milk samples of Simmental cows. Contents of milk protein fractions and genetic variants were quantified by reversed-phase HPLC. The most accurate predictions were those obtained for total protein, casein (CN), α_S1-CN, β-lactoglobulin (LG), glycosylated κ-CN, and whey protein content, which exhibited coefficients of determination between predicted and measured values in cross-validation (1-VR) ranging from 0.61 to 0.78. Less favorable were results for β-CN (1-VR = 0.53), α_S2-CN, and κ-CN (1-VR = 0.49). Neither the content of α-LA nor that of γ-CN was accurately predicted by MIR. Predicting the content of the most common milk protein genetic variants (κ-CN A and B; β-CN A¹, A², and B; and β-LG A and B) was unfeasible (1-VR <0.15 for the content of κ-CN genetic variants and 1-VR <0.01 for the content of β-CN variants). The best predictions were obtained for β-LG A and β-LG B contents (1-VR of 0.60 and 0.44, respectively). Results indicated that MIR is not applicable for predicting individual milk protein composition with high accuracy. However, MIR spectroscopy predictions may play a role as indicator traits in selective breeding to enhance milk protein composition. The genetic correlation between MIR spectroscopy predictions and measures of milk protein composition needs to be investigated, as it affects the suitability of MIR spectroscopy predictions as indicator traits in selective breeding.     

Genetic parameters for αS1-casein and αS2-casein phosphorylation isoforms in Dutch Holstein Friesian

Relative concentrations of α_S1-casein and α_S2-casein (α_S1-CN and α_S2-CN) phosphorylation isoforms vary considerably among milk of individual cows. We estimated heritabilities for α_S2-CN phosphorylation isoforms, determined by capillary zone electrophoresis from 1,857 morning milk samples, and genetic correlations among α_S2-CN phosphorylation isoforms in Dutch Holstein Friesian. To investigate if phosphorylation of α_S1-CN and α_S2-CN are due to the same genetic mechanism, we also estimated genetic correlations between α_S1-CN and α_S2-CN phosphorylation isoforms as well as the genetic correlations between the phosphorylation degrees (PD) of α_S1-CN and α_S2-CN defined as the proportion of isoforms with higher degrees of phosphorylation in total α_S1-CN and α_S2-CN, respectively. The intra-herd heritabilities for the relative concentrations of α_S2-CN phosphorylation isoforms were high and ranged from 0.54 for α_S2-CN-10P to 0.89 for α_S2-CN-12P. Furthermore, the high intra-herd heritabilities of α_S1-CN PD and α_S2-CN PD imply a strong genetic control of the phosphorylation process, which is independent of casein production. The genetic correlations between α_S2-CN phosphorylation isoforms are positive and moderate to high (0.33–0.90). Furthermore, the strong positive genetic correlation (0.94) between α_S1-CN PD and α_S2-CN PD suggests that the phosphorylation processes of α_S1-CN and α_S2-CN are related. This study shows the possibility of breeding for specific α_S1-CN and α_S2-CN phosphorylation isoforms, and relations between the phosphorylation degrees of α_S1-CN and α_S2-CN and technological properties of milk need to be further investigated to identify potential benefits for the dairy industry.     

Impact of somatic cell count combined with differential somatic cell count on milk protein fractions in Holstein cattle

Udder health in dairy herds is a very important issue given its implications for animal welfare and the production of high-quality milk. Somatic cell count (SCC) is the most widely used means of assessing udder health status. However, differential somatic cell count (DSCC) has recently been proposed as a new and more effective means of evaluating intramammary infection dynamics. Differential SCC represents the combined percentage of polymorphonuclear neutrophils and lymphocytes (PMN-LYM) in the total SCC, with macrophages (MAC) accounting for the remaining proportion. The aim of this study was to evaluate the association between SCC and DSCC and the detailed milk protein profile in a population of 1,482 Holstein cows. A validated reversed-phase HPLC method was used to quantify 4 caseins (CN), namely α_S1-CN, α_S2-CN, κ-CN, and β-CN, and 3 whey protein fractions, namely β-lactoglobulin, α-lactalbumin, and lactoferrin, which were expressed both quantitatively (g/L) and qualitatively (as a percentage of the total milk nitrogen content, %N). A linear mixed model was fitted to explore the associations between somatic cell score (SCS) combined with DSCC and the protein fractions expressed quantitatively and qualitatively. We ran an additional model that included DSCC expressed as PMN-LYM and MAC counts, obtained by multiplying the percentages of PMN-LYM and MAC by SCC for each cow in the data set. When the protein fractions were expressed as grams per liter, SCS was significantly negatively associated with almost all the casein fractions and positively associated with the whey protein α-lactalbumin, while DSCC was significantly associated with α_S1-CN, β-CN, and α-lactalbumin, but in the opposite direction to SCS. We observed the same pattern with the qualitative data (i.e., %N), confirming opposite effects of SCS and DSCC on milk protein fractions. The PMN-LYM count was only slightly associated with the traits of concern, although the pattern observed was the same as when both SCS and DSCC were included in the model. The MAC count, however, generally had a greater impact on many casein fractions, in particular decreasing both β-CN content (g/L) and proportion (%N), and exhibited the opposite pattern to the PMN-LYM count. Our results show that information obtained from both SCS and DSCC may be useful in assessing milk quality and protein fractions. They also demonstrate the potential of MAC count as a novel udder health trait.