Interrelationships among physical and chemical traits of cheese: Explanatory latent factors and clustering of 37 categories of cheeses.

Cheese presents extensive variability in physical, chemical, and sensory characteristics according to the variety of processing methods and conditions used to create it. Relationships between the many characteristics of cheeses are known for single cheese types or by comparing a few of them, but not for a large number of cheese types. This case study used the properties recorded on 1,050 different cheeses from 107 producers grouped into 37 categories to analyze and quantify the interrelationships among the chemical and physical properties of many cheese types. The 15 cheese traits considered were ripening length, weight, firmness, adhesiveness, 6 different chemical characteristics, and 5 different color traits. As the 105 correlations between the 15 cheese traits were highly variable, a multivariate analysis was carried out. Four latent explanatory factors were extracted, representing 86% of the covariance matrix: the first factor (38% of covariance) was named Solids because it is mainly linked positively to fat, protein, water-soluble nitrogen, ash, firmness, adhesiveness, and ripening length, and negatively to moisture and lightness; the second factor (24%) was named Hue because it is linked positively to redness/blueness, yellowness/greenness, and chroma, and negatively to hue; the third factor (17%) was named Size because it is linked positively to weight, ripening length, firmness, and protein; and the fourth factor (7%) was named Basicity because it is linked positively to pH. The 37 cheese categories were grouped into 8 clusters and described using the latent factors: the Grana Padano cluster (characterized mainly by high Size scores); hard mountain cheeses (mainly high Solids scores); very soft cheeses (low Solids scores); blue cheeses (high Basicity scores), yellowish cheeses (high Hue scores), and 3 other clusters (soft cheeses, pasta filata and treated rind, and firm mountain cheeses) according to specific combinations of intermediate latent factors and cheese traits. In this case study, the high variability and interdependence of 15 major cheese traits can be substantially explained by only 4 latent factors, allowing us to identify and characterize 8 cheese type clusters.

Lactation modeling and the effects of rotational crossbreeding on milk production traits and milk-spectra-predicted enteric methane emissions.

Rotational crossbreeding has not been widely studied in relation to the enteric methane emissions of dairy cows, nor has the variation in emissions during lactation been modeled. Milk infrared spectra could be used to predict proxies of methane emissions in dairy cows. Therefore, the objective of this work was to study the effects of crossbreeding on the predicted infrared proxies of methane emissions and the variation in the latter during lactation. Milk samples were taken once from 1,059 cows reared in 2 herds, and infrared spectra of the milk were used to predict milk fat (mean ± SD; 3.79 ± 0.81%) and protein (3.68 ± 0.36%) concentrations, yield (21.4 ± 1.5 g/kg dry matter intake), methane intensity (14.2 ± 2.0 g/kg corrected milk), and daily methane production (358 ± 108 g/d). Of these cows, 620 were obtained from a 3-breed (Holstein, Montbéliarde, and Viking Red) rotational mating system, and the rest were purebred Holsteins. Milk production data and methane traits were analyzed using a nonlinear model that included the fixed effects of herd, genetic group, and parity, and the 4 parameters (a, b, c, and k) of a lactation curve modeled using the Wilmink function. Milk infrared spectra were found to be useful for direct prediction of qualitative proxies, such as methane yield and intensity, but not quantitative traits, such as daily methane production, which appears to be better estimated (450 ± 125 g/d) by multiplying a measured daily milk yield by infrared-predicted methane intensity. Lactation modeling of methane traits showed daily methane production to have a zenith curve, similar to that of milk yield but with a delayed peak (53 vs. 37 d in milk), whereas methane intensity is characterized by an upward curve that increases rapidly during the first third of lactation and then slowly till the end of lactation (10.5 g/kg at 1 d in milk to 15.2 g/kg at 300 d in milk). However, lactation modeling was not useful in explaining methane yield, which is almost constant during lactation. Lastly, the methane yield and intensity of cows from 3-breed rotational crossbreeding are not greater, and their methane production is lower than that of purebred Holsteins (452 vs. 477 g/d). Given the greater longevity of crossbred cows, and their lower replacement rate, rotational crossbreeding could be a way of mitigating the environmental impact of milk production.

Detailed mineral profile of milk, whey, and cheese from cows, buffaloes, goats, ewes and dromedary camels, and efficiency of recovery of minerals in their cheese.

Milk and dairy products are important in the human diet not only for the macro nutrients, such as proteins and fats, that they provide, but also for the supply of essential micronutrients, such as minerals. Minerals are present in milk in soluble form in the aqueous phase and in colloidal form associated with the macronutrients of the milk. These 2 forms affect the nutritional functions of the minerals and their contribution to the technological properties of milk during cheese-making. The aim of the present work was to study and compare the detailed mineral profiles of dairy foods (milk, whey, and cheese) obtained from cows, buffaloes, goats, ewes and dromedary camels, and to analyze the recovery in the curd of the individual minerals according to a model cheese-making procedure applied to the milk of these 5 dairy species. The detailed mineral profile of the milk samples was obtained by inductively coupled plasma - optical emission spectroscopy (ICP - OES). We divided the 21 minerals identified in the 3 different matrices into essential macro- and micro-minerals, and environmental micro-minerals, and calculated the recovery of the individual minerals in the cheeses. The complete mineral profiles and the recoveries in the cheeses were then analyzed using a linear mixed model with Species and Food, and their interaction included as fixed effects, and Sample within Species as a random effect. The mineral profiles of each food matrix were then analyzed separately with a general linear model in which only the fixed effect of Species was included. The results showed that the species could be divided into 2 groups: those producing a more diluted milk characterized by a higher content of soluble minerals (in particular K), and those with a more concentrated milk with a higher colloidal mineral content in the skim of the milk (such as Ca and P). The recoveries of the minerals in the curd were in line with the initial content in the milk, and also highlighted the fact that the influence of the brine was not limited to the Na content but to its whole mineral makeup. These results provide valuable information for the evaluation of the nutritional and technological properties of milk, and for the uses made of the byproducts of cheese making from the milk of different species.

Body and milk production traits as indicators of energy requirements and efficiency of purebred Holstein and 3-breed rotational crossbred cows from Viking Red, Montbéliarde, and Holstein sires.

This study aimed to compare rotational 3-breed crossbred cows of Viking Red, Montbéliarde, and Holstein breeds with purebred Holstein cows for a range of body measurements, as well as different metrics of the cows' productivity and production efficiency. The study involved 791 cows (440 crossbreds and 351 purebreds), that were managed across 2 herds. Within each herd, crossbreds and purebreds were reared and milked together, fed the same diets, and managed as one group. The heart girth, height at withers, and body length were measured, and body condition score (BCS) was determined on all the cows on a single test day. The body weight (BW) of 225 cows were used to develop an equation to predict BW from body size traits, parity, and days in milk, which was then used to estimate the BW of all the cows. Equations from the literature were used to estimate body protein and lipid contents using the predicted BW and BCS. Evidence suggests that maintenance energy requirements may be closely related to body protein mass, and Holstein and crossbred cows may be different in body composition. Therefore, we computed the requirements of net energy for maintenance (NEM) on the basis either of the metabolic weight (NEM-MW: 0.418 MJ/kg of metabolic BW) or of the estimated body protein mass according to a coefficient (NEM-PM: 0.631 MJ/kg body protein mass) computed on the subset comprising the purebred Holstein. On the same day when body measurements were collected, individual test-day milk yield and fat and protein contents were retrieved once from the official Italian milk recording system, and milk was sampled to determine fresh cheese yield. Measures of NEM were used to scale the production traits. Statistical analyses of all variables included the fixed effects of herd, days in milk, parity, and genetic group (purebred Holstein and crossbred), and the herd × genetic group interaction. External validation of the equation predicting BW yielded a correlation coefficient of 0.94 and an average bias of -4.95 ± 36.81 kg. The crossbreds had similar predicted BW and NEM-MW compared with the Holsteins. However, NEM-PM of crossbreds was 3.8% lower than that of the Holsteins, due to their 11% greater BCS and different estimated body composition. The crossbred cows yielded 4.8% less milk and 3.4% less milk energy than the purebred Holsteins. However, the differences between genetic groups were no longer significant when the production traits were scaled on NEM-PM, suggesting that the crossbreds and purebreds have the same productive ability and efficiency per unit of body protein mass. In conclusion, measures of productivity and efficiency that combine the cows' production capability with traits related to body composition and the energy cost of production seem to be more effective criteria for comparing crossbred and purebred Holstein cows than just milk, fat, and protein yields.

Cheese yield and nutrients recovery in the curd predicted by Fourier-transform spectra from individual sheep milk samples.

The objectives of this study were to explore the use of Fourier-transform infrared (FTIR) spectroscopy on individual sheep milk samples for predicting cheese-making traits, and to test the effect of the farm variability on their prediction accuracy. For each of 121 ewes from 4 farms, a laboratory model cheese was produced, and 3 actual cheese yield traits (fresh cheese, cheese solids, and cheese water) and 4 milk nutrient recovery traits (fat, protein, total solids, and energy) in the curd were measured. Calibration equations were developed using a Bayesian approach with 2 different scenarios: (1) a random cross-validation (80% calibration; 20% validation set), and (2) a leave-one-out validation (3 farms used as calibration, and the remaining one as validation set) to assess the accuracy of prediction of samples from external farms, not included in calibration set. The best performance was obtained for predicting the yield and recovery of total solids, justifying for the practical application of the method at sheep population and dairy industry levels. Performances for the remaining traits were lower, but still useful for the monitoring of the milk processing in the case of fresh curd and recovery of energy. Insufficient accuracies were found for the recovery of protein and fat, highlighting the complex nature of the relationships among the milk nutrients and their recovery in the curd. The leave-one-out validation procedure, as expected, showed lower prediction accuracies, as a result of the characteristics of the farming systems, which were different between calibration and validation sets. In this regard, the inclusion of information related to the farm could help to improve the prediction accuracy of these traits. Overall, a large contribution to the prediction of the cheese-making traits came from the areas known as "water" and "fingerprint" regions. These findings suggest that, according to the traits studied, the inclusion of water regions for the development of the prediction equation models is fundamental to maintain a high prediction accuracy. However, further studies are necessary to better understand the role of specific absorbance peaks and their contribution to the prediction of cheese-making traits, to offer reliable tools applicable along the dairy ovine chain.

Milk metagenomics and cheese-making properties as affected by indoor farming and summer highland grazing.

The study of the complex relationships between milk metagenomics and milk composition and cheese-making efficiency as affected by indoor farming and summer highland grazing was the aim of the present work. The experimental design considered monthly sampling (over 5 mo) of the milk produced by 12 Brown Swiss cows divided into 2 groups: the first remained on a lowland indoor farm from June to October, and the second was moved to highland pastures in July and then returned to the lowland farm in September. The resulting 60 milk samples (2 kg each) were used to analyze milk composition, milk coagulation, curd firming, and syneresis processes, and to make individual model cheeses to measure cheese yields and nutrient recoveries in the cheese. After DNA extraction and Illumina Miseq sequencing, milk microbiota amplicons were also processed by means of an open-source pipeline called Quantitative Insights Into Microbial Ecology (Qiime2, version 2018.2; https://qiime2.org). Out of a total of 44 taxa analyzed, 13 bacterial taxa were considered important for the dairy industry (lactic acid bacteria, LAB, 5 taxa; and spoilage bacteria, 4) and for human (other probiotics, 2) and animal health (pathogenic bacteria, 2). The results revealed the transhumant group of cows transferred to summer highland pastures showed an increase in almost all the LAB taxa, bifidobacteria, and propionibacteria, and a reduction in spoilage taxa. All the metagenomic changes disappeared when the transhumant cows were moved back to the permanent indoor farm. The relationships between 17 microbial traits and 30 compositional and technological milk traits were investigated through analysis of correlation and latent explanatory factor analysis. Eight latent factors were identified, explaining 75.3% of the total variance, 2 of which were mainly based on microbial traits: pro-dairy bacteria (14% of total variance, improving during summer pasturing) and pathogenic bacteria (6.0% of total variance). Some bacterial traits contributed to other compositional-technological latent factors (gelation, udder health, and caseins).

Effects of breed, farm intensiveness, and cow productivity on infrared predicted milk urea.

Milk urea content is receiving growing interest from science and industry as a tool to infer the protein adequacy of dairy cows' diets, nitrogen excretion and its environmental impact, and efficiency of animals' protein metabolism. Fourier-transform infrared (FTIR) prediction is a high-throughput method for rapidly and cheaply evaluating milk urea content at the population level. Existing knowledge of the major sources of variation (e.g., year, season, farming system, individual herd, and the cow's breed, parity, stage of lactation, and productive potential) is fragmentary. The objective of this work was to study at the population level the simultaneous effects of all the major sources of variation and their most important interactions. Milk urea content in 1,759,706 test day milk samples collected from 291,129 lactations of 115,819 cows from 6,430 herds over 8 yr was predicted by FTIR. The milk urea content data (and also milk protein percentage, for comparison) were analyzed using a linear model that included the effects of parity, days in milk (DIM) class, year, month, herd intensiveness level, cow productivity level, breed, and herd intensiveness and cow productivity levels within breed. All sources of variation of milk urea content proved highly significant, the most important in terms of F-value being breed > year > herd intensiveness level > parity. The ranking for milk protein was very different (DIM class > herd intensiveness level > parity > breed). The patterns of the least squares means for urea and protein contents of milk were also very different and sometimes contrasting. The seasonal variation in urea was sinusoidal. Urea content increased during the first 4 mo of lactation and then remained almost stable before decreasing after 11 mo. Specialized dairy breeds had lower average milk urea content than dual-purpose breeds; in the former case it was lower in Holsteins than in Brown Swiss, and in the latter it was lower in Simmentals than in Alpine Greys. The effect of herd intensiveness level was much stronger than the effect of cow productivity level; the increase in milk urea with increasing herd average daily milk yield was almost linear in the case of dairy breeds but curvilinear in dual-purpose breeds. The large differences in breed and the modest relationships with the cow's productive potential require further analysis at the genetic level to obtain information of potential use in genetic improvement of the dairy cow populations.

Invited review: A comprehensive review of visible and near-infrared spectroscopy for predicting the chemical composition of cheese.

Substantial research has been carried out on rapid, nondestructive, and inexpensive techniques for predicting cheese composition using spectroscopy in the visible and near-infrared radiation range. Moreover, in recent years, new portable and handheld spectrometers have been used to predict chemical composition from spectra captured directly on the cheese surface in dairies, storage facilities, and food plants, removing the need to collect, transport, and process cheese samples. For this review, we selected 71 papers (mainly dealing with prediction of the chemical composition of cheese) and summarized their results, focusing our attention on the major sources of variation in prediction accuracy related to cheese variability, spectrometer and spectra characteristics, and chemometrics techniques. The average coefficient of determination obtained from the validation samples ranged from 86 to 90% for predicting the moisture, fat, and protein contents of cheese, but was lower for predicting NaCl content and cheese pH (79 and 56%, respectively). There was wide variability with respect to all traits in the results of the various studies (standard deviation: 9-30%). This review draws attention to the need for more robust equations for predicting cheese composition in different situations; the calibration data set should consist of representative cheese samples to avoid bias due to an overly specific field of application and ensure the results are not biased for a particular category of cheese. Different spectrometers have different accuracies, which do not seem to depend on the spectrum extension. Furthermore, specific areas of the spectrum-the visible, infrared-A, or infrared-B range-may yield similar results to broad-range spectra; this is because several signals related to cheese composition are distributed along the spectrum. Small, portable instruments have been shown to be viable alternatives to large bench-top instruments. Last, chemometrics (spectra pre-treatment and prediction models) play an important role, especially with regard to difficult-to-predict traits. A proper, fully independent, validation strategy is essential to avoid overoptimistic results.

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.

Composition, coagulation properties, and predicted cheesemaking traits of bulk goat milk from different farming systems, breeds, and stages of production.

At the global level, the quantity of goat milk produced and its gross production value have increased considerably over the last 2 decades. Although many scientific papers on this topic have been published, few studies have been carried out on bulk goat milk samples. The aim of the present study was to investigate in the field the effects of farming system, breed type, individual flock, and stage of production on the composition, coagulation properties (MCP), curd firming over time parameters (CFt), predicted cheese yield (CY%), and nutrient recovery traits (REC) of 432 bulk milk samples from 161 commercial goat farms in Sardinia, Italy. We found that the variance due to individual flock was of the same order as the residual variance for almost all composition and cheesemaking traits. With regard to the fixed effects, the effect of farming system on bulk milk variability was not highly significant for the majority of traits (it was lower than individual flock), whereas the effects of breed type and stage of production were much higher. More specifically, the intensive farms produced milk with the best concentrations of almost all constituents, whereas extensive farms exhibited faster rennet coagulation times, a slower rate of curd firming, lower potential curd firmness, and lower percentages of fat and energy recoveries in the fresh curd. Farms rearing the local breed, Sarda, alone or together with the Maltese breed, produced milk with the best concentrations of fat and protein, superior curd firmness, and better predicted percentage of fresh curd (CYCURD) and recovery traits. The results show the potential of both types of breed, either for their quantitative (specialized breeds) or their qualitative (local breeds) attributes. As expected, the concentrations of fat, protein fractions, and lactose were influenced by the stage of production, with samples collected in the early stage of production (in February and March) having a greater quantity of the main constituents. Somatic cells reached the highest levels in the late stage of production, which corresponds to the goats' advanced stage of lactation (June-July), although no differences were present in the logarithmic bacterial counts between the early and late stages. Regarding cheesemaking potential, bulk milk samples of the late stage were characterized by delayed rennet coagulation and curd firming times, the lowest values of curd firmness, and a general reduction in CY%, and REC traits. In conclusion, we highlight several issues regarding the effects of the most important sources of variation on bulk goat milk, and point to some critical factors relevant for improving dairy goat farming and milk production.

Composition and aptitude for cheese-making of milk from cows, buffaloes, goats, sheep, dromedary camels, and donkeys.

Bovines produce about 83% of the milk and dairy products consumed by humans worldwide, the rest represented by bubaline, caprine, ovine, camelid, and equine species, which are particularly important in areas of extensive pastoralism. Although milk is increasingly used for cheese production, the cheese-making efficiency of milk from the different species is not well known. This study compares the cheese-making ability of milk sampled from lactating females of the 6 dairy species in terms of milk composition, coagulation properties (using lactodynamography), curd-firming modeling, nutrients recovered in the curd, and cheese yield (through laboratory model-cheese production). Equine (donkey) milk had the lowest fat and protein content and did not coagulate after rennet addition. Buffalo and ewe milk yielded more fresh cheese (25.5 and 22.9%, respectively) than cow, goat, and dromedary milk (15.4, 11.9, and 13.8%, respectively). This was due to the greater fat and protein contents of the former species with respect to the latter, but also to the greater recovery of fat in the curd of bubaline (88.2%) than in the curd of camelid milk (55.0%) and consequent differences in the recoveries of milk total solids and energy in the curd; protein recovery, however, was much more similar across species (from 74.7% in dromedaries to 83.7% in bovine milk). Compared with bovine milk, the milk from the other Artiodactyla species coagulated more rapidly, reached curd firmness more quickly (especially ovine milk), had a more pronounced syneresis (especially caprine milk), had a greater potential asymptotical curd firmness (except dromedary and goat milk), and reached earlier maximum curd firmness (especially caprine and ovine milk). The maximum measured curd firmness was greater for bubaline and ovine milk, intermediate for bovine and caprine milk, and lower for camelid milk. The milk of all ruminant species can be used to make cheese, but, to improve efficiency, cheese-making procedures need to be optimized to take into account the large differences in their coagulation, curd-firming, and syneresis properties.

Integrating genomic and infrared spectral data improves the prediction of milk protein composition in dairy cattle.

BACKGROUND: Over the past decade, Fourier transform infrared (FTIR) spectroscopy has been used to predict novel milk protein phenotypes. Genomic data might help predict these phenotypes when integrated with milk FTIR spectra. The objective of this study was to investigate prediction accuracy for milk protein phenotypes when heterogeneous on-farm, genomic, and pedigree data were integrated with the spectra. To this end, we used the records of 966 Italian Brown Swiss cows with milk FTIR spectra, on-farm information, medium-density genetic markers, and pedigree data. True and total whey protein, and five casein, and two whey protein traits were analyzed. Multiple kernel learning constructed from spectral and genomic (pedigree) relationship matrices and multilayer BayesB assigning separate priors for FTIR and markers were benchmarked against a baseline partial least squares (PLS) regression. Seven combinations of covariates were considered, and their predictive abilities were evaluated by repeated random sub-sampling and herd cross-validations (CV). RESULTS: Addition of the on-farm effects such as herd, days in milk, and parity to spectral data improved predictions as compared to those obtained using the spectra alone. Integrating genomics and/or the top three markers with a large effect further enhanced the predictions. Pedigree data also improved prediction, but to a lesser extent than genomic data. Multiple kernel learning and multilayer BayesB increased predictive performance, whereas PLS did not. Overall, multilayer BayesB provided better predictions than multiple kernel learning, and lower prediction performance was observed in herd CV compared to repeated random sub-sampling CV. CONCLUSIONS: Integration of genomic information with milk FTIR spectral can enhance milk protein trait predictions by 25% and 7% on average for repeated random sub-sampling and herd CV, respectively. Multiple kernel learning and multilayer BayesB outperformed PLS when used to integrate heterogeneous data for phenotypic predictions.

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.

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 R2 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. 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.

Purebreeding with sexed semen and crossbreeding with semen from double-muscled sires to improve beef production from dairy herds: Live and slaughter performances of crossbred calves.

The use of sexed semen to produce purebred replacement heifers allows a large proportion of dairy cows to be mated to beef sires, and quantitative and qualitative improvements to be made to beef production from dairy herds. The major dairy and beef breeds are undergoing rapid genetic improvement as a result of more efficient selection methods, prompting a need to evaluate the meat production of crossbred beef × dairy cattle produced using current genetics. As part of a large project involving 125 commercial dairy farms, we evaluated the combined use of purebreeding with sexed semen and crossbreeding with semen from beef sires, particularly double-muscled breeds. A survey of 1,530 crossbred calves revealed that, whereas purebred dairy calves are destined almost exclusively for veal production, beef × dairy crossbred calves are also destined for beef production after fattening on either the dairy farm of birth or by specialized fatteners. In veal production, compared with Belgian Blue-sired calves (taken as the reference), double-muscled INRA 95-sired calves had a lighter slaughter weight (303 vs. 346 kg), but a greater dressing percent (62.3 vs. 58.4%). Limousin (also known as Limousine)-sired calves had a smaller average daily gain (1.26 vs. 1.34 kg/d), and lighter slaughter (314 vs. 346 kg) and carcass weights (182 vs. 201 kg). Last, Simmental-sired calves had a similar growth rate, but lighter carcass weight (177 vs. 201 kg), smaller dressing percentage (55.3 vs. 58.4%), and smaller muscularity scores (3.25 vs. 3.72). In the case of young bulls and heifers fattened on the dairy farm of birth, Belgian Blue-, Piemontese (also known as Piedmontese)-, and Limousin-sired calves performed similarly; the only exception was that Piemontese-sired calves had a greater dressing percentage. Belgian Blue- and Limousin-sired calves performed similarly when fattened by specialized beef producers. In both veal and beef production, the effects of dam breed were less important than sire breed. Considering the entire project, we can conclude that the combined use of sexed semen for purebreeding and conventional beef semen for terminal crossbreeding improves meat production from dairy herds, especially when the sires are double-muscled beef breeds.

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.

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 (CFt) 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 CFt 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.

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

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

Shift in the cow milk microbiota during alpine pasture as analyzed by culture dependent and high-throughput sequencing techniques.

In the present study, two groups of cows from a permanent lowland farm (PF) were divided during summer and reared in the PF or in a temporary alpine farm (ALP), respectively. Microbiological analyses were performed with the objective to investigate the microbial evolution of milk before, during, and after summer transhumance comparing, in particular, the two groups of cows to determine whether the alpine pasture could directly influence the milk microbiota. A significant increase of all microbial groups was registered in milk samples collected in the ALP. Interestingly, many strains belonging to species with well reported technological and probiotic activities were isolated from Alpine milk (20% Lactococcus lactis subsp. lactis/cremoris, 18% Lactobacillus paracasei, 14% Bifidobacterium crudilactis and 18% Propionibacterium sp.), whereas only 16% of strains isolated from the permanent farm milk belonged to the species Lactococcus lactis subsp. lactis/cremoris, 6% to Lactobacillus paracasei, 2% to Bifidobacterium crudilactis and 5% to Propionibacterium sp. The MiSeq Illumina data showed that Alpine milk presented a significant reduction of Pseudomonas and an increase of Lactococcus, Bifidobacterium and Lactobacillus genera. These data confirmed the practice of Alpine pasture as one of the main drivers affecting the milk microbiota. All the microbial changes disappeared when cows were delivered back from Alpine pasture to the indoor farm.

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

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