Prediction accuracies of cheese-making traits using Fourier-transform infrared spectra in goat milk.

The objectives of this study were to explore the use of Fourier-transform infrared (FITR) spectroscopy on 458 goat milk samples for predicting cheese-making traits, and to test the effect of the farm variability on their prediction accuracy. Calibration equations were developed using a Bayesian approach with three different scenarios: i) a random cross-validation (CV) [80% calibration (CAL); 20% validation (VAL) set], ii) a stratified CV [(SCV), 13 farms used as CAL, and the remaining one as VAL set], and iii) a SCV where 20% of the goats randomly selected from the VAL farm were included in the CAL set (SCV80). The best prediction performance was obtained for cheese yield solids, justifying for its practical application at population level. Overall results were similar to or outperformed those reported for bovine milk. Our results suggest considering specific procedures for calibration development to propose reliable tools applicable along the dairy goat chain.

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

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

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 (%CYCURD), total solids (%CYSOLIDS), and water retained (%CYWATER) 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 (R2VAL) and the root mean square error of validation (RMSEVAL). In CrV, the formula with the best prediction accuracy for all %CY traits included fat, casein, and UHI (R2VAL = 0.65, 0.96, and 0.23 for %CYCURD, %CYSOLIDS, and %CYWATER, respectively). The WAV procedure showed R2VAL 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 %CYCURD and %CYWATER 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.

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 (RCTeq), 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 (SCV80). To assess model performance, coefficient of determination (R2VAL) and the root mean squared error of validation were recorded. The R2VAL varied between 0.14 and 0.45 (instant curd-firming rate constant and RCTeq, 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 SCV80, the maximum R2VAL increased at about 15% across traits, with the highest observed for time to curd firmness of 20 mm (20%) and the lowest for RCTeq (6%). Further investigation evidenced important variability among farms, with R2VAL 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.

Test positivity for Maedi-Visna virus and Mycobacterium avium ssp. paratuberculosis in Sarda ewes: Effects on milk composition and coagulation traits and heritability estimates for susceptibility.

Maedi-Visna virus (MVV) and Mycobacterium avium ssp. paratuberculosis (MAP) are two pathogens that cause chronic, production-limiting diseases in dairy sheep. Although they are present worldwide, there are no detailed reports on their actual effects on milk traits in the literature. This study was designed to investigate the effects of test positivity to MVV and MAP on ovine milk yield, composition and coagulation properties, and curd-firming over time (CFt) variables in clinically healthy animals at the field level. The additive genetic variation and heritabilities of MVV and MAP positivity were also estimated. Milk samples were collected from 1,079 Sarda sheep kept on 23 farms, and pedigree information was obtained from the flock book. Milk yield was also recorded on the sampling date. Positivity for MVV and MAP was determined from milk samples using indirect ELISA test kits. Milk composition traits were measured by spectroscopy, milk coagulation properties were measured with a Formagraph (Foss Italia, Padua, Italy), and CFt traits were calculated using the data from the Formagraph diagram. The effects of MVV and MAP positivity on milk traits were determined through a set of mixed linear models, which took into account various sources of variation, such as days in milk, parity, and flock effects, and included the effects (positive or negative) of the 2 pathogens. A Bayesian threshold sire model with sire relationship was used to estimate genetic variation and heritability. The overall animal prevalence of MVV-positive ewes was 43.6%; on only 1 farm of the 23 tested were all sampled ewes negative. An overall animal prevalence of 10.6% was recorded for MAP, with 4 farms at 0%. Positivity for MVV significantly affected the logarithmic score of the bacterial count, curd firmness after 30 min and 45 min, and the curd-firming instant rate constant. We found significant effects of MAP infection on milk composition, pH, and rennet coagulation time. The mean of the posterior distributions of heritability estimates on the liability scale was 0.15 for MAP and 0.07 for MVV. Our results demonstrate that only a few traits are negatively affected by MVV and MAP positivity, and that there is exploitable genetic variation in MVV and MAP susceptibility in dairy sheep.

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

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 (k20), and curd firmness at 30 min (a30)], and modeled curd-firming over time (CFt) parameters [estimated rennet coagulation time (RCTeq), curd-firming instant rate constant (kCF), and potential curd firmness (CFP)] 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 CFt 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 CFt 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 (%CYWATER). We found an optimum value of milk k20, associated with the highest recovery of components and cheese yield in solids (%CYSOLIDS). In addition, a lower level of curd water retention and an increased fresh curd yield (%CYCURD) 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 RCTeq on %REC traits and %CYSOLIDS, which had not previously been revealed for traditional RCT. Moreover, higher kCF values were associated with increased %CYCURD and %CYSOLIDS. Given that CFt 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.

Association between the GHR, GHRHR and IGF1 gene polymorphisms and milk coagulation properties in Sarda sheep.

We investigated whether variation of the sheep Growth Hormone Receptor (GHR), Growth Hormone Releasing Hormone Receptor (GHRHR) and Insulin-Like Growth Factor 1 (IGF1) genes were associated with milk coagulation properties (MCP) in sheep. The GHR, GHRHR and IGF1 genes are part of the GH system, which is known to modulate metabolism, growth and reproduction as well as mammogenesis and galactopoiesis in dairy species. A total of 380 dairy Sarda sheep were genotyped for 36 SNPs mapping to these three genes. Traditional MCP were measured as rennet coagulation time (RCT), curd-firming time (k20) and curd firmness at 30 m (a30). Modeling of curd firming over time (CFt) was based on a 60 m lactodynamographic test, generating a total of 240 records of curd firmness (mm) for each milk sample. The model parameters obtained included: the rennet coagulation time as a result of modeling all data available (RCTeq, min); the asymptotic potential value of curd firmness (CFP, mm) at an infinite time; the CF instant rate constant (kCF, %/min); the syneresis instant rate constant (kSR, %/min); the maximum value of CF (CFmax, mm) and the time at achievement of CFmax (tmax, min). Statistical analysis revealed that variation of the GHR gene was significantly associated with RCT, kSR and CFP (P < 0.05). No other significant associations were detected. These findings may be useful for the dairy industry, as well as for selection programs.

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.

Association between the GHR, GHRHR, and IGF1 gene polymorphisms and milk yield and quality traits in Sarda sheep.

The growth hormone receptor (GHR), the growth hormone releasing hormone receptor (GHRHR), and the insulin-like growth factor 1 (IGF1) genes are known to modulate growth, reproduction, and lactation traits in livestock. The aim of the current work was to investigate if the variation of the sheep GHR, GHRHR, and IGF1 genes is associated with milk yield and quality traits. Three hundred eighty dairy Sarda sheep were genotyped for 36 single nucleotide polymorphisms (SNP) mapping to these 3 loci, and records for milk yield and daily fat and protein yield, as well as for fat, protein, casein, lactose, and milk urea contents, pH, somatic cell score, logarithmic bacterial count, and milk energy were obtained. The linkage disequilibrium analysis was performed only for GHR, as both GHRHR and IGF1 had only 1 polymorphic SNP. Haplotype analysis revealed the existence of 7 haplotype blocks in GHR. Two haplotype blocks, including part of the intron 1 and the upstream region, were clearly separated from the remaining 5 blocks by SNP rs412986330, which may be a recombination hotspot. The latter 5 blocks were contiguous, spanning from intron 2 to exon 10. Statistical analysis revealed that the GHR polymorphism is significantly associated with milk traits for daily fat and protein yield and fat, milk urea, and lactose content. Moreover, variation in IGF1 was associated with milk protein and casein content. Data generated in this research provide new insights into the allelic effects of the ovine GHRHR, GHR, and IGF1 genes on milk production and quality traits, information that may be useful in gene-assisted selection programs.

Cheese yield, cheesemaking efficiency, and daily production of 6 breeds of goats.

Little is known about the complex process of cheesemaking at the individual level of dairy goats because of the difficulties of producing a high number of model cheeses. The objectives of this work were (1) to study the cheesemaking ability of goat milk; (2) to investigate the variability of cheesemaking-related traits among different farms; (3) to assess the effects of stage of lactation and parity; and (4) to compare 6 breeds of goat (Saanen and Camosciata delle Alpi for the Alpine type; Murciano-Granadina, Maltese, Sarda and Sarda Primitiva for the Mediterranean type) for their cheesemaking ability. For each goat (n = 560) we studied (1) 8 milk quality traits (fat, protein, total solids, casein, lactose, pH, somatic cell score, and bacterial count); (2) 4 milk nutrient recovery traits (fat, protein, total solids, and energy) in curd; (3) 3 actual cheese yield traits (fresh cheese, cheese solids, and cheese water); (4) 2 theoretical cheese yield values (fresh cheese and cheese solids) and the related cheesemaking efficiencies; and (5) daily milk yield and 3 daily cheese yield traits (fresh cheese, cheese solids, and water retained in the curd). With respect to individual animal factors, farm was not particularly important for recovery traits or actual and theoretical cheese yield and estimates of efficiency, whereas it highly influenced daily productions. Parity of goats influenced daily cheese production, whereas DIM slightly affected recovery as well as percent and daily cheese yield traits. Breed was the most important source of variation for almost all cheesemaking traits. Compared with those of Alpine type, the 4 Mediterranean breeds had, on average, lower daily milk and cheese productions, greater actual and theoretical cheese yield, and higher recovery of nutrients in the curd. Among Alpine type, Camosciata delle Alpi was characterized by greater nutrients recovery than Saanen. Within the 4 Mediterranean types, the 3 Italians produced much less milk per day, with much more fat and protein and greater recovery traits than the Murciano-Granadina, resulting in greater actual cheese yield. Within the Italian breeds, milk from Sarda and Sarda Primitiva was characterized by lower daily yields, higher protein and fat content, and greater recoveries of nutrients than Maltese goats. These results confirmed the potential of goat milk for cheese production and could be useful to give new possibilities and direction in breeding programs.

Effect of composition on coagulation, curd firming, and syneresis of goat milk.

The present study investigated the effect of different levels of fat, protein, and casein on (1) traditional milk coagulation properties, and (2) curd firming over time parameters of 1,272 goat milk samples. Relationships between fat, protein, and casein and some indicators of udder health status (lactose, pH, somatic cells, bacterial count, and NaCl) were also investigated. Traditional milk coagulation properties and modeled curd-firming parameters were analyzed using a mixed model that considered the effect of days in milk, parity, farm, breed, the pendulum of the instrument, and different levels of fat, protein, and casein. Fat, protein, and casein were also tested with the same model but one at a time. Information provided by this model demonstrated the effect of one component alone, without contemporarily considering that of the others. The results allowed us to clarify the effect of the major milk nutrients on coagulation, curd firming, and syneresis ability of goat milk. In particular, milk rich in fat was associated with better coagulation properties, whereas milk rich in protein was associated with delayed coagulation. The high correlation of fat with protein and casein contents suggests that the effect of fat on the cheese-making process is also attributable to the effects of protein and casein. When only protein or only casein was included in the statistical model, the pattern of coagulation, curd firming, and syneresis was almost indistinguishable. The contemporary inclusion of protein and casein in the statistical model did not generate computing problems and allowed us to better characterize the role of protein and casein. Consequently, given their strong association, we also tested the effect of casein-to-protein ratio (i.e., casein number). Higher values of casein number led to a general improvement in the coagulation ability of milk, suggesting that casein-to-protein ratio, not just protein or casein, should be considered when milk is destined for cheese making. These results are especially useful for dairy farmers who want to increase their profits by improving the technological quality of the milk produced.

Modeling of coagulation, curd firming, and syneresis of goat milk from 6 breeds.

Traditional milk coagulation properties are used to predict the suitability of milk for cheese-making. In bovine and ovine species, the introduction of the concept of curd firming over time, continuously recorded by a lactodynamograph during prolonged tests, provides additional information about milk coagulation, curd-firming, and syneresis processes. The aims of present study were (1) to test the adaptability of a 4-parameter curd-firming model in the assessment of goat milk (also comparing published data of other species); (2) to describe variability of coagulation, curd firming, and syneresis processes among individual goat milk samples; (3) to quantify the effects of farm and animal factors (breed, parity, and stage of lactation); and (4) to compare 6 goat breeds for their model parameters. Milk samples from 1,272 goats reared in 35 farms were collected. Goats were of 6 breeds: Saanen and Camosciata delle Alpi for the Alpine type; and Murciano-Granadina, Maltese, Sarda, and Sarda Primitiva for the Mediterranean type. During a lactodynamographic analysis (60 min), 240 measures of curd firmness (mm) were recorded for each milk sample. The modeling of curd firming allowed us to achieve the rennet coagulation time estimated on the basis of all the data points (min); the curd firming and the curd syneresis instant rate constants; the asymptotical potential value of curd firming; the actual maximum curd firmness; and the time at which the curd firming maximum level is attained. Modeling parameter data were analyzed using a linear mixed model. Comparison with other dairy species showed several differences: goat milk coagulated later than sheep but earlier than bovine, and curd firming and curd syneresis instant rate constants were greater in small ruminants. Modeling parameters of goat milk were mostly affected by the farm effect (37% of the total variance, on average) compared with the results found for bovine and ovine samples, and this was probably attributable to the marked differences among goat farming systems. Small differences were demonstrated between Alpine and Mediterranean breeds, but the time of maximum curd firmness was lower in Murciano-Granadina compared with Maltese, Sarda, and Sarda Primitiva. Sarda and Sarda Primitiva were very similar and exhibited the most favorable coagulation properties of milk. For almost all the model parameters, the direct effect of breed was increased after correction for milk yield and composition. In conclusion, this approach allowed us to fully depict the effects of the different factors on coagulation of goat milk, and clarified the different renneting pattern among goat breeds, and with other species. Results could be used for the valorization of goat dairy products, also when these are linked to particular local breeds, and to stimulate further studies about relationships between coagulation and cheese-making traits.

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 (k20, min), curd firmness at 30, 45, and 60 min after rennet addition (a30, a45, and a60, 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 a30). The proportion of variance regarding dMY, dFPY, and milk composition due to the farm was moderate, whereas for MCP it was low, except for a60, 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 k20 decreased rapidly and a30 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 k20 and a greater a30. 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 k20, a30, and a45. 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.

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 (CFt, mm) using a nonlinear procedure. The most accurate and informative results were observed using a modified 4-parameter model, structured as follows: [Formula: see text] where t is time, RCTeq (min) is the gelation time, CFP (mm) is the potential asymptotical CF at an infinite time, kCF (%/min) is the curd firming rate constant, and kSR (%/min) is the curd syneresis rate constant. To avoid nonconvergence and computational problems due to interrelations among the equation parameters, CFP was preliminarily defined as a function of maximum observed curd firmness (CFmax, 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.

Phenotypic and genetic relationships between indicators of the mammary gland health status and milk composition, coagulation, and curd firming in dairy sheep.

The present study investigated the effect of somatic cell count, lactose, and pH on sheep milk composition, coagulation properties (MCP), and curd firming (CF) parameters. Individual milk samples were collected from 1,114 Sarda ewes reared in 23 farms. Milk composition, somatic cell count, single point MCP (rennet coagulation time, RCT; curd firming time, k20; and curd firmness, a30, a45, and a60), and CF model parameters were achieved. Phenotypic traits were statistically analyzed using a mixed model to estimate the effects of the different levels of milk somatic cell score (SCS), lactose, and pH, respectively. Additive genetic, herd, and residual correlations among these 3 traits, and with milk composition, MCP and CF parameters, were inferred using a Bayesian approach. From a phenotypic point of view, higher SCS levels caused a delayed gelification of milk. Lactose concentration and pH were significant for many milk quality traits, with a very intense effect on both coagulation times and curd firming. These traits (RCT, RCT estimated using the curd firming over time equation, and k20) showed an unfavorable increase of about 20% from the highest to the lowest level of lactose. Milk samples with pH values lower than 6.56 versus higher than 6.78 were characterized by an increase of RCT (from 6.00 to 14.3 min) and k20 (from 1.65 to 2.65 min) and a decrease of all the 3 curd firmness traits. From a genetic point of view, the marginal posterior distribution of heritability estimates evidenced a large and exploitable variability for all 3 phenotypes. The mean intra-farm heritability estimates were 0.173 for SCS, 0.418 for lactose content, and 0.206 for pH. Lactose (favorably), and SCS and pH (unfavorably), at phenotypic and genetic levels, were correlated mainly with RCT and RCT estimated using the curd firming over time equation and scarcely with the other curd firming traits. The SCS, lactose, and pH were significantly correlated with each other's. In conclusion, results reported in the present study suggest that SCS, pH, and lactose affect, contemporarily and independently, milk quality and MCP. These phenotypes, easily available during milk recording schemes measured by infrared spectra prediction, could be used as potential indicators traits for improving cheese-making ability of ovine milk.

Relationships between milk characteristics and somatic cell score in milk from primiparous browsing goats.

To determine milk yield and composition, total microbic count (TMC) and somatic cell count (SCC) of browsing goats throughout the first lactation, 100 goats of Sarda breed, equally distributed in four flocks (F1, F2, F3 and F4), were selected. They were exclusively fed pasture and hand-milked once daily. Individual milk samples and daily milk yield were taken from each goat at monthly intervals, from March to July. Milk samples were analyzed for: total protein, fat, lactose, urea, freezing point (FP), pH, TMC and SCC. The data was subjected to analysis of variance and to correlation matrix. On the whole, in all the flocks, milk yield showed the highest production in April and May. Fat content increased (P < 0.01) throughout the lactation. Protein content showed the lowest value (P < 0.01) in June (4.15%). Urea and pH values were fluctuating. FP was lower (P < 0.01) at the start of lactation (-0.562 Hortvet degrees). TMC log10 values were low, considering the hand milking and inadequacy of facilities on the farms. SCC increased (P < 0.01) throughout the lactation and, on the whole, SCC and TMC were not correlated.