Variation of milk technological properties in sheep milk: Relationships among composition,coagulation and cheese-making traits

The relationships between milk composition, coagulation properties and cheese-making traits in sheep milk were characterised. Ten traits related to milk coagulation (RCT_eq, k_CF, CF_p), cheese yield (%CY_CURD, %CY_SOLIDS, %CY_WATER), and curd nutrients recovery or whey loss (%REC_FAT, %REC_PROTEIN, %REC_SOLIDS, %REC_ENERGY) were recorded. To obtain a measure of the efficiency in terms of %CY, the ratio between the observed and the theoretical %CY (Ef-%CY_CURD, Ef-%CY_SOLIDS) was calculated. Sheep milk showed good qualities for coagulation and cheese production; milk lactose appeared to be the component most linked to gelation, curd firming time and water retained in the curd. In the case of milk protein, an opposite relationship with gelation time was observed. Milk fat and protein positively affected total solids recovery and yield inducing higher %CY_CURD. Relationships with CF_t parameters were limited; curd firming instant rate seems to be the most informative trait to assess the efficiency of the cheese-making process.     

Associations between pathogen-specific cases of subclinical mastitis and milk yield,quality, protein composition,and cheese-making traits in dairy cows

The aim of this study was to investigate associations between pathogen-specific cases of subclinical mastitis and milk yield, quality, protein composition, and cheese-making traits. Forty-one multibreed herds were selected for the study, and composite milk samples were collected from 1,508 cows belonging to 3 specialized dairy breeds (Holstein Friesian, Brown Swiss, and Jersey) and 3 dual-purpose breeds of Alpine origin (Simmental, Rendena, and Grey Alpine). Milk composition [i.e., fat, protein, casein, lactose, pH, urea, and somatic cell count (SCC)] was analyzed, and separation of protein fractions was performed by reversed-phase high performance liquid chromatography. Eleven coagulation traits were measured: 5 traditional milk coagulation properties [time from rennet addition to milk gelation (RCT, min), curd-firming rate as the time to a curd firmness (CF) of 20 mm (k₂₀, min), and CF at 30, 45, and 60 min from rennet addition (a₃₀, a₄₅, and a₆₀, mm)], and 6 new curd firming and syneresis traits [potential asymptotical CF at an infinite time (CF_P, mm), curd-firming instant rate constant (k_CF, % × min^?1), curd syneresis instant rate constant (k_SR, % × min^?1), modeled RCT (RCT_eq, min), maximum CF value (CF_max, mm), and time at CF_max (t_max, min)]. We also measured 3 cheese yield traits, expressing the weights of total fresh curd (%CY_CURD), dry matter (%CY_SOLIDS), and water (%CY_WATER) in the curd as percentages of the weight of the processed milk, and 4 nutrient recovery traits (REC_PROTEIN, REC_FAT, REC_SOLIDS, and REC_ENERGY), representing the percentage ratio between each nutrient in the curd and milk. Milk samples with SCC > 100,000 cells/mL were subjected to bacteriological examination. All samples were divided into 7 clusters of udder health (UH) status: healthy (cows with milk SCC < 100,000 cells/mL and uncultured); culture-negative samples with low, medium, or high SCC; and culture-positive samples divided into contagious, environmental, and opportunistic intramammary infection (IMI). Data were analyzed using a linear mixed model. Significant variations in the casein to protein ratio and lactose content were observed in all culture-positive samples and in culture-negative samples with medium to high SCC compared to normal milk. No differences were observed among contagious, environmental, and opportunistic pathogens, suggesting an effect of inflammation rather than infection. The greatest impairment in milk quantity and composition, clotting ability, and cheese production was observed in the 2 UH status groups with the highest milk SCC (i.e., contagious IMI and culture-negative samples with high SCC), revealing a discrepancy between the bacteriological results and inflammatory status, and thus confirming the importance of SCC as an indicator of udder health and milk quality.     

Genetic parameters of different measures of cheese yield and milk nutrient recovery from an individual model cheese-manufacturing process

Cheese yield (CY) is an important technological trait in the dairy industry, and the objective of this study was to estimate the genetic parameters of cheese yield in a dairy cattle population using an individual model-cheese production procedure. A total of 1,167 Brown Swiss cows belonging to 85 herds were sampled once (a maximum of 15 cows were sampled per herd on a single test day, 1 or 2 herds per week). From each cow, 1,500 mL of milk was processed according to the following steps: milk sampling and heating, culture addition, rennet addition, gelation-time recording, curd cutting, whey draining and sampling, wheel formation, pressing, salting in brine, weighing, and cheese sampling. The compositions of individual milk, whey, and curd samples were determined. Three measures of percentage cheese yield (%CY) were calculated: %CY_CURD, %CY_SOLIDS, and %CY_WATER, which represented the ratios between the weight of fresh curd, the total solids of the curd, and the water content of the curd, respectively, and the weight of the milk processed. In addition, 3 measures of daily cheese yield (dCY, kg/d) were defined, considering the daily milk yield. Three measures of nutrient recovery (REC) were computed: REC_FAT, REC_PROTEIN, and REC_SOLIDS, which represented the ratio between the weights of the fat, protein, and total solids in the curd, respectively, and the corresponding nutrient in the milk. Energy recovery, REC_ENERGY, represented the energy content of the cheese versus that in the milk. For statistical analysis, a Bayesian animal model was implemented via Gibbs sampling. The effects of parity (1 to ≥4), days in milk (6 classes), and laboratory vat (15 vats) were assigned flat priors; those of herd-test-date, animal, and residual were given Gaussian prior distributions. Intra-herd heritability estimates of %CY_CURD, _%CY_SOLIDS, and %CY_WATER ranged from 0.224 to 0.267; these were larger than the estimates obtained for milk yield (0.182) and milk fat content (0.122), and similar to that for protein content (0.275). Daily cheese yields showed heritability estimates similar to those of daily milk yield. The trait %CY_WATER showed a highly positive genetic correlation with %CY_SOLIDS (0.87), whereas their phenotypic correlation was moderate (0.37), and the fat and protein contents of milk showed high genetic correlations with %CY traits. The heritability estimates of REC_PROTEIN and REC_FAT were larger (0.490 and 0.208, respectively) than those obtained for the protein and fat contents of milk, and the genetic relationships between REC_PROTEIN and REC_FAT with milk protein and fat content were low or moderate; REC_PROTEIN and REC_FAT were moderately correlated with the %CY traits and highly correlated with REC_SOLIDS and REC_ENERGY. Both REC_SOLIDS and REC_ENERGY were heritable (0.274 and 0.232), and showed high correlations with each other (0.96) and with the %CY traits (0.83 to 0.97). Together, these findings demonstrate the existence of economically important, genetically determined variability in cheese yield that does not depend solely upon the fat and protein contents of milk, but also relies on the ability of the coagulum to retain the highest possible proportions of the available protein, fat, and water. Exploitation of this interesting genetic variation does not seem to be feasible through direct measurement of the phenotype in cows at the population level. Instead, further research is warranted to examine possible means for indirect prediction, such as through assessing the mid-infrared spectra of milk samples.     

Genome-wide association study for cheese yield and curd nutrient recovery in dairy cows

Cheese production and consumption are increasing in many countries worldwide. As a result, interest has increased in strategies for genetic selection of individuals for technological traits of milk related to cheese yield (CY) in dairy cattle breeding. However, little is known about the genetic background of a cow's ability to produce cheese. Recently, a relatively large panel (1,264 cows) of different measures of individual cow CY and milk nutrient and energy recoveries in the cheese (REC) became available. Genetic analyses showed considerable variation for CY and for aptitude to retain high proportions of fat, protein, and water in the coagulum. For the dairy industry, these characteristics are of major economic importance. Nevertheless, use of this knowledge in dairy breeding is hampered by high costs, intense labor requirement, and lack of appropriate technology. However, in the era of genomics, new possibilities are available for animal breeding and genetic improvement. For example, identification of genomic regions involved in cow CY might provide potential for marker-assisted selection. The objective of this study was to perform genome-wide association studies on different CY and REC measures. Milk and DNA samples from 1,152 Italian Brown Swiss cows were used. Three CY traits expressing the weight (wt) of fresh curd (%CY_CURD), curd solids (%CY_SOLIDS), and curd moisture (%CY_WATER) as a percentage of weight of milk processed, and 4 REC (REC_FAT, REC_PROTEIN, REC_SOLIDS, and REC_ENERGY, calculated as the % ratio between the nutrient in curd and the corresponding nutrient in processed milk) were analyzed. Animals were genotyped with the Illumina BovineSNP50 Bead Chip v.2. Single marker regressions were fitted using the GenABEL R package (genome-wide association using mixed model and regression–genomic control). In total, 103 significant associations (88 single nucleotide polymorphisms) were identified in 10 chromosomes (2, 6, 9, 11, 12, 14, 18, 19, 27, 28). For REC_FAT and REC_PROTEIN, high significance peaks were identified in Bos taurus autosome (BTA) 6 and BTA11, respectively. Marker ARS-BFGL-NGS-104610 (～104.3 Mbp) was highly associated with REC_PROTEIN and Hapmap52348-rs29024684 (～87.4 Mbp), closely located to the casein genes on BTA6, with REC_FAT. Genomic regions identified may enhance marker-assisted selection in bovine cheese breeding beyond the use of protein (casein) and fat contents, whereas new knowledge will help to unravel the genomic background of a cow's ability for cheese production.     

Genetic analysis of coagulation properties,curd firming modeling,milk yield,composition, and acidity in Sarda dairy sheep

Sheep milk is an important source of food, especially in Mediterranean countries, and is used in large part for cheese production. Milk technological traits are important for the sheep dairy industry, but research is lacking into the genetic variation of such traits. Therefore the aim of this study was to estimate the heritability of traditional milk coagulation properties and curd firmness modeled on time t (CF_t) parameters, and their genetic relationships with test-day milk yield, composition (fat, protein, and casein content), and acidity in Sarda dairy sheep. Milk samples from 1,121 Sarda ewes from 23 flocks were analyzed for 5 traditional coagulation properties by lactodynamographic tests conducted for up to 60 min: rennet coagulation time (min), curd-firming time (k₂₀, min), and 3 measures of curd firmness (a₃₀, a₄₅, and a₆₀, mm). The 240 curd firmness observations (1 every 15 s) from each milk sample were recorded, and 4 parameters for each individual sample equation were estimated: rennet coagulation time estimated from the equation (RCT_eq), the asymptotic potential curd firmness (CF_P), the curd firming instant rate constant (k_CF), and the syneresis instant rate constant (k_SR). Two other derived traits were also calculated (CF_max, the maximum curd firmness value; and t_max, the attainment time). Multivariate analyses using Bayesian methodology were performed to estimate the genetic relationships of milk coagulation properties and CF_t with the other traits; statistical inference was based on the marginal posterior distributions of the parameters of concern. The marginal posterior distribution of heritability estimates of milk yield (0.16 ± 0.07) and composition (0.21 ± 0.11 to 0.28 ± 0.10) of Sarda ewes was similar to those often obtained for bovine species. The heritability of rennet coagulation time as a single point trait was also similar to that frequently obtained for cow milk (0.19 ± 0.09), whereas the same trait calculated as an individual equation parameter exhibited larger genetic variation and a higher heritability estimate (0.32 ± 0.11). The other curd firming and syneresis traits, whether as traditional single point observations or as individual equation parameters and derived traits, were characterized by heritability estimates lower than for coagulation time and for the corresponding bovine milk traits (0.06 to 0.14). Phenotypic and additive genetic correlations among the 11 technological traits contribute to describing the interdependencies and meanings of different traits. The additive genetic relationships of these technological traits with the single test-day milk yield and composition were variable and showed milk yield to have unfavorable effects on all measures of curd firmness (a₃₀, a₄₅, a₆₀, CF_P, and CF_max) and t_max, but favorable effects on both instant rate constants (k_CF and k_SR). Milk fat content had a positive effect on curd firmness traits, especially on those obtained from CF_t equations, whereas the negative effects on both coagulation time traits were attributed to the milk protein and casein contents. Finally, in view of the estimated heritabilities and additive genetic correlations, enhancement of technological traits of sheep milk through selective breeding could be feasible in this population.     

Factors affecting variation of different measures of cheese yield and milk nutrient recovery from an individual model cheese-manufacturing process

Cheese yield (CY) is the most important technological trait of milk, because cheese-making uses a very high proportion of the milk produced worldwide. Few studies have been carried out at the level of individual milk-producing animals due to a scarcity of appropriate procedures for model-cheese production, the complexity of cheese-making, and the frequent use of the fat and protein (or casein) contents of milk as a proxy for cheese yield. Here, we report a high-throughput cheese manufacturing process that mimics all phases of cheese-making, uses 1.5-L samples of milk from individual animals, and allows the simultaneous processing of 15 samples per run. Milk samples were heated (35°C for 40 min), inoculated with starter culture (90 min), mixed with rennet (51.2 international milk-clotting units/L of milk), and recorded for gelation time. Curds were cut twice (10 and 15 min after gelation), separated from the whey, drained (for 30 min), pressed (3 times, 20 min each, with the wheel turned each time), salted in brine (for 60 min), weighed, and sampled. Whey was collected, weighed, and sampled. Milk, curd, and whey samples were analyzed for pH, total solids, fat content, and protein content, and energy content was estimated. Three measures of percentage cheese yield (%CY) were calculated: %CY_CURD, %CY_SOLIDS, and %CY_WATER, representing the ratios between the weight of fresh curd, the total solids of the curd, and the water content of the curd, respectively, and the weight of the milk processed. In addition, 3 measures of daily cheese yield (dCY, kg/d) were defined, considering the daily milk yield. Three measures of nutrient recovery (REC) were computed: REC_FAT, REC_PROTEIN, and REC_SOLIDS, which represented the ratio between the weights of the fat, protein, and total solids in the curd, respectively, and the corresponding components in the milk. Energy recovery, REC_ENERGY, represented the energy content of the cheese compared with that in the milk. This procedure was used to process individual milk samples obtained from 1,167 Brown Swiss cows reared in 85 herds of the province of Trento (Italy). The assessed traits exhibited almost normal distributions, with the exception of REC_FAT. The average values (± SD) were as follows: %CY_CURD = 14.97 ± 1.86, %CY_SOLIDS = 7.18 ± 0.92, %CY_WATER = 7.77 ± 1.27, dCY_CURD = 3.63 ± 1.17, dCY_SOLIDS = 1.74 ± 0.57, dCY_WATER = 1.88 ± 0.63, REC_FAT = 89.79 ± 3.55, REC_PROTEIN = 78.08 ± 2.43, REC_SOLIDS = 51.88 ± 3.52, and REC_ENERGY = 67.19 ± 3.29. All traits were highly influenced by herd-test-date and days in milk of the cow, moderately influenced by parity, and weakly influenced by the utilized vat. Both %CY_CURD and dCY_CURD depended not only on the fat and protein (casein) contents of the milk, but also on their proportions retained in the curd; the water trapped in curd presented an higher variability than that of %CY_SOLIDS. All REC traits were variable and affected by days in milk and parity of the cows. The described model cheese-making procedure and the results obtained provided new insight into the phenotypic variation of cheese yield and recovery traits at the individual level.     

Breed of cow and herd productivity affect milk composition and modeling of coagulation,curd firming,and syneresis

Milk coagulation properties (MCP) have been widely investigated in the past using milk collected from different cattle breeds and herds. However, to our knowledge, no previous studies have assessed MCP in individual milk samples from several multi-breed herds characterized by either high or low milk productivity, thereby allowing the effects of herd and cow breed to be evaluated independently. Multi-breed herds (n = 41) were classified into 2 categories based on milk productivity (high vs. low), defined according to the average milk net energy yielded daily by lactating cows. Milk samples were taken from 1,508 cows of 6 different breeds: 3 specialized dairy (Holstein-Friesian, Brown Swiss, Jersey) and 3 dual-purpose (Simmental, Rendena, Alpine Grey) breeds, and analyzed in duplicate (3,016 tests) using 2 lactodynamographs to obtain 240 curd firming (CF) measurements over 60 min (1 every 15 s) for each duplicate. The 5 traditional single-point MCP (RCT, k₂₀, a₃₀, a₄₅, and a₆₀) were yielded directly by the instrument from the available CF measures. All 240 CF measures of each replicate were also used to estimate 4 individual equation parameters: RCT estimated according to curd firm change over time modeling (RCT_eq), asymptotic potential curd firmness (CF_P), curd firming instant rate constant (k_CF), and syneresis instant rate constant (k_SR) and 2 derived traits: maximum curd firmness achieved within 45 min (CF_max) and time at achievement of CF_max (t_max) by curvilinear regression using a nonlinear procedure. Results showed that the effect of herd-date on traditional and modeled MCP was modest, ranging from 6.1% of total variance for k₂₀ to 10.7% for RCT, whereas individual animal variance was the highest, ranging from 32.0% for t_max to 82.5% for RCT_eq. The repeatability of MCP was high (>80%) for all traits except those associated with the last part of the lactodynamographic curve (i.e., a₆₀, k_SR, k_CF, and t_max: 57 to 71%). Reproducibility, taking into account the effect of instrument, was equal to or slightly lower than repeatability. Milk samples collected in farms characterized by high productivity exhibited delayed coagulation (RCT_eq: 18.6 vs. 16.3 min) but greater potential curd firmness (CF_P: 76.8 vs. 71.9 mm) compared with milk samples collected from low-productivity herds. Parity and days in milk influenced almost all MCP. Large differences in all MCP traits were observed among breeds, both between specialized and dual-purpose breeds and within these 2 groups of breeds, even after adjusting for milk quality and yield. Milk quality and MCP of samples from Jersey cows, and coagulation time of samples from Rendena cows were better than in milk from Holstein-Friesian cows, and intermediate results were found with the other breeds of Alpine origin. The results of this study, taking into account the intrinsic limitation of this technique, show that the effects of breed on traditional and modeled MCP are much greater than the effects of herd productivity class, parity, and DIM. Moreover, the variance in individual animals is much greater than the variance in individual herds within herd productivity class. It seems that improvement in MCP depends more on genetics (e.g., breed, selection) than on environmental and management factors.     

Milk protein fractions strongly affect the patterns of coagulation,curd firming,and syneresis

The aim of this study was to assess the role of milk protein fractions in the coagulation, curd firming, and syneresis of bovine milk. Analyses were performed on 1,271 individual milk samples from Brown Swiss cows reared in 85 herds classified into 4 types of farming systems, from the very traditional (tied cows, feed manually distributed, summer highland pasture) to the most modern (loose cows, use of total mixed rations with or without silage). Fractions α_S1-casein (CN), α_S2-CN, β-CN, κ-CN, β-lactoglobulin (LG), and α-lactalbumin (LA) and genotypes at CSN2, CSN3, and BLG were obtained by reversed-phase HPLC. The following milk coagulation properties were measured with a lactodynamograph, with the testing time extended to 60 min: rennet coagulation time (RCT, min), curd firming time (min), and curd firmness at 30 and 45 min (mm). All the curd firmness measures recorded over time (total of 240 observations/sample) were used in a 4-parameter nonlinear model to obtain parameters of coagulation, curd firming, and syneresis: RCT estimated from the equation (min), asymptotic potential curd firmness (mm), the curd firming and syneresis instant rate constants (%/min), and the maximum curd firmness value (CF_max, mm) and the time taken to reach it (min). All the aforementioned traits were analyzed with 2 linear mixed models, which tested the effects of the protein fractions expressed in different ways: in the first, quantitative model, each protein fraction was expressed as content in milk; in the second, qualitative model, each protein fraction was expressed as a percentage of total casein content. Besides proteins, additional nuisance parameters were herd (included as a random effect), daily milk production (only for the quantitative model), casein content (only for the qualitative model), dairy system, parity, days in milk, the pendulum of the lactodynamograph, and the CSN2, CSN3, and BLG genotypes. Both α_S1-CN and β-CN showed a clear and favorable effect on CF_max, where the former effect was almost double the latter. Milk coagulation ability was favorably affected by κ-CN, which reduced both the RCT and RCT estimated from the equation, increased the curd firming and syneresis instant rate constants, and allowed a higher CF_max to be reached. In contrast, α_S2-CN delayed gelation time and β-LG worsened curd firming, both resulting in a low CF_max. The results of this study suggest that modification of the relative contents of specific protein fractions can have an enormous effect on the technological behavior of bovine milk.     

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.     

Modeling rennet coagulation time and curd firmness of milk

Milk coagulation properties (MCP) are traditionally expressed using rennet coagulation time (RCT), time to curd firmness (CF) of 20 mm (k₂₀), and CF 30 min after enzyme addition (a₃₀) values, all of which are single-point measures taken from the output of computerized renneting meters, such as the Formagraph. Thus, traditional MCP use only some of the available information. Moreover, because of the worldwide spreading of breeds such as the Holstein-Friesian, characterized by late-coagulating milk, it happens often that some samples do not coagulate at all, that a₃₀ is strongly and negatively related to RCT, and that k₂₀ is not measurable. The aim of the present work was to model CF as a function of time (CF_t, mm) over a 30-min interval. The model tested was CFt=CFP×(1−e−kCF×(t−RCT)), where CF_P (mm) is the potential asymptotical CF at an infinite time, k_CF (min⁻¹) is the curd firming rate constant, and RCT is measured in minutes. The CF_t model was initially applied to data of milk of each of 105 Brown Swiss cows from 7 herds, each sampled once (trial 1). Four samples did not coagulate within 30 min. Eighty-seven of the 101 individual equations obtained fit the CF data of milk samples very well, even though the samples differed in composition, and were produced by cows of different ages and days in milk, reared on different farms (coefficient of determination >0.99; average residual standard deviation = 0.21 mm). Samples with a very late RCT (slowly coagulating samples) yielded so few observational data points that curve parameters could not be precisely estimated. The repeatability of CF_t equation parameters was estimated using data obtained from 5 replicates of each of 2 samples of bulk milk from 5 Holstein-Friesian cows analyzed every day for 5 consecutive days (trial 2). Repeatability of RCT was better than that of the other 2 parameters. Moreover, traditional MCP values (RCT, a₃₀, and k₂₀) can be obtained from the individual CF_t equations, using all available information. The MCP estimated from equations were very similar to the single-point measures yielded by the computerized renneting meter (coefficient of determination >0.97), but repeatability was slightly better. The model allowed the estimation of k₂₀ for samples with a very late coagulation or with very slow curd firming. Finally, the 3 novel parameters used to assess different milk samples were less interdependent than are the traditional measures, and their practical and scientific utility requires further study.     

Invited review: Genetics and modeling of milk coagulation properties

Milk coagulation properties (MCP) are conventionally measured using computerized renneting meters, mechanical or optical devices that record curd firmness over time (CF_t). The traditional MCP are rennet coagulation time (RCT, min), curd firmness (a₃₀, mm), and curd-firming time (k₂₀, min). The milk of different ruminant species varies in terms of CF_t pattern. Milk from Holstein-Friesian and some Scandinavian cattle breeds yields higher proportions of noncoagulating samples, samples with longer RCT and lower a₃₀, and samples for which k₂₀ is not estimable, than does milk from Brown Swiss, Simmental, and other local Alpine breeds. The amount, proportion, and genetic variants (especially κ-casein) of milk protein fractions strongly influence MCP and explain variable proportions of the observed differences among breeds and among individuals of the same breed. In addition, other major genes have been shown to affect MCP. Individual repeatability of MCP is high, whereas any herd effect is low; thus, the improvement of MCP should be based principally on selection. Exploitable additive genetic variation in MCP exists and has been assessed using different breeds in various countries. Several models have been formulated that either handle noncoagulating samples or not. The heritability of MCP is similar to that of other milk quality traits and is higher than the heritability of milk yield. Rennet coagulation time and a₃₀ are highly correlated, both phenotypically and genetically. This means that the use of a₃₀ data does not add valuable information to that obtainable from RCT; both traits are genetically correlated mainly with milk acidity. Moreover, a₃₀ is correlated with casein content. The major limitations of traditional MCP can be overcome by prolonging the observation period and by using a novel CF_t modeling, which uses all available information provided by computerized renneting meters and allows the estimation of RCT, the potential asymptotic curd firmness, the curd-firming rate, and the syneresis rate. Direct measurements of RCT obtained from both mechanical and optical devices show similar heritabilities and exhibit high phenotypic and genetic correlations. Moreover, mid-infrared reflectance spectroscopy can predict MCP. The heritabilities of predicted MCP are higher than those of measured MCP, and the 2 sets of values are strongly correlated. Therefore, mid-infrared reflectance spectroscopy is a reliable and cheap method whereby MCP can be improved at the population level; this is because such spectra are already routinely acquired from the milk of cows enrolled in milk recording schemes.     

Variations in milk protein fractions affect the efficiency of the cheese-making process

The aim of this study was to assess the influence of the amounts of the α_S1-, α_S2-, β-, and κ-casein (CN) and the α-lactalbumin and β-lactoglobulin protein fractions on the efficiency of the cheese-making process independently of their genetic polymorphisms. The study was carried out on milk samples from 1,271 Brown Swiss cows from 85 herds classified into 4 categories according to management, feeding, and housing characteristics (traditional and modern systems). To assess the efficiency of the cheese-making process, we processed the milk samples according to a laboratory cheese-making procedure (1,500 mL/sample) and obtained the following measures: (1) 3 percentage cheese yields (%CY_curd, %CY_solids, %CY_water), (2) 2 daily cheese yields obtained by multiplying %CY (curd and total solids) by daily milk yields (dCY_curd, dCY_solids), (3) 4 measures of nutrient recovery in the curd (REC_fat, REC_protein, REC_solids, REC_energy), and (4) 2 measures of cheese-making efficiency in terms of the ratio between the observed and theoretical %CY (Ef-%CY_curd, Ef-%CY_solids). All the aforementioned traits were analyzed by fitting 2 linear mixed models with protein fractions as fixed effects expressed as percentage in the milk (model M-%milk) and as percentage of the total casein content (model M-%cas) together with the effects of total casein content (only in model M-%cas), daily milk yield (only in model M-%milk; not for dCY traits), dairy system, herd (random effect), days in milk, parity, and vat. The efficiency of overall cheese yield (Ef-%CY_curd) was mostly positively associated with β-CN content in the milk, whereas Ef-%CY_solids was greater with higher amounts of κ-CN and α_S1-CN (M-%milk) due to the strong influence of both fractions on the recovery rate of milk components in the curd (fat and total solids, protein with α_S1-CN only) when expressed as percentage of milk and of total casein; only β-CN was more important for REC_protein. In contrast, we found β-lactoglobulin to be highly negatively related to all the traits related to the cheese-making process and to the daily cheese yield per cow, whereas α-lactalbumin was positively associated with the latter traits. Additional research on this topic is needed, with particular focus on the genetic and genomic aspects of the role of protein fractions in the cheese-making process and on the associations between genetic polymorphisms in milk protein and milk nutrient recovery in the curd.     

Goat farm variability affects milk Fourier-transform infrared spectra used for predicting coagulation properties

Driven by the large amount of goat milk destined for cheese production, and to pioneer the goat cheese industry, the objective of this study was to assess the effect of farm in predicting goat milk-coagulation and curd-firmness traits via Fourier-transform infrared spectroscopy. Spectra from 452 Sarda goats belonging to 14 farms in central and southeast Sardinia (Italy) were collected. A Bayesian linear regression model was used, estimating all spectral wavelengths' effects simultaneously. Three traditional milk-coagulation properties [rennet coagulation time (min), time to curd firmness of 20 mm (min), and curd firmness 30 min after rennet addition (mm)] and 3 curd-firmness measures modeled over time [rennet coagulation time estimated according to curd firmness change over time (RCT_eq), instant curd-firming rate constant, and asymptotical curd firmness] were considered. A stratified cross validation (SCV) was assigned, evaluating each farm separately (validation set; VAL) and keeping the remaining farms to train (calibration set) the statistical model. Moreover, a SCV, where 20% of the goats randomly taken (10 replicates per farm) from the VAL farm entered the calibration set, was also considered (SCV₈₀). To assess model performance, coefficient of determination (R²_VAL) and the root mean squared error of validation were recorded. The R²_VAL varied between 0.14 and 0.45 (instant curd-firming rate constant and RCT_eq, respectively), albeit the standard deviation was approximating half of the mean for all the traits. Although average results of the 2 SCV procedures were similar, in SCV₈₀, the maximum R²_VAL increased at about 15% across traits, with the highest observed for time to curd firmness of 20 mm (20%) and the lowest for RCT_eq (6%). Further investigation evidenced important variability among farms, with R²_VAL for some of them being close to 0. Our work outlined the importance of considering the effect of farm when developing Fourier-transform infrared spectroscopy prediction equations for coagulation and curd-firmness traits in goats.     

Technical note: Improving modeling of coagulation,curd firming,and syneresis of sheep milk

The importance of milk coagulation properties for milk processing, cheese yield, and quality is widely recognized. The use of traditional coagulation traits presents several limitations for testing bovine milk and even more for sheep milk, due to its rapid coagulation and curd firming, and early syneresis of coagulum. The aim of this technical note is to test and improve model fitting for assessing coagulation, curd firming, and syneresis of sheep milk. Using milk samples from 87 Sarda ewes, we performed in duplicate lactodynamographic testing. On each of the 174 analyzed milk aliquots, using 180 observations from each aliquot (one every 15 s for 45 min after rennet addition), we compared 4 different curd firming models as a function of time (CF_t, mm) using a nonlinear procedure. The most accurate and informative results were observed using a modified 4-parameter model, structured as follows: ${\text{CF}}_t={\text{CF}}_{\text{P}}\times \left(1?{\text{e}}^{?{\text{k}}_{\text{CF}}\left({\text{RCT}}_{\text{eq}}\right)}\right)\times {\text{e}}^{{\text{k}}_{\text{SR}}\times \left({\text{t-RCT}}_{\text{eq}}\right)}$ where t is time, RCT_eq (min) is the gelation time, CF_P (mm) is the potential asymptotical CF at an infinite time, k_CF (%/min) is the curd firming rate constant, and k_SR (%/min) is the curd syneresis rate constant. To avoid nonconvergence and computational problems due to interrelations among the equation parameters, CF_P was preliminarily defined as a function of maximum observed curd firmness (CF_max, mm) recorded during the analysis. For this model, all the modeling equations of individual sheep milk aliquots were converging, with a negligible standard error of the estimates (coefficient of determination >0.99 for all individual sample equations). Repeatability of the modeled parameters was acceptable, also in the presence of curd syneresis during the lactodynamographic analysis.     

Milk yield,quality, and coagulation properties of 6 breeds of goats: Environmental and individual variability

Goat milk and cheese production is continuously increasing and milk composition and coagulation properties (MCP) are useful tools to predict cheesemaking aptitude. The present study was planned to investigate the extension of lactodynamographic analysis up to 60 min in goat milk, to measure the farm and individual factors, and to investigate differences among 6 goat breeds. Daily milk yield (dMY) was recorded and milk samples collected from 1,272 goats reared in 35 farms. Goats were of 6 different breeds: Saanen and Camosciata delle Alpi for the Alpine type, and Murciano-Granadina, Maltese, Sarda, and Sarda Primitiva for the Mediterranean type. Milk composition (fat, protein, lactose, pH; somatic cell score; logarithmic bacterial count) and MCP [rennet coagulation time (RCT, min), curd-firming time (k₂₀, min), curd firmness at 30, 45, and 60 min after rennet addition (a₃₀, a₄₅, and a₆₀, mm)] were recorded, and daily fat and protein yield (dFPY g/d) was calculated as the sum of fat and protein concentration multiplied by the dMY. Data were analyzed using different statistical models to measure the effects of farm, parity_, stage of lactation and breed; lastly, the direct and the indirect effect of breed were quantified by comparing the variance of breed from models with or without the inclusion of linear regression of fat, protein, lactose, pH, bacterial, somatic cell counts, and dMY. Orthogonal contrasts were performed to compare least squares means. Almost all traits exhibited high variability, with coefficients of variation between 32 (for RCT) and 63% (for a₃₀). The proportion of variance regarding dMY, dFPY, and milk composition due to the farm was moderate, whereas for MCP it was low, except for a₆₀, at 69%. Parity affected both yield and quality traits of milk, with least squares means of dMY and dFPY showing an increase and RCT and curd firmness traits a decrease from the first to the last parity class. All milk quality traits, excluding fat, were affected by the stage of lactation; RCT and k₂₀ decreased rapidly and a₃₀ was higher from the first to the last part of lactation. Alpine breeds showed the highest dMY and dFPY but Mediterranean the best percentage of protein, fat, and lactose and a shorter k₂₀ and a greater a₃₀. Among the Mediterranean goats, Murciano-Granadina goats had the highest milk yield, fat, and protein contents, whereas Maltese, Sarda, and Sarda Primitiva were characterized by much more favorable technological properties in terms of k₂₀, a₃₀, and a₄₅. In conclusion, as both the farm and individual factors highly influenced milk composition and MCP traits, improvements of these traits should be based both on modifying management and individual goat factors. As expected, several differences were attributable to the breed effect, with the best milk production for the Alpines and milk quality and coagulation for the Mediterranean goats.     

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.     

Genetic parameters for milk coagulation properties in Estonian Holstein cows

The objective of this study was to estimate heritabilities and repeatabilities for milk coagulation traits [milk coagulation time (RCT) and curd firmness (E₃₀)] and genetic and phenotypic correlations between milk yield and composition traits (milk fat percentage and protein percentage, urea, somatic cell count, pH) in first-lactation Estonian Holstein dairy cattle. A total of 17,577 test-day records from 4,191 Estonian Holstein cows in 73 herds across the country were collected during routine milk recordings. Measurements of RCT and E₃₀ determined with the Optigraph (Ysebaert, Frepillon, France) are based on an optical signal in the near-infrared region. The cows had at least 3 measurements taken during the period from April 2005 to January 2009. Data were analyzed using a repeatability animal model. There was substantial variation in milk coagulation traits with a coefficient of variation of 27% for E₃₀ and 9% for the log-transformed RCT. The percentage of variation explained by herd was 3% for E₃₀ and 4% for RCT, suggesting that milk coagulation traits are not strongly affected by herd conditions (e.g., feeding). Heritability was 0.28 for RCT and 0.41 for E₃₀, and repeatability estimates were 0.45 and 0.50, respectively. Genetic correlation between both milk coagulation traits was negligible, suggesting that RCT and E₃₀ have genetically different foundations. Milk coagulation time had a moderately high positive genetic (0.69) and phenotypic (0.61) correlation with milk pH indicating that a high pH is related to a less favorable RCT. Curd firmness had a moderate positive genetic (0.48) and phenotypic (0.45) correlation with the protein percentage. Therefore, a high protein percentage is associated with favorable curd firmness. All reported genetic parameters were statistically significantly different from zero. Additional univariate random regression analysis for milk coagulation traits yielded slightly higher average heritabilities of 0.38 and 0.47 for RCT and E₃₀ compared with the heritabilities of the repeatability model.     

Milk quality,coagulation properties,and curd firmness modeling of purebred Holsteins and first- and second-generation crossbred cows from Swedish Red,Montbéliarde,and Brown Swiss bulls

The objective of the present study was to investigate how the crossbreeding of Holstein (HO) cows with bulls from Nordic and Alpine European breeds affect milk quality traits, traditional milk coagulation properties (MCP), and curd firmness modeling obtained from individual milk samples. A total of 506 individual milk samples were collected from evening milking at 3 commercial farms located in Northern Italy. Over the past decade, the 3 farms have followed crossbreeding programs in part of their herds, whereas the remainder of the animals consisted of purebred HO. The basic scheme was a 3-breed rotation based on the use of Swedish Red (SR) semen on HO cows (SR × HO), the use of Montbéliarde (MO) semen on first-cross cows [MO × (SR × HO)], and the use of HO semen in the third cross. In all herds, a smaller proportion of purebred HO were mated to M and Brown Swiss (BS) bulls, and these first crosses were mated to SR and MO bulls, respectively. Milk samples were analyzed for milk composition and MCP, and parameters for curd firmness were modeled. Compared with purebred HO, crossbred cows produced less milk with lower lactose content, higher fat and protein content, and a tendency for higher casein content. Crossbred cows generally produced milk with a more favorable curd-firming rate (k₂₀) and curd firmness 30 min after rennet addition, among traditional MCP, and better trends of curd firmness measures as shown by model parameters: estimated rennet coagulation time, asymptotical potential value of curd firmness, and curd-firming instant rate constant. Among crossbred cows, SR × HO presented longer rennet coagulation time compared with MO × HO and BS × HO cows, and MO × HO showed shorter k₂₀ compared with BS × HO cows. Among second-generation cows, those sired by SR bulls showed a lower incidence of noncoagulated samples, higher curd firmness 30 min after rennet addition and asymptotical potential value of curd firmness, and faster curd-firming instant rate constant compared with animals sired by MO bulls. Our results revealed that different sire breeds were characterized by specific technological aptitudes, but that these were not strictly related to other milk quality traits. Furthermore, the favorable characteristics (in terms of the quality and technological properties of milk) could be maintained in the third generation of 3-way crosses without negative effects on milk yield, even though the HO heritage had been reduced from 50 to 25%. Our findings, therefore, suggest that different types of sires can be chosen (depending on the intended use of the milk) to ensure the optimization of farm crossbreeding programs.     

Milk coagulation ability of five dairy cattle breeds

Samples of herd milk (506) were analyzed to assess sources of variation for milk coagulation properties (MCP) for 5 different dairy cattle breeds. Data were recorded in 55 single-breed dairy herds in the Trento province, a mountain area in northeast Italy. The 5 cattle breeds were Holstein-Friesian (8 herds), Brown Swiss (16 herds), Simmental (10 herds), Rendena (13 herds), and Alpine Gray (8 herds). Herd milk samples were analyzed for the MCP traits, milk rennet coagulation time (RCT), curd-firming time, and curd firmness (a₃₀), as well as protein and fat percentages, somatic cell count, Soxhlet-Henkel acidity, and bacterial count. An ANOVA was performed to study the effect of breed, herd within breed, DIM, month of lactation, protein and fat percentages, somatic cell score, titratable acidity, and log bacterial count within breed on MCP. Breed was the most important source of variation. In particular, the Rendena breed showed the best MCP traits at 13.5 min and 27.0 mm for RCT and a₃₀, respectively. The Holstein-Friesian breed had the worst coagulation properties at 18.0 min and 17.5 mm for RCT and a₃₀, respectively. The other 3 breeds showed intermediate coagulation properties. The RCT values were better at the beginning of lactation, whereas RCT and a₃₀ values were better in September and October (14.3 min and 25.7 mm, respectively). Among the composition traits, only the titratable acidity affected MCP traits of herd milk positively.     

Phenotypic factors affecting coagulation properties of milk from Sarda ewes

In this study, milk-coagulation properties (MCP) were characterized in the Sarda sheep breed. Milk composition and MCP [rennet-coagulation time (RCT), curd-firming time [time to reach a curd firmness of 20 mm (k₂₀)], and curd firmness (a₃₀), (a₄₅), and (a₆₀)] were obtained extending the lactodynamographic analysis from 30 to 60 min from a population of 1,121 ewes from 23 different farms. Managerial characteristics of farms and parity, individual daily milk yields and stage of lactation of ewes were recorded. Data were analyzed using a mixed-model procedure with fixed effects of days in milk, parity, daily milk yield, and flock size and the random effect of the flock/test day nested within flock size. Sampled farms were classified as small (<300 ewes) and medium (300 to 600 ewes), and these were kept by family operations, or as large (>600 ewes), often operated through hired workers. Daily milk yield was, on average, 1.58 ± 0.79 L/d and variability for this trait was very high. The average content of fat, protein, and casein was respectively 6.41, 5.39, and 4.20%. The class of flock size had a significant effect only on curd firmness, whereas days in milk affected RCT and k₂₀. The flock test day, parity, and daily milk yield were important sources of variation for all MCP. The mean value of RCT (8.6 min) and the low occurrence of noncoagulating samples (0.44%) confirmed the excellent coagulation ability of sheep milk compared with cattle milk. A more rapid coagulation was observed in mid-lactating, primiparous, and high-yielding ewes. The k₂₀ was usually reached in less than 2 min after gelation, with the most favorable values at mid lactation. The mean value of curd firmness 30 min after rennet addition (a₃₀) was, on average, 50 mm and decreased to 46 and 42 mm respectively after 45 (a₄₅) and 60 min (a₆₀). The decreasing value of curd-firmness traits was likely to be caused by curd syneresis and whey expulsion. The correlation between RCT and a₃₀ was much lower than in dairy cows and about null for a₄₅ and a₆₀. This means that curd firmness in dairy ewes is almost independent of gelation time and this can provide specific information for this species. In conclusion, this study showed that milk from Sarda sheep is characterized by an earlier gelation, a faster increase in curd firmness with time, and greater curd firmness after 30 min compared with dairy cows. Furthermore, correlations between MCP in sheep are much lower than in bovines and some of the assumptions and interpretations related to cows cannot be applied to sheep.