Milk composition and its technological properties are traits of interest for the dairy sheep industry because almost all milk produced is processed into cheese. However, several variables define milk technological properties and a complex correlation pattern exists among them. In the present work, we measured milk composition, coagulation properties, and individual cheese yields in a sample of 991 Sarda breed ewes in 47 flocks. The work aimed to study the correlation pattern among measured variables and to obtain new synthetic indicators of milk composition and cheese-making properties. Multivariate factor analysis was carried out on individual measures of milk coagulation parameters; cheese yield; fat, protein, and lactose percentages; somatic cell score; casein percentage; NaCl content; pH; and freezing point. Four factors that were able to explain about 76% of the original variance were extracted. They were clearly interpretable: the first was associated with composition and cheese yield, the second with udder health status, the third with coagulation, and the fourth with curd characteristics. Factor scores were then analyzed by using a mixed linear model that included the fixed effect of parity, lambing month, and lactation stage, and the random effect of flock-test date. The patterns of factor scores along lactation stages were coherent with their technical meaning. A relevant effect of flock-test date was detected, especially on the 2 factors related to milk coagulation properties. Results of the present study suggest the existence of a simpler latent structure that regulates relationships between variables defining milk composition and coagulation properties in sheep. Heritability estimates for the 4 extracted factors were from low to moderate, suggesting potential use of these new variables as breeding goals.

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.

The suitability of multivariate factor analysis (MFA) to extract a small number of latent variables\ud able to explain the correlation pattern among fatty acids (FA) in buffalo milk was evaluated. FA\ud profile of milk samples from 214 Italian water buffaloes was analysed by gas chromatography.\ud MFA, performed on the correlation matrix of 52 FA, was able to extract 10 latent factors with\ud specific biological meaning related to a common metabolic origin for FA associated with the same\ud factor. Scores of the factors were treated as new quantitative phenotypes to evaluate the effect of\ud age, month of calving and lactation stage. MFA approach was effective in describing the FA profile\ud of buffalo milk by using a low number of new latent variables that clustered FA having similar\ud metabolic origin and function. The new variables were also useful to test the effect of\ud environmental and individual animal factors on milk FA composition

Fatty acid (FA) composition is a key component of sheep milk nutritional quality. However, breeding for FA composition in dairy sheep is hampered by the logistic and phenotyping costs. This study was aimed at estimating genetic parameters for sheep milk FA and to test the feasibility of their routine measurement by using Fourier-transform IR (FTIR) spectroscopy. Milk FA composition of 989 Sarda ewes farmed in 48 flocks was measured by gas chromatography (FAGC). Moreover, FTIR spectrum was collected for each sample, and it was used to predict FA composition (FAFTIR) by partial least squares regression: 700 ewes were used for estimating model parameters, whereas the remaining 289 animals were used to validate the model. One hundred replicates were performed by randomly assigning animals to estimation and validation data set, respectively. Variance components for both measured and predicted traits were estimated with an animal model that included the fixed effects of parity, days in milking interval, lambing month, province, altitude of flock location, the random effects of flock-test-date and animal genetic additive. Genetic correlations among FAGC, and between corresponding FAGC and FAFTIR were estimated by bivariate analysis. Coefficients of determination between FAGC and FAFTIR ranged from moderate (about 0.50 for odd- and branched-chain FA) to high (about 0.90 for de novo FA) values. Low-to-moderate heritabilities were observed for individual FA (ranging from 0.01 to 0.47). The largest value was observed for GC measured C16:0. Low-to-moderate heritabilities were estimated for FA groups. In most of cases, heritabilites were slightly larger for FAGC than FAFTIR. Estimates of genetic correlations among FAGC showed a large variability in magnitude and sign. The genetic correlation between FAFTIR and FAGC was higher than 60% for all investigated traits. Results of the present study confirm the existence of genetic variability of the FA composition in sheep and suggest the feasibility of using FAFTIR as proxies for these traits in large scale breeding programs.

Simple SummaryConsidered that all sheep milk in Italy is destined for cheese processing, traits describing rennet coagulation aptitude should be among the most important selection goals for dairy breeds. To reduce the costs and logistics related to the large-scale recording of these traits, mid-infrared (MIR) spectroscopy could be conveniently used to generate reliable predictions without any additional cost. The aims of this research were to predict the milk coagulation properties (MCP) and individual cheese yield (ILCY) in sheep by MIR spectrometry using partial least squares regression (PLS), and to compare different data pre-treatment procedures. The prediction results observed in the present study, although moderate, suggest the possibility of adding novel phenotypes (e.g., MCP and ILCY) in breeding schemes for dairy sheep breeds. Mid-infrared spectroscopy coupled with PLS regression could allow the prediction of phenotypes at the population level without additional costs.AbstractThe objectives of this study were (i) the prediction of sheep milk coagulation properties (MCP) and individual laboratory cheese yield (ILCY) from mid-infrared (MIR) spectra by using partial least squares (PLS) regression, and (ii) the comparison of different data pre-treatments on prediction accuracy. Individual milk samples of 970 Sarda breed ewes were analyzed for rennet coagulation time (RCT), curd-firming time (k20), and curd firmness (a30) using the Formagraph instrument; ILCY was measured by micro-manufacturing assays. An Furier-transform Infrared (FTIR) milk-analyzer was used for the estimation of the milk gross composition and the recording of MIR spectrum. The dataset (n = 859, after the exclusion of 111 noncoagulating samples) was divided into two sub-datasets: the data of 700 ewes were used to estimate prediction model parameters, and the data of 159 ewes were used to validate the model. Four prediction scenarios were compared in the validation, differing for the use of whole or reduced MIR spectrum and the use of raw or corrected data (locally weighted scatterplot smoothing). PLS prediction statistics were moderate. The use of the reduced MIR spectrum yielded the best results for the considered traits, whereas the data correction improved the prediction ability only when the whole MIR spectrum was used. In conclusion, PLS achieves good accuracy of prediction, in particular for ILCY and RCT, and it may enable increasing the number of traits to be included in breeding programs for dairy sheep without additional costs and logistics.