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.     

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.     

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.     

Short communication: Factors affecting coagulation properties of Mediterranean buffalo milk

The aim of this study was to investigate sources of variation of milk coagulation properties (MCP) of buffalo cows. Individual milk samples were collected from 200 animals in 5 herds located in northern Italy from January to March 2010. Rennet coagulation time (RCT, min) and curd firmness after 30 min from rennet addition (a₃₀, mm) were measured using the Formagraph instrument (Foss Electric, Hillerød, Denmark). In addition to MCP, information on milk yield, fat, protein, and casein contents, pH, and somatic cell count (SCC) was available. Sources of variation of RCT and a₃₀ were investigated using a linear model that included fixed effects of herd, days in milk (DIM), parity, fat content, casein content (only for a₃₀), and pH. The coefficient of determination was 51% for RCT and 48% for a₃₀. The most important sources of variation of MCP were the herd and pH effects, followed by DIM and fat content for RCT, and casein content for a₃₀. The relevance of acidity in explaining the variation of both RCT and a₃₀, and of casein content in explaining that of a₃₀, confirmed previous studies on dairy cows. Although future research is needed to investigate the effect of these sources of variation on cheese yield, findings from the present study suggest that casein content and acidity may be used as indicator traits to improve technological properties of buffalo milk.     

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.     

Prediction of coagulation properties, titratable acidity, and pH of bovine milk using mid-infrared spectroscopy

This study investigated the potential application of mid-infrared spectroscopy (MIR 4,000-900 cm⁻¹) for the determination of milk coagulation properties (MCP), titratable acidity (TA), and pH in Brown Swiss milk samples (n = 1,064). Because MCP directly influence the efficiency of the cheese-making process, there is strong industrial interest in developing a rapid method for their assessment. Currently, the determination of MCP involves time-consuming laboratory-based measurements, and it is not feasible to carry out these measurements on the large numbers of milk samples associated with milk recording programs. Mid-infrared spectroscopy is an objective and nondestructive technique providing rapid real-time analysis of food compositional and quality parameters. Analysis of milk rennet coagulation time (RCT, min), curd firmness (a₃₀, mm), TA (SH°/50 mL; SH° = Soxhlet-Henkel degree), and pH was carried out, and MIR data were recorded over the spectral range of 4,000 to 900 cm⁻¹. Models were developed by partial least squares regression using untreated and pretreated spectra. The MCP, TA, and pH prediction models were improved by using the combined spectral ranges of 1,600 to 900 cm⁻¹, 3,040 to 1,700 cm⁻¹, and 4,000 to 3,470 cm⁻¹. The root mean square errors of cross-validation for the developed models were 2.36 min (RCT, range 24.9 min), 6.86 mm (a₃₀, range 58 mm), 0.25 SH°/50 mL (TA, range 3.58 SH°/50 mL), and 0.07 (pH, range 1.15). The most successfully predicted attributes were TA, RCT, and pH. The model for the prediction of TA provided approximate prediction (R² = 0.66), whereas the predictive models developed for RCT and pH could discriminate between high and low values (R² = 0.59 to 0.62). It was concluded that, although the models require further development to improve their accuracy before their application in industry, MIR spectroscopy has potential application for the assessment of RCT, TA, and pH during routine milk analysis in the dairy industry. The implementation of such models could be a means of improving MCP through phenotypic-based selection programs and to amend milk payment systems to incorporate MCP into their payment criteria.     

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.     

Selection for milk coagulation properties predicted by Fourier transform infrared spectroscopy in the Italian Holstein-Friesian breed

Milk coagulation is based on a series of physicochemical changes at the casein micelle level, resulting in formation of a gel. Milk coagulation properties (MCP) are relevant for cheese quality and yield, important factors for the dairy industry. They are also evaluated in herd bulk milk to reward or penalize producers of Protected Designation of Origin cheeses. The economic importance of improving MCP justifies the need to account for this trait in the selection process. A pilot study was carried out to determine the feasibility of including MCP in the selection schemes of the Italian Holstein. The MCP were predicted in 1,055 individual milk samples collected in 16 herds (66 ± 24 cows per herd) located in Brescia province (northeastern Italy) by means of Fourier transform infrared (FTIR) spectroscopy. The coefficient of determination of prediction models indicated moderate predictions for milk rennet coagulation time (RCT = 0.65) and curd firmness (a₃₀ = 0.68), and poor predictions for curd-firming time (k₂₀ = 0.49), whereas the range error ratio (8.9, 6.9, and 9.5 for RCT, k₂₀, and a₃₀, respectively) indicated good practical utility of the predictive models for all parameters. Milk proteins were genotyped and casein haplotypes (α_S1-, β-, α_S2-, and κ-casein) were reconstructed. Data from 51 half-sib families (19.9 ± 16.4 daughters per sire) were analyzed by an animal model to estimate (1) the genetic parameters of predicted RCT, k₂₀, and a₃₀; (2) the breeding values for these predicted clotting variables; and (3) the effect of milk protein genotypes and casein haplotypes on predicted MCP (pMCP). This is the first study to estimate both genetic parameters and breeding values of pMCP, together with the effects of milk protein genotypes and casein haplotypes, that also considered k₂₀, probably the most important parameter for the dairy industry (because it indicates the time for the beginning of curd-cutting). Heritability of predicted RCT (0.26) and k₂₀ (0.31) were close to the average heritability described in literature, whereas the heritability of a₃₀ was higher (0.52 vs. 0.27). The effects of milk proteins were statistically significant and similar to those obtained on measured MCP. In particular, haplotypes including uncommon variants showed positive (B-I-A-B) or negative (B-A¹-A-E) effects. Based on these findings, FTIR spectroscopy-pMCP is proposed as a potential selection criterion for the Italian Holstein.     

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.     

The mineral profile affects the coagulation pattern and cheese-making efficiency of bovine milk

Natural variations in milk minerals, their relationships, and their associations with the coagulation process and cheese-making traits present an opportunity for the differentiation of milk destined for high-quality natural products, such as traditional specialties or Protected Designation of Origin (PDO) cheeses. The aim of this study was to quantify the effects of the native contents of Ca, P, Na, K, and Mg on 18 traits describing traditional milk coagulation properties (MCP), curd firming over time (CF_t) equation parameters, cheese yield (CY) measures, and nutrient recoveries in the curd (REC) using models that either included or omitted the simultaneous effects of milk fat and casein contents. The results showed that, by including milk fat and casein and the minerals in the statistical model, we were able to determine the specific effects of each mineral on coagulation and cheese-making efficiency. In general, about two-thirds of the apparent effects of the minerals on MCP and the CF_t equation parameters are actually mediated by their association with milk composition, especially casein content, whereas only one-third of the effects are direct and independent of milk composition. In the case of cheese-making traits, the effects of the minerals were mediated only negligibly by their association with milk composition. High Ca content had a positive effect on the coagulation pattern and cheese-making traits, favoring water retention in the curd in particular. Phosphorus positively affected the cheese-making traits in that it was associated with an increase in CY in terms of curd solids, and in all the nutrient recovery traits. However, a very high P content in milk was associated with lower fat recovery in the curd. The variation in the Na content in milk only mildly affected coagulation, whereas with regard to cheese-making, protein recovery was negatively associated with high concentrations of this mineral. Potassium seemed not to be actively involved in coagulation and the cheese-making process. Magnesium content tended to slow coagulation and reduce CY measures. Further studies on the relationships of minerals with casein and protein fractions could deepen our knowledge of the role of all minerals in coagulation and the cheese-making process.     

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.     

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.     

Reproducibility and repeatability of measures of milk coagulation properties and predictive ability of mid-infrared reflectance spectroscopy

The objectives of the study were to estimate the reproducibility and repeatability of milk coagulation properties (MCP) measured by a computerized renneting meter (CRM) and to evaluate the predictive ability of mid-infrared spectroscopy (MIRS) as an innovative technology for the assessment of rennet coagulation time (RCT, min) and curd firmness (a₃₀, mm). Four samples without addition of preservative (NP) and 4 samples with Bronopol addition (PS) were collected from each of 83 Holstein-Friesian cows. Six hours after collection, 2 replicated measures of MCP were obtained with CRM using 1 NP and 1 PS sample from each cow. Mid-infrared spectra of the remaining NP and PS samples from each animal were recorded after 6 h, 4 d, and 8 d after sampling. Two groups of calibration equations were developed using MIRS spectra and CRM measures of MCP as reference data obtained from analysis of NP and PS, respectively. Reproducibility and repeatability of CRM measures were obtained from REML estimation of variance components on the basis of a linear model including the fixed effects of herd and days in milk class and the random effects of cows, sample treatment (addition or no addition of preservative), and the interaction between cow and sample treatment. Coefficient of reproducibility is an indicator of the agreement between 2 measurements of MCP for the same milk sample preserved with or without addition of Bronopol. Coefficient of repeatability is an indicator of the agreement between repeated measures of MCP. Pearson correlations between MCP measures for NP and PS were 0.97 and 0.83 for RCT and a₃₀, respectively. Reproducibility of CRM measures under different preserving conditions of milk was 93.5% for RCT and 64.6% for a₃₀. Repeatabilities of RCT and a₃₀ measures were 95.7 and 77.3%, respectively. Based on the estimated cross-validation standard errors and coefficients of determination and ratios of standard errors of cross-validation to standard deviation of reference data, the predictive ability of MIRS calibration equations was moderate for RCT and unsatisfactory for a_30. Predictive ability of equations based on spectra and MCP measures of PS was greater than that of equations based on data of NP. The study did not provide conclusive evidence on the effectiveness of MIRS as a predictive tool for MCP and it requires an enlargement of the variability of milk sampling circumstances. Because the relevance of MIRS predictions in relation to breeding programs for MCP based on indicator traits relies on the genetic variation of MIRS predictions and on phenotypic and genetic correlations between MIRS predictions and MCP measures, additional specific investigations on these topics are needed.     

Relationships between milk coagulation property traits analyzed with different methodologies

Milk coagulation properties (MCP) analysis is performed using a wide range of methodologies in different countries and laboratories, using different instruments, coagulant activity in the milk, and type of coagulant. This makes it difficult to compare results and data from different research. The aims of this study were to propose a method for the transformation of values of rennet coagulation time (RCT) and curd firmness (a₃₀) and to predict the noncoagulation (NC) probability of milk samples analyzed using different methodologies. Individual milk samples were collected during the morning milking in October 2010 from each of 165 Holstein-Friesian dairy cows in 2 freestall barns in Italy, and sent to 3 laboratories for MCP analysis. For each laboratory, MCP analysis was performed using a different methodology: A, with a computerized renneting meter instrument using 0.051 international milk clotting units (IMCU)/mL of coagulant activity; B, with a Lattodinamografo (Foss-Italia, Padova, Italy) using 0.051 IMCU/mL of coagulant activity; and C, with an Optigraph (Ysebaert, Frépillon, France) using 0.120 IMCU/mL of coagulant activity. The relationships between MCP traits were analyzed with correlation and regression analyses for each pair of methodologies. For each MCP trait, 2 regression models were applied: model 1 was a single regression model, where the dependent and independent variables were the same MCP trait determined by 2 different methodologies; in model 2, both a₃₀ and RCT were included as independent variables. The NC probabilities for laboratories with the highest number of NC samples were predicted based on the RCT and a₃₀ values measured in the laboratories with lower number of NC samples using logistic regression and receiver operating characteristic analysis. The percentages of NC samples were 4.2, 11.5, and 0.6% for A, B, and C, respectively. The transformation of MCP traits was more precise with model 1 for RCT (R²: 0.77-0.82) than for a₃₀ (R²: 0.28-0.63). The application of model 2 was needed when the C measurements were transformed into the other scales. The analyses of NC probabilities of milk samples showed that NC samples from one methodology were well distinguishable (with an accuracy of 0.972-0.996) based on the rennet coagulation time measured with the other methodology. A standard definition for MCP traits analysis is needed to enable reliable comparisons between MCP traits recorded in different laboratories and in different animal populations and breeds.     

Genetic and environmental parameters for milk production, udder health, and fertility traits in Mexican Holstein cows

Genetic and phenotypic parameters for Mexican Holstein cows were estimated for first- to third-parity cows with records from 1998 to 2003 (n = 2,971-15,927) for 305-d mature equivalent milk production (MEM), fat production (MEF), and protein production (MEP), somatic cell score (SCS), subsequent calving interval (CAI), and age at first calving (AFC). Genetic parameters were obtained by average information matrix-REML methodology using 6-trait (first-parity data) and 5-trait (second- and third-parity data) animal models. Heritability estimates for production traits were between 0.17 ± 0.02 and 0.23 ± 0.02 for first- and second-parity cows and between 0.12 ± 0.03 and 0.13 ± 0.03 for third-parity cows. Heritability estimates for SCS were similar for all parities (0.10 ± 0.02 to 0.11 ± 0.03). For CAI, estimates of heritability were 0.01 ± 0.05 for third-parity cows and 0.02 ± 0.02 for second-parity cows. The heritability for AFC was moderate (0.28 ± 0.03). No unfavorable estimates of correlations were found among MEM, MEF, MEP, CAI, and SCS. Estimates of environmental and phenotypic correlations were large and positive among production traits; favorable between SCS and CAI; slightly favorable between MEM, MEF, and MEP and SCS, between AFC and SCS, and between SCS and CAI; and small but unfavorable between production traits and CAI. Estimates of genetic variation and heritability indicate that selection would result in genetic improvement of production traits, AFC, and SCS. Estimates of both heritability and genetic variation for CAI were small, which indicates that genetic improvement would be difficult.     

Factors associated with herd bulk milk composition and technological traits in the Italian dairy industry

The aim of the present study was to investigate sources of variation of milk composition and technological characteristics routinely collected in field conditions in the Italian dairy industry. A total of 40,896 bulk milk records from 620 herds and 10 regions across Italy were analyzed. Composition traits were fat, protein, and casein percentages, urea content, and somatic cell score; and technological characteristics were rennet coagulation time, curd firming time, curd firmness 30 min after rennet addition to milk, and titratable acidity. Data of herd bulk milks were analyzed using a model that included fixed effects of region, herd nested within region, and season of milk analysis. An average good milk quality was reported in the dairy industry (especially concerning fat, protein, and casein percentages), and moderate to high correlations between composition and technological traits were observed. All factors included in the statistical model were significant in explaining the variation of the studied traits except for region effect in the analysis of casein and somatic cell score. Northeast and central-southern Italian regions showed the best performance for composition and technological features, respectively. Traits varied greatly across regions, which could reflect differences in herd management and strategies. Overall, less suitable milk for dairy processing was observed in summer. Results of the present study suggested that a constant monitoring of technological traits in the dairy industry is necessary to improve production quality at herd level and it may be a way to segregate milk according to its processing characteristics.     

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.     

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.     

Short communication: Genetic correlation and heritability of milk coagulation traits within and across lactations in Holstein cows using multiple-lactation random regression animal models

Genetic parameters of milk rennet coagulation time (RCT) and curd firmness (a₃₀) among the first 3 lactations in Holstein cows were estimated. The data set included 39,960 test-day records from 5,216 Estonian Holstein cows (the progeny of 306 sires), which were recorded from April 2005 to May 2010 in 98 herds across the country. A multiple-lactation random regression animal model was used. Individual milk samples from each cow were collected during routine milk recording. These samples were analyzed for milk composition and coagulation traits with intervals of 2 to 3 mo in each lactation (7 to 305 DIM) and from first to third lactation. Mean heritabilities were 0.36, 0.32, and 0.28 for log-transformed RCT [ln(RCT)] and 0.47, 0.40, and 0.62 for a₃₀ for parities 1, 2, and 3, respectively. Mean repeatabilities for ln(RCT) were 0.53, 0.55, and 0.56, but 0.59, 0.61, and 0.68 for a₃₀ for parities 1, 2 and 3, respectively. Mean genetic correlations between ln(RCT) and a₃₀ were −0.19, −0.14, and 0.02 for parities 1, 2, and 3, respectively. Mean genetic correlations were 0.91, 0.79, and 0.99 for ln(RCT), and 0.95, 0.94, and 0.94 for a₃₀ between parities 1 and 2, 1 and 3, and 2 and 3, respectively. Due to these high genetic correlations, we concluded that for a proper genetic evaluation of milk coagulation properties it is sufficient to record RCT and a₃₀ only in the first lactation.