Estimation of variation in concentration, phosphorylation and genetic polymorphism of milk proteins using capillary zone electrophoresis

Advances in separation and identification of biologically important milk proteins and peptides

Milk is a rich source of biologically important proteins and peptides. In addition, milk contains a variety of extracellular vesicles (EVs), including exosomes, that carry their own proteome cargo. EVs are essential for cell-cell communication and modulation of biological processes. They act as nature carriers of bioactive proteins/peptides in targeted delivery during various physiological and pathological conditions. Identification of the proteins and protein-derived peptides in milk and EVs and recognition of their biological activities and functions had a tremendous impact on food industry, medicine research, and clinical applications. Advanced separation methods, mass spectrometry (MS)-based proteomic approaches and innovative biostatistical procedures allowed for characterization of milk protein isoforms, genetic/splice variants, posttranslational modifications and their key roles, and contributed to novel discoveries. This review article discusses recently published developments in separation and identification of bioactive proteins/peptides from milk and milk EVs, including MS-based proteomic approaches.

The effect of calcium removal from skim milk by ion exchange on the properties of the ultrafiltration retentate

New processes are needed to produce concentrated milk feedstocks with tailored calcium content, due to the direct link between calcium concentration and final product texture and functionality. Skim milk treatment with cation exchange resin 1% (w/v) or 2% (w/v) prior to ultrafiltration to a volumetric concentration factor (VCF) of 2.5 or 5 successfully decreased the calcium concentration by 20-30% and produced concentrates with solids content at ∼22-24 g 100 g-1 at a VCF of 5. Calcium reduction partially solubilized the casein micelles, increasing the concentration of soluble protein and individual caseins, leading to decreased turbidity but increased protein hydration and hydrophobicity. Decalcification (2% (w/v) resin treatment) reduced thermal stability, significantly decreasing the denaturation temperature of α-lactalbumin and β-lactoglobulin in the milk by ∼3 °C and ∼1 °C respectively. Filtration was also altered, reducing permeation flux and the gel concentration and increased filtration time. When combined, calcium reduction and filtration altered functional properties including soluble calcium, soluble protein and sedimentable solids, with increased milk protein hydration also contributing to increased viscosity. This study provides a route to produce calcium-reduced milk concentrates with potential for use in retentate-based dairy products with tailored functionality.

Polymorphism at the CSN1S1 Locus and Energy Intake Level Affect Milk Traits and Casein Profiles in Rossa Mediterranea Goats

A total of twenty-seven Rossa Mediterranea lactating goats, consisting of nine homozygous for strong alleles (AA), twelve heterozygous (AF) and six homozygous for weak alleles (FF) at the CSN1S1 locus, were used to evaluate the effect of genotype, diet and genotype × diet interaction on goat milk traits and casein profile. The goats were used in a 3 × 3 factorial arrangement of treatments, with three genotypes (AA, AF and FF) and three different energy intake levels: high (H), medium (M) and low (L). The diets supplied a complete pelleted feed containing 65% of alfalfa hay, respectively, at 150%, 100% and 70% of the total energy requirements. Milk yield was significantly affected by the genotype and diet: Lower levels were found in FF goats than in AA and AF genotypes (673.7 vs. 934.5 and 879.8 d/g, respectively; p = 0.002) as well as in goats fed with the L diet (651.5 vs. 1041 and 852.9 g/d for H and M diet, respectively, p < 0.001). The genotype influenced the casein profile. Specifically, AA goat milk exhibited higher concentrations of total casein and α_s1-casein compared to AF and FF genotypes (for total casein and αs1-casein, respectively: 24.9 vs. 20.4 and 19.8 g/kg, p = 0.001; 7.2 vs. 3.7 and 0.7 g/kg, p < 0.001), while the FF genotype showed higher values for α_s2-casein concentrations compared to homozygous AA and heterozygous AF goats (3.1 vs. 2.4 and 2.5 g/kg, respectively, p < 0.001). A significant genotype x diet interaction occurred for αs2-casein levels (g/kg) (p = 0.034) and αs₁-casein yields (p = 0.027): The αs₂-casein level was not affected by the diet in AA goats, whereas it increased with energy intake in AF and FF genotypes. Conversely, the αs₁-casein yield gradually increased with energy intake in AA and AF groups, whereas the diet in FF goats did not modify it. The results demonstrated that high energy input, as well as the strong allele at the CSN1S1 locus, enhanced milk production and casein concentrations. Furthermore, they confirmed the existence of an interaction between α_s1-casein polymorphism and diets, influencing the milk casein composition and yield.

Separation and characterization of bovine milk proteins by capillary electrophoresis-mass spectrometry

Protein profiling of major bovine milk proteins (i.e., whey and casein proteins) is of great interest in food science and technology. This complex set of protein proteoforms may vary with breed, genetics, lactation stage, health, and nutritional status of the animal. Current routine methods for bovine milk protein profiling at the intact level are typically based on capillary electrophoresis-ultraviolet, which does not allow confirming unequivocally the identity of the separated proteins. As an alternative, in this study, we describe for the first time a novel and simple capillary electrophoresis-mass spectrometry method in positive electrospray ionization mode. Under the optimized conditions, capillary electrophoresis-mass spectrometry allowed the separation and identification at the intact level of major bovine milk whey and casein proteins in less than 15 min. Furthermore, high-resolution mass spectrometry confirmed its importance in the reliable characterization of bovine milk protein proteoforms, especially those with slight molecular mass differences, such as β-casein A1 and A2, which are relevant to unequivocally identify milk with specific β-casein compositions (e.g., A2A2 milk, which is widely known as A2 milk). This differentiation was not possible by matrix-assisted laser desorption/ionization mass spectrometry, which provided rapidly and easily a rich but less accurate fingerprint of bovine milk proteins due to the lower mass resolution.

Short communication: Caseins and α-lactalbumin content of camel milk (Camelus dromedarius) determined by capillary electrophoresis

Camel milk has unique physical, nutritional, and technological properties when compared with other milks, especially bovine. Because proteins confer many of the properties of milk and its products, this study aimed to determine the proteins of camel milk, their correlations, and relative distribution. Raw milk samples were collected from 103 dromedary camels in the morning and evening. Capillary electrophoresis results showed wide variation in the concentrations (g/L) of proteins between samples as follows: α-lactalbumin, 0.3 to 2.9; α_S1-casein, 2.4 to 10.3; α_S2-casein, 0.3 to 3.9; β-casein, 5.5 to 29.0; κ-casein, 0.1 to 2.4; unknown casein protein 1, 0.0 to 3.4; and unknown casein protein 2, 0.0 to 4.6. The range in percent composition of the 4 caseins were as follows: α_S1, 12.7 to 35.3; α_S2, 1.8 to 20.8; β, 42.3 to 77.4; and κ, 0.6 to 17.4. The relative proportion of α_S1-, α_S2-, β-, and κ-caseins in camel milk (26:4:67:3, wt/wt) differed from that of bovine milk (38:10:36:12, wt/wt). This difference might explain the dissimilarity between the 2 milks with respect to technical and nutritional properties.

Quantitative and qualitative detailed milk protein profiles of 6 cattle breeds: Sources of variation and contribution of protein genetic variants

Different fractions of milk nitrogenous compounds (not only caseins) have different effects on the nutritional value of milk, its coagulation and curd firming properties, and its cheese-making efficiency. To assess different sources of variation, especially the cows' breed and genetic variants of the main protein fractions, milk samples were collected from 1,504 cows belonging to 3 dairy breeds (Holstein-Friesian, Brown Swiss, and Jersey) and 3 dual-purpose breeds (Simmental, Rendena, and Alpine Grey) reared in 41 multibreed herds. Beyond crude protein, casein (CN), and urea, 7 protein fractions were analyzed using HPLC, and 5 other N fraction traits were calculated. All 15 traits were measured qualitatively (% of milk N) and quantitatively (g/L of milk). The HPLC technique allowed us to discriminate between the main genetic variants of β-CN, κ-CN, and β-lactoglobulin and thus to genotype the cows for the CSN2, CSN3, and BLG genes, respectively. Data were analyzed using 2 mixed models, both including the effects of herd-date, breed, parity, and lactation stage, and only one also including the effects of the genotypes of the milk proteins. Breed of cow explained 2 to 36% of phenotypic variability for all the N fractions, with the exception of the urea and total casein contents of milk and the urea and β-CN proportions of total milk N. Lactation stage had a considerable influence on the amount (g/L) of almost all the protein fractions in milk, but neither the nonprotein N fractions nor the percentage of milk N protein profile were affected. The inclusion of the CSN2, CSN3, and BLG genotypes in the model explained a large part of the total variability in all the milk protein and nonprotein fractions except urea. It also reduced the variance explained by breed and residual factors. An exception was shown by the proportion of α_S1-CN variance explained by breed that moved from 13 to 28%. Similarly, for amount (g/L) of β-CN, the effect of breed became significant (12%), whereas it was almost null before inclusion of genotypes. In terms of percentage of milk N, the genotypes of CSN3 notably affected all the casein fractions, whereas the BLG genotypes had a much greater influence on most noncasein traits. The genotypes of the CSN2 gene exerted an appreciable effect on α_S2-CN and not β-CN, as expected. Comparing the 2 models, we were also able to discriminate the effect of the breed on a milk N fraction, both quantitatively and qualitatively, in 2 quotas: the first due to the milk protein polymorphisms (major genes) and the second due to other genetic factors (polygene), after correcting for the effect of herd-date of sampling, parity, and lactation stage. The knowledge about the detailed milk protein profile of different cattle breeds provided by this study could be of great benefit for the dairy industry, providing new tools for the enhancement of milk payment systems and breeding program designs.

An improved LC-MS method to profile molecular diversity and quantify the six main bovine milk proteins, including genetic and splicing variants as well as post-translationally modified isoforms

•

Identification of the main milk proteins isoforms are inferred from a mass database.

•

Quantification of the six main milk proteins is achieved from corrected UV at 214 nm.

•

Multiple isoforms can be quantified from mass signal, within each protein family.

•

Glycosylation isoforms of κ-casein which impact micelle stability are quantified.

Here we describe a method based on Liquid Chromatography coupled with Mass Spectrometry (LC-MS) that provides an accurate determination of the six main bovine milk proteins, including allelic and splicing variants, as well as isoforms resulting from post-translational modifications, with an unprecedented level of resolution. Proteins are identified from observed molecular masses in comparison with theoretical masses of intact proteins indexed in an “in-house” database that includes nearly 3000 entries. Quantification was performed either from UV (214 nm) or mass signals. Thus, up to one hundred molecules, derived from the six major milk proteins, can be identified and quantified from an individual milk sample. This powerful and reliable method, initially developed as an anchoring method to estimate the composition of the six main bovine milk proteins from MIR spectra, is transferable to several mammalian species, including small ruminants, camels, equines, rabbits, etc., for which specific mass databases are available.

Ripening of Hard Cheese Produced from Milk Concentrated by Reverse Osmosis

The application of reverse osmosis (RO) for preconcentration of milk (RO-milk) on farms can decrease the overall transportation costs of milk, increase the capacity of cheese production, and may be highly attractive from the cheese manufacturer’s viewpoint. In this study, an attempt was made to produce a hard cheese from RO-milk with a concentration factor of 1.9 (RO-cheese). Proteolysis, volatile profiles, and sensory properties were evaluated throughout six months of RO-cheese ripening. Moderate primary proteolysis took place during RO-cheese ripening: about 70% of α_s1-casein and 45% of β-casein were hydrolyzed by the end of cheese maturation. The total content of free amino acids (FAA) increased from 4.3 to 149.9 mmol kg⁻¹, with Lys, Pro, Glu, Leu, and γ-aminobutyric acid dominating in ripened cheese. In total, 42 volatile compounds were identified at different stages of maturation of RO-cheese; these compounds have previously been found in traditional Gouda-type and hard-type cheeses of prolonged maturation. Fresh RO-cheese was characterized by a milky and buttery flavor, whereas sweetness, saltiness, and umami flavor increased during ripening. Current results prove the feasibility of RO-milk for the production of hard cheese with acceptable sensory characteristics and may encourage further research and implementation of RO technology in cheese manufacture.

Bacterial proteolysis of casein leading to UHT milk gelation: An applicative study

Heat-stable peptidases released in refrigerated raw milk by psychrotrophic bacteria are responsible for UHT milk gelation. K-casein-derived caseinomacropeptides, identified by mass spectrometry, were constantly detected in gelled milk by capillary electrophoresis. Strains of Pseudomonas fluorescens, Ps. poae and Chryseobacterium joostei, selected among aprX-positive strains from raw milk, were incubated in milk up to 6 days at 4 °C before sterilization (98 °C/4 min). Samples were then stored at 25 or 40 °C, visually observed for gelation, and analysed for presence of caseinomacropeptides throughout 90 days of storage. Depending on cold pre-incubation time, caseinomacropeptides accumulated well before gelation onset in milk stored at 25 °C. Caseinomacropeptides were successively degraded, especially in milk stored at 40 °C, due to extensive proteolysis, and an abundant sediment developed instead of a gel. The caseinomacropeptides are here presented as an early indicator of UHT milk gelation and a mechanism explaining this phenomenon is proposed.

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.

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.

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.

Production of Cow's Milk Free from Beta-Casein A1 and Its Application in the Manufacturing of Specialized Foods for Early Infant Nutrition

Beta-casein (BC) is frequently expressed as BC A2 and BC A1 in cow’s milk. Gastrointestinal digestion of BC A1 results in the release of the opioid peptide beta-casomorphin 7 (BCM7) which is less likely to occur from BC A2. This work was aimed to produce milk containing BC A2 with no BC A1 (BC A2 milk) using genetically selected CSN2 A2A2 Jersey cows. Additionally, we aimed to develop an infant formula (IF) suitable for healthy full-term infants during the first six months of life based on BC A2 milk. The concentration of BCM7 released from BC A2 IF, from commercially available IFs as well as from human milk and raw cow’s milk was evaluated after simulated gastrointestinal digestion (SGID). BC A2 IF presented the lowest mean relative abundance of BC A1 (IF 1 = 0.136 ± 0.010), compared with three commercially available IFs (IF 2 = 0.597 ± 0.020; IF 3 = 0.441 ± 0.014; IF 4 = 0.503 ± 0.011). Accordingly, SGID of whole casein fraction from BC A2 IF resulted in a significantly lower release of BCM7 (IF 1 = 0.860 ± 0.014 µg/100 mL) compared to commercially available IFs (IF 2 = 2.625 ± 0.042 µg/100 mL; IF 3 = 1.693 ± 0.012 µg/100 mL; IF 4 = 1.962 ± 0.067 µg/100 mL). Nevertheless, BCM7 levels from BC A2 IF were significantly higher than those found in SGID hydrolysates of BC A2 raw milk (0.742 ± 0.008 µg/100 mL). Interestingly, results showed that BCM7 was also present in human milk in significantly lower amounts (0.697 ± 0.007 µg/100 mL) than those observed in IF 1 and BC A2 milk. This work demonstrates that using BC A2 milk in IF formulation significantly reduces BCM7 formation during SGID. Clinical implications of BC A2 IF on early infant health and development need further investigations.

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.

Characterizing a region on BTA11 affecting β-lactoglobulin content of milk using high-density genotyping and haplotype grouping

Background

Milk β-lactoglobulin (β-LG) content is of interest as it is associated with nutritional and manufacturing properties. It is known that milk β-LG content is strongly affected by genetic factors. In cattle, most of the genetic differences are associated with a chromosomal region on BTA11, which contains the β-LG gene. The aim of this study was to characterize this region using 777 k SNP data (BovineHDbeadChip) and perform a haplotype-based association study. A statistical approach was developed to build haplotypes that capture the genetic variation associated with this genomic region.

Results

The SNP with the most significant effect on β-lactoglobulin content was one of the 2 causal mutations responsible for the β-lactoglobulin protein variants A/B. Haplotypes based on 2 to 5 selected lead SNP were clustered in groups with different effects on β-lactoglobulin content. Four different groups were identified suggesting that β-lactoglobulin variant A and B can be further refined in A₁, A₂, B₁ and B₂.

Conclusions

This study showed that β-lactoglobulin protein variants A/B do not explain all genetic variation associated with the tail part of BTA11 but this region contains more than one mutation with an effect on β-lactoglobulin content. These findings can be used for selection of cows with higher cheese yield, which is desirable for the dairy industry.

Quantitation and Identification of Intact Major Milk Proteins for High-Throughput LC-ESI-Q-TOF MS Analyses

Cow's milk is an important source of proteins in human nutrition. On average, cow's milk contains 3.5% protein. The most abundant proteins in bovine milk are caseins and some of the whey proteins, namely beta-lactoglobulin, alpha-lactalbumin, and serum albumin. A number of allelic variants and post-translationally modified forms of these proteins have been identified. Their occurrence varies with breed, individuality, stage of lactation, and health and nutritional status of the animal. It is therefore essential to have reliable methods of detection and quantitation of these proteins. Traditionally, major milk proteins are quantified using liquid chromatography (LC) and ultra violet detection method. However, as these protein variants co-elute to some degree, another dimension of separation is beneficial to accurately measure their amounts. Mass spectrometry (MS) offers such a tool. In this study, we tested several RP-HPLC and MS parameters to optimise the analysis of intact bovine proteins from milk. From our tests, we developed an optimum method that includes a 20-28-40% phase B gradient with 0.02% TFA in both mobile phases, at 0.2 mL/min flow rate, using 75°C for the C8 column temperature, scanning every 3 sec over a 600-3000 m/z window. The optimisations were performed using external standards commercially purchased for which ionisation efficiency, linearity of calibration, LOD, LOQ, sensitivity, selectivity, precision, reproducibility, and mass accuracy were demonstrated. From the MS analysis, we can use extracted ion chromatograms (EICs) of specific ion series of known proteins and integrate peaks at defined retention time (RT) window for quantitation purposes. This optimum quantitative method was successfully applied to two bulk milk samples from different breeds, Holstein-Friesian and Jersey, to assess differences in protein variant levels.

The relationships among bovine αS-casein phosphorylation isoforms suggest different phosphorylation pathways

Casein (CN) phosphorylation is an important posttranslational modification and is one of the key factors responsible for constructing and stabilizing casein micelles. Variation in phosphorylation degree of αS-CN is of great interest because it is suggested to affect milk technological properties. This study aimed to investigate the variation in phosphorylation degree of αS-CN among milk of individual cows and to explore relationships among different phosphorylation isoforms of αS-CN. For this purpose, we analyzed morning milk samples from 529 French Montbéliarde cows using liquid chromatography coupled with electrospray ionization mass spectrometry. We detected 3 new phosphorylation isoforms: αS2-CN-9P, αS2-CN-14P, and αS2-CN-15P in bovine milk, in addition to the known isoforms αS1-CN-8P, αS1-CN-9P, αS2-CN-10P, αS2-CN-11P, αS2-CN-12P, and αS2-CN-13P. The relative concentrations of each αS-CN phosphorylation isoform varied considerably among individual cows. Furthermore, the phenotypic correlations and hierarchical clustering suggest at least 2 regulatory systems for phosphorylation of αS-CN: one responsible for isoforms with lower levels of phosphorylation (αS1-CN-8P, αS2-CN-10P, and αS2-CN-11P), and another responsible for isoforms with higher levels of phosphorylation (αS1-CN-9P, αS2-CN-12P, αS2-CN-13P, and αS2-CN-14P). Identifying all phosphorylation sites of αS2-CN and investigating the genetic background of different αS2-CN phosphorylation isoforms may provide further insight into the phosphorylation mechanism of caseins.

Phosphorylation of αS1-casein is regulated by different genes

Casein phosphorylation is a posttranslational modification catalyzed by kinase enzymes that attach phosphate groups to specific AA in the protein sequence. This modification is one of the key factors responsible for the stabilization of calcium phosphate nanoclusters in casein micelles and for the internal structure of the casein micelles. α(S1)-Casein (α(s1)-CN) is of special interest because it constitutes up to 40% of the total casein fraction in milk, and it has 2 common phosphorylation states, with 8 (α(S1)-CN-8P) and 9 (α(S1)-CN-9P) phosphorylated serine residues. Factors affecting this variation in the degree of phosphorylation are not currently known. The objective of this research was to determine the genetic background of α(S1)-CN-8P and α(S1)-CN-9P. The genetic and phenotypic correlation between α(S1)-CN-8P and α(S1)-CN-9P was low (0.18 and 0.19, respectively). This low genetic correlation suggests a different genetic background. These differences were further investigated by means of a genome-wide association study, which showed that both α(S1)-CN-8P and α(S1)-CN-9P were affected by a region on Bos taurus autosome (BTA) 6, but only α(S1)-CN-8P was affected by a region on BTA11 that contains the gene that encodes for β-lactoglobulin (β-LG), and only α(S1)-CN-9P was affected by a region on BTA14 that contains the diacylglycerol acyltransferase 1 (DGAT1) gene. Estimated effects of β-LG protein genotypes showed that only α(S1)-CN-8P was associated with the β-LG A/B polymorphism (g.1772G>A and g.3054C>T); the AA genotype of β-LG was associated with a lower concentration of α(S1)-CN-8P (-0.32% wt/wt) than the BB genotype (+0.41% wt/wt). Estimated effects of DGAT1 K232A genotypes showed that only α(S1)-CN-9P was associated with the DGAT1 gene polymorphism; DGAT1 AA genotype was associated with a higher α(S1)-CN-9P concentration (+0.53% wt/wt) than the DGAT1 KK genotype (-0.44% wt/wt). The results give insight in phosphorylation of α(S1)-CN-8P and α(S1)-CN-9P, which seem to be regulated by a different set of genes.

Casein breakdown in bovine milk by a field strain of Staphylococcus aureus

The objectives of this study were to establish the proteolytic effects of Staphylococcus aureus during mastitis on economically important milk proteins. Concentrations of milk proteins were determined by capillary electrophoresis in an experimental model using a field strain of S. aureus. The pathogen was inoculated into bacteria-free control milk to imitate proteolysis caused by the pathogen in the mammary gland between milkings. Milk content of caseins (CN) α(S1), α(S2), κ, β(A1), and β(A2) and whey proteins α-lactalbumin and β-lactoglobulin were analyzed initially and after 6 h of incubation. After 6 h, the overall CN content was significantly reduced (21%) in milk inoculated with S. aureus compared with the bacteria-free control milk. S. aureus significantly lowered concentration of α(S1)-CN (2.5%), β(A1)-CN (3%), and β(A2)-CN (5%). S. aureus also hydrolyzed κ-CN into para-κ-CN, with significant reduction of κ-CN (7.4%) as a consequence.

The occurrence of noncoagulating milk and the association of bovine milk coagulation properties with genetic variants of the caseins in 3 Scandinavian dairy breeds

Substantial variation in milk coagulation properties has been observed among dairy cows. Consequently, raw milk from individual cows and breeds exhibits distinct coagulation capacities that potentially affect the technological properties and milk processing into cheese. This variation is largely influenced by protein composition, which is in turn affected by underlying genetic polymorphisms in the major milk proteins. In this study, we conducted a large screening on 3 major Scandinavian breeds to resolve the variation in milk coagulation traits and the frequency of milk with impaired coagulation properties (noncoagulation). In total, individual coagulation properties were measured on morning milk collected from 1,299 Danish Holstein (DH), Danish Jersey (DJ), and Swedish Red (SR) cows. The 3 breeds demonstrated notable interbreed differences in coagulation properties, with DJ cows exhibiting superior coagulation compared with the other 2 breeds. In addition, milk samples from 2% of DH and 16% of SR cows were classified as noncoagulating. Furthermore, the cows were genotyped for major genetic variants in the αS1- (CSN1S1), β- (CSN2), and κ-casein (CSN3) genes, revealing distinct differences in variant frequencies among breeds. Allele I of CSN2, which had not formerly been screened in such a high number of cows in these Scandinavian breeds, showed a frequency around 7% in DH and DJ, but was not detected in SR. Genetic polymorphisms were significantly associated with curd firming rate and rennet coagulation time. Thus, CSN1S1 C, CSN2 B, and CSN3 B positively affected milk coagulation, whereas CSN2 A(2), in particular, had a negative effect. In addition to the influence of individual casein genes, the effects of CSN1S1-CSN2-CSN3 composite genotypes were also examined, and revealed strong associations in all breeds, which more or less reflected the single gene results. Overall, milk coagulation is under the influence of additive genetic variation. Optimal milk for future cheese production can be ensured by monitoring the frequency of unfavorable variants and thus preventing an increase in the number of cows producing milk with impaired coagulation. Selective breeding for variants associated with superior milk coagulation can potentially increase raw milk quality and cheese yield in all 3 Scandinavian breeds.

Newly identified mutations at the CSN1S1 gene in Ethiopian goats affect casein content and coagulation properties of their milk

Very high casein content and good coagulation properties previously observed in some Ethiopian goat breeds led to investigating the αs1-casein (CSN1S1) gene in these breeds. Selected regions of the CSN1S1 gene were sequenced in 115 goats from 5 breeds (2 indigenous: Arsi-Bale and Somali, 1 exotic: Boer, and 2 crossbreeds: Boer × Arsi-Bale and Boer × Somali). The DNA analysis resulted in 35 new mutations: 3 in exons, 3 in the 5' untranslated region (UTR), and 29 in the introns. The mutations in exons that resulted in an amino acid shift were then picked to evaluate their influence on individual casein content (αs1-, αs2-, β-, and κ-CN), micellar size, and coagulation properties in the milk from the 5 goat breeds. A mutation at nucleotide 10657 (exon 10) involved a transversion: CAG→CCG, resulting in an amino acid exchange Gln77→Pro77. This mutation was associated with the indigenous breeds only. Two new mutations, at nucleotide 6072 (exon 4) and 12165 (exon 12), revealed synonymous transitions: GTC→GTT in Val15 and AGA→AGG in Arg100 of the mature protein. Transitions G→A and C→T at nucleotides 1374 and 1866, respectively, occurred in the 5' UTR, whereas the third mutation involved a transversion T→G at nucleotide location 1592. The goats were grouped into homozygote new (CC), homozygote reference (AA), and heterozygote (CA) based on the nucleotide that involved the transversion. The content of αs1-CN (15.32g/kg) in milk samples of goats homozygous (CC) for this newly identified mutation, Gln77→Pro77 was significantly higher than in milks of heterozygous (CA; 9.05g/kg) and reference (AA; 7.61g/kg) genotype animals. The αs2-, β-, and κ-CN contents showed a similar pattern. Milk from goats with a homozygous new mutation had significantly lower micellar size. Milk from both homozygote and heterozygote new-mutation goats had significantly shorter coagulation rate and stronger gel than the reference genotype. Except the transversion, the sequence corresponded to allele A and presumably derived from it. Therefore, this allele is denoted by A3. All goats from the reference genotype (AA) were homozygous for the allele at nucleotide position 1374 and 1866, whereas all mutations in the 5' UTR existed in a heterozygous form in both heterozygous (CA) and the new mutation (CC) genotype. The newly identified mutation (CC) detected in some of the goat breeds is, therefore, important in selection for genetic improvement and high-quality milk for the emerging goat cheese-producing industries. The finding will also benefit farmers raising these goat breeds due to the increased selling price of goats. Further studies should investigate the effect of this amino acid exchange on the secondary and tertiary structure of the αs1-CN molecule and on the susceptibility of peptide hydrolysis by digestive enzymes.

Short communication: A new bovine milk-protein variant: α-lactalbumin variant D

Capillary zone electrophoresis of 1,948 Holstein-Friesian cows suggested the presence of an unknown protein variant of α-lactalbumin (α-LA) in the milk of 1 cow. Sequencing genomic DNA of this cow showed a polymorphism in the α-LA gene (LAA) that appears to be responsible for this protein variant. This single nucleotide polymorphism g.600G > T was located in exon 2 of LAA and causes the amino acid change 65Gln > His in the α-LA protein. This α-LA protein variant is a new protein variant and should be called α-LA protein variant D. This amino acid change is not expected to affect protein function. Genomic DNA of 156 bulls of various dairy cattle breeds was screened to examine the presence of the new α-LA protein variant D. Single nucleotide polymorphism g.600G > T, responsible for α-LA protein variant D, was not found in any of the 156 bulls. However, 10 other polymorphisms in the coding and promoter regions of LAA were detected that were used to construct haplotypes.

Protein degradation in bovine milk caused byStreptococcus agalactiae

Streptococcus(Str.)agalactiaeis a contagious mastitis bacterium, often associated with cases of subclinical mastitis. Different mastitis bacteria have been evaluated previously from a diagnostic point of view, but there is a lack of knowledge concerning their effect on milk composition. Protein composition is important in achieving optimal yield and texture when milk is processed to fermented products, such as cheese and yoghurt, and is thus of great economic value. The aim of thisin vitrostudy was to evaluate protein degradation mainly caused by exogenous proteases originating from naturally occurringStr. agalactiae. The samples were incubated at 37°C to imitate degradation caused by the bacteria in the udder. Protein degradation caused by different strains ofStr. agalactiaewas also investigated. Protein degradation was observed to occur whenStr. agalactiaewas added to milk, but there were variations between strains of the bacteria. Caseins, the most economically important proteins in milk, were degraded up to 75% in milk inoculated withStr. agalactiaein relation to sterile ultra-high temperature (UHT) milk, used as control milk. The major whey proteins, α-lactalbumin and β-lactoglobulin, were degraded up to 21% in relation to the sterile control milk. These results suggest that different mastitis bacteria but also different strains of mastitis bacteria should be evaluated from a milk quality perspective to gain knowledge about their ability to degrade the economically important proteins in milk.