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.

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

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.

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

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.