Effects of genetic variants and sialylation on in vitro digestibility of purified κ-casein.

Milk with different κ-casein (CN) phenotypes has previously been found to influence its gastric digestion rate. Therefore, the aim of the present study is to disentangle contributions of genetic variation and its related sialylation on the in vitro digestion process of κ-CN. Accordingly, κ-CN was purified from milk representing homozygous cows with κ-CN phenotypes AA, BB, or EE and used as substrate molecules in model studies using the INFOGEST 2.0 in vitro static digestion model. Furthermore, the effect of removal of the terminal sialic acids present on the O-linked oligosaccharides of the purified κ-CN A, B, and E protein variants were studied by desialylation enzymatic assays. The κ-CN proteins were purified by reducing anion exchange chromatography with purities of variants A, B, and E of 93.0, 97.1, and 90.0%, respectively. Protein degradations of native and desialylated κ-CN isolates in gastric and intestinal phases were investigated by sodium dodecyl sulfate-PAGE, degree of hydrolysis (DH), and liquid chromatography electrospray ionization mass spectrometry. It was shown that after purification, the κ-CN molecules reassembled into multimer states, which then constituted the basis for the digestion studies. As assessed by DH, purified variants A and E were found to exhibit faster in vitro digestion rates in both gastric and intestinal phases compared with variant B. Desialylation increased both gastric and intestinal digestion rates for all variants, as measured by DH. In the gastric phase, desialylation promoted digestion of variant B at a rate comparable with native variants A and E, whereas in the intestinal phase, desialylation of variant B promoted better digestion than native A or E. Taken together, the results confirm that low glycosylation degree of purified κ-CN promotes faster in vitro digestion rates, and that desialylation of the O-linked oligosaccharides further promotes digestion. This finding could be applied to produce dairy products with enhanced digestibility.

Phosphorylation and glycosylation isoforms of bovine κ-casein variant E in homozygous Swedish Red cow milk detected by liquid chromatography-electrospray ionization mass spectrometry.

Variations in the phosphorylation and glycosylation patterns of the common κ-casein (CN) variants A and B have been explored, whereas studies on variant E heterogeneity are scarce. This study reports for the first time the detailed phosphorylation and glycosylation pattern of the κ-CN variant E in comparison with variants A and B. Individual cow milk samples representing κ-CN genotype EE (n = 12) were obtained from Swedish Red cows, and the natural posttranslational modifications of its κ-CN were identified and quantified by liquid chromatography-electrospray mass spectrometry. In total, 12 unique isoform masses of κ-CN variant E were identified. In comparison, AA and BB milk consisted of 14 and 17 unique isoform masses, respectively. The most abundant κ-CN E isoform detected in the EE milk was the monophosphorylated, unglycosylated [1P 0G, ∼70%; where P indicates phosphorylation from single to triple phosphorylation (1-3P), and G indicates glycosylation from single to triple glycosylation (1-3G)] form, followed by diphosphorylated, unglycosylated (2P 0G, ∼12%) form, resembling known patterns from variants A and B. However, a clear distinction was the presence of the rare triphosphorylated, nonglycosylated (3P 0G, ∼0.05%) κ-CN isoform in the EE milk. All isoforms detected in variant E were phosphorylated, giving a phosphorylation degree of 100%. This is comparable with the phosphorylation degree of variants A and B, being also almost 100%, though with very small amounts of nonphosphorylated, glycosylated isoforms detected. The glycosylation degree of variant E was found to be around 17%, a bit higher than observed for variant B (around 14%), and higher than variant A (around 7%). Among glycosylation, the glycan e was the most common type identified for all 3 variants, followed by c/d (straight and branched chain trisaccharides, respectively), and b. In contrast to κ-CN variants A and B, no glycan of type a was found in variant E. Taken together, this study shows that the posttranslational modification pattern of variant E resembles that of known variants to a large extent, but with subtle differences.

Comparison of bovine milk oligosaccharides in native North European cattle breeds.

Milk oligosaccharides are of high interest due to their bioactive properties. This study is the first to characterise milk oligosaccharides from native North European cattle breeds, as represented by 80 milk samples collected from eight native breeds originated from Norway (Norwegian Doela cattle and Norwegian Telemark cattle), Sweden (Swedish Mountain cattle), Denmark (Danish Red anno 1970), Iceland (Icelandic cattle), Lithuania (native Lithuanian Black and White) and Finland (Western Finncattle and Eastern Finncattle). Using high-performance liquid-chromatography chip/quadrupole time-of-flight mass-spectrometry, 18 unique monosaccharide compositions and a multitude of isomers were identified. No N-glycolylneuraminic acid was identified among these breeds. Western Finncattle milk was most abundant in neutral, acidic and fucosylated oligosaccharides. Further, Eastern Finncattle milk was significantly higher in acidic oligosaccharides and Icelandic cattle milk significantly higher in fucosylated oligosaccharides, compared to the mean. This study highlights specific native breeds of particular interest for future exploitation of milk oligosaccharides and breeding strategies.

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 A2 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 A2, κ-CN E and A variants, as well as ß-LG B, and the predominant composite genotype of ß- and κ-CN in the NC group was A2A2/AA. Generally, the A2A2/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 A2A2/AA through selective breeding could have a positive effect on milk coagulation properties.

Association between the bovine milk metabolome and rennet-induced coagulation properties of milk.

The milk metabolomes of 407 individual Swedish Red dairy cows were analyzed by nuclear magnetic resonance spectroscopy as part of the Danish-Swedish Milk Genomics Initiative. By relating these metabolite profiles to total milk protein concentration and rheological measurements of rennet-induced milk coagulation together using multivariate data analysis techniques, we were able to identify several different associations of the milk metabolome to technological properties of milk. Several novel correlations of milk metabolites to protein content and rennet-induced coagulation properties were demonstrated. Metabolites associated with the prediction of total protein content included choline, N-acetyl hexosamines, creatinine, glycerophosphocholine, glutamate, glucose 1-phosphate, galactose 1-phosphate, and orotate. In addition, levels of lactate, acetate, glutamate, creatinine, choline, carnitine, galactose 1-phosphate, and glycerophosphocholine were significantly different when comparing noncoagulating and well-coagulating milks. These findings suggest that the mentioned metabolites are associated with milk protein content and rennet-induced coagulation properties and may act as quality markers for cheese milk.

Re-sequencing of the casein genes in Swedish Red cattle giving milk with diverse protein profiles and extreme rennet coagulation properties.

Impaired rennet coagulation properties in milk could lead to prolonged processing times and production losses. Heritability for milk coagulation has previously been estimated to be 0.28 to 0.45, indicating that genetic selection can be used to manipulate this trait. The CN proteins are expressed by the genes CSN1S1, CSN2, CSN1S2, and CSN3 and are located on bovine chromosome 6. To better understand the effect of genetic variation in the CN genes on milk coagulation, blood and milk samples from 30 Swedish Red Dairy Cattle (RDC) with divergent coagulation properties were investigated. DNA from the 30 cows was sequenced for the CN genes to determine the theoretical AA sequence and to look for genetic variation in the untranslated regions. The aim is to confirm the protein genetic variants previously reported, while searching for additional genetic variation in the CN genes of 30 RDC. We observed genetic variation in 116 SNPs in the known CN genes where 10% of the SNPs are exon variants and the remaining 90% are intron variants. A total of 2.5% of the SNPs are found in the 5'- or 3'-untranslated region (UTR) regions of the exons; 2% are synonymous variants and 6% are missense variants that concurred with the known protein variants for CSN1S1, CSN2, and CSN3. Furthermore, 6% of the SNPs are splice polypyrimidine tract intron variants. The 2 genetic variants in the 5'- and 3'-UTR in CSN1S1 and CSN3 are found with protein variants CSN1S1C and CSN3B. Because both UTR variants are associated with gain and loss of micro RNA and transcription factors, this could explain differences in expression of the genetic protein variants. Preliminary chi-squared analysis and comparison with previous GWAS studies showed potential connections between the identified SNPs and coagulation properties of milk. By advancing the knowledge of the connection between the DNA sequence and the functional properties of the CN proteins, we hope to learn more about the cheese coagulation properties of milk from RDC.

Genomic selection in relation to bovine milk composition and processability.

Genomic selection is a new technology in which selection decisions are based on direct genomic values (DGVs) or genomic enhanced breeding values (GEBVs). The objective of this study was to evaluate the relations between DGVs and several milk traits important for both the nutritional value and processability of milk. This is a new approach and can be used to increase the knowledge on how genomic selection can be used in practice. Morning milk samples from Swedish Holstein cows were analyzed for milk composition and technological properties. DGVs were received for each cow for milk, protein and fat yield, milk index, udder health, Nordic total merit and a quota was calculated between fat and milk yield as well as protein and milk yield. The results show that linear correlations exist (P<0·10) between the studied DGVs and contents and yields of parameters in the protein (P=0·002-0·097), fat (P=0·024-0·055) and mineral profiles (P=0·001-0·099) as well as for cheese characteristics (P=0·004-0·065), thus making it possible to obtain detailed information on milk traits that are not registered in the milk recording scheme. Hence, genomic selection will be an efficient tool for breeding and dairy industry to select cows early in life for targeted milk production.

Effect of polymorphisms in the leptin, leptin receptor and acyl-CoA:diacylglycerol acyltransferase 1 (DGAT1) genes and genetic polymorphism of milk proteins on bovine milk composition.

The relations between cow genetics and milk composition have gained a lot of attention during the past years, however, generally only a few compositional traits have been examined. The aim of this study was to determine if polymorphisms in the leptin (LEP), leptin receptor (LEPR) and acyl-CoA:diacylglycerol acyltransferase 1 (DGAT1) genes as well as genetic polymorphism of ß-casein (ß-CN), κ-CN and ß-lactoglobulin (ß-LG) impact several bovine milk composition traits. Individual milk samples from the Swedish Red and Swedish Holstein breeds were analyzed for components in the protein, lipid, carbohydrate and mineral profiles. Cow alleles were determined on the following SNP: A1457G, A252T, A59V and C963T on the LEP gene, T945M on the LEPR gene and Nt984+8(A-G) on the DGAT1 gene. Additionally, genetic variants of ß-CN, κ-CN and ß-LG were determined. For both the breeds, the same tendency of minor allele frequency was found for all SNPs and protein genes, except on LEPA1457G and LEPC963T. This study indicated significant (P<0·05) associations between the studied SNPs and several compositional parameters. Protein content was influenced by LEPA1457G (G>A) and LEPC963T (T>C), whereas total Ca, ionic Ca concentration and milk pH were affected by LEPA1457G, LEPA59V, LEPC963T and LEPRT945M. However, yields of milk, protein, CN, lactose, total Ca and P were mainly affected by ß-CN (A2>A1) and κ-CN (A>B>E). ß-LG was mainly associated with whey protein yield and ionic Ca concentration (A>B). Thus, this study shows possibilities of using these polymorphisms as markers within genetic selection programs to improve and adjust several compositional parameters.

Naturally Occurring Glycosidases in Milk from Native Cattle Breeds: Activity and Consequences on Free and Protein Bound-Glycans.

Little is known about the extent of variation and activity of naturally occurring milk glycosidases and their potential to degrade milk glycans. A multi-omics approach was used to investigate the relationship between glycosidases and important bioactive compounds such as free oligosaccharides and O-linked glycans in bovine milk. Using 4-methylumbelliferone (4-MU) assays activities of eight indigenous glycosidases were determined, and by mass spectrometry and 1H NMR spectroscopy various substrates and metabolite products were quantified in a subset of milk samples from eight native North European cattle breeds. The results showed a clear variation in glycosidase activities among the native breeds. Interestingly, negative correlations between some glycosidases including ß-galactosidase, N-acetyl-ß-d-glucosaminidase, certain oligosaccharide isomers as well as O-linked glycans of κ-casein were revealed. Further, a positive correlation was found for free fucose content and α-fucosidase activity (r = 0.37, p-value < 0.001) indicating cleavage of fucosylated glycans in milk at room temperature. The results obtained suggest that milk glycosidases might partially degrade valuable glycans, which would result in lower recovery of glycans and thus represent a loss for the dairy ingredients industry if these activities are pronounced.

Revealing the size, conformation, and shape of casein micelles and aggregates with asymmetrical flow field-flow fractionation and multiangle light scattering.

Casein (CN) micelles are naturally occurring colloidal protein aggregates present in a dispersed state in milk. In this paper we aim to obtain a detailed description of physicochemical properties of CN micelles over the entire size distribution using asymmetrical flow field-flow fractionation (AsFlFFF) connected to multiangle light scattering (MALS) and refractive index (RI) detection. Conclusions are drawn on the colloidal level regarding shape and conformation by comparison with models of colloidal particles. By using AsFlFFF-MALS-RI, it is concluded that the CN micelles are highly polydisperse with an average rms radius and hydrodynamic radius of 177 and 116 nm, respectively. The results show that the majority of CN micelles have a spherical shape, whereas a low concentration exists of larger and elongated aggregates. By comparison with models of aggregates of colloidal particles, the aggregates are shown to be anisotropic, e.g., aggregating linearly (threadlike) or in a sheet, rather than forming randomly spherical clusters. The results show that the characterization of colloidal dispersions with AsFlFFF-MALS-RI and the comparison with theoretical models are of a general character and, thus, of fundamental importance for colloidal dispersions.

Identification of QTL on Chromosome 18 Associated with Non-Coagulating Milk in Swedish Red Cows.

Non-coagulating (NC) milk, defined as milk not coagulating within 40 min after rennet-addition, can have a negative influence on cheese production. Its prevalence is estimated at 18% in the Swedish Red (SR) cow population. Our study aimed at identifying genomic regions and causal variants associated with NC milk in SR cows, by doing a GWAS using 777k SNP genotypes and using imputed sequences to fine map the most promising genomic region. Phenotypes were available from 382 SR cows belonging to 21 herds in the south of Sweden, from which individual morning milk was sampled. NC milk was treated as a binary trait, receiving a score of one in case of non-coagulation within 40 min. For all 382 SR cows, 777k SNP genotypes were available as well as the combined genotypes of the genetic variants of αs1-ß-κ-caseins. In addition, whole-genome sequences from the 1000 Bull Genome Consortium (Run 3) were available for 429 animals of 15 different breeds. From these sequences, 33 sequences belonged to SR and Finish Ayrshire bulls with a large impact in the SR cow population. Single-marker analyses were run in ASReml using an animal model. After fitting the casein loci, 14 associations at -Log10(P-value) > 6 identified a promising region located on BTA18. We imputed sequences to the 382 genotyped SR cows using Beagle 4 for half of BTA18, and ran a region-wide association study with imputed sequences. In a seven mega base-pairs region on BTA18, our strongest association with NC milk explained almost 34% of the genetic variation in NC milk. Since it is possible that multiple QTL are in strong LD in this region, 59 haplotypes were built, genetically differentiated by means of a phylogenetic tree, and tested in phenotype-genotype association studies. Haplotype analyses support the existence of one QTL underlying NC milk in SR cows. A candidate gene of interest is the VPS35 gene, for which one of our strongest association is an intron SNP in this gene. The VPS35 gene belongs to the mammary gene sets of pre-parturient and of lactating cows.