Noncoagulation of milk in Finnish Ayrshire and Holstein-Friesian cows and effect of herds on milk coagulation ability

The objectives of this study were to compare milk coagulation ability (MCA) and the prevalence of noncoagulation of milk within the main Finnish dairy breeds, Finnish Ayrshire (FA) and Holstein-Friesian (HOL), as well as to study the herd effect on MCA. Data used in the statistical analyses consisted of individual milk samples of 959 FA, 399 HOL, and 50 crossbred cows from 84 herds. Data were collected before the grazing season in the spring 1999. Milk samples were analyzed for the milk coagulation traits (milk renneting time, R and curd firmness, E(30)) and pH. In addition, information on the 305-d milk production traits from the year 1999, and background information about feeding and management regimes of the herds were obtained. Variance components for the random herd and animal effects were estimated using REML methodology and an animal model. Breed, parity, lactation stage (for R, E(30) and pH only), and a measuring unit (for R and E(30) only) were included as fixed effects in the model. When the effects of concentrate feeding frequency and type of concentrate were studied, the random effect of herd was excluded from the model. A relationship matrix included parents, grandparents, and great grandparents of the cows with observations. The HOL cows were superior to FA cows in MCA when both the proportion of poorly coagulating (PC) and noncoagulating (NC) milk, and the differences in curd firmness were considered. About 30% of the FA cows and 12% of the HOL cows produced PC milk. Only 1.3% of the HOL cows and 8.6% of the FA cows produced NC milk. Herd effect explained only a minor part of the variation in MCA (8%) compared with that in 305-d milk production traits (about 43%). Frequent feeding of the concentrate was associated with good MCA as well as for the high milk, protein and fat yields, but it was not associated with the prevalence of the NC milk.     

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.     

Short communication: milk fat composition of 4 cattle breeds in the Netherlands

Milk fatty acid (FA) composition was compared among 4 cattle breeds in the Netherlands: Dutch Friesian (DF; 47 animals/3 farms), Meuse-Rhine-Yssel (MRY; 52/3), Groningen White Headed (GWH; 45/3), and Jersey (JER; 46/3). Each cow was sampled once between December 2008 and March 2009 during the indoor housing season, and samples were analyzed using gas chromatography. Significant breed differences were found for all traits including fat and protein contents, 13 major individual FA, 9 groups of FA, and 5 indices. The saturated fatty acid proportion, which is supposed to be unfavorable for human health, was smaller for GWH (68.9%) compared with DF (74.1%), MRY (72.3%), and JER (74.3%) breeds. The proportion of conjugated linoleic acid and the unsaturation index, which are associated positively with human health, were both highest for GWH. Differences in milk fat composition can be used in strategies to breed for milk with a FA profile more favorable for human health. Our results support the relevance of safeguarding the local Dutch breeds.     

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.     

Lead excretion in milk of accidentally exposed dairy cattle

Lead (Pb) exposure in dairy cattle is associated with economic losses due to mortality and treatment costs, but with production animals there is also risk to the human food chain. The first objective of this study was to quantify the Pb concentration in milk from Pb-exposed cattle. The second objective was to correlate blood and milk Pb concentrations from individual cows. The third objective was long-term monitoring to determine the duration of milk contamination after exposure ceased. A dairy herd of more than 100 cows was accidentally exposed to Pb-contaminated feed. Milk and blood were collected for Pb analysis. Serial collection of milk samples continued for 2.5 years. The initial concentration of Pb in bulk tank milk was 0.0999 mg l^–1. The highest milk Pb concentration from an individual cow was 0.4657 mg l^–1 and the highest blood Pb concentration was 1.216 mg l^–1. One milk sample collected at the end of the study (day 922) contained 0.0117 mg Pb l^–1 of Pb. The calculated relationship between milk (y) and blood (x) Pb concentration was ln(y) = 3.4(x) – 2.21 (R² = 0.98).     

Genetic analysis of herd life in Canadian dairy cattle on a lactation basis using a Weibull proportional hazards model

The objectives of this study were to identify the most important factors that influence functional survival and to estimate the genetic parameters of functional survival for Canadian dairy cattle. Data were obtained from lactation records extracted for the May 2002 genetic evaluation of Holstein, Jersey, and Ayrshire breeds that calved between July 1, 1985 and April 5, 2002. Analysis was performed using a Weibull proportional hazard model, and the baseline hazard function was defined on a lactation basis instead of the traditional analysis of the whole length of life. The statistical model included the effects of stage of lactation; season of production; the annual change in herd size; type of milk recording supervision; age at first calving; effects of milk, fat, and protein yields calculated within herd-year-parity deviations; and the random effects of herd-year-season of calving and sire. All effects fitted in the model had a significant effect on functional survival of cows in all breeds. Milk yield was by far the most important factor influencing survival, and the hazard increased as the milk production of the cows decreased. The hazard also increased as the fat content increased compared with the average group. Heifers that were older at calving were at higher risk of being culled, and expanding herds were at lower risk of being culled compared with stable herds. More culling was found in unsupervised herds than in supervised herds. The heritability values obtained were 0.14, 0.10, and 0.09 for Holstein, Jersey, and Ayrshire, respectively. Rank correlation between estimated breeding values (EBV) obtained from the current national genetic evaluation of direct herd life and the survival kit used in this study ranged from 0.65 to 0.87, depending on the number of daughters per sire. Estimated genetic trend obtained using the survival kit was overestimated.     

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.     

Selection with inbreeding control in simulated young bull schemes for local dairy cattle breeds

Local breeds are rarely subject to modern selection techniques; however, selection programs will be required if local breeds are to remain a viable livelihood option for farmers. Selection in small populations needs to take into account accurate inbreeding control. Optimum contribution selection (OCS) is efficient in controlling inbreeding and maximizes genetic gain. The current paper investigates genetic progress in simulated dairy cattle populations from 500 to 6,000 cows undergoing young bull selection schemes with OCS compared with truncation selection (TS) at an annual inbreeding rate of 0.003. Selection is carried out for a dairy trait with a base heritability of 0.3. A young bull selection scheme was used because of its simplicity in implementation. With TS, annual genetic gain from 0.111 standard deviation units with 500 cows increases rapidly to 0.145 standard deviation units with 4,000 cows. Then, genetic gain increases more slowly up to 6,000 cows. At the same inbreeding rate, OCS produces higher genetic progress than TS. Differences in genetic gain between OCS and TS vary from to 2 to 6.3%. Genetic gain is also improved by increasing the number of years that males can be used as sires of sires. When comparing OCS versus TS at different heritabilities, we observe an advantage of OCS only at high heritability, up to 8% with heritability of 0.9. By increasing the constraint on inbreeding, the difference of genetic gain between the 2 selection methods increases in favor of OCS, and the advantage at the inbreeding rate of 0.001 per generation is 6 times more than at the inbreeding rate of 0.003. Opportunities exist for selection even in dairy cattle populations of a few hundred females. In any case, selection in local breeds will most often require specific investments in infrastructure and manpower, including systems for accurate data recording and selection skills and the presence of artificial insemination and breeders organizations. A cost-benefit analysis is therefore advisable before considering the implementation of selection schemes in local dairy cattle breeds.     

Estimation of heritability for milk urea and genetic correlations with milk production traits in 3 Danish dairy breeds

The aim of this study was to estimate genetic parameters for milk urea (MU) content in 3 main Danish dairy breeds. As a part of the Danish milk recording system, milk samples from cows on commercial farms were analyzed for MU concentration (mmol/L) and the percentages of fat and protein. There were 323,800 Danish Holstein, 70,634 Danish Jersey, and 27,870 Danish Red cows sampled with a total of 1,436,580, 368,251, and 133,922 test-day records per breed, respectively, included in the data set. Heritabilities for MU were low to moderate (0.22, 0.18, and 0.24 for the Holstein, Jersey, and Red breeds, respectively). The genetic correlation was close to zero between MU and milk yield in Jersey and Red, and −0.14 for Holstein. The genetic correlations between MU and fat and protein percentages, respectively, were positive for all 3 dairy breeds. Herd-test-day explained 51%, 54%, and 49% of the variation in MU in Holstein, Jersey, and Red, respectively. This indicates that MU levels in milk can be reduced by farm management. The current study shows that there are possibilities to influence MU by genetic selection as well as by farm management.     

Efficiency of multi-breed genomic selection for dairy cattle breeds with different sizes of reference population

Single-breed genomic selection (GS) based on medium single nucleotide polymorphism (SNP) density (~50,000; 50K) is now routinely implemented in several large cattle breeds. However, building large enough reference populations remains a challenge for many medium or small breeds. The high-density BovineHD BeadChip (HD chip; Illumina Inc., San Diego, CA) containing 777,609 SNP developed in 2010 is characterized by short-distance linkage disequilibrium expected to be maintained across breeds. Therefore, combining reference populations can be envisioned. A population of 1,869 influential ancestors from 3 dairy breeds (Holstein, Montbéliarde, and Normande) was genotyped with the HD chip. Using this sample, 50K genotypes were imputed within breed to high-density genotypes, leading to a large HD reference population. This population was used to develop a multi-breed genomic evaluation. The goal of this paper was to investigate the gain of multi-breed genomic evaluation for a small breed. The advantage of using a large breed (Normande in the present study) to mimic a small breed is the large potential validation population to compare alternative genomic selection approaches more reliably. In the Normande breed, 3 training sets were defined with 1,597, 404, and 198 bulls, and a unique validation set included the 394 youngest bulls. For each training set, estimated breeding values (EBV) were computed using pedigree-based BLUP, single-breed BayesC, or multi-breed BayesC for which the reference population was formed by any of the Normande training data sets and 4,989 Holstein and 1,788 Montbéliarde bulls. Phenotypes were standardized by within-breed genetic standard deviation, the proportion of polygenic variance was set to 30%, and the estimated number of SNP with a nonzero effect was about 7,000. The 2 genomic selection (GS) approaches were performed using either the 50K or HD genotypes. The correlations between EBV and observed daughter yield deviations (DYD) were computed for 6 traits and using the different prediction approaches. Compared with pedigree-based BLUP, the average gain in accuracy with GS in small populations was 0.057 for the single-breed and 0.086 for multi-breed approach. This gain was up to 0.193 and 0.209, respectively, with the large reference population. Improvement of EBV prediction due to the multi-breed evaluation was higher for animals not closely related to the reference population. In the case of a breed with a small reference population size, the increase in correlation due to multi-breed GS was 0.141 for bulls without their sire in reference population compared with 0.016 for bulls with their sire in reference population. These results demonstrate that multi-breed GS can contribute to increase genomic evaluation accuracy in small breeds.     

Comparing local and commercial breeds on functional traits and profitability: the case of Reggiana dairy cattle

The objective of this study was to compare fertility, longevity, milkability, and profitability of cows from the Reggiana and Holstein breeds in northern Italy. Profitability was gauged for each breed, with consideration of economic incentive programs and alternative milk pricing scenarios. Calving to first service interval, days open, and calving interval were significantly shorter in Reggiana than in Holstein cows. Reggiana cows conceived approximately one estrus cycle before Holstein and had a calving interval 33 d shorter. Holstein cows released a significantly higher quantity of milk per unit of time (1.81 vs. 1.28 kg/min). Reggiana cows had longer expected total and productive lives than Holstein cows, by 5.8 and 10.0 mo, respectively. Replacement rate was 26% higher in the Holstein. Standard 305-d milk production was 5,360 and 7,870 kg in Reggiana and Holstein, respectively. Comparing breeds on annual milk and meat production, instead of standard 305-d milk yield, changed marginally the difference in annual profitability between the Reggiana and Holstein, from −€696 to −€679 per cow per year. Including feeding, milking, replacement, and insemination costs reduced the gap between breeds by 32%, from −€679, measured on annual milk and meat production, to −€460. These differences in profitability assumed a pricing scenario referring to milk sold to the dairy industry where protein and fat contents are valued but not the breed origin of milk. Incentive payments to farmers of endangered cattle compensated partially (22%) the lower income from Reggiana cows. When Reggiana milk production was sold as branded Parmigiano Reggiano cheese, Reggiana cows were more profitable than Holstein cows by €1,953 per cow per year.     

Joint analysis of quantitative trait loci for clinical mastitis and somatic cell score on five chromosomes in three Nordic dairy cattle breeds

Five chromosomes were selected for joint quantitative trait loci (QTL) analyses for clinical mastitis (CM) and somatic cell score (SCS) in 3 breeds: Finnish Ayrshire (FA), Swedish Red and White (SRB), and Danish Red (DR). In total, 19 grandsires and 672 sons in FA, 19 grandsires and 499 sons in SRB, and 8 grandsires and 258 sons in DR were used in the study. These individuals were genotyped with the 61 microsatellite markers used in any of the previous QTL scans on the selected chromosomes. Within-family QTL analyses based on linear regression models were carried out for CM and SCS to identify the segregating sires for each region. On the segregating families, joint single-trait and 2-trait analyses were performed using variance components models. The analyses confirmed that QTL affecting CM or SCS, or both, segregate on Bos taurus autosomes (BTA) 9, 11, 14, and 18, whereas a QTL on BTA29 could not be confirmed. Our results indicate that there may be at least 2 linked QTL on BTA9, one that primarily affects CM and a second that primarily affects SCS. On chromosomes BTA11, 14, and 18, the joint analyses were only significant for SCS.     

Alterations in milk metabolome and coagulation ability during the lactation of dairy cows

Milk composition has been known to change during lactation. To help understand the changes in metabolic profile throughout the whole lactation, liquid chromatography mass-spectrometry was used to analyze 306 milk samples from 82 primi- and multiparous dairy cows. Changes in metabolic profile common to all cows throughout lactation were ascertained based on principal component and general linear model analysis. Sets of specific markers; for instance, 225, 397, and 641–642 m/z (positive mode), and 186, 241, and 601–604 (negative mode), with at least a 1.5-fold higher intensity during the first 60 d compared with the last 60 d of lactation were observed. The metabolome was affected by parity and milking time. Markers, identified as peptides differentiating parity, were observed. A significant increase for citrate was observed in evening milk. Milk coagulation traits were strongly animal specific. The curd firmness values were influenced by milking time. Sets of markers were associated with curd firmness in positive (197 m/z) and negative (612, 737, 835, 836, 902, 1000, 1038, and 1079 m/z) ion mode.     

Random herd curves in a test-day model for milk, fat, and protein production of dairy cattle in The Netherlands

With random regression models, genetic parameters of test-day milk production records of dairy cattle can be estimated directly from the data. However, several researchers that used this method have reported unrealistically high variances at the borders of the lactation trajectory and low genetic correlations between beginning and end of lactation. Recently, it has been proposed to include herd-specific regression curves in the random regression model. The objective was to study the effect of including random herd curves on estimated genetic parameters. Genetic parameters were estimated with 2 models; both included random regressions for the additive genetic and permanent environmental effect, whereas the second model also included a random regression effect for herd x 2-yr period of calving. All random regressions were modeled with fourth-order Legendre polynomials. Bayesian techniques with Gibbs sampling were used to estimate all parameters. The data set comprised 857,255 test-day milk, fat, and protein records from lactations 1, 2, and 3 of 43,990 Holstein cows from 544 herds. Genetic variances estimated by the second model were lower in the first 100 d and at the end of the lactation, especially in lactations 2 and 3. Genetic correlations between d 50 and the end of lactation were around 0.25 higher in the second model and were consistent with studies where lactation stages are modeled as different traits. Subsequently, estimated heritabilities for persistency were up to 0.14 lower in the second model. It is suggested to include herd curves in a random regression model when estimating genetic parameters of test-day production traits in dairy cattle.     

Extending the postpartum anovulatory period in dairy cattle with estradiol cypionate

Estradiol cypionate (ECP), a long-acting estrogen, has been used therapeutically in early postpartum (PP) dairy cows. In experiment 1, effects of ECP on circulating reproductive hormones, cyclicity, and ovarian function in early PP dairy cows were investigated. Lactating Holsteins received 10 mg of ECP (ECP; n = 17) or placebo (CON; n = 16) on d 7 PP. Serum and ultrasound data were acquired from 5 to 90 d of lactation. Compared to CON, ECP cows had greater serum estradiol for 10 d and lower serum FSH for 15 d posttreatment. After ECP, the appearance of follicles > or = 10 mm and time to first ovulation were delayed. Nevertheless, by 90 d PP, normal estrous cycles were found in only 50% of CON versus 88% of ECP cows. Primiparous, but not multiparous, cows receiving ECP had higher milk yields. Experiment 2 investigated effects of ECP on reproduction and milk production on a commercial dairy. Cows were blocked by parity and randomized to three treatments: 0 (n = 85), 4 (n = 85), or 10 (n = 86) mg of ECP on d 5 to 8 PP. Data included cycling status (two ultrasound examinations at 30 to 33 d PP and 7 d later), individual reproduction records, and daily milk yields from 10 to 90 d PP. In primiparous cows, ECP had no significant effects on ovulatory status or milk yields. By 40 d PP, a greater percentage of multiparous cows receiving 10 mg of ECP remained anovulatory compared with those receiving 0 or 4 mg. Milk yields were highest for multiparous cows receiving 4 mg of ECP, intermediate for the 10-mg dose, and lowest for controls. Lower conception was observed in multiparous cows receiving 4 mg of ECP. In summary, ECP delayed time to first ovulation particularly in multiparous cows, an effect associated with observed inhibition of circulating FSH. Milk yield responses to ECP were inconsistent within parity groups across the two experiments.     

Survival analysis as a statistical methodology for analyzing factors that affect milk coagulation time in Holstein-Friesian and Brown Swiss cows

The aim of this study was to propose the survival analysis technique as a statistical approach for the analysis of rennet coagulation time (RCT) able to make use of coagulating and noncoagulating (NC) milk information in order to estimate potential sources of variations that affects RCT. A total of 1,025 Italian Holstein-Friesian cows (HF; progeny of 54 sires) and 1,234 Brown Swiss cows (BS; progeny of 58 sires) reared in 34 and 38 herds, respectively, were milk-sampled once. Rennet coagulation time was analyzed with a semiparametric proportional hazard model (i.e., a Cox model), with the NC samples considered as censored records. Furthermore, a different censoring scenario, with a new end point at 18 min, was considered after the rearrangement of the time space originally used for the observation of RCT. The percentage of NC samples was almost 10% for HF and 3.5% for BS cows in in the 31 min set, whereas it increased to 44 and 24.9%, respectively, in the 18 min set. Estimated hazard ratios indicated that the most important factors affecting the coagulation process were herd, days in milk, casein number, and milk acidity (expressed in terms of titratable acidity) for both HF and BS, whereas the SCS was relevant only for BS. The survival model seems to be particularly suitable for this analysis, as it can properly account for censored and uncensored records and appropriately use all available information. Moreover, this methodology allows us to rearrange the time space used for the observation of RCT and to define alternative traits (i.e., RCT with an end point at 18 min). Our restriction of the time space and the increased percentage of censored records did not highlight any substantial differences in terms of the risk of coagulating with respect to the traditional 31 min testing time. Although further research is needed to investigate the effect of these sources of variation on cheese yield, our results indicate that casein number, acidity, and SCS may be used as indicator traits for enhancing the technological properties of bovine milk.     

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 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: Breed differences affecting dairy cattle welfare in traditional alpine tie-stall husbandry systems

The aim of this investigation was to compare the prevalence of indicators of poor welfare among 5 Italian cattle breeds (Italian Holstein-Friesian, Italian Bruna, Pezzata Rossa Italiana, Grigia Alpina, and Pezzata Rossa d’Oropa) kept in tie-stalls in the Italian Alps under similar housing and management conditions. We recorded the presence of integument alterations (hairless patch areas, lesion/swollen areas, or overgrown claws) and lameness in 612 cows. Additionally, we checked 834 cows for the presence of physical malformations (“open” shoulders). In general, the prevalence of welfare problems showed a decreasing trend from the more productive to the less productive breeds. Local breeds (Grigia Alpina and Pezzata Rossa d’Oropa) showed a significantly lower prevalence of welfare problems compared with the other 3 breeds, whereas Italian Holstein-Friesian usually had the highest percentage of individuals with problems. No differences were found between Pezzata Rossa Italiana and Italian Bruna, both of which showed fewer problems than Italian Holstein-Friesian. The effect of the breed significantly affected the welfare of dairy cows in tie-stalls in alpine traditional husbandry systems. The prevalence of the negative welfare indicators studied was lower in local breeds, which are better adapted to local breeding conditions. Our results indicate an urgent need to promote changes in the criteria used for genetic selection in the dairy industry and underline the importance of maintaining the diversity of local breeds, which should be carefully chosen for each specific environmental condition.     

Genetic parameters for milk coagulation properties in Estonian Holstein cows

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