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

FAM20C: A key protein kinase in multiple diseases

Family with sequence similarity 20 C (FAM20C) is a Golgi protein kinase that phosphorylates the serine residue in the S-x-E/pS motif of target proteins. FAM20C phosphorylates most secreted proteins, which play important roles in multiple biological processes, including cancer progression, biomineralization, and lipid homeostasis. Numerous studies have documented the potential contribution of FAM20C to the growth, invasion, and metastasis of glioma, breast cancer, and other cancers, as well as to the mineralization process of teeth and bone. In addition, FAM20C has been found to be associated with the occurrence and development of certain cardiovascular diseases and endocrine metabolism disorders. It raises hopes that understanding the disease-specific mechanisms of FAM20C may hold the key to developing new strategies for these diseases. This review comprehensively covers the existing literature to provide a summary of the structure and biological functions of FAM20C, with a particular focus on its roles in the disease context.

Profiling the differential phosphoproteome between breast milk and infant formula through a titanium (IV)-immobilized magnetic nanoplatform

Breast milk (BM) fulfills the nutritional needs of infants and sets the standard for infant formula (IF). However, profiling the differential phosphoproteome between BM and IF remains unclear. Herein, a titanium (IV) (Ti4+)-immobilized magnetic nanoplatform (Fe3O4@GO@PDA-Ti4+) was constructed by self-assembly polymerization of dopamine on magnetic graphene oxide, followed by immobilizing Ti4+ through chelation for phosphopeptide enrichment. Fe3O4@GO@PDA-Ti4+ possessed outstanding selectivity (1/1000, a molar ratio of β-casein digests to bovine serum albumin digests) and favorable sensitivity (2.5 fmol/μL), along with rapid magnetic separation. Excellent phosphopeptide capture efficiencies were obtained for BM and IF using Fe3O4@GO@PDA-Ti4+ as an adsorbent coupled with liquid chromatography-mass spectrometry/mass spectrometry. There were 191 and 239 phosphopeptides found in BM and IF, respectively, with 36 phosphoproteins identified in both. However, BM and IF shared only 17 phosphopeptides and 4 phosphoproteins. The variation in the phosphoproteome between BM and IF provides valuable insights into the optimization of IF humanization.

Revealing the diversity of endogenous peptides and parent proteins in human colostrum and mature milk through peptidomics analysis

Endogenous peptides and their parent proteins are important nutritional components with diverse biological functions. The objective of this study was to analyze and compare endogenous peptides and parent proteins found in human colostrum (HC) and human mature milk (HM) using a 4D label-free technique. In total, 5162 and 940 endogenous peptides derived from 258 parent proteins were identified in human milk by database (DB) search and de novo, respectively. Among these peptides, 2446 differentially expressed endogenous peptides with various bioactivities were identified. The Gene Ontology analysis unveiled the cellular components, biological processes, and molecular functions associated with these parent proteins. Metabolic pathway analysis suggested that neutrophil extracellular trap formation had the greatest significance with 24 parent proteins. These findings will offer a fresh perspective on the development of infant formula powder, highlighting the potential for incorporating these changes to enhance its nutritional composition and benefits.

Sequencing and Characterization of αs2-Casein Gene (CSN1S2) in the Old-World Camels Have Proven Genetic Variations Useful for the Understanding of Species Diversification

The CSN1S2 gene encodes αs2-casein, the third most abundant protein in camel milk. Despite its importance in foals, human nutrition, and dairy processing, the CSN1S2 gene in camels has received little attention. This study presents the first complete characterization of the CSN1S2 gene sequence in Old-World camels (Camelus bactrianus and Camelus dromedarius). Additionally, the gene promoter, consisting of 752 bp upstream of exon 1, was analyzed. The entire gene comprises 17 exons, ranging in length from 24 bp (exons 4, 8, 11, and 13) to 280 bp (exon 17). Interesting was the identification of the exon 12 in both species. The promoter analysis revealed 24 putative binding sites in the Bactrian camel and 22 in dromedary camel. Most of these sites were typical elements associated with milk protein, such as C/EBP-α, C/EBP-β, Oct-1, and AP1. The SNP discovery showed relatively high genetic diversity compared to other camel casein genes (CSN1S1, CSN2, and CSN3), with a total of 34 polymorphic sites across the two species. Particularly noteworthy is the transition g.311G>A in the CSN1S2 promoter, creating a new putative consensus binding site for a C/EBP-β in the Bactrian camel. At the exon level, two novel variants were found. One was detected in exon 6 of the Bactrian camel (g.3639C>G), resulting in an amino acid replacement, p.36Ile>Met. The second variant was found in noncoding exon 17 of dromedary CSN1S2 (g.1511G>T). Although this mutation occurs in the 3'-UnTranslated Region, it represents the first example of exonic polymorphism in the CSN1S2 for this species. This SNP also affects the binding sites of different microRNAs, including the seed sequence of the miRNA 4662a-3p, highlighting its role as a regulatory factor for CSN1S2 gene. A PCR-RFLP was set up for genotyping a dromedary Tunisian population (n = 157), and the minor allele frequency was found to be 0.27 for the G allele, indicating a potential yield improvement margin. The interspersed elements (INEs) analysis revealed 10 INEs covering 7.34% and 8.14% of the CSN1S2 sequence in the Bactrian and dromedary camels, respectively. Furthermore, six elements (A, B, F, H, I, and L) are shared among cattle and camels and are partially found in other ruminants, suggesting a common ancestral origin of these retrotransposons. Conversely, elements C, D, E, and G are specific to camels.

Milk Osteopontin and Human Health

Osteopontin (OPN) is a multifunctional protein found in all vertebrates. OPN is expressed in many different cell types, and is consequently found in most tissues and physiological secretions. OPN is involved in a multitude of biological processes, such as activation and regulation of the immune system; biomineralization; tissue-transformative processes, including growth and development of the gut and brain; interaction with bacteria; and many more. OPN is found in the highest concentrations in milk, where it is believed to initiate and regulate developmental, immunological and physiological processes in infants who consume milk. Processes for the isolation of bovine OPN for use in infant formula have been developed, and in recent years, many studies have investigated the effects of the intake of milk OPN. The purpose of this article is to review and compare existing knowledge about the structure and function of milk OPN, with a particular focus on the effects of milk OPN on human health and disease.

Caseins: Versatility of Their Micellar Organization in Relation to the Functional and Nutritional Properties of Milk

The milk of mammals is a complex fluid mixture of various proteins, minerals, lipids, and other micronutrients that play a critical role in providing nutrition and immunity to newborns. Casein proteins together with calcium phosphate form large colloidal particles, called casein micelles. Caseins and their micelles have received great scientific interest, but their versatility and role in the functional and nutritional properties of milk from different animal species are not fully understood. Caseins belong to a class of proteins that exhibit open and flexible conformations. Here, we discuss the key features that maintain the structures of the protein sequences in four selected animal species: cow, camel, human, and African elephant. The primary sequences of these proteins and their posttranslational modifications (phosphorylation and glycosylation) that determine their secondary structures have distinctively evolved in these different animal species, leading to differences in their structural, functional, and nutritional properties. The variability in the structures of milk caseins influence the properties of their dairy products, such as cheese and yogurt, as well as their digestibility and allergic properties. Such differences are beneficial to the development of different functionally improved casein molecules with variable biological and industrial utilities.

Probing of alpha, beta and kappa-caseins polymorphic variants in Gangatiri cow milk with the use of polyacrylamide gel electrophoresis and HRAMS

This research communication aimed to probe the genetic polymorphisms of alpha, beta, and kappa caseins in Gangatiri cows (an indigenous Indian cattle). Detection of variants has received considerable research interest in the dairy industry and animal breeding in recent years as a source of good quality protein, but also of bioactive peptides that may be linked to health implications. The polymorphic nature of casein fractions and their association with milk production traits, composition, and quality also attracted several efforts in evaluating the allelic distribution of different casein locus as a potential dairy trait marker. Molecular techniques of polyacrylamide gel electrophoresis and high-resolution accurate mass-spectrometry have been applied to this probe. Sequence analysis of different casein types in the cows showed the presence of four specific variants.

Quantitative Phosphoproteome of Infant Formula: New Insights into the Difference of Phosphorylation in Milk Proteins between Bovine and Goat Species

Phosphorylation is a broad post-translational protein modification, and the level of phosphorylation of milk proteins is associated with lactation, coagulation properties, and digestibility. However, phosphoproteins in bovine milk-based and goat milk-based infant formula have not been systematically explored. Here, we have analyzed six bovine and six goat milk-based infant formula using a quantitative phosphoproteomics approach, from which we identified 200 phosphoproteins with 276 phosphorylation sites and 156 phosphorylation sites from 75 phosphoproteins, respectively. Of these, 99 phosphorylation sites from 26 shared phosphoproteins were differentially expressed between bovine and goat milk-based infant formula. Especially, CSN1S1 was the most phosphoprotein with 25 quantified phosphorylation sites. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that the identified phosphoproteins not only provide nutrition to the infant but also have anti-inflammatory, antipathogenic, and other biological functions. Our results shed light on the composition, phosphorylation sites, and biological functions of phosphoproteins in bovine milk and goat milk-based infant formula, which provide new insights into the key role of protein modifications during infant development. It also helps us to better understand the differences in digestibility of infant formula from different animal milk sources and thus guides the choice of milk source for infant formula.

Milk proteins as a feed restriction signature indicating the metabolic adaptation of dairy cows

Milk production in dairy cows is affected by numerous factors, including diet. Feed restriction is known to have little impact on milk total protein content but its effect on the fine protein composition is still poorly documented. The objective of this study was to describe the effects of two feed restriction trials of different intensities on the milk protein composition of Holstein cows. One restriction trial was of high intensity (H: 8 mid-lactation Holstein cows) and the second of moderate intensity (M: 19 peak lactation Holstein cows). Feed restriction decreased the milk protein yield for caseins under the M trial and of all six major milk proteins under the H trial. These decreased yields lead to lower concentrations of αs1-, αs2- and β-caseins during the H trial. The milk proteome, analyzed on 32 milk samples, was affected as a function of restriction intensity. Among the 345 proteins identified eight varied under the M trial and 160 under the H trial. Ontology analyses revealed their implication in carbohydrate, lipid and protein metabolisms as well as in the immune system. These proteins reflected adaptations of the animal and mammary gland physiology to feed restriction and constituted a signature of this change.

Electrophoresis of Phosphoproteins on a Tandem Polymerized Gel

A substrate with n phosphorylated sites may have 2ⁿ phosphor-forms for temporal-spatial regulation of biological events. Because phosphates do not significantly change molecular masses but net charges of proteins, those isoforms cannot be separated by regular mass-based sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS PAGE). A tandem polymerized gel was developed to resolve phosphor-isoforms with different masses, charges, and posttranslational modifications. Without the usage of SDS, the electrophoresis was primarily performed on three adjacent acidic polyacrylamide gels. After being concentrated on a stacking gel, protonated proteins were then separated on the Zr⁴⁺ immobilized gel through the coordination of metal ions with phosphates followed by further charge and mass (z/m)-based electrophoretic separation on a TiO₂ containing gel. The presence of TiO₂ nanoparticles in the third gel is aimed for the initiation of the polymerization of acrylamide in acidic conditions upon ultraviolet irradiation. Distinct isoforms of α-S1-casein, α-S2-casein, β-casein, and κ casein model proteins located on 11, 8, 8, and 7 different bands of the tandem gel were unambiguously identified, respectively. With the tandem polymerized gel electrophoresis, new phosphorylation events that may occur simultaneously or sequentially were discovered in not only model proteins but also complex biological samples including human saliva, chicken egg, and sprouting maize. This provides a new tool to dissect complex biological processes that are triggered by dynamic phosphorylation events.

Selecting processing robust markers using high resolution mass spectrometry for the detection of milk in food products

BACKGROUND

Cow's milk allergy is one of the most reported food allergies in Europe. To help patients suffering from food allergies it is important to be able to detect milk in different foods. An analytical method that is gaining interest in the field of allergen detection is Ultra-High Performance Liquid Chromatography-tandem Mass Spectrometry, where the analyte is a target peptide. When these peptide biomarkers are selected the effect of food processing should be taken into account to allow a robust detection method.

OBJECTIVE

This works aims at identifying such processing stable peptide markers for milk for the Ultra-High Performance Liquid Chromatography-tandem Mass Spectrometry based detection of food allergens in different food products.

METHODS

Milk-incurred food materials that underwent several processing techniques were produced. This was followed by establishing tryptic peptide profiles from each matrix using Ultra-High Performance Liquid Chromatography-High Resolution Mass Spectrometry .

RESULT

A careful comparison of peptide profiles/intensities and the use of specific exclusion criteria resulted in the selection of 8 peptide biomarkers suitable for application in Ultra-High Performance Liquid Chromatography-tandem Mass Spectrometry based milk detection methods. One of these markers is a α-lactalbumin specific peptide, which has been determined to be stable in different incurred materials for the first time.

CONCLUSION

To our knowledge, this is the first systematic and experimentally based approach for the selection of suitable milk peptide biomarkers robust towards multiple, often applied food processing techniques for milk. Ensuring the exact knowledge of the food processing circumstances by starting from well-defined raw material and using fully controlled settings to produce incurred test material allowed the construction of a peptide database with robust markers. These robust markers can be used for the development of a robust detection method for milk in different food matrices.

Partial-, Double-Enzymatic Dephosphorylation and EndoGluC Hydrolysis as an Original Approach to Enhancing Identification of Casein Phosphopeptides (CPPs) by Mass Spectrometry

The identification of phosphopeptides is currently a challenge when they are part of a complex matrix of peptides, such as a milk protein enzymatic hydrolysate. This challenge increases with both the number of phosphorylation sites on the phosphopeptides and their amino acid length. Here, this paper reports a four-phase strategy from an enzymatic casein hydrolysate before a mass spectrometry analysis in order to enhance the identification of phosphopeptides and phosphosites: (i) the control protein hydrolysate, (ii) a two-step enzymatic dephosphorylation of the latter, allowing for the almost total dephosphorylation of peptides, (iii) a one-step enzymatic dephosphorylation, allowing for the partial dephosphorylation of the peptides and (iv) an additional endoGluC enzymatic hydrolysis, allowing for the cleavage of long-size peptides into shorter ones. The reverse-phase high-pressure liquid chromatography–tandem mass spectrometry (RP-HPLC-MS/MS) analyses of hydrolysates that underwent this four-phase strategy allowed for the identification of 28 phosphorylation sites (90%) out of the 31 referenced in UniprotKB/Swiss-Prot (1 June 2021), compared to 17 sites (54%) without the latter. The alpha-S2 casein phosphosites, referenced by their similarity in the UniProt database, were experimentally identified, whereas pSer148, pThr166 and pSer187 from a multiphosphorylated long-size kappa-casein were not. Data are available via ProteomeXchange with identifier PXD027132.

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.

The ABCs of the atypical Fam20 secretory pathway kinases

The study of extracellular phosphorylation was initiated in late 19th century when the secreted milk protein, casein, and egg-yolk protein, phosvitin, were shown to be phosphorylated. However, it took more than a century to identify Fam20C, which phosphorylates both casein and phosvitin under physiological conditions. This kinase, along with its family members Fam20A and Fam20B, defined a new family with altered amino acid sequences highly atypical from the canonical 540 kinases comprising the kinome. Fam20B is a glycan kinase that phosphorylates xylose residues and triggers peptidoglycan biosynthesis, a role conserved from sponges to human. The protein kinase, Fam20C, conserved from nematodes to humans, phosphorylates well over 100 substrates in the secretory pathway with overall functions postulated to encompass endoplasmic reticulum homeostasis, nutrition, cardiac function, coagulation, and biomineralization. The preferred phosphorylation motif of Fam20C is SxE/pS, and structural studies revealed that related member Fam20A allosterically activates Fam20C by forming a heterodimeric/tetrameric complex. Fam20A, a pseudokinase, is observed only in vertebrates. Loss-of-function genetic alterations in the Fam20 family lead to human diseases such as amelogenesis imperfecta, nephrocalcinosis, lethal and nonlethal forms of Raine syndrome with major skeletal defects, and altered phosphate homeostasis. Together, these three members of the Fam20 family modulate a diverse network of secretory pathway components playing crucial roles in health and disease. The overarching theme of this review is to highlight the progress that has been made in the emerging field of extracellular phosphorylation and the key roles secretory pathway kinases play in an ever-expanding number of cellular processes.

Identification of Genomic Regions Associated with Concentrations of Milk Fat, Protein, Urea and Efficiency of Crude Protein Utilization in Grazing Dairy Cows

The objective of this study was to identify genomic regions associated with milk fat percentage (FP), crude protein percentage (CPP), urea concentration (MU) and efficiency of crude protein utilization (ECPU: ratio between crude protein yield in milk and dietary crude protein intake) using grazing, mixed-breed, dairy cows in New Zealand. Phenotypes from 634 Holstein Friesian, Jersey or crossbred cows were obtained from two herds at Massey University. A subset of 490 of these cows was genotyped using Bovine Illumina 50K SNP-chips. Two genome-wise association approaches were used, a single-locus model fitted to data from 490 cows and a single-step Bayes C model fitted to data from all 634 cows. The single-locus analysis was performed with the Efficient Mixed-Model Association eXpedited model as implemented in the SVS package. Single nucleotide polymorphisms (SNPs) with genome-wide association p-values ≤ 1.11 × 10⁻⁶ were considered as putative quantitative trait loci (QTL). The Bayes C analysis was performed with the JWAS package and 1-Mb genomic windows containing SNPs that explained > 0.37% of the genetic variance were considered as putative QTL. Candidate genes within 100 kb from the identified SNPs in single-locus GWAS or the 1-Mb windows were identified using gene ontology, as implemented in the Ensembl Genome Browser. The genes detected in association with FP (MGST1, DGAT1, CEBPD, SLC52A2, GPAT4, and ACOX3) and CPP (DGAT1, CSN1S1, GOSR2, HERC6, and IGF1R) were identified as candidates. Gene ontology revealed six novel candidate genes (GMDS, E2F7, SIAH1, SLC24A4, LGMN, and ASS1) significantly associated with MU whose functions were in protein catabolism, urea cycle, ion transportation and N excretion. One novel candidate gene was identified in association with ECPU (MAP3K1) that is involved in post-transcriptional modification of proteins. The findings should be validated using a larger population of New Zealand grazing dairy cows.

Post translational modifications of milk proteins in geographically diverse goat breeds

Goat milk is a source of nutrition in difficult areas and has lesser allerginicity than cow milk. It is leading in the area for nutraceutical formulation and drug development using goat mammary gland as a bioreactor. Post translational modifications of a protein regulate protein function, biological activity, stabilization and interactions. The protein variants of goat milk from 10 breeds were studied for the post translational modifications by combining highly sensitive 2DE and Q-Exactive LC-MS/MS. Here we observed high levels of post translational modifications in 201 peptides of 120 goat milk proteins. The phosphosites observed for CSN2, CSN1S1, CSN1S2, CSN3 were 11P, 13P, 17P and 6P, respectively in 105 casein phosphopeptides. Whey proteins BLG and LALBA showed 19 and 4 phosphosites respectively. Post translational modification was observed in 45 low abundant non-casein milk proteins mainly associated with signal transduction, immune system, developmental biology and metabolism pathways. P_asp is reported for the first time in 47 sites. The rare conserved peptide sequence of (SSSEE) was observed in αS1 and αS2 casein. The functional roles of identified phosphopeptides included anti-microbial, DPP-IV inhibitory, anti-inflammatory and ACE inhibitory. This is first report from tropics, investigating post translational modifications in casein and non-casein goat milk proteins and studies their interactions.

LC-MS Quantification of Site-Specific Phosphorylation Degree by Stable-Isotope Dimethyl Labeling Coupled with Phosphatase Dephosphorylation

Protein phosphorylation is a crucial post-translational modification that plays an important role in the regulation of cellular signaling processes. Site-specific quantitation of phosphorylation levels can help decipher the physiological functions of phosphorylation modifications under diverse physiological statuses. However, quantitative analysis of protein phosphorylation degrees is still a challenging task due to its dynamic nature and the lack of an internal standard simultaneously available for the samples differently prepared for various phosphorylation extents. In this study, stable-isotope dimethyl labeling coupled with phosphatase dephosphorylation (DM + deP) was tried to determine the site-specific degrees of phosphorylation in proteins. Firstly, quantitation accuracy of the (DM + deP) approach was confirmed using synthetic peptides of various simulated phosphorylation degrees. Afterwards, it was applied to evaluate the phosphorylation stoichiometry of milk caseins. The phosphorylation degree of Ser130 on α-S1-casein was also validated by absolute quantification with the corresponding synthetic phosphorylated and nonphosphorylated peptides under a selected reaction monitoring (SRM) mode. Moreover, this (DM + deP) method was used to detect the phosphorylation degree change of Ser82 on the Hsp27 protein of HepG2 cells caused by tert-butyl hydroperoxide (t-BHP) treatment. The results showed that the absolute phosphorylation degree obtained from the (DM + deP) approach was comparable with the relative quantitation resulting from stable-isotope dimethyl labeling coupled with TiO₂ enrichment. This study suggested that the (DM + deP) approach is promising for absolute quantification of site-specific degrees of phosphorylation in proteins, and it may provide more convincing information than the relative quantification method.

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.

Comparison of Bayesian and partial least squares regression methods for mid-infrared prediction of cheese-making properties in Montbéliarde cows

Assessing the cheese-making properties (CMP) of milks with a rapid and cost-effective method is of particular interest for the Protected Designation of Origin cheese sector. The aims of this study were to evaluate the potential of mid-infrared (MIR) spectra to estimate coagulation and acidification properties, as well as curd yield (CY) traits of Montbéliarde cow milk. Samples from 250 cows were collected in 216 commercial herds in Franche-Comté with the objectives to maximize the genetic diversity as well as the variation in milk composition. All coagulation and CY traits showed high variability (10 to 43%). Reference analyses performed for soft (SC) and pressed cooked (PCC) cheese technology were matched with MIR spectra. Prediction models were built on 446 informative wavelengths not tainted by the water absorbance, using different approaches such as partial least squares (PLS), uninformative variable elimination PLS, random forest PLS, Bayes A, Bayes B, Bayes C, and Bayes RR. We assessed equation performances for a set of 20 CMP traits (coagulation: 5 for SC and 4 for PCC; acidification: 5 for SC and 3 for PCC; laboratory CY: 3) by comparing prediction accuracies based on cross-validation. Overall, variable selection before PLS did not significantly improve the performances of the PLS regression, the prediction differences between Bayesian methods were negligible, and PLS models always outperformed Bayesian models. This was likely a result of the prior use of informative wavelengths of the MIR spectra. The best accuracies were obtained for curd yields expressed in dry matter (CY_DM) or fresh (CY_FRESH) and for coagulation traits (curd firmness for PCC and SC) using the PLS regression. Prediction models of other CMP traits were moderately to poorly accurate. Whatever the prediction methodology, the best results were always obtained for CY traits, probably because these traits are closely related to milk composition. The CY_DM predictions showed coefficient of determination (R²) values up to 0.92 and 0.87, and RS_y,x values of 3 and 4% for PLS and Bayes regressions, respectively. Finally, we divided the data set into calibration (2/3) and validation (1/3) sets and developed prediction models in external validation using PLS regression only. In conclusion, we confirmed, in the validation set, an excellent prediction for CY_DM [R² = 0.91, ratio of performance to deviation (RPD) = 3.39] and a very good prediction for CY_FRESH (R² = 0.84, RPD = 2.49), adequate for analytical purposes. We also obtained good results for both PCC and SC curd firmness traits (R² ≥ 0.70, RPD ≥1.8), which enable quantitative prediction.

Alternative splicing events expand molecular diversity of camel CSN1S2 increasing its ability to generate potentially bioactive peptides

In a previous study on camel milk from Kazakhstan, we reported the occurrence of two unknown proteins (UP1 and UP2) with different levels of phosphorylation. Here we show that UP1 and UP2 are isoforms of camel αs2-CN (αs2-CNsv1 and αs2-CNsv2, respectively) arising from alternative splicing events. First described as a 178 amino-acids long protein carrying eight phosphate groups, the major camel αs2-CN isoform (called here αs2-CN) has a molecular mass of 21,906 Da. αs2-CNsv1, a rather frequent (35%) isoform displaying a higher molecular mass (+1,033 Da), is present at four phosphorylation levels (8P to 11P). Using cDNA-sequencing, αs2-CNsv1 was shown to be a variant arising from the splicing-in of an in-frame 27-nucleotide sequence encoding the nonapeptide ENSKKTVDM, for which the presence at the genome level was confirmed. αs2-CNsv2, which appeared to be present at 8P to 12P, was shown to include an additional decapeptide (VKAYQIIPNL) revealed by LC-MS/MS, encoded by a 3'-extension of exon 16. Since milk proteins represent a reservoir of biologically active peptides, the molecular diversity generated by differential splicing might increase its content. To evaluate this possibility, we searched for bioactive peptides encrypted in the different camel αs2-CN isoforms, using an in silico approach. Several peptides, putatively released from the C-terminal part of camel αs2-CN isoforms after in silico digestion by proteases from the digestive tract, were predicted to display anti-bacterial and antihypertensive activities.

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.

Combining different proteomic approaches to resolve complexity of the milk protein fraction of dromedary, Bactrian camels and hybrids, from different regions of Kazakhstan

Nutritional suitability of milk is not only related to gross composition, but is also strongly affected by the microheterogeniety of the protein fraction. Hence, to go further into the evaluation of the potential suitability of non-bovine milks in human/infant nutrition it is necessary to have a detailed characterization of their protein components. Combining proven proteomic approaches (SDS-PAGE, LC-MS/MS and LC-ESI-MS) and cDNA sequencing, we provide here in depth characterization of the milk protein fraction of dromedary and Bactrian camels, and their hybrids, from different regions of Kazakhstan. A total 391 functional groups of proteins were identified from 8 camel milk samples. A detailed characterization of 50 protein molecules, relating to genetic variants and isoforms arising from post-translational modifications and alternative splicing events, belonging to nine protein families (κ-, α_s1-, α_s2-, β-; and γ-CN, WAP, α-LAC, PGRP, CSA/LPO) was achieved by LC-ESI-MS. The presence of two unknown proteins UP1 (22,939 Da) and UP2 (23,046 Da) was also reported as well as the existence of a β-CN short isoform (946 Da lighter than the full-length β-CN), arising very likely in both genetic variants (A and B) from proteolysis by plasmin. In addition, we report, for the first time to our knowledge, the occurrence of a α_s2-CN phosphorylation isoform with 12P groups within two recognition motifs, suggesting thereby the existence of two kinase systems involved in the phosphorylation of caseins in the mammary gland. Finally, we demonstrate that genetic variants, which hitherto seemed to be species- specific (e.g. β-CN A for Bactrian and β-CN B for dromedary), are in fact present both in Camel dromedarius and C. bactrianus.

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.

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.