The use of milk Fourier-transform infrared spectra for predicting cheesemaking traits in Grana Padano Protected Designation of Origin cheese

The prediction of the cheese yield (%CY) traits for curd, solids, and retained water and the amount of fat, protein, solids, and energy recovered from the milk into the curd (%REC) by Bayesian models, using Fourier-transform infrared spectroscopy (FTIR), can be of significant economic interest to the dairy industry and can contribute to the improvement of the cheese process efficiency. The yields give a quantitative measure of the ratio between weights of the input and output of the process, whereas the nutrient recovery allows to assess the quantitative transfer of a component from milk to cheese (expressed in % of the initial weight). The aims of this study were: (1) to investigate the feasibility of using bulk milk spectra to predict %CY and %REC traits, and (2) to quantify the effect of the dairy industry and the contribution of single-spectrum wavelengths on the prediction accuracy of these traits using vat milk samples destined to the production of Grana Padano Protected Designation of Origin cheese. Information from 72 cheesemaking days (in total, 216 vats) from 3 dairy industries were collected. For each vat, the milk was weighed and analyzed for composition (total solids [TS], lactose, protein, and fat). After 48 h from cheesemaking, each cheese was weighed, and the resulting whey was sampled for composition as well (TS, lactose, protein, and fat). Two spectra from each milk sample were collected in the range between 5,011 and 925 cm-1 and averaged before the data analysis. The calibration models were developed via a Bayesian approach by using the BGLR (Bayesian Generalized Linear Regression) package of R software. The performance of the models was assessed by the coefficient of determination (R2VAL) and the root mean squared error (RMSEVAL) of validation. Random cross-validation (CVL) was applied [80% calibration and 20% validation set] with 10 replicates. Then, a stratified cross-validation (SCV) was performed to assess the effect of the dairy industry on prediction accuracy. The study was repeated using a selection of informative wavelengths to assess the necessity of using whole spectra to optimize prediction accuracy. Results showed the feasibility of using FTIR spectra and Bayesian models to predict cheesemaking traits. The R2VAL values obtained with the CVL procedure were promising in particular for the %CY and %REC for protein, ranging from 0.44 to 0.66 with very low RMSEVAL (from 0.16 to 0.53). Prediction accuracy obtained with the SCV was strongly influenced by the dairy factory industry. The general low values gained with the SCV do not permit a practical application of this approach, but they highlight the importance of building calibration models with a dataset covering the largest possible sample variability. This study also demonstrated that the use of the full FTIR spectra may be redundant for the prediction of the cheesemaking traits and that a specific selection of the most informative wavelengths led to improved prediction accuracy. This could lead to the development of dedicated spectrometers using selected wavelengths with built-in calibrations for the online prediction of these innovative traits.

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.

Effects of Supplementation of a Mycotoxin Mitigation Feed Additive in Lactating Dairy Cows Fed Fusarium Mycotoxin-Contaminated Diet for an Extended Period

Fusarium mycotoxins are inactivated by rumen flora; however, a certain amount can pass the rumen and reticulum or be converted into biological active metabolites. Limited scientific evidence is available on the impact and mitigation of Fusarium mycotoxins on dairy cows' performance and health, particularly when cows are exposed for an extended period (more than 2 months). The available information related to these mycotoxin effects on milk cheese-making parameters is also very poor. The objective of this study was to evaluate a commercially available mycotoxin mitigation product (MMP, i.e., TOXO® HP-R, Selko, Tilburg, The Netherlands) in lactating dairy cows fed a Fusarium mycotoxin-contaminated diet, and the repercussions on the dry matter intake, milk yield, milk quality, cheese-making traits and health status of cows. The MMP contains smectite clays, yeast cell walls and antioxidants. In the study, 36 lactating Holstein cows were grouped based on the number of days of producing milk, milk yield, body condition score and those randomly assigned to specific treatments. The study ran over 2 periods (March/May-May/July 2022). In each period, six animals/treatment were considered. The experimental periods consisted of 9 days of adaptation and 54 days of exposure. The physical activity, rumination time, daily milk production and milk quality were measured. The cows were fed once daily with the same total mixed ration (TMR) composition. The experimental groups consisted of a control (CTR) diet, with a TMR with low contamination, high moisture corn (HMC), and beet pulp; a mycotoxins (MTX) diet, with a TMR with highly contaminated HMC, and beet pulp; and an MTX diet supplemented with 100 g/cow/day of the mycotoxin mitigation product (MMP). The trial has shown that the use of MMP reduced the mycotoxin's negative effects on the milk yield and quality (protein, casein and lactose). The MTX diet had a lower milk yield and feed efficiency than the CTR and MMP HP-R diets. The MMP limited the negative effect of mycotoxin contamination on clotting parameters, mitigating the variations on some coagulation properties; however, the MMP inclusion tended to decrease the protein and apparent starch digestibility of the diet. These results provide a better understanding of mycotoxin risk on dairy cows' performances and milk quality. The inclusion of an MMP product mitigated some negative effects of the Fusarium mycotoxin contamination during this trial. The major effects were on the milk yield and quality in both the experimental periods. These results provide better insight on the effects of mycotoxins on the performance and quality of milk, as well as the cheese-making traits. Further analyses should be carried out to evaluate MMP's outcome on immune-metabolic responses and diet digestibility.

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 (NE_M) on the basis either of the metabolic weight (NE_M-MW: 0.418 MJ/kg of metabolic BW) or of the estimated body protein mass according to a coefficient (NE_M-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 NE_M 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 NE_M-MW compared with the Holsteins. However, NE_M-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 NE_M-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.

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.

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 (CF_t), 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 (CY_CURD) 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.

Predictive formulas for different measures of cheese yield using milk composition from individual goat samples

The objective of this study was to develop formulas based on milk composition of individual goat samples for predicting cheese yield (%CY) traits (fresh curd, milk solids, and water retained in the curd). The specific aims were to assess and quantify (1) the contribution of major milk components (fat, protein, and casein) and udder health indicators (lactose, somatic cell count, pH, and bacterial count) on %CY traits (fresh curd, milk solids, and water retained in the curd); (2) the cheese-making method; and (3) goat breed effects on prediction accuracy of the %CY formulas. The %CY traits were analyzed in duplicate from 600 goats, using an individual laboratory cheese-making procedure (9-MilCA method; 9 mL of milk per observation) for a total of 1,200 observations. Goats were reared in 36 herds and belonged to 6 breeds (Saanen, Murciano-Granadina, Camosciata delle Alpi, Maltese, Sarda, and Sarda Primitiva). Fresh %CY (%CY_CURD), total solids (%CY_SOLIDS), and water retained (%CY_WATER) in the curd were used as response variables. Single and multiple linear regression models were tested via different combinations of standard milk components (fat, protein, casein) and indirect udder health indicators (UHI; lactose, somatic cell count, pH, and bacterial count). The 2 %CY observations within animal were averaged, and a cross-validation (CrV) scheme was adopted, in which 80% of observations were randomly assigned to the calibration (CAL) set and 20% to the validation (VAL) set. The procedure was repeated 10 times to account for sampling variability. Further, the model presenting the best prediction accuracy in CrV (i.e., comprehensive formula) was used in a secondary analysis to assess the accuracy of the %CY predictive formulas as part of the laboratory cheese-making procedure (within-animal validation, WAV), in which the first %CY observation within animal was assigned to CAL, and the second to the VAL set. Finally, a stratified CrV (SCrV) was adopted to assess the %CY traits prediction accuracy across goat breeds, again using the best model, in which 5 breeds were included in CAL and the remaining one in the VAL set. Fitting statistics of the formulas were assessed by coefficient of determination of validation (R²_VAL) and the root mean square error of validation (RMSE_VAL). In CrV, the formula with the best prediction accuracy for all %CY traits included fat, casein, and UHI (R²_VAL = 0.65, 0.96, and 0.23 for %CY_CURD, %CY_SOLIDS, and %CY_WATER, respectively). The WAV procedure showed R²_VAL higher than those obtained in CrV, evidencing a low effect of the 9-MilCA method and, indirectly, its high repeatability. In the SCrV, large differences for %CY_CURD and %CY_WATER among breeds evidenced that the breed is a fundamental factor to consider in %CY predictive formulas. These results may be useful to monitor milk composition and quantify the influence of milk traits in the composite selection indices of specific breeds, and for the direct genetic improvement of cheese production.

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

Bovines account for 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.

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.

Genetic parameters of differential somatic cell count, milk composition, and cheese-making traits measured and predicted using spectral data in Holstein cows

Mastitis is one of the most prevalent diseases in dairy cattle and is the cause of considerable economic losses. Alongside somatic cell count (SCC), differential somatic cell count (DSCC) has been recently introduced as a new indicator of intramammary infection. The DSCC is expressed as a count or a proportion (%) of polymorphonuclear neutrophils plus lymphocytes (PMN-LYM) in milk somatic cells. These numbers are complemented to total somatic cell count or to 100 by macrophages (MAC). The aim of this study was to investigate the genetic variation and heritability of DSCC, and its correlation with milk composition, udder health indicators, milk composition, and technological traits in Holstein cattle. Data used in the analysis consisted in single test-day records from 2,488 Holstein cows reared in 36 herds located in northern Italy. Fourier-transform infrared (FTIR) spectroscopy was used to predict missing information for some milk coagulation and cheese-making traits, to increase sample size and improve estimation of the genetic parameters. Bayesian animal models were implemented via Gibbs sampling. Marginal posterior means of the heritability estimates were 0.13 for somatic cell score (SCS); 0.11 for DSCC, MAC proportion, and MAC count; and 0.10 for PMN-LYM count. Posterior means of additive genetic correlations between SCS and milk composition and udder health were low to moderate and unfavorable. All the relevant genetic correlations between the SCC traits considered and the milk traits (composition, coagulation, cheese yield and nutrients recovery) were unfavorable. The SCS showed genetic correlations of -0.30 with the milk protein proportion, -0.56 with the lactose proportion and -0.52 with the casein index. In the case of milk technological traits, SCS showed genetic correlations of 0.38 with curd firming rate (k₂₀), 0.45 with rennet coagulation time estimated using the curd firming over time equation (RCT_eq), -0.39 with asymptotic potential curd firmness, -0.26 with maximum curd firmness (CF_max), and of -0.31 with protein recovery in the curd. Differential somatic cell count expressed as proportion was correlated with SCS (0.60) but had only 2 moderate genetic correlations with milk traits: with lactose (-0.32) and CF_max (-0.33). The SCS was highly correlated with the log PMN-LYM count (0.79) and with the log MAC count (0.69). The 2 latter traits were correlated with several milk traits: fat (-0.38 and -0.43 with PMN-LYM and MAC counts, respectively), lactose percentage (-0.40 and -0.46), RCT_eq (0.53 and 0.41), t_max (0.38 and 0.48). Log MAC count was correlated with k₂₀ (+0.34), and log PMN-LYM count was correlated with CF_max (-0.26) and weight of water curd as percentage of weight of milk processed (-0.26). The results obtained offer new insights into the relationships between the indicators of udder health and the milk technological traits in Holstein cows.

Coagulation process in Manchega sheep milk from Spain: A path analysis approach

Characteristics of sheep milk are of great interest for the dairy industry, as almost the totality of production is intended for cheesemaking. However, the existing relationships between these variables are complex. This study assessed composition, hygienic quality, coagulation properties, and curd yield of 1,200 individual Manchega sheep milk samples. The aim was to compare the effect of composition and hygienic quality on coagulation and curdling, and to evaluate the relationship between curd yields and the coagulation process and the effect of other features by using path analysis methodologies. Outcomes proved path analysis to be a useful and effective tool to assess these relationships through direct and indirect paths within the same model. Results showed that the factors that had a direct influence on milk coagulation were lactose concentration, casein content, and initial pH of milk. Contrastingly, somatic cells did not seem to have any effect (direct or indirect) on the coagulation process. Factors that directly affected curd yield were fat content, lactose concentration, casein content, and curd moisture. However, technological parameters showed little effect over curd yield.

Associations between differential somatic cell count and milk yield, quality, and technological characteristics in Holstein cows

The aim of this study was to investigate the associations between differential somatic cell count (DSCC) and milk quality and udder health traits, and for the first time, between DSCC and milk coagulation properties and cheesemaking traits in a population of 1,264 Holstein cows reared in northern Italy. Differential somatic cell count represents the combined proportions of polymorphonuclear neutrophils plus lymphocytes (PMN-LYM) in the total somatic cell count (SCC), with macrophages (MAC) making up the remaining proportion. The milk traits investigated in this study were milk yield (MY), 8 traits related to milk composition and quality (fat, protein, casein, casein index, lactose, urea, pH, and milk conductivity), 9 milk coagulation traits [3 milk coagulation properties (MCP) and 6 curd firming (CF) traits], 7 cheesemaking traits, 3 cheese yield (CY) traits, and 4 milk nutrient recovery in the curd (REC) traits. A linear mixed model was fitted to explore the associations between SCS combined with DSCC and the aforementioned milk traits. An additional model was run, which included DSCC expressed as the PMN-LYM and MAC counts, obtained by multiplying the percentage of PMN-LYM and MAC by SCC in the milk for each cow in the data set. The unfavorable association between SCS and milk quality and technological traits was confirmed. Increased DSCC was instead associated with a linear increase in MY, casein index, and lactose proportion and a linear decrease in milk fat and milk conductivity. Accordingly, DSCC was favorably associated with all MCP and CF traits (with the exception of the time needed to achieve maximum, CF), particularly with rennet coagulation time, and it always displayed linear relationships. Differential somatic cell count was also positively associated with the recovery of milk nutrients in the curd (protein, fat, and energy), which increased linearly with increasing DSCC. The PMN-LYM count was rarely associated with milk traits, even though the pattern observed confirmed the results obtained when both SCS and DSCC were included in the model. The MAC count, however, showed the opposite pattern: MY, casein index, and lactose percentage decreased and milk conductivity increased with an increasing MAC count. No significant association was found between PMN-LYM count and MCP, CF, CY, and REC traits, whereas MAC count was unfavorably associated with MCP, CF traits, some CY traits, and all REC traits. Our results showed that the combined information derived from SCS and DSCC might be useful to monitor milk quality and cheesemaking-related traits.

Exploring the Genotype at CSN3 Gene, Milk Composition, Coagulation and Cheese-Yield Traits of the Sardo-Modicana, an Autochthonous Cattle Breed from the Sardinia Region, Italy

Simple Summary

The Sardo-Modicana is a local cattle breed from Sardinia, Italy. It originated from the crossing of Sardinian local cows with Modicana bulls from Sicily, imported in the late 19th century. In the 1950s, approximately 60,000 heads were present, but nowadays the total population has decreased to about 1800 animals. It is a multipurpose breed, and animals are farmed using extensive methods. Traditionally, cows are hand-milked and milk is destined to produce a traditional pasta filata cheese. In the literature, the information about the dairy potential of this breed is scarce. The present study evidenced the favorable genetic patterns, milk composition and coagulation traits of the Sardo-Modicana cattle breed; such information will be useful for the preservation and enhancement of the breed.

Abstract

The Sardo-Modicana is a local cattle breed from Sardinia, Italy. No information about its dairy potential is available in the literature. This study investigated the genotype at the CSN3 gene and milk traits of the Sardo-Modicana cattle breed. Fifty-four cows were sampled for DNA extraction and genotyping at the κ-casein gene locus, CSN3. Forty individual milk samples were analyzed for milk composition, milk coagulation properties and cheese yield (CY%). All the Sardo-Modicana cows were BB homozygotes at CSN3. Hence, the results were compared with the other two local Sardinian breeds. Eighty-three Sarda and 21 Sardo-Bruna cows were genotyped, and the A allele was found (at frequencies of 0.416 and 0.405, respectively). As regards milk traits, the mean protein value was 3.74 g/100 mL, and the mean casein value was 2.98 g/100 mL. Total bacterial and somatic cell counts showed excellent levels of hygiene considering the extensive farming and hand milking. In addition, milk produced by Sardo-Modicana cows was characterized by favorable values of coagulation properties and cheese yield. This information may represent a starting point for the conservation and enhancement of this breed.

Characterization of milk composition, coagulation properties, and cheese-making ability of goats reared in extensive farms

The aims of this study were to explore the variability of milk composition, coagulation properties, and cheese-making traits of the Sarda goat breed, and to investigate the effects of animal and farm factors, and the geographic area (Central-East vs. South-West) of an insular region of Italy, Sardinia. A total of 570 Sarda goats reared in 21 farms were milk-sampled during morning milking. Individual milk samples were analyzed for composition, traditional milk coagulation properties (MCP), modeled curd-firming over time parameters (CF_t), and cheese-making traits (cheese yield, %CY; recovery of nutrients, %REC; daily cheese yield, dCY). Farms were classified into 2 categories based on milk energy level (MEL; high or low), defined according to the average net energy of milk daily produced by the lactating goats. Milk yield and composition were analyzed using a mixed model including the fixed effects of MEL, geographic area, days in milk, and parity, and the random effect of farm within MEL and geographic area. Data about MCP, CF_t, and the cheese-making process were analyzed using the same model, with the inclusion of the effects of animal and pendulum of the lactodynamograph instrument, allowing the measure of repeatability of these traits. Results showed that animal had greater influence on coagulation and cheese-making traits compared with farm effect. Days in milk influenced milk composition, whose changes partly reflected the modifications of %CY traits. Moreover, large differences were observed between primiparous and multiparous goats: primiparous goats produced less milk of better quality (higher fat, lower somatic cell and bacterial counts) and less cheese, but with higher recovery of fat and protein in the curd, compared with multiparous goats. The repeatability was very high, for both coagulation (84.0 to 98.8%) and cheese-making traits (89.7 to 99.9%). The effect of MEL was significant for daily productions of milk and cheese, coagulation time, and recovery of protein in the curd, which were better in high-MEL farms. As regards geographic area, milk composition and percentage cheese yield were superior in the Central-East area, whereas daily milk and cheese production and MCP were better in the South-West. This result was explainable by the phenomenon of crossbreeding Sarda goats with Maltese bucks, which occurred with greater intensity in the South-West than in the Central-East area of the island. The results provided by this study could be of great interest for the goat dairy sector. Indeed, the methods described in the present study could be applicable for other farming methods, goat breeds, and geographic areas. The collection of a wide range of phenotypes at individual animal level is fundamental for the characterization of local populations and can be used to guarantee breed conservation and the persistence of traditional farming systems, and to increase farmers' profit.

Effect of Feeding Adaptation of Italian Simmental Cows before Summer Grazing on Animal Behavior and Milk Characteristics

Simple Summary

The traditional transhumant system of rearing dairy cows in mountain areas expects animals to remain indoors in the valley during the cold season, whereas during the summer they are moved to pastures at progressively higher altitudes. The animals transferred from the valley farm to the alpine pasture must adapt to various management changes. This study aimed to evaluate whether a gradual inclusion of fresh grass in the diet of dairy cows in the valley farm can improve the performance and milk characteristics during summer grazing. Three groups of six animals each were considered: one group was kept in the stable, one was transferred from the valley to the summer farm without adaptation, and the other was progressively adapted to grazing with a feeding adaptation period. Compared to animals kept indoors, grazing animals had similar performance and milk characteristics, higher rumination time and, with respect to volatile compounds in milk, higher concentrations of alcohols, aldehydes, hydrocarbons, and ketones but lower concentrations of organic acids, phenolic compounds, and dimethyl sulfone, regardless of the feeding adaptation. In conclusion, the gradual inclusion of fresh grass in the diet in the valley farm did not improve the performance and milk characteristics during summer grazing.

Abstract

According to the alpine transhumance system, dairy cows are moved from indoor feeding with conserved forage to fresh herbage feeding on pasture. The aim of this study was to assess, as a feeding adaptation technique, the effect of a gradual inclusion of fresh herbage in the diet of Italian Simmental dairy cows before their transfer to alpine pasture on performance, behavior, and milk characteristics. Eighteen cows were assigned to three groups: animals transferred to alpine pasture with a 10-d feeding adaptation period consisting in gradual access to a pasture close to the valley farm (GT), animals transferred to alpine pasture without a feeding adaptation period (AT), and animals kept in the valley farm (IND). During the first two weeks of summer grazing, GT and AT showed higher rumination time and different concentrations of ketones, hydrocarbons, organic acids, toluene, alcohols, phenols, and dimethyl sulfone in milk as compared to IND, whereas no differences were found in milk yield, composition, or coagulation properties. No differences between GT and AT were evident for the studied variables. The feeding adaptation technique used in this study did not influence the performance and milk characteristics of Italian Simmental dairy cows grazing on alpine pasture.

Goat cheese yield and recovery of fat, protein, and total solids in curd are affected by milk coagulation properties

The aims of the present research were to quantify the effects of each coagulation trait, traditional milk coagulation properties [MCP: rennet coagulation time (RCT), curd-firming time (k₂₀), and curd firmness at 30 min (a₃₀)], and modeled curd-firming over time (CF_t) parameters [estimated rennet coagulation time (RCT_eq), curd-firming instant rate constant (k_CF), and potential curd firmness (CF_P)] directly on the following: (1) recovery of 3 milk components in the curd (%REC), (2) 3 measures of cheese yield (%CY), and (3) 3 daily cheese yield traits (dCY) from goat milk. Cheese-making traits were analyzed using 2 mixed different models, the first to test MCP and the second to test CF_t parameters. Pearson correlations were also calculated. Significant and favorable relationships (negative for time intervals and positive for CF measures) were found between the traditional MCP and the CF_t parameters and %REC and %CY traits. The effects of milk fat and protein contents were particularly important on all cheese-making traits, with the only exception being the effect of fat content on water retention in cheese (%CY_WATER). We found an optimum value of milk k₂₀, associated with the highest recovery of components and cheese yield in solids (%CY_SOLIDS). In addition, a lower level of curd water retention and an increased fresh curd yield (%CY_CURD) were associated with greater recovery of fat. The collection of all available information during the process of milk coagulation and curd-firming allowed us to discover the effect of RCT_eq on %REC traits and %CY_SOLIDS, which had not previously been revealed for traditional RCT. Moreover, higher k_CF values were associated with increased %CY_CURD and %CY_SOLIDS. Given that CF_t parameters showed a high level of independence from one another, these can also be easily used and characterized in future applications at the industry level. Information provided by traditional and modeled coagulation properties could efficiently support the goat dairy industry and lay the foundations for a quality payment scheme for goat milk.

Coagulation Traits of Sheep and Goat Milk

Milk production from sheep and goat species is continuously growing worldwide, and its main use is for cheesemaking. Given that the final quality of cheese is linked to the traits of raw milk cheese yield at dairy plants, it is often calculated by using predictive formulas based on fat and protein content. Predictive formulas have been studied for bovine milk and are very effective but not appropriate for sheep and goat milk. Several methods, which simulate the actual coagulation processes, are available at the laboratories. This article reviews the available literature about rennet coagulation and cheese yield traits from sheep and goat milk and the methods used at the laboratory level. In general, if compared to cow milk, sheep and goat milk are characterized by shorter rennet coagulation times and a very limited amount of non-coagulating samples. Curd firmness of sheep milk is almost independent from the rennet coagulation time, and some coagulation traits can be predicted by infrared spectra. In addition, coagulation traits are characterized by appropriate values of heritability to be considered in selective breeding plans. With regard to goat milk, rennet coagulation time and cheese yield are strongly influenced by the breed effect.

Effects of indirect indicators of udder health on nutrient recovery and cheese yield traits in goat milk

In dairy goats, very little is known about the effect of the 2 most important indirect indicators of udder health [somatic cell count (SCC) and total bacterial count (TBC)] on milk composition and cheese yield, and no information is available regarding the effects of lactose levels, pH, and NaCl content on the recovery of nutrients in the curd, cheese yield traits, and daily cheese yields. Because large differences exist among dairy species, conclusions from the most studied species (i.e., bovine) cannot be drawn for all types of dairy-producing animals. The aims of this study were to quantify, using milk samples from 560 dairy goats, the contemporary effects of a pool of udder health indirect indicators (lactose level, pH, SCC, TBC, and NaCl content) on the recovery of nutrients in the curd (%REC), cheese yield (%CY), and daily cheese yields (dCY). Cheese-making traits were analyzed using a mixed model, with parity, days in milk (DIM), lactose level, pH, SCC, TBC, and NaCl content as fixed effects, and farm, breed, glass tube, and animal as random effects. Results indicated that high levels of milk lactose were associated with reduced total solids recovery in the curd and lower cheese yields, because of the lower milk fat and protein contents in samples rich in lactose. Higher pH correlated with higher recovery of nutrients in the curd and higher cheese yield traits. These results may be explained by the positive correlation between pH and milk fat, protein, and casein in goat milk. High SCC were associated with higher recovery of solids and energy in the curd but lower recovery of protein. The higher cheese yield obtained from milk with high SCC was due to both increased recovery of lactose in the curd and water retention. Bacterial count proved to be the least important factor affecting cheese-making traits, but it decreased daily cheese yields, suggesting that, even if below the legal limits, TBC should be considered in order to monitor flock management and avoid economic losses. The effect of NaCl content on milk composition was linked with lower recovery of all nutrients in the curd during cheese-making. In addition, high milk NaCl content led to reductions in fresh cheese yield and cheese solids. The indirect indicators of the present study significantly affected the cheese-making process. Such information should be considered, to adjust the milk-to-cheese economic value and the milk payment system.

Effect of goat milk composition on cheesemaking traits and daily cheese production

Cheese yield is strongly influenced by the composition of milk, especially fat and protein contents, and by the efficiency of the recovery of each milk component in the curd. The real effect of milk composition on cheesemaking ability of goat milk is still unknown. The aims of this study were to quantify the effects of milk composition; namely, fat, protein, and casein contents, on milk nutrient recovery in the curd, cheese yield, and average daily yield. Individual milk samples were collected from 560 goats of 6 different breeds. Each sample was analyzed in duplicate using the 9-laboratory milk cheesemaking assessment, a laboratory method that mimicked cheesemaking procedures, with milk heating, rennet addition, coagulation, curd cutting, and draining. Data were submitted to statistical analysis; results showed that the increase of milk fat content was associated with a large improvement of cheese yield because of the higher recovery of all milk nutrients in the curd, and thus a higher individual daily cheese yield. The increase of milk protein content affected the recovery of fat, total solids, and energy in the curd. Casein number, calculated as casein-to-protein ratio, did not affect protein recovery but strongly influenced the recovery of fat, showing a curvilinear pattern and the most favorable data for the intermediate values of casein number. In conclusion, increased fat and protein contents in the milk had an effect on cheese yield not only for the greater quantity of nutrients available but also for the improved efficiency of the recovery in the curd of all nutrients. These results are useful to improve knowledge on cheesemaking processes in the caprine dairy industry.

A Study on the Effects of Rumen Acidity on Rumination Time and Yield, Composition, and Technological Properties of Milk from Early Lactating Holstein Cows

Simple Summary

The increase in milk yield achieved in recent decades by the dairy sector has been sustained by feeding dairy cows with more concentrates and less forage. This leads to increasing rumen acidity, a status widespread in high-producing dairy cows that may affect feed intake, impair ruminal digestion, and cause diarrhea, laminitis, inflammation, and liver abscesses. The effects of rumen acidity on milk yield and composition are controversial, while those on milk coagulation properties and cheese yield have not yet been explored. This study investigated whether the rumen acidity status affects rumination time, and the production, composition, coagulation properties and cheese yield of milk obtained by 100 early-lactating Holstein cows. The variation in rumen acidity was associated with changes in the cows’ rumen fluid composition and circadian pattern of rumination time. Moreover, daily milk yield linearly decreased as the rumen acidity increased. Conversely, the composition and technological properties of milk were unaffected, even when there were differences in rumen acidity, suggesting that variation in rumen acidity has little impact on cheese-making traits.

Abstract

The use of high grain rations in dairy cows is related to an increase in rumen acidity. This study investigated whether the rumen acidity status affects rumination time (RT), and the production, composition, coagulation properties (MCPs) and cheese yield (CY) of milk. One hundred early-lactating Holstein cows with no clinical signs of disease and fed total mixed rations were used. Rumen fluid was collected once from each cow by rumenocentesis to determine pH and volatile fatty acid (VFA) content. The cows were classified according to the quartile of rumen acidity (QRA), a factor defined by multivariate analysis and associated with VFA and pH. Rumen fluid pH averaged 5.61 in the first quartile and 6.42 in the fourth, and total VFA content increased linearly with increasing rumen acidity. In addition, RT increased as rumen acidity increased, but only in the daily time interval from 08:00 to 12:00. Milk yield linearly decreased as rumen acidity increased, whereas QRA did not affect pH, fat or protein contents of milk. Furthermore, the MCPs, assessed by lactodynamograph, and CY were unaffected by QRA. It is suggested that differences in rumen acidity have little influence on the nutrient content, coagulation properties and CY of milk.