Multivariate analysis of the milk coagulation process in ovine breeds from Spain

In Spain, ewe milk is mainly used for cheesemaking, and farming systems have traditionally been based on the use of autochthonous breeds. However, in recent years, the progressive introduction of highly productive foreign breeds in Spanish farms has led to an increasing interest in the characterization of dairy sheep breeds to evaluate whether genetic selection schemes should focus on productivity or milk technological aptitude. The purpose of this work was to explore milk composition and coagulation to classify 4 of the main dairy sheep breeds used in Spain. This study included 832 individual ewe milk samples from the breeds Manchega, Assaf, Merino de Grazalema, and Merino de Los Pedroches. Samples were analyzed for native pH, composition (fat, protein, lactose, and total solids), coagulation properties, and individual laboratory curd yield. An indicator of coagulation efficiency was also determined. Canonical discriminant analysis was performed to establish differences and similarities among breeds based on the measured variables. In addition, cluster analysis was performed to study and quantify the concrete relationships among the discriminated groups. Discriminant analysis proved to be a powerful tool to accurately draw distinctions between breeds. In all cases, discrimination among breeds was evident and the 4 breeds could be easily differentiated. Cluster analysis showed greater similarity between Merino de Grazalema and Assaf compared with the other breeds, and F-statistics indicated a higher discriminating ability for the variables related to milk composition. However, Merino de Grazalema and Manchega were difficult to separate according to milk composition, but the coagulation process differenced them clearly. Coagulation also evidenced similarities between Manchega and Merino de Los Pedroches, although the latter was revealed to be the most different breed of all 4, which could lay the ground for its differentiation as an independent breed in the Official Catalogue of Spanish Livestock Breeds.     

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.     

Effects of seasonal changes in feeding management under part-time grazing on the evolution of the composition and coagulation properties of raw milk from ewes

Ewe raw milk composition, rennet coagulation parameters, and curd texture were monitored throughout the milk production season in 11 commercial flocks reared under a part-time grazing system. Milking season lasted from February to July. During that period, the diet of the animals shifted from indoor feeding, consisting of concentrate and forage, to an outdoor grazing diet. Lean dry matter, fat, protein, calcium, and magnesium contents increased throughout the milking season, as did rennet coagulation time, curd firmness, and curd resistance to compression. However, lean dry matter, protein content, and curd resistance to compression stabilized when sheep started to graze. Principal component analysis correlated curd resistance to compression and proteins, whereas curd firmness was highly correlated with fat content and minerals. Discriminant analysis distributed milk samples according to the feeding management. Curd firmness, fat, and magnesium turned out to be discriminant variables. Those variables reflected the evolution of the composition and coagulation parameters when fresh pasture prevailed over other feeds in the diet of the flocks. The present study shows that seasonal changes associated with feeding management influence milk technological quality and that milk of good processing quality can be obtained under part-time grazing.     

Effect of microparticulated whey proteins on milk coagulation properties

The enhancement of milk coagulation properties (MCP) and the reuse of whey produced by the dairy industry are of great interest to improve the efficiency of the cheese-making process. Native whey proteins (WP) can be aggregated and denatured to obtain colloidal microparticulated WP (MWP). The objective of this study was to assess the effect of MWP on MCP; namely, rennet coagulation time (RCT), curd-firming time, and curd firmness 30 min after rennet addition. Six concentrations of MWP (vol/vol; 1.5, 3.0, 4.5, 6.0, 7.5, and 9.0%) were added to 3 bulk milk samples (collected and analyzed during 3 d), and a sample without MWP was used as control. Within each day of analysis, 6 replicates of MCP for each treatment were obtained, changing the position of the treatment in the rack. For control samples, 2 replicates per day were performed. In addition to MCP, WP fractions were measured on each treatment during the 3 d of analysis. Milk coagulation properties were measured on 144 samples by using a Formagraph (Foss Electric, Hillerød, Denmark). Increasing the amount of MWP added to milk led to a longer RCT. In particular, significant differences were found between RCT of the control samples (13.5 min) and RCT of samples with 3.0% (14.6 min) or more MWP. A similar trend was observed for curd-firming time, which was shortest in the control samples and longest in samples with 9.0% MWP (21.4 min). No significant differences were detected for curd firmness at 30 min across concentrations of MWP. Adjustments in cheese processing should be made when recycling MWP, in particular during the coagulation process, by prolonging the time of rennet activity before cutting the curd.     

Genetic analysis of detailed milk protein composition and coagulation properties in Simmental cattle

The objective of this study was to estimate genetic parameters for milk protein fraction contents, milk protein composition, and milk coagulation properties (MCP). Contents of α_S1-, α_S2-, β-, γ-, and κ-casein (CN), β-lactoglobulin (β-LG), and α-lactalbumin (α-LA) were measured by reversed-phase HPLC in individual milk samples of 2,167 Simmental cows. Milk protein composition was measured as percentage of each CN fraction in CN (α_S1-CN%, α_S2-CN%, β-CN%, γ-CN%, and κ-CN%) and as percentage of β-LG in whey protein (β-LG%). Rennet clotting time (RCT) and curd firmness (a₃₀) were measured by a computerized renneting meter. Heritabilities for contents of milk proteins ranged from 0.11 (α-LA) to 0.52 (κ-CN). Heritabilities for α_S1-CN%, κ-CN%, and β-CN% were similar and ranged from 0.63 to 0.69, whereas heritability of α_S2-CN%, γ-CN%, and β-LG% were 0.28, 0.18, and 0.34, respectively. Effects of CSN2-CSN3 haplotype and BLG genotype accounted for more than 80% of the genetic variance of α_S1-CN%, β-CN%, and κ-CN% and 50% of the genetic variance of β-LG%. The genetic correlations among the contents of CN fractions and between CN and whey protein fractions contents were generally low. When the data were adjusted for milk protein gene effects, the magnitude of the genetic correlations among the contents of milk protein fractions markedly increased, indicating that they undergo a common regulation. The proportion of β-CN in CN correlated negatively with κ-CN% (r = −0.44). The genetic relationships between CN and whey protein composition were trivial. Low milk pH correlated with favorable MCP. Genetically, contents and proportions of α_S1- and α_S2-CN in CN were positively correlated with RCT. The relative proportion of β-CN in CN exhibited a genetic correlation with RCT of −0.26. Both the content and the relative proportion of κ-CN in CN did not correlate with RCT. Weak curds were genetically associated with increased proportions in CN of α_S1- and α_S2-CN, decreased contents of β-CN and κ-CN, and decreased proportion of κ-CN in CN. Negligible effects on the estimated correlations between a₃₀ and κ-CN contents or proportion in CN were observed when the model accounted for milk protein gene effects. Increasing β-CN and κ-CN contents and relative proportions in CN and decreasing the content and proportions of α_S1-CN and α_S2-CN and milk pH through selective breeding exert favorable effects on MCP.     

Use of multivariate factor analysis to define new indicator variables for milk composition and coagulation properties in Brown Swiss cows

The aim of this study was to elucidate the structure of relationships between milk yield, composition, and coagulation properties of Brown Swiss cattle. Multivariate factor analysis was used to derive new synthetic variables that can be used for selection purposes. For this reason, genetic parameters of these new variables were estimated. Individual records on milk yield, fat and protein percentages, casein content, lactose percentage, somatic cell count, titratable acidity, and pH were taken on 1,200 Italian Brown Swiss cows located in 38 herds. Factor analysis was able to extract 4 latent variables with an associated communality equal to 70% of the total original variance. The 4 latent factors were interpreted as indicators of milk composition, coagulation, acidity, and mammary gland health, respectively. Factor scores calculated for each animal exhibited coherent patterns along the lactation and across different parities. Estimation of genetic parameters of factor scores carried out with a multiple-trait Bayesian hierarchical model showed moderate to low heritabilities (raging from 0.10 to 0.23) and genetic correlations (from −0.15 to 0.46). Results of the present study support the hypothesis of a simpler structure that controls, at least in part, the covariance of milk composition and coagulation properties. Moreover, extracted variables may be useful for both breeding and management purposes, being able to represent, with a single value for each animal, complex traits such as milk coagulation properties or health status of the mammary gland.     

Association between the bovine milk metabolome and rennet-induced coagulation properties of milk

The milk metabolomes of 407 individual Swedish Red dairy cows were analyzed by nuclear magnetic resonance spectroscopy as part of the Danish-Swedish Milk Genomics Initiative. By relating these metabolite profiles to total milk protein concentration and rheological measurements of rennet-induced milk coagulation together using multivariate data analysis techniques, we were able to identify several different associations of the milk metabolome to technological properties of milk. Several novel correlations of milk metabolites to protein content and rennet-induced coagulation properties were demonstrated. Metabolites associated with the prediction of total protein content included choline, N-acetyl hexosamines, creatinine, glycerophosphocholine, glutamate, glucose 1-phosphate, galactose 1-phosphate, and orotate. In addition, levels of lactate, acetate, glutamate, creatinine, choline, carnitine, galactose 1-phosphate, and glycerophosphocholine were significantly different when comparing noncoagulating and well-coagulating milks. These findings suggest that the mentioned metabolites are associated with milk protein content and rennet-induced coagulation properties and may act as quality markers for cheese milk.     

Genetic and phenotypic correlations between milk coagulation properties, milk production traits, somatic cell count, casein content, and pH of milk

Genetic and phenotypic correlations between milk coagulation properties (MCP: coagulation time and curd firmness), milk yield, fat content, protein content, ln(somatic cell count) (SCS), casein content, and pH of milk and heritability of these traits were estimated from data consisting of milk samples of 4664 Finnish Ayrshire cows sired by 91 bulls. In addition, differences in average estimated breeding values (EBV) for the above traits between the cows with noncoagulating (NC) milk and those with milk that coagulated (CO samples) were examined. The estimations were carried out to study the possibilities of indirect genetic improvement of MCP by use of the above characteristics. The genetic and phenotypic correlations between MCP and the milk production traits were low or negligible. The genetic associations between desirable MCP and low SCS were rather strong (-0.45 to 0.29). Desirable MCP correlated both genetically and phenotypically with low pH of milk (-0.51 to 0.50). The rather high heritability estimates for curd firmness in different forms (0.22 to 0.39), and the wide variation in the proportion of daughters producing NC milk between the sires (0 to 47%) suggested that noncoagulation of milk is partly caused by additive genetic factors. Based on the genetic correlations between curd firmness and SCS and the high EBV for SCS obtained for the cows with NC-milk, it is possible that the loci causing noncoagulation of milk and increasing somatic cell count of milk are closely linked or partly the same. One means to genetically improve MCP and to reduce the occurrence of NC milk could thus be selection for low somatic cell count of milk.     

Composition and properties of bovine milk: A study from dairy farms in northern Sweden; Part I. Effect of dairy farming system

This study was part of a larger project that aimed to understand the causes for increasing variation in cheese ripening in a cheese-producing region in northern Sweden. The influence of different on-farm factors on raw milk composition and properties was investigated and is described in this paper, whereas the monthly variation in the milk quality traits during 1 yr is described in our companion paper. The dairy farming systems on a total of 42 dairy farms were characterized through a questionnaire and farm visits. Milk from farm tanks was sampled monthly over 1 yr and analyzed for quality attributes important for cheese making. On applying principal component analyses to evaluate the variation in on-farm factors, different types of farms were distinguished. Farms with loose housing and automatic milking system (AMS) or milking parlor had a higher number of lactating cows, and predominantly Swedish Holstein (SH) breed. Farms associated with tiestalls had a lower number of lactating cows and breeds other than SH. Applying principal component analyses to study the variation in composition and properties of tank milk samples from farms revealed a tendency for the formation of 2 clusters: milk from farms with AMS or a milking parlor, and milk from farms with tiestall milking. The interaction between the milking system, housing system, and breed probably contributed to this grouping. Other factors that were used in the characterization of the farming systems only showed a minor influence on raw milk quality. Despite the interaction, milk from tiestall farms with various cow breeds had higher concentrations (g/100 g of milk) of fat (4.74) and protein (3.63), and lower lactose concentrations (4.67) than milk from farms with predominantly SH cows and AMS (4.32, 3.47, and 4.74 g/100 g of milk, respectively) or a milking parlor (4.47, 3.54, and 4.79 g/100 g of milk, respectively). Higher somatic cell count (195 × 10³/mL) and lower free fatty acid concentration (0.75 mmol/100 g of fat) were observed in milk from farms with AMS than in milk from tiestall systems (150 × 10³/mL and 0.83 mmol/100 g of fat, respectively). Type of farm influenced milk gel strength, with milk from farms with predominantly SH cows showing the lowest gel strength (65.0 Pa), but not a longer rennet coagulation time. Effects of dairy farming system (e.g., dominant breed, milking system, housing, and herd size) on milk quality attributes indicate a need for further studies to evaluate the in-depth effects of farm-related factors on milk quality attributes.     

Genetic parameters of milk coagulation properties and their relationships with milk yield and quality traits in Italian Holstein cows

Milk coagulation properties (MCP) are an important aspect in assessing cheese-making ability. Several studies showed that favorable conditions of milk reactivity with rennet, curd formation rate, and curd strength, as well as curd syneresis, have a positive effect on the entire cheese-making process and subsequently on the ripening of cheese. Moreover, MCP were found to be heritable, but little scientific literature is available about their genetic aspects. The aims of this study were to estimate heritability of MCP and genetic correlations among MCP and milk production and quality traits. A total of 1,071 Italian Holstein cows (progeny of 54 sires) reared in 34 herds in Northern Italy were sampled from January to July 2004. Individual milk samples were collected during the morning milking and analyzed for coagulation time (RCT), curd firmness (a₃₀), pH, titratable acidity, fat, protein, and casein contents, and somatic cell count. About 10% of individual milk samples did not coagulate in 31 min, so they were removed from the analyses. Estimates of heritability for RCT and a₃₀ were 0.25 ± 0.04 and 0.15 ± 0.03, respectively. Estimates of genetic correlations between MCP traits and milk production traits were negligible except for a₃₀ with protein and casein contents (0.44 ± 0.10 and 0.53 ± 0.09, respectively). Estimates of genetic correlations between MCP traits and somatic cell score were strong and favorable, as well as those between MCP and pH and titratable acidity. Selecting for high casein content, milk acidity, and low somatic cell count might be an indirect way to improve MCP without reducing milk yield and quality traits.     

Mid-infrared spectroscopy predictions as indicator traits in breeding programs for enhanced coagulation properties of milk

The aims of this study were to investigate variation of milk coagulation property (MCP) measures and their predictions obtained by mid-infrared spectroscopy (MIR), to investigate the genetic relationship between measures of MCP and MIR predictions, and to estimate the expected response from a breeding program focusing on the enhancement of MCP using MIR predictions as indicator traits. Individual milk samples were collected from 1,200 Brown Swiss cows (progeny of 50 artificial insemination sires) reared in 30 herds located in northern Italy. Rennet coagulation time (RCT, min) and curd firmness (a₃₀, mm) were measured using a computerized renneting meter. The MIR data were recorded over the spectral range of 4,000 to 900 cm⁻¹. Prediction models for RCT and a₃₀ based on MIR spectra were developed using partial least squares regression. A cross-validation procedure was carried out. The procedure involved the partition of available data into 2 subsets: a calibration subset and a test subset. The calibration subset was used to develop a calibration equation able to predict individual MCP phenotypes using MIR spectra. The test subset was used to validate the calibration equation and to estimate heritabilities and genetic correlations for measured MCP and their predictions obtained from MIR spectra and the calibration equation. Point estimates of heritability ranged from 0.30 to 0.34 and from 0.22 to 0.24 for RCT and a₃₀, respectively. Heritability estimates for MCP predictions were larger than those obtained for measured MCP. Estimated genetic correlations between measures and predictions of RCT were very high and ranged from 0.91 to 0.96. Estimates of the genetic correlation between measures and predictions of a₃₀ were large and ranged from 0.71 to 0.87. Predictions of MCP provided by MIR techniques can be proposed as indicator traits for the genetic enhancement of MCP. The expected response of RCT and a₃₀ ensured by the selection using MIR predictions as indicator traits was equal to or slightly less than the response achievable through a single measurement of these traits. Breeding strategies for the enhancement of MCP based on MIR predictions as indicator traits could be easily and immediately implemented for dairy cattle populations where routine acquisition of spectra from individual milk samples is already performed.     

Genetic parameters for rennet- and acid-induced coagulation properties in milk from Swedish Red dairy cows

Milk coagulation is an important processing trait, being the basis for production of both cheese and fermented products. There is interest in including technological properties of these products in the breeding goal for dairy cattle. The aim of the present study was therefore to estimate genetic parameters for milk coagulation properties, including both rennet- and acid-induced coagulation, in Swedish Red dairy cattle using genomic relationships. Morning milk samples and blood samples were collected from 395 Swedish Red cows that were selected to be as genetically unrelated as possible. Using a rheometer, milk samples were analyzed for rennet- and acid-induced coagulation properties, including gel strength (G′), coagulation time, and yield stress (YS). In addition to the technological traits, milk composition was analyzed. A binary trait was created to reflect that milk samples that had not coagulated 40 min after rennet addition were considered noncoagulating milk. The cows were genotyped by using the Illumina BovineHD BeadChip (Illumina Inc., San Diego, CA). Almost 600,000 markers remained after quality control and were used to construct a matrix of genomic relationships among the cows. Multivariate models including fixed effects of herd, lactation stage, and parity were fitted using the ASReml software to obtain estimates of heritabilities and genetic and phenotypic correlations. Heritability estimates (h²) for G′ and YS in rennet and acid gels were found to be high (h² = 0.38–0.62) and the genetic correlations between rennet-induced and acid-induced coagulation properties were weak but favorable, with the exception of YS_rennet with G′_acid and YS_acid, both of which were strong. The high heritability (h² = 0.45) for milk coagulating ability expressed as a binary trait suggests that noncoagulation could be eliminated through breeding. Additionally, the results indicated that the current breeding objective could increase the frequency of noncoagulating milk and lead to deterioration of acid-induced coagulation through unfavorable genetic associations with protein content (0.38) and milk yield (−0.61 to −0.71), respectively. The outcome of this study suggests that by including more detailed compositional traits genetically associated with milk coagulation or by including milk coagulation properties directly within the breeding goal, it appears possible to breed cows that produce milk better suited for production of cheese and fermented products.     

Genetic parameters of coagulation properties, milk yield, quality, and acidity estimated using coagulating and noncoagulating milk information in Brown Swiss and Holstein-Friesian cows

The aim of this study was to estimate heritabilities of rennet coagulation time (RCT) and curd firmness (a₃₀) and their genetic correlations with test-day milk yield, composition (fat, protein, and casein content), somatic cell score, and acidity (pH and titratable acidity) using coagulating and noncoagulating (NC) milk information. Data were from 1,025 Holstein-Friesian (HF) and 1,234 Brown Swiss (BS) cows, which were progeny of 54 HF and 58 BS artificial insemination sires, respectively. Milk coagulation properties (MCP) of each cow were measured once using a computerized renneting meter and samples not exhibiting coagulation within 31 min after rennet addition were classified as NC milk. For NC samples, RCT was unobserved. Multivariate analyses, using Bayesian methodology, were performed to estimate the genetic relationships of RCT or a₃₀ with the other traits and statistical inference was based on the marginal posterior distributions of parameters of concern. For analyses involving RCT, a right-censored Gaussian linear model was used and records of NC milk samples, being censored records, were included as unknown parameters in the model implementing a data augmentation procedure. Rennet coagulation time was more heritable [heritability (h²) = 0.240 and h² = 0.210 for HF and BS, respectively] than a₃₀ (h² = 0.148 and h² = 0.168 for HF and BS, respectively). Milk coagulation properties were more heritable than a single test-day milk yield (h² = 0.103 and h² = 0.097 for HF and BS, respectively) and less heritable than milk composition traits whose heritability ranged from 0.275 to 0.275, with the only exception of fat content of BS milk (h² = 0.108). A negative genetic correlation, lower than −0.85, was estimated between RCT and a₃₀ for both breeds. Genetic relationships of MCP with yield and composition were low or moderate and favorable. The genetic correlation of somatic cell score with RCT in BS cows was large and positive and even more positive were those of RCT with pH and titratable acidity in both breeds, ranging from 0.80 to 0.94. Including NC milk information in the data affected the estimated correlations and decreased the uncertainty associated with the estimation process. On the basis of the estimated heritabilities and genetic correlations, enhancement of MCP through selective breeding with no detrimental effects on yield and composition seems feasible in both breeds. Milk acidity may play a role as an indicator trait for indirect enhancement of MCP.     

Short communication: Fourier-transform mid-infrared spectroscopy to predict coagulation and acidity traits of sheep bulk milk

Sheep milk is mainly transformed into cheese; thus, the dairy industry seeks more rapid and cost-effective methods of analysis to determine milk coagulation and acidity traits. This study aimed to assess the feasibility of Fourier-transform mid-infrared spectroscopy to determine milk coagulation and acidity traits of sheep bulk milk and to classify milk samples according to their renneting capacity. A total of 465 bulk milk samples collected in 140 single-breed flocks of Comisana (84 samples, 24 flocks) and Sarda (381 samples, 116 flocks) breeds located in Central Italy were analyzed for coagulation properties (rennet coagulation time, curd firming time, and curd firmness) and acidity traits (pH and titratable acidity) using standard laboratory procedures. Fourier-transform mid-infrared spectroscopy prediction models for these traits were built using partial least squares regression analysis and were externally validated by randomly dividing the full data set into a calibration set (75%) and a validation set (25%). The discriminant capacity of the rennet coagulation time prediction model was determined using partial least squares discriminant analysis. Prediction models were more accurate for acidity traits than for milk coagulation properties, and the ratio of prediction to deviation ranged from 1.01 (curd firmness) to 2.14 (pH). Moreover, the discriminant analysis led to an overall accuracy of 74 and 66% for the calibration and validation sets, respectively, with greater sensitivity for samples that coagulated between 10 and 20 min and greater specificity to detect early-coagulating (<10 min) and late-coagulating (20–30 min) samples. Results suggest that Fourier-transform mid-infrared spectroscopy has the potential to help the dairy sheep industry identify milk with better coagulation ability for cheese production and thus improve milk transformation efficiency. However, further research is needed before this information can be exploited at the industry level.     

Short communication: Prediction of milk coagulation and acidity traits in Mediterranean buffalo milk using Fourier-transform mid-infrared spectroscopy

Milk coagulation and acidity traits are important factors to inform the cheesemaking process. Those traits have been deeply studied in bovine milk, whereas scarce information is available for buffalo milk. However, the dairy industry is interested in a method to determine milk coagulation and acidity features quickly and in a cost-effective manner, which could be provided by Fourier-transform mid-infrared (FT-MIR) spectroscopy. The aim of this study was to evaluate the potential of FT-MIR to predict coagulation and acidity traits of Mediterranean buffalo milk. A total of 654 records from 36 herds located in central Italy with information on milk yield, somatic cell score, milk chemical composition, milk acidity [pH, titratable acidity (TA)], and milk coagulation properties (rennet coagulation time, curd firming time, and curd firmness) were available for statistical analysis. Reference measures of milk acidity and coagulation properties were matched with milk spectral information, and FT-MIR prediction models were built using partial least squares regression. The data set was divided into a calibration set (75%) and a validation set (25%). The capacity of FT-MIR spectroscopy to correctly classify milk samples based on their renneting ability was evaluated by a canonical discriminant analysis. Average values for milk coagulation traits were 13.32 min, 3.24 min, and 39.27 mm for rennet coagulation time, curd firming time, and curd firmness, respectively. Milk acidity traits averaged 6.66 (pH) and 7.22 Soxhlet-Henkel degrees/100 mL (TA). All milk coagulation and acidity traits, except for pH, had high variability (17 to 46%). Prediction models of coagulation traits were moderately to scarcely accurate, whereas the coefficients of determination of external validation were 0.76 and 0.66 for pH and TA, respectively. Canonical discriminant analysis indicated that information on milk coagulating ability is present in the MIR spectra, and the model correctly classified as noncoagulating the 91.57 and 67.86% of milk samples in the calibration and validation sets, respectively. In conclusion, our results can be relevant to the dairy industry to classify buffalo milk samples before processing.     

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.     

Short communication: Influence of composite casein genotypes on additive genetic variation of milk production traits and coagulation properties in Holstein-Friesian cows

The aim of the study was to quantify the effects of composite β- and κ-casein (CN) genotypes on genetic variation of milk coagulation properties (MCP); milk yield; fat, protein, and CN contents; somatic cell score; pH; and titratable acidity (TA) in 1,042 Italian Holstein-Friesian cows. Milk coagulation properties were defined as rennet coagulation time (RCT) and curd firmness (a₃₀). Variance components were estimated using 2 animal models: model 1 included herd, days in milk, and parity as fixed effects and animal and residual as random effects, and model 2 was model 1 with the addition of composite β- and κ-CN genotype as a fixed effect. Genetic correlations between RCT and a₃₀ and between these traits and milk production traits were obtained with bivariate analyses, based on the same models. The inclusion of casein genotypes led to a decrease of 47, 68, 18, and 23% in the genetic variance for RCT, a₃₀, pH, and TA, respectively, and less than 6% for other traits. Heritability of RCT and a₃₀ decreased from 0.248 to 0.143 and from 0.123 to 0.043, respectively. A moderate reduction was found for pH and TA, whereas negligible changes were detected for other milk traits. Estimates of genetic correlations were comparable between the 2 models. Results show that composite β- and κ-CN genotypes are important for RCT and a₃₀ but cannot replace the recording of MCP themselves.     

Effects of composite beta- and kappa-casein genotypes on milk coagulation, quality, and yield traits in Italian Holstein cows

The aim of the study was to estimate the effect of the composite CSN2 and CSN3 genotypes on milk coagulation, quality, and yield traits in Italian Holstein cows. A total of 1,042 multiparous Holstein cows reared on 34 commercial dairy herds were sampled once, concurrently with monthly herd milk recording. The data included the following traits: milk coagulation time; curd firmness; pH and titratable acidity; fat, protein, and casein contents; somatic cell score; and daily milk, fat, and protein yields. A single-trait animal model was assumed with fixed effects of herd, days in milk, parity, composite casein genotype of CSN2 and CSN3 (CSN2-CSN3), and random additive genetic effect of an animal. The composite genotype of CSN2-CSN3 showed a strong effect on both milk coagulation traits and milk and protein yields, but not on fat and protein contents and other milk quality traits. For coagulation time, the best CSN2-CSN3 genotypes were those with at least one B allele in both the CSN2 and CSN3 loci. The CSN3 locus was associated more strongly with milk coagulation traits, whereas the CSN2 locus was associated more with milk and protein yields. However, because of the tight linkage between the 2 loci, the composite genotypes, or haplotypes, are more appropriate than the single-locus genotypes if they were considered for use in selection.     

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.     

Economic weights for major milk constituents of Manchega dairy ewes

The objective of this study was to estimate economic weights of major components of milk (fat and protein) for the Manchega dairy sheep breed. An economic study was carried out and the profit associated with fat and protein yields of milk was calculated as the difference between incomes and costs. Incomes were obtained from milk sales to cheese industry and a reference marked price was used. Costs were calculated considering the energy necessary to produce each of the components of milk, and the price per milk forage unit was calculated as the total expense in feedstuff divided by the total (theoretical) needs of the flock. Economic values were defined as partial derivatives of the profit function with respect to each trait. Economic weights for fat and protein yields were similar, being slightly greater for protein in all cases. For carrier, economic weights were close to zero and negative because an increase in carrier production without changes in fat and protein composition leads to an increase in energy demands, holding the sale price of milk constant. When genetic standard deviations were taken into account and standardized economic values were calculated, an increase in economic value of protein and a decrease in economic value of fat yields were observed. The consequences that different changes in production system conditions have on the estimated economic weights were also studied. In general, economic weights were relatively insensitive to changes in production levels and market prices except for changes on milk price. Given the economic importance shown for fat and protein, milk components should be taken into consideration when breeding objectives for dairy sheep are established.