Milk test-day records of 5728 lactations of Italian Simmental cows were analyzed with multivariate factor analysis in order to extract 2 common factors, whose scores were used as quantitative measures of 2 main features of lactation curve shape-i.e., the increasing rate of yield in the first part of lactation and the rate of decline of milk yield after the lactation peak. The 2 indices, objectively derived from the correlation matrix of original test-day records, showed a high discriminant power in separating lactation curves with different shapes. The weak correlation between the 2 factors (0.11), together with the high correlation of factors and the total 305-d yield (about 0.70), suggests that an increase in lactation yield could be achieved by acting only on one of the 2 factors related to lactation-curve shape, with the other kept constant at a medium or low value. The suitability of the 2 factors as descriptors of lactation patterns has been confirmed by the relationships found between factor scores and the main environmental effects known to affect the shape of the lactation curve, such as parity and season of calving.

Heat stress represents a key factor that negatively affects the productive and reproductive performance of farm animals. In the present work, a new measure of tolerance to heat stress for dairy cattle was developed using principal component analysis. Data were from 590,174 test-day records for milk yield, fat and protein percentages, and somatic cell score of 39,261 Italian Holstein cows. Test-day records adjusted for main systematic factors were grouped into 11 temperature-humidity index (THI) classes. Daughter trait deviations (DTD) were calculated for 1,540 bulls as means of the adjusted test-day records for each THI class. Principal component analysis was performed on the DTD for each bull. The first 2 principal components (PC) explained 42 to 51% of the total variance of the system across the 4 traits. The first PC, a measure of the level at which the curve is located, was interpreted as a measure of the level at which the DTD curve was located. The second PC, which shows the slope of increasing or decreases DTD curves, synthesized the behavior of the DTD pattern. Heritability of the 2 component scores was moderate to high for level across all traits (range = 0.23-0.82) and low to moderate for slope (range = 0.16-0.28). For each trait, phenotypic and genetic correlations between level and slope were equal to zero. A genome-wide association analysis was carried out on a subsample of 423 bulls genotyped with the Illumina 50K bovine bead chip (Illumina, San Diego, CA). Two single nucleotide polymorphisms were significantly associated with slope for milk yield, 4 with level for fat percentage, and 2 with level and slope of protein percentage, respectively. The gene discovery was carried out considering windows of 0.5 Mb surrounding the significant markers and highlighted some interesting candidate genes. Some of them have been already associated with the mechanism of heat tolerance as the heat shock transcription factor (HSF1) and the malonyl-CoA-acyl carrier protein transacylase (MCAT). The 2 PC were able to describe the overall level and the slope of response of milk production traits across increasing levels of THI index. Moreover, they exhibited genetic variability and were genetically uncorrelated. These features suggest their use as measures of thermotolerance in dairy cattle breeding schemes.

Objective of this study was to estimate genetic parameters of milk coagulation properties (MCPs) and individual laboratory cheese yield (ILCY) in a sample of 1018 Sarda breed ewes farmed in 47 flocks. Rennet coagulation time (RCT), curd-firming time (k 20) and curd firmness (a 30) were measured using Formagraph instrument, whereas ILCY were determined by a micromanufacturing protocol. About 10% of the milk samples did not coagulate within 30 min and 13% had zero value for k 20. The average ILCY was 36%. (Co)variance components of considered traits were estimated by fitting both single- and multiple-trait animal models. Flock-test date explained from 13% to 28% of the phenotypic variance for MCPs and 26% for ILCY, respectively. The largest value of heritability was estimated for RCT (0.23±0.10), whereas it was about 0.15 for the other traits. Negative genetic correlations between RCT and a 30 (-0.80±0.12), a 30 and k 20 (-0.91±0.09), and a 30 and ILCY (-0.67±0.08) were observed. Interesting genetic correlations between MCPs and milk composition (r G>0.40) were estimated for pH, NaCl and casein. Results of the present study suggest to use only one out of three MCPs to measure milk renneting ability, due to high genetic correlations among them. Moreover, negative correlations between ILCY and MCPs suggest that great care should be taken when using these methods to estimate cheese yield from small milk samples.