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: statistical models for the analysis of coagulation traits using coagulating and noncoagulating milk information

Aims of this study were to propose statistical models for the analysis of rennet coagulation time (RCT) suitable for making use of coagulating and noncoagulating (NC) milk information, to estimate heritabilities and to obtain rank correlations for sire merit. A total of 1,025 Holstein cows (progeny of 54 sires) reared in 34 herds were milk-sampled once. Data were analyzed using 4 alternative models: a standard linear (SLM), a right-censored linear Gaussian (CLM), a survival (SUM), and a threshold (THM) model. Model SLM analyzed coagulated milk records only, whereas analysis with CLM or SUM considered information of NC samples as censored records. Model THM analyzed occurrence of milk coagulation as a dichotomous trait. An artificial censoring scenario with an endpoint at 18 min (SET18) was considered after the rearrangement of the timeframe originally used for the observation of RCT (SET31). Heritabilities ranged from 0.12 to 0.25. Correlations of sire rankings ranged from 0.23 to 0.92. Differences in sire rankings between SLM and CLM or SUM increased when the proportion of NC records increased. Correlations between sire rankings obtained for SET31 and SET18 were high for CLM and SUM, indicating that rankings provided by these models tended to be stable even when a large fraction of samples with observed RCT was re-classified as NC milk. Results indicate that CLM and SUM are more suitable than SLM and THM for the analysis of coagulation ability when data contain NC milk information.     

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.     

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.     

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.     

Milk, fat, protein, somatic cell score, and days open among Holstein, Brown Swiss, and their crosses

The objectives of this study were to compare Holstein (HO), Brown Swiss (BS), and their crosses for milk, fat, and protein yields, somatic cell score (SCS), days open (DO), and age at first calving (AFC), and to estimate the effects of heterosis and recombination. First through fifth lactation records were obtained from 19 herds milking crosses among BS and HO. The edited data set included 6,534 lactation records from 3,473 cows of the following breed combinations: 2,125 pure HO, 926 pure BS, 256 BS sire × HO dam (SH), 105 backcrosses to BS (SX), 18 HO sire × BS dam, and 43 backcrosses to HO. Least squares means for daily milk, fat, and protein yields, mature-equivalent milk, fat, and protein yields, SCS, DO, and AFC were calculated for breed combinations with a model that included fixed effects of age within parity (except for AFC), days in milk for daily yield and SCS, herd-year-season of calving, and breed combination. Cow and error were random effects. Breed combination was replaced with regressions on coefficients for heterosis and recombination in a second analysis. Last, data were analyzed with a 5-trait animal model that included a single pedigree file for both breeds and coefficients for heterosis and recombination. The least squares means for fat production were 1.21, 1.15, 1.27, and 1.16 kg for HO, BS, SH, and SX, respectively, which corresponds to a heterosis estimate of 7.30% and a recombination estimate of −3.76%. Heterosis and recombination estimates for protein production were 5.63% and −3.31%, respectively. Heterosis estimates increased for fat yield (10.38%) and protein yield (7.07%) when maternal grandsire identification from a known artificial insemination sire was required. Regression coefficients indicated an 11.44-d reduction in DO due to heterosis. Heterosis estimates for SCS were inconsistent. Regression on heterosis for SCS was significant and favorable (−0.22) when the breed of sire was BS, but nonsignificant and unfavorable when sire breed was HO (0.43). Heterosis estimates were favorable for all traits, whereas recombination effects tended to be unfavorable for yield traits. Reduced performance of future generations did not appear to be the result of inseminating crossbred cows with inferior sires. Results indicated that first-generation crosses among BS and HO compared favorably with HO. Yield in subsequent generations was somewhat below expectations, perhaps due to recombination loss in HO.     

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

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

Composition and effect of blending of noncoagulating, poorly coagulating, and well-coagulating bovine milk from individual Danish Holstein cows

The aim of the present investigation was to study the underlying causes of noncoagulating (NC) milk. Based on an initial screening in a herd of 53 Danish Holstein-Friesians, 20 individual Holstein-Friesian cows were selected for good and poor chymosin-induced coagulation properties; that is, the 10 cows producing milk with the poorest and best coagulating properties, respectively. These 20 selected cows were followed and resampled on several occasions to evaluate possible changes in coagulation properties. In the follow-up study, we found that among the 10 cows with the poorest coagulating properties, 4 cows consistently produced poorly coagulating (PC) or NC milk, corresponding to a frequency of 7%. Noncoagulating milk was defined as milk that failed to form a coagulum, defined as increase in the storage modulus (G′) in oscillatory rheometry, within 45 min after addition of chymosin. Poorly coagulating milk was characterized by forming a weak coagulum of low G′. Milk proteomic profiling and contents of different casein variants, ionic contents of Ca, P and Mg, κ-casein (CN) genotypes, casein micelle size, and coagulation properties of the 4 NC or PC samples were compared with milk samples of 4 cows producing milk with good coagulation properties. The studies included determination of production of caseinomacropeptide to ascertain whether noncoagulation could be ascribed to the first or second phase of chymosin-induced coagulation. Caseinomacropeptide was formed in all 8 milk samples after addition of chymosin, indicating that the first step (cleavage of κ-CN) was not the cause of inability to coagulate. Furthermore, the effect of mixing noncoagulating and well-coagulating milk was studied. By gradually blending NC with well-coagulating milk, the coagulation properties of the well-coagulating samples were compromised in a manner similar to titration. Milk samples from cows that consistently produced NC milk were further studied at the udder quarter level. The coagulation properties of the quarter milk samples were not significantly different from those of the composite milk sample, showing that poor coagulation traits and noncoagulation traits of the composite milk were not caused by the milk quality of a single quarter. The milk samples exhibiting PC or NC properties were all of the κ-CN variant AA genotype, and contained casein micelles with a larger mean diameter and a lower fraction of κ-CN relative to total CN than milk with good coagulation properties. Interestingly, the relative proportions of different phosphorylation forms of α-CN differed between well-coagulating milk and PC or NC milk samples. The PC and NC milk samples contained a lower proportion of the 2 less-phosphorylated variants of α-CN (α_S1-CN-8P and α_S2-CN-11P) compared with samples of milk that coagulated well.     

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.     

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.     

Recursive relationships between milk yield and somatic cell score of Canadian Holsteins from finite mixture random regression models

Finite mixture, multiple-trait, random regression animal models with recursive links between phenotypes for milk yield and somatic cell score (SCS) on the same test-day were applied to first lactation Canadian Holstein data. All models included fixed herd-test-day effects and fixed regressions within region-age at calving-season of calving classes, and animal additive genetic and permanent environmental regressions with random coefficients. Causal links between phenotypes for milk yield and SCS were fitted separately for records from healthy cows and cows with a putative, subclinical form of mastitis. Bayesian methods via Gibbs sampling were used for the estimation of model parameters. Bayes factors indicated superiority of the model with recursive link from milk to SCS over the reciprocal recursive model and the standard multiple-trait model. Differences between models measured by other, single-trait model comparison criteria (i.e., weighted mean squared error, squared bias, and correlation between observed and expected data) were negligible. Approximately 20% of test-day records were classified as originating from cows with mastitis in recursive mixture models. The proportion of records from cows infected with mastitis was largest at the beginning of lactation. Recursive mixture models exhibited different distributions of data from healthy and infected cows in different parts of lactation. A negative effect of milk to SCS (up to −0.15 score points for every kilogram of milk for healthy cows from 5 to 45 d in milk) was estimated for both mixture components (healthy and infected) in all stages of lactation for the most plausible model. The magnitude of this effect was stronger for healthy cows than for cows infected with mastitis. Different patterns of genetic and environmental correlations between milk and SCS for healthy and infected records were revealed, due to heterogeneity of structural coefficients between mixture components. Estimated breeding values for SCS from the best fitting model for sires of infected daughters were more related to estimated breeding values for the same trait from the regular multiple-trait model than evaluations for sires of mastitis-free cows.     

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.     

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.     

A comparison between Holstein-Friesian and Jersey dairy cows and their F1 cross with regard to milk yield, somatic cell score, mastitis, and milking characteristics under grazing conditions

The key objectives of this study were to investigate differences in milking characteristics and udder health between Holstein-Friesian (HF), Jersey (J), and Jersey × Holstein-Friesian (F₁) cows and to determine possible associations between milking characteristics and udder health. Records were available from 329 lactations (162 cows): 65 HF, 48 J, and 49 F₁. Data included lactation mean milk yield, somatic cell score (SCS), incidence of mastitis, average milk flow (AMF), peak milk flow (PMF), and milking duration (MD). Breed group had a significant effect on milk yield and was higher with the HF cows (18.0 kg/d) compared with the J cows (14.2 kg/d). Udder health (SCS and incidence of mastitis at least once during lactation) were similar across the breed groups. Average milk flow was greater with the HF cows (1.36 kg/min) compared with the J cows (1.09 kg/min). Peak milk flow also tended to greater with the HF cows. No difference in MD was observed between the breed groups. The performance of the F₁ cows tended to be similar to the mid-parent (breed) mean for udder health and MD, but heterosis was evident for milk yield, AMF, and PMF. Correlations examined showed that phenotypic milk yield was negatively associated with SCS. Increased milk yield was synonymous with increased AMF, PMF, and MD. Correlations between SCS and milking characteristics were weak. Correlations also showed that cows with low AMF and PMF had extended MD. Therefore, no difference in udder health was observed between HF, J, or F₁ cows. The fact that higher yielding animals exhibit faster milking speeds was confirmed; however, no difference in MD was observed between the breed groups. Such findings indicate that regularity in the milking process will be maintained within mixed-breed herds.     

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.     

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.     

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.     

Analysis of the relationship between somatic cell score and functional longevity in Canadian dairy cattle

The aim of this study was to assess the level of somatic cell count (SCC) and to explore the impact of somatic cell score (SCS) on the functional longevity of Canadian dairy cattle by using a Weibull proportional hazards model. Data consisted of 1,911,428 cows from 15,970 herds sired by 7,826 sires for Holsteins, 80,977 cows in 2,036 herds from 1,153 sires for Ayrshires, and 53,114 cows in 1,372 herds from 1,758 sires for Jerseys. Functional longevity was defined as the number of days from the first calving to culling, death, or censoring. The test-day SCC was transformed to a linear score, and the resulting SCS were averaged within each lactation. The average SCS were grouped into 10 classes. The statistical model included the effects of stage of lactation; season of production; annual change in herd size; type of milk recording supervision; age at first calving; effects of milk, fat, and protein yields, calculated as within-herd-year-parity deviations; herd-year-season of calving; SCS class; and sire. The relative culling rate was calculated for animals in each SCS class after accounting for the aforementioned effects. The overall average SCC for Holsteins was 167,000 cells/mL, for Ayrshires was 155,000 cells/mL, and for the Jerseys was 212,000 cells/mL. In all breeds there were no appreciable differences in the relative risk of culling among classes of SCS breed averages (i.e., up to a SCS of 5). However, as the SCS increased beyond the breed average, the relative risk of cows being culled increased considerably. For instance, Holstein, Ayrshire, and Jersey cows with the highest classes of SCS had, respectively, a 4.95, 6.73, and 6.62 times greater risk of being culled than cows with average SCS.     

Effect of Holstein Friesian and Brown Swiss breeds on quality of milk and cheese

In Italy, more than 75% of milk is used for cheese making. For this reason, milk composition and coagulation traits and cheese quality represent the most important tools for the economic development of the dairy sector. In particular, cheese quality varies in relation to cheese-making technology and breed of cow. The aim of this study was to investigate the effect of 3 types of milk, originating from Holstein-Friesian (HF), Brown Swiss (BS), and mixed of both breeds, on vat milk characteristics, cheese yield, and quality in 3 different typical Italian cheese-making conditions (Casolet, Vezzena, and Grana Trentino). One hundred forty-four cows (66 HF and 78 BS) were involved, and a total of 24 vats of milk were evaluated. At maturity, 30, 21, and 16 wheels of Casolet, Vezzena, and Grana Trentino cheese were analyzed. Brown Swiss cows yielded 9% less milk per day than HF cows, but milk showed greater contents of protein, casein, titratable acidity, and better rennet coagulation time and curd firmness than HF milk. The chemical composition and cholesterol content of the 3 types of cheese were similar between breeds, whereas the cheese made with BS milk showed greater contents of monounsaturated and polyunsaturated fatty acids. Cheese made with BS milk had greater b* (yellow component) than HF. Cheese yield, recorded at different ripening times, demonstrated that BS milk yielded more cheese than HF. Mixed milk showed values, on average, intermediate to HF and BS milk characteristics, and this trend was confirmed in cheese yield at different ripening times.     

Effects of breed, feeding system, and parity on udder health and milking characteristics

The objective of this study was to quantify differences in udder health and milking characteristics among the Holstein-Friesian (HF), Montbéliarde (MB), Normande (NM), Norwegian Red (NRF), Montbéliarde × Holstein-Friesian (MBX), and Normande × Holstein-Friesian (NMX) genotypes, while considering the effect of feeding system and parity. A total of 749 lactations were available for inclusion in the analysis from 309 cows in 1 research herd over 5 yr. Somatic cell score (SCS; i.e., natural logarithm of somatic cell count) was used as an indicator of udder health. Milking duration (seconds/d) was defined as the sum of the milking duration in the a.m. and milking duration in the p.m. Average daily milk flow (AMF; kg/min) was defined as total daily milk yield divided by total daily milking duration. Peak milk flow (kg/min) was defined as the maximum rate of milk flow achieved in the daily milking process. The SCS of the NRF (10.31 units) and MB (10.47 units) breeds was less than that of the HF (10.96 SCS units), whereas that of the NM (10.88 SCS units), MBX (10.93 SCS units), and NMX (10.84 SCS units) breeds was similar to that of the HF. The MBX and NMX had the greatest AMF (1.56 and 1.54 kg/min, respectively) and the NM had the lowest (1.33 kg/min). Animals offered a high concentrate diet had greater AMF, peak milk flow, and milking duration. The differences expressed by the divergent breeds may reflect differences in the past breeding goals among the breeds, namely the inclusion of traits aimed at maintaining or improving udder health.