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

Chromosome-level dairy goat genome reveals the regulatory landscape of lactation

Goat milk is rich in various nutrients that are beneficial for human health. However, the genomic evolution and genetic basis underlying the nutritional value and unique flavor formation in dairy goats remain poorly understood. In the present study, we generate a chromosome-level genome assembly for dairy goats comprising 2.63 Gb with a contig N50 of 43 Mb and a scaffold N50 of 101 Mb. Genome quality comparisons revealed that the dairy goat genome has higher integrity and continuity than the published goat and sheep genomes. The identification of genes under positive selection in dairy goats highlights potential candidates to explain their high milk production. Comparative genomic analysis elucidates the adaptive evolutionary mechanisms of dairy goats such as strong disease resistance, broad adaptability, and unique milk flavor. Moreover, we demonstrate the conservation of the lactation gene network and identify new potential regulators associated with lipid metabolism. Additionally, we establish the regulatory landscape of lactation for the first time in dairy goats, revealing its unique gene regulatory characteristics. Hence, our study not only provides the first chromosome-level reference genome for dairy goat, but also offers potential research directions for dairy production and genetic improvement.

Spectral Profiling (Fourier Transform Infrared Spectroscopy) and Machine Learning for the Recognition of Milk from Different Bovine Breeds

The Podolica cattle breed is widespread in southern Italy, and its productivity is characterized by low yields and an extraordinary quality of milk and meats. Most of the milk produced is transformed into "Caciocavallo Podolico" cheese, which is made with 100% Podolica milk. Fourier Transform Infrared Spectroscopy (FTIR) is the technique that, in this research work, was applied together with machine learning to discriminate 100% Podolica milk from contamination of other Calabrian cattle breeds. The analysis on the test set produced a misclassification percentage of 6.7%. Among the 15 non-Podolica samples in the test set, 2 were misclassified and recognized as Podolica milk even though the milk was from other species. The correct classification rate improved to 100% when the same method was applied to the recognition of Podolica and Pezzata Rossa milk produced by the same farm. Furthermore, this technique was tested for the recognition of Podolica milk mixed with milk from other bovine species. The multivariate model and the respective confusion matrices obtained showed that all the 14 Podolica samples (test set) mixed with 40% non-Podolica milk were correctly classified. In addition, Pezzata Rossa milk produced by the same farm was detected as a contaminant in Podolica milk from the same farm down to concentrations as little as 5% with a 100% correct classification rate in the test set. The method described yielded higher accuracy values when applied to the discrimination of milks from different breeds belonging to the same farm. One of the reasons for this phenomenon could be linked to the elimination of the environmental variable. However, the results obtained in this work demonstrate the possibility of using FTIR to discriminate between milks from different breeds.

Effect of intramammary infection and inflammation on milk protein profile assessed at the quarter level in Holstein cows

In this study we wanted to investigate the associations between naturally occurring subclinical intramammary infection (IMI) caused by different etiological agents (i.e., Staphylococcus aureus, Streptococcus agalactiae, Streptococcus uberis, and Prototheca spp.), in combination with somatic cell count (SCC), on the detailed milk protein profile measured at the individual mammary gland quarter. An initial bacteriological screening (time 0; T0) conducted on individual composite milk from 450 Holstein cows reared in 3 herds, was performed to identify cows with subclinical IMI. We identified 78 infected animals which were followed up at the quarter level at 2 different sampling times: T1 and T2, 2 and 6 wk after T0, respectively. A total of 529 quarter samples belonging to the previously selected animals were collected at the 2 sampling points and analyzed with a reversed phase HPLC (RP-HPLC) validated method. Specifically, we identified and quantified 4 caseins (CN), namely αS1-CN, αS2-CN, κ-CN, and β-CN, and 3 whey protein fractions, namely β-lactoglobulin, α-lactalbumin, and lactoferrin (LF), which were later expressed both quantitatively (g/L) and qualitatively (as a percentage of the total milk nitrogen content, % N). Data were analyzed with a hierarchical linear mixed model with the following fixed effects: days in milk (DIM), parity, herd, SCC, bacteriological status (BACT), and the SCC × BACT interaction. The random effect of individual cow, nested within herd, DIM and parity was used as the error term for the latter effects. Both IMI (i.e., BACT) and SCC significantly reduced the proportion of β-CN and αS1-CN, ascribed to the increased activity of both milk endogenous and microbial proteases. Less evident alterations were found for whey proteins, except for LF, which being a glycoprotein with direct and undirect antimicrobial activity, increased both with IMI and SCC, suggesting its involvement in the modulation of both the innate and adaptive immune response. Finally, increasing SCC in the positive samples was associated with a more marked reduction of total caseins at T1, and αS1-CN at T2, suggesting a synergic effect of infection and inflammation, more evident at high SCC. In conclusion, our work helps clarify the behavior of protein fractions at quarter level in animals having subclinical IMI. The inflammation status driven by the increase in SCC, rather the infection, was associated with the most significant changes, suggesting that the activity of endogenous proteolytic enzymes related to the onset of inflammation might have a pivotal role in directing the alteration of the milk protein profile.

CSN1S1 and LALBA Polymorphisms and Other Factors Influencing Yield, Composition, Somatic Cell Score, and Technological Properties of Cow's Milk

We evaluated the influence of CSN1S1 and LALBA polymorphisms on cow's milk yield and quality. The analysis was done on Czech Simmental and Holstein cows. Non-genetic factors were included as well. CSN1S1 did not influence the milk performance in the first lactation. In the second lactation, cows with the BB genotype had significantly higher milk, protein, and fat yields than BC. The differences between LALBA genotypes were non-significant in the first lactation, while in the second lactation, the fat percentage was significantly higher in BB than in AB. The farm significantly influenced milk, protein, and fat yields in both the first and second lactations and fat percentage in the first lactation. The effect of CSN1S1 and LALBA genotypes on the milk technological quality was non-significant. Breed did not influence any of the evaluated technological traits and SCS. The ethanol test was not influenced by farm, season, lactation phase, protein percentage, breed, or non-fat solids percentage. Farm, season, and protein percentage significantly influenced milk fermentation ability, renneting, and SCS. The lactose content is a good indicator of udder health.

Herd and animal factors affect the variability of total and differential somatic cell count in bovine milk

The aim of this study was to quantify some environmental (individual herds, herd productivity, milking system, and season) and animal factors [individual animals, breed, days in milk (DIM) and parity] on the variability of the log-10 transformation of somatic cell count (LSCC) and differential somatic cell count (DSCC) on individual bovine milk. A total of 159,360 test-day records related to milk production and composition were extracted from 12,849 Holstein-Friesian and 9,275 Simmental cows distributed across 223 herds. Herds were classified into high and low productivity, defined according to the average daily milk net energy output (DMEO) yielded by the cows. Data included daily milk yield (DYM; kg/d), milk fat, protein, lactose, SCC, and DSCC, and information on herds (i.e., productivity, milking system). The daily production of total and differential somatic cells in milk was calculated and then log-10 transformed, obtaining DLSCC and DLDSCC, respectively. Data were analyzed using a mixed model including the effects of individual herd, animal, repeated measurements intra animal as random, and herd productivity, milking system, season, breed, DIM, parity, DIM × parity, breed × season, DIM × milking system and parity × milking system as fixed factors. Herds with a high DMEO were characterized by a lower content of LSCC and DSCC, and higher DLSCC and DLDSCC, compared to the low DMEO herds. The association between milking system and somatic cell traits suggested that the use of the automatic milking systems would not allow for a rapid intervention on the cow, as evidenced by the higher content of all somatic cell traits compared to the other milking systems. Season was an important source of variation, as evidenced by high LSCC and DSCC content in milk during summer. Breed of cow had a large influence, with Holstein-Friesian having greater LSCC, DSCC, DLSCC, and DLDSCC compared to Simmental. With regard to DIM, the variability of LSCC was mostly related to that of DSCC, showing an increase from calving to the end of lactation, and suggesting the higher occurrence of chronic mastitis in cows toward the end of lactation. All the somatic cell traits increased across number of parities, possibly because older cows may have increased susceptibility to intramammary infections.

Obtaining of Recombinant Camel Chymosin and Testing Its Milk-Clotting Activity on Cow's, Goat's, Ewes', Camel's and Mare's Milk

In the cheese-making industry, commonly chymosin is used as the main milk-clotting enzyme. Bactrian camel (Camelus bactrianus) chymosin (BacChym) has a milk-clotting activity higher than that of calf chymosin for cow's, goat's, ewes', mare's and camel's milk. A procedure for obtaining milk-clotting reagent based on recombinant camel chymosin is proposed here. Submerged fermentation by a recombinant yeast (Pichia pastoris GS115/pGAPZαA/ProchymCB) was implemented in a 50 L bioreactor, and the recombinant camel chymosin was prepared successfully. The activity of BacChym in yeast culture was 174.5 U/mL. The chymosin was concentrated 5.6-fold by cross-flow ultrafiltration and was purified by ion exchange chromatography. The activity of the purified BacChym was 4700 U/mL. By sublimation-drying with casein peptone, the BacChym powder was obtained with an activity of 36,000 U/g. By means of this chymosin, cheese was prepared from cow's, goat's, ewes', camel's and mare's milk with a yield of 18%, 17.3%, 15.9%, 10.4% and 3%, respectively. Thus, the proposed procedure for obtaining a milk-clotting reagent based on BacChym via submerged fermentation by a recombinant yeast has some prospects for biotechnological applications. BacChym could be a prospective milk-clotting enzyme for different types of milk and their mixtures.

Prediction of fresh and ripened cheese yield using detailed milk composition and udder health indicators from individual Brown Swiss cows

The composition of raw milk is of major importance for dairy products, especially fat, protein, and casein (CN) contents, which are used worldwide in breeding programs for dairy species because of their role in human nutrition and in determining cheese yield (%CY). The aim of the study was to develop formulas based on detailed milk composition to disentangle the role of each milk component on %CY traits. To this end, 1,271 individual milk samples (1.5 L/cow) from Brown Swiss cows were processed according to a laboratory model cheese-making procedure. Fresh %CY (%CY_CURD), total solids and water retained in the fresh cheese (%CY_SOLIDS and %CY_WATER), and 60-days ripened cheese (%CY_RIPENED) were the reference traits and were used as response variables. Training-testing linear regression modeling was performed: 80% of observations were randomly assigned to the training set, 20% to the validation set, and the procedure was repeated 10 times. Four groups of predictive equations were identified, in which different combinations of predictors were tested separately to predict %CY traits: (i) basic composition, i.e., fat, protein, and CN, tested individually and in combination; (ii) udder health indicators (UHI), i.e., fat + protein or CN + lactose and/or somatic cell score (SCS); (iii) detailed protein profile, i.e., fat + protein fractions [CN fractions, whey proteins, and nonprotein nitrogen (NPN) compounds]; (iv) detailed protein profile + UHI, i.e., fat + protein fractions + NPN compounds and/or UHI. Aside from the positive effect of fat, protein, and total casein on %CY, our results allowed us to disentangle the role of each casein fraction and whey protein, confirming the central role of β-CN and κ-CN, but also showing α-lactalbumin (α-LA) to have a favorable effect, and β-lactoglobulin (β-LG) a negative effect. Replacing protein or casein with individual milk protein and NPN fractions in the statistical models appreciably increased the validation accuracy of the equations. The cheese industry would benefit from an improvement, through genetic selection, of traits related to cheese yield and this study offers new insights into the quantification of the influence of milk components in composite selection indices with the aim of directly enhancing cheese production.

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

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

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

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

Effects of breed, farm intensiveness, and cow productivity level on cheese-making ability predicted using infrared spectral data at the population level

Fourier-transform infrared (FTIR) spectra collected during milk recording schemes at population level can be used for predicting novel traits of interest for farm management, cows' genetic improvement, and milk payment systems. The aims of this study were as follows. (1) To predict cheese yield traits using FTIR spectra from routine milk recordings exploiting previously developed calibration equations. (2) To compare the predicted cheese-making abilities of different dairy and dual-purpose breeds. (3) To analyze the effects of herds' level of intensiveness (HL) and of the cow's level of productivity (CL). (4) To compare the patterns of predicted cheese yields with the patterns of milk composition in different breeds to discern the drivers of cheese-making efficiency. The major sources of variation of FTIR predictions of cheese yield ability (fresh cheese or cheese solids produced per unit milk) of individual milk samples were studied on 115,819 cows of 4 breeds (2 specialized dairy breeds, Holstein and Brown Swiss, and 2 dual-purpose breeds, Simmental and Alpine Grey) from 6,430 herds and exploiting 1,759,706 FTIR test-day spectra collected over 7 yr of milk sampling. Calibration equations used were previously developed on 1,264 individual laboratory model cheese procedures (cross-validation R² 0.85 and 0.95 for fresh and solids cheese yields, respectively). The linear model used for statistical analysis included the effects of parity, lactation stage, year of calving, month of sampling, HL, CL, breed of cow, and the interactions breed × HL and breed × CL. The HL and CL stratifications (5 classes each) were based on average daily secretion of milk net energy per cow. All effects were highly significant (P < 0.001). The major conclusions were as follows. (1) The FTIR-based prediction of cheese yield of milk goes beyond the knowledge of fat and protein content, partially explaining differences in cheese-making ability in different cows, breeds and herds. (2) Differences in cheese yields of different breeds are only partially explained by milk fat and protein composition, and less productive breeds are characterized by a higher milk nutrient content as well as a higher recovery of nutrients in the cheese. (3) High-intensive herds not only produce much more milk, but the milk has a higher nutrient content and a higher cheese yield, whereas within herds, compared with less productive cows, the more productive cows have a much greater milk yield, milk with a greater content of fat but not of protein, and a moderate improvement in cheese yield, differing little from expectations based on milk composition. Finally, (4) the effects of HL and CL on milk quality and cheese-making ability are similar but not identical in different breeds, the less productive ones having some advantage in terms of cheese-making ability. We can obtain FTIR-based prediction of cheese yield from individual milk samples retrospectively at population level, which seems to go beyond the simple knowledge of milk composition, incorporating information on nutrient retention ability in cheese, with possible advantages for management of farms, genetic improvement of dairy cows, and milk payment systems.

Nonlinear modeling to describe the pattern of 15 milk protein and nonprotein compounds over lactation in dairy cows

The protein profile of milk includes several caseins, whey proteins, and nonprotein nitrogen compounds, which influence milk's value for human nutrition and its cheesemaking properties for the dairy industry. To fill in the gap in current knowledge of the patterns of these individual nitrogenous compounds throughout lactation, we tested the ability of a parametric nonlinear lactation model to describe the pattern of each N compound expressed qualitatively (as % of total milk N), quantitatively (in g/L milk), and as daily yield (in g/d). The lactation model was tested on a data set of detailed milk nitrogenous compound profiles (15 fractions-12 protein traits and 3 nonproteins-for each expression mode: 45 traits) obtained from 1,342 cows reared in 41 multibreed herds. Our model was a modified version of Wilmink's model, often used for describing milk yield during lactation because of its reliability and ease of parameter interpretation from a biological point of view. We allowed the sign of the persistency coefficient (parameter c) that explained the variation in the long-term milk component (parameter a) to be positive or negative. We also allowed the short-term milk component (parameter b) to be positive or negative, and we estimated a specific speed of adaptation parameter (parameter k) for each trait rather than assumed a value a priori, as in the original model (k = 0.05). These 4 parameters were included in a nonlinear mixed model with cow breed and parity order as fixed effects, and herd-date as random. Combinations of the positive and negative signs of the b and c parameters allowed us to identify 4 differently shaped lactation curves, all found among the patterns exhibited by the nitrogenous fractions as follows: the "zenith" curve (with a maximum peak; for milk yield and 10 other N traits), the "nadir" curve (with a minimum point; for 20 traits, including almost all those expressed in g/L of milk), the "downward" curve (continuously decreasing; for 14 traits, including almost all those in g/d), and the "upward" curve (continuously increasing; only for κ-casein, in % N). Direct estimation of the k parameters specific to each trait showed the large variability in the adaptation speed of fresh cows and greatly increased the model's flexibility. The results indicated that nonlinear parametric mathematical models can effectively describe the different and complex patterns exhibited by individual nitrogenous fractions during lactation; therefore, they could be useful tools for interpreting milk composition variations during lactation.

The mineral profile affects the coagulation pattern and cheese-making efficiency of bovine milk

Natural variations in milk minerals, their relationships, and their associations with the coagulation process and cheese-making traits present an opportunity for the differentiation of milk destined for high-quality natural products, such as traditional specialties or Protected Designation of Origin (PDO) cheeses. The aim of this study was to quantify the effects of the native contents of Ca, P, Na, K, and Mg on 18 traits describing traditional milk coagulation properties (MCP), curd firming over time (CF_t) equation parameters, cheese yield (CY) measures, and nutrient recoveries in the curd (REC) using models that either included or omitted the simultaneous effects of milk fat and casein contents. The results showed that, by including milk fat and casein and the minerals in the statistical model, we were able to determine the specific effects of each mineral on coagulation and cheese-making efficiency. In general, about two-thirds of the apparent effects of the minerals on MCP and the CF_t equation parameters are actually mediated by their association with milk composition, especially casein content, whereas only one-third of the effects are direct and independent of milk composition. In the case of cheese-making traits, the effects of the minerals were mediated only negligibly by their association with milk composition. High Ca content had a positive effect on the coagulation pattern and cheese-making traits, favoring water retention in the curd in particular. Phosphorus positively affected the cheese-making traits in that it was associated with an increase in CY in terms of curd solids, and in all the nutrient recovery traits. However, a very high P content in milk was associated with lower fat recovery in the curd. The variation in the Na content in milk only mildly affected coagulation, whereas with regard to cheese-making, protein recovery was negatively associated with high concentrations of this mineral. Potassium seemed not to be actively involved in coagulation and the cheese-making process. Magnesium content tended to slow coagulation and reduce CY measures. Further studies on the relationships of minerals with casein and protein fractions could deepen our knowledge of the role of all minerals in coagulation and the cheese-making process.

Evaluating the performance of machine learning methods and variable selection methods for predicting difficult-to-measure traits in Holstein dairy cattle using milk infrared spectral data

Fourier-transform infrared (FTIR) spectroscopy is a powerful high-throughput phenotyping tool for predicting traits that are expensive and difficult to measure in dairy cattle. Calibration equations are often developed using standard methods, such as partial least squares (PLS) regression. Methods that employ penalization, rank-reduction, and variable selection, as well as being able to model the nonlinear relations between phenotype and FTIR, might offer improvements in predictive ability and model robustness. This study aimed to compare the predictive ability of 2 machine learning methods, namely random forest (RF) and gradient boosting machine (GBM), and penalized regression against PLS regression for predicting 3 phenotypes differing in terms of biological meaning and relationships with milk composition (i.e., phenotypes measurable directly and not directly in milk, reflecting different biological processes which can be captured using milk spectra) in Holstein-Friesian cattle under 2 cross-validation scenarios. The data set comprised phenotypic information from 471 Holstein-Friesian cows, and 3 target phenotypes were evaluated: (1) body condition score (BCS), (2) blood β-hydroxybutyrate (BHB, mmol/L), and (3) κ-casein expressed as a percentage of nitrogen (κ-CN, % N). The data set was split considering 2 cross-validation scenarios: samples-out random in which the population was randomly split into 10-folds (8-folds for training and 1-fold for validation and testing); and herd/date-out in which the population was randomly assigned to training (70% herd), validation (10%), and testing (20% herd) based on the herd and date in which the samples were collected. The random grid search was performed using the training subset for the hyperparameter optimization and the validation set was used for the generalization of prediction error. The trained model was then used to assess the final prediction in the testing subset. The grid search for penalized regression evidenced that the elastic net (EN) was the best regularization with increase in predictive ability of 5%. The performance of PLS (standard model) was compared against 2 machine learning techniques and penalized regression using 2 cross-validation scenarios. Machine learning methods showed a greater predictive ability for BCS (0.63 for GBM and 0.61 for RF), BHB (0.80 for GBM and 0.79 for RF), and κ-CN (0.81 for GBM and 0.80 for RF) in samples-out cross-validation. Considering a herd/date-out cross-validation these values were 0.58 (GBM and RF) for BCS, 0.73 (GBM and RF) for BHB, and 0.77 (GBM and RF) for κ-CN. The GBM model tended to outperform other methods in predictive ability around 4%, 1%, and 7% for EN, RF, and PLS, respectively. The prediction accuracies of the GBM and RF models were similar, and differed statistically from the PLS model in samples-out random cross-validation. Although, machine learning techniques outperformed PLS in herd/date-out cross-validation, no significant differences were observed in terms of predictive ability due to the large standard deviation observed for predictions. Overall, GBM achieved the highest accuracy of FTIR-based prediction of the different phenotypic traits across the cross-validation scenarios. These results indicate that GBM is a promising method for obtaining more accurate FTIR-based predictions for different phenotypes in dairy cattle.

A2 Bovine Milk and Caprine Milk as a Means of Remedy for Milk Protein Allergy

A new type of cow’s milk, called A2 milk, has appeared in the dairy aisles of supermarkets in recent years. Cows’ milk generally contains two major types of beta-casein as A1 and A2 types, although there are 13 genetic variants of β-casein: A1, A2, A3, A4, B, C, D, E, F, H1, H2, I and G. Studies have shown that A1 β-casein may be harmful, and A2 β-casein is a safer choice for human health especially in infant nutrition and health. The A2 cow milk is reportedly easier to digest and better absorb than A1 or other types of milk. The structure of A2 cow’s milk protein is more comparable to human breast milk, as well as milk from goats, sheep and buffalo. Digestion of A1 type milk produces a peptide called β-casomorphin-7 (BCM-7), which is implicated with adverse gastrointestinal effects on milk consumption. In addition, bovine milk contains predominantly αs1-casein and low levels or even absent in αs2-casein, whereby caprine milk has been recommended as an ideal substitute for patients suffering from allergies against cow milk protein or other food sources. Since goat milk contains relatively low levels of αs1-casein or negligible its content, and αs2-casein levels are high in the milk of most dairy goat breeds, it is logical to assume that children with a high milk sensitivity to αs1-casein should tolerate goat milk well. Cow milk protein allergy (CMPA) is considered a common milk digestive and metabolic disorder or allergic disease with various levels of prevalence from 2.5% in children during the first 3 years of life to 12–30% in infants less than 3 months old, and it can go up to even as high as 20% in some countries. CMPA is an IgE-mediated allergy where the body starts to produce IgE antibodies against certain protein (allergens) such as A1 milk and αs1-casein in bovine milk. Studies have shown that ingestion of β-casein A1 milk can cause ischemic heart disease, type-1 diabetes, arteriosclerosis, sudden infant death syndrome, autism, schizophrenia, etc. The knowledge of bovine A2 milk and caprine αs2-casein has been utilized to rescue CMPA patients and other potential disease problems. This knowledge has been genetically applied to milk production in cows or goats or even whole herds of the two species. This practice has happened in California and Ohio, as well as in New Zealand, where this A2 cow milk has been now advanced commercially. In the USA, there have been even promotions of bulls, whose daughters have been tested homozygous for the A2 β-casein protein.