Search

This study aimed to assess the environmental impact (cradle-to-farm gate LCA) of milk production of cereal-based dairy farms in northern Italy and to analyse which traits related to farm management (MAN_F), lactating cows’ ration composition (DIET_F), and animal responses (ANI_F) could be identified as main determinants to address specific mitigation strategies. Data originated from 28 dairy farms. The functional unit was 1 kg of fat- and protein-corrected milk (FPCM). Impact categories were global warming (GWP), acidification (AP) and eutrophication (EP) potentials and land occupation (LO). Impact values were analysed with a linear model to test separately the effects of MAN_F, DIET_F and ANI_F variables. One kg FPCM was associated with an average impact of 1.10 kg CO2-eq (GWP), 17.1 g SO2-eq (AP), 6.0 g PO4-eq (EP), and 1.3 m2/y (LO). Regarding MAN_F, GWP and LO linearly decreased with increasing values of stocking rate and feed self-sufficiency (FSS), whereas minimum EP was found for intermediate FSS values (50% DM). Regarding DIET_F, GWP linearly decreased at increasing levels of dietary inclusion of cereal silages, whereas AP and EP linearly decreased at decreasing dietary crude protein (CP) content values. Among ANI_F, GWP and EP were quadratically associated with milk yield (MY), evidencing the lowest values at intermediate MY (9100 ± 477 kg FPCM/cow per 305d-lactation). The GWP, AP and EP decreased with decreasing age at first calving. These results can contribute at drawing good practices for farmers and consultants to promote more environmentally sustainable dairy production while supporting farm functioning and farmers’ income.HIGHLIGHTSFarm and animal determinants of Life Cycle Assessment of cereal-based dairy farms were studied.Mitigation effects can be achieved by operating on farm organisation, cow’s diet, and animal responses.Good practices for farmers and consultants can be proposed for a better dairy production sustainability. Farm and animal determinants of Life Cycle Assessment of cereal-based dairy farms were studied. Mitigation effects can be achieved by operating on farm organisation, cow’s diet, and animal responses. Good practices for farmers and consultants can be proposed for a better dairy production sustainability.

To devise better selection strategies in dairy cattle breeding programs, a deeper knowledge of the role of the major genes encoding for milk protein fractions is required. The aim of the present study was to assess the effect of the CSN2, CSN3, and BLG genotypes on individual protein fractions (α

Background Blood metabolic profiles can be used to assess metabolic disorders and to evaluate the health status of dairy cows. Given that these analyses are time-consuming, expensive, and stressful for the cows, there has been increased interest in Fourier transform infrared (FTIR) spectroscopy of milk samples as a rapid, cost-effective alternative for predicting metabolic disturbances. The integration of FTIR data with other layers of information such as genomic and on-farm data (days in milk (DIM) and parity) has been proposed to further enhance the predictive ability of statistical methods. Here, we developed a phenotype prediction approach for a panel of blood metabolites based on a combination of milk FTIR data, on-farm data, and genomic information recorded on 1150 Holstein cows, using BayesB and gradient boosting machine (GBM) models, with tenfold, batch-out and herd-out cross-validation (CV) scenarios. Results The predictive ability of these approaches was measured by the coefficient of determination (R2). The results show that, compared to the model that includes only FTIR data, integration of both on-farm (DIM and parity) and genomic information with FTIR data improves the R2 for blood metabolites across the three CV scenarios, especially with the herd-out CV: R2 values ranged from 5.9 to 17.8% for BayesB, from 8.2 to 16.9% for GBM with the tenfold random CV, from 3.8 to 13.5% for BayesB and from 8.6 to 17.5% for GBM with the batch-out CV, and from 8.4 to 23.0% for BayesB and from 8.1 to 23.8% for GBM with the herd-out CV. Overall, with the model that includes the three sources of data, GBM was more accurate than BayesB with accuracies across the CV scenarios increasing by 7.1% for energy-related metabolites, 10.7% for liver function/hepatic damage, 9.6% for oxidative stress, 6.1% for inflammation/innate immunity, and 11.4% for mineral indicators. Conclusions Our results show that, compared to using only milk FTIR data, a model integrating milk FTIR spectra with on-farm and genomic information improves the prediction of blood metabolic traits in Holstein cattle and that GBM is more accurate in predicting blood metabolites than BayesB, especially for the batch-out CV and herd-out CV scenarios.

Dairy bull fertility is traditionally evaluated using semen production and quality traits; however, these attributes explain only part of the differences observed in fertility among bulls. Alternatively, bull fertility can be directly evaluated using cow field data. The main objective of this study was to investigate bull fertility in the Italian Brown Swiss dairy cattle population using confirmed pregnancy records. The data set included a total of 397,926 breeding records from 1,228 bulls and 129,858 lactating cows between first and fifth lactation from 2000 to 2019. We first evaluated cow pregnancy success, including factors related to the bull under evaluation, such as bull age, bull inbreeding, and AI organization, and factors associated with the cow that receives the dose of semen, including herd-year-season, cow age, parity, and milk yield. We then estimated sire conception rate using only factors related to the bull. Model predictive ability was evaluated using 10-fold cross-validation with 10 replicates. Interestingly, our analyses revealed that there is a substantial variation in conception rate among Brown Swiss bulls, with more than 20% conception rate difference between high-fertility and low-fertility bulls. We also showed that the prediction of bull fertility is feasible as our cross-validation analyses achieved predictive correlations equal to 0.30 for sire conception rate. Improving reproduction performance is one of the major challenges of the dairy industry worldwide, and for this, it is essential to have accurate predictions of service sire fertility. This study represents the foundation for the development of novel tools that will allow dairy producers, breeders, and artificial insemination companies to make enhanced management and selection decisions on Brown Swiss male fertility.

The protein profile of milk includes several caseins, whey proteins, and nonprotein nitrogen compounds, which influence milk's value for human nutrition and its cheesemaking properties for the dairy industry. To fill in the gap in current knowledge of the patterns of these individual nitrogenous compounds throughout lactation, we tested the ability of a parametric nonlinear lactation model to describe the pattern of each N compound expressed qualitatively (as % of total milk N), quantitatively (in g/L milk), and as daily yield (in g/d). The lactation model was tested on a data set of detailed milk nitrogenous compound profiles (15 fractions-12 protein traits and 3 nonproteins-for each expression mode: 45 traits) obtained from 1,342 cows reared in 41 multibreed herds. Our model was a modified version of Wilmink's model, often used for describing milk yield during lactation because of its reliability and ease of parameter interpretation from a biological point of view. We allowed the sign of the persistency coefficient (parameter c) that explained the variation in the long-term milk component (parameter a) to be positive or negative. We also allowed the short-term milk component (parameter b) to be positive or negative, and we estimated a specific speed of adaptation parameter (parameter k) for each trait rather than assumed a value a priori, as in the original model (k = 0.05). These 4 parameters were included in a nonlinear mixed model with cow breed and parity order as fixed effects, and herd-date as random. Combinations of the positive and negative signs of the b and c parameters allowed us to identify 4 differently shaped lactation curves, all found among the patterns exhibited by the nitrogenous fractions as follows: the "zenith" curve (with a maximum peak; for milk yield and 10 other N traits), the "nadir" curve (with a minimum point; for 20 traits, including almost all those expressed in g/L of milk), the "downward" curve (continuously decreasing; for 14 traits, including almost all those in g/d), and the "upward" curve (continuously increasing; only for κ-casein, in % N). Direct estimation of the k parameters specific to each trait showed the large variability in the adaptation speed of fresh cows and greatly increased the model's flexibility. The results indicated that nonlinear parametric mathematical models can effectively describe the different and complex patterns exhibited by individual nitrogenous fractions during lactation; therefore, they could be useful tools for interpreting milk composition variations during lactation.

Fourier-transform infrared (FTIR) spectroscopy is a powerful high-throughput phenotyping tool for predicting traits that are expensive and difficult to measure in dairy cattle. Calibration equations are often developed using standard methods, such as partial least squares (PLS) regression. Methods that employ penalization, rank-reduction, and variable selection, as well as being able to model the nonlinear relations between phenotype and FTIR, might offer improvements in predictive ability and model robustness. This study aimed to compare the predictive ability of 2 machine learning methods, namely random forest (RF) and gradient boosting machine (GBM), and penalized regression against PLS regression for predicting 3 phenotypes differing in terms of biological meaning and relationships with milk composition (i.e., phenotypes measurable directly and not directly in milk, reflecting different biological processes which can be captured using milk spectra) in Holstein-Friesian cattle under 2 cross-validation scenarios. The data set comprised phenotypic information from 471 Holstein-Friesian cows, and 3 target phenotypes were evaluated: (1) body condition score (BCS), (2) blood β-hydroxybutyrate (BHB, mmol/L), and (3) κ-casein expressed as a percentage of nitrogen (κ-CN, % N). The data set was split considering 2 cross-validation scenarios: samples-out random in which the population was randomly split into 10-folds (8-folds for training and 1-fold for validation and testing); and herd/date-out in which the population was randomly assigned to training (70% herd), validation (10%), and testing (20% herd) based on the herd and date in which the samples were collected. The random grid search was performed using the training subset for the hyperparameter optimization and the validation set was used for the generalization of prediction error. The trained model was then used to assess the final prediction in the testing subset. The grid search for penalized regression evidenced that the elastic net (EN) was the best regularization with increase in predictive ability of 5%. The performance of PLS (standard model) was compared against 2 machine learning techniques and penalized regression using 2 cross-validation scenarios. Machine learning methods showed a greater predictive ability for BCS (0.63 for GBM and 0.61 for RF), BHB (0.80 for GBM and 0.79 for RF), and κ-CN (0.81 for GBM and 0.80 for RF) in samples-out cross-validation. Considering a herd/date-out cross-validation these values were 0.58 (GBM and RF) for BCS, 0.73 (GBM and RF) for BHB, and 0.77 (GBM and RF) for κ-CN. The GBM model tended to outperform other methods in predictive ability around 4%, 1%, and 7% for EN, RF, and PLS, respectively. The prediction accuracies of the GBM and RF models were similar, and differed statistically from the PLS model in samples-out random cross-validation. Although, machine learning techniques outperformed PLS in herd/date-out cross-validation, no significant differences were observed in terms of predictive ability due to the large standard deviation observed for predictions. Overall, GBM achieved the highest accuracy of FTIR-based prediction of the different phenotypic traits across the cross-validation scenarios. These results indicate that GBM is a promising method for obtaining more accurate FTIR-based predictions for different phenotypes in dairy cattle.

The aim of this study was to investigate the associations between differential somatic cell count (DSCC) and milk quality and udder health traits, and for the first time, between DSCC and milk coagulation properties and cheesemaking traits in a population of 1,264 Holstein cows reared in northern Italy. Differential somatic cell count represents the combined proportions of polymorphonuclear neutrophils plus lymphocytes (PMN-LYM) in the total somatic cell count (SCC), with macrophages (MAC) making up the remaining proportion. The milk traits investigated in this study were milk yield (MY), 8 traits related to milk composition and quality (fat, protein, casein, casein index, lactose, urea, pH, and milk conductivity), 9 milk coagulation traits [3 milk coagulation properties (MCP) and 6 curd firming (CF) traits], 7 cheesemaking traits, 3 cheese yield (CY) traits, and 4 milk nutrient recovery in the curd (REC) traits. A linear mixed model was fitted to explore the associations between SCS combined with DSCC and the aforementioned milk traits. An additional model was run, which included DSCC expressed as the PMN-LYM and MAC counts, obtained by multiplying the percentage of PMN-LYM and MAC by SCC in the milk for each cow in the data set. The unfavorable association between SCS and milk quality and technological traits was confirmed. Increased DSCC was instead associated with a linear increase in MY, casein index, and lactose proportion and a linear decrease in milk fat and milk conductivity. Accordingly, DSCC was favorably associated with all MCP and CF traits (with the exception of the time needed to achieve maximum, CF), particularly with rennet coagulation time, and it always displayed linear relationships. Differential somatic cell count was also positively associated with the recovery of milk nutrients in the curd (protein, fat, and energy), which increased linearly with increasing DSCC. The PMN-LYM count was rarely associated with milk traits, even though the pattern observed confirmed the results obtained when both SCS and DSCC were included in the model. The MAC count, however, showed the opposite pattern: MY, casein index, and lactose percentage decreased and milk conductivity increased with an increasing MAC count. No significant association was found between PMN-LYM count and MCP, CF, CY, and REC traits, whereas MAC count was unfavorably associated with MCP, CF traits, some CY traits, and all REC traits. Our results showed that the combined information derived from SCS and DSCC might be useful to monitor milk quality and cheesemaking-related traits.

BackgroundOver the past decade, Fourier transform infrared (FTIR) spectroscopy has been used to predict novel milk protein phenotypes. Genomic data might help predict these phenotypes when integrated with milk FTIR spectra. The objective of this study was to investigate prediction accuracy for milk protein phenotypes when heterogeneous on-farm, genomic, and pedigree data were integrated with the spectra. To this end, we used the records of 966 Italian Brown Swiss cows with milk FTIR spectra, on-farm information, medium-density genetic markers, and pedigree data. True and total whey protein, and five casein, and two whey protein traits were analyzed. Multiple kernel learning constructed from spectral and genomic (pedigree) relationship matrices and multilayer BayesB assigning separate priors for FTIR and markers were benchmarked against a baseline partial least squares (PLS) regression. Seven combinations of covariates were considered, and their predictive abilities were evaluated by repeated random sub-sampling and herd cross-validations (CV).ResultsAddition of the on-farm effects such as herd, days in milk, and parity to spectral data improved predictions as compared to those obtained using the spectra alone. Integrating genomics and/or the top three markers with a large effect further enhanced the predictions. Pedigree data also improved prediction, but to a lesser extent than genomic data. Multiple kernel learning and multilayer BayesB increased predictive performance, whereas PLS did not. Overall, multilayer BayesB provided better predictions than multiple kernel learning, and lower prediction performance was observed in herd CV compared to repeated random sub-sampling CV.ConclusionsIntegration of genomic information with milk FTIR spectral can enhance milk protein trait predictions by 25% and 7% on average for repeated random sub-sampling and herd CV, respectively. Multiple kernel learning and multilayer BayesB outperformed PLS when used to integrate heterogeneous data for phenotypic predictions.

The use of sexed semen to produce purebred replacement heifers allows a large proportion of dairy cows to be mated to beef sires, and quantitative and qualitative improvements to be made to beef production from dairy herds. The major dairy and beef breeds are undergoing rapid genetic improvement as a result of more efficient selection methods, prompting a need to evaluate the meat production of crossbred beef × dairy cattle produced using current genetics. As part of a large project involving 125 commercial dairy farms, we evaluated the combined use of purebreeding with sexed semen and crossbreeding with semen from beef sires, particularly double-muscled breeds. A survey of 1,530 crossbred calves revealed that, whereas purebred dairy calves are destined almost exclusively for veal production, beef × dairy crossbred calves are also destined for beef production after fattening on either the dairy farm of birth or by specialized fatteners. In veal production, compared with Belgian Blue-sired calves (taken as the reference), double-muscled INRA 95-sired calves had a lighter slaughter weight (303 vs. 346 kg), but a greater dressing percent (62.3 vs. 58.4%). Limousin (also known as Limousine)-sired calves had a smaller average daily gain (1.26 vs. 1.34 kg/d), and lighter slaughter (314 vs. 346 kg) and carcass weights (182 vs. 201 kg). Last, Simmental-sired calves had a similar growth rate, but lighter carcass weight (177 vs. 201 kg), smaller dressing percentage (55.3 vs. 58.4%), and smaller muscularity scores (3.25 vs. 3.72). In the case of young bulls and heifers fattened on the dairy farm of birth, Belgian Blue-, Piemontese (also known as Piedmontese)-, and Limousin-sired calves performed similarly; the only exception was that Piemontese-sired calves had a greater dressing percentage. Belgian Blue- and Limousin-sired calves performed similarly when fattened by specialized beef producers. In both veal and beef production, the effects of dam breed were less important than sire breed. Considering the entire project, we can conclude that the combined use of sexed semen for purebreeding and conventional beef semen for terminal crossbreeding improves meat production from dairy herds, especially when the sires are double-muscled beef breeds.

To promote local food products could be a valuable strategy for preserving autochthonous sheep breeds. A total of 115 lambs and ewes of local sheep breeds (Alpagota, Bogna, Foza and Lamon) of the Veneto region were used to evaluate the effect of different feeding systems on the production and quality of lambs and mutton meat. We carried out three fattening trials involving 24 ewes and their 31 suckling lambs, 24 weaned lambs, and 24 heavy lambs, fed indoor diets with or without a rumen-protected conjugated linoleic acids (rp-CLA) supplement (C18:2c9,t11 and C18:2t10,c12) to improve the nutritional value of the meat. Controls were a group of 12 pasture-fed weaned lambs. The weaned lambs at pasture had good fattening performances (33 kg live-weight at slaughter, 45% dressing percentage), whereas lambs fattened indoors on hay and some concentrates had lower weight (28 kg) and dressing percentages (41%), and require improvements to their feed. Feeding the animals on total mixed diets resulted in good in vivo and post mortem performances of culled ewes (64 kg slaughter weight, 48% dressing percentage), suckling lambs (21 kg, 50%) and heavy lambs (61 kg, 48%), showing that valuable supplies of lamb and mutton meat can be obtained across different seasons. Supplementing indoor diets with rp-CLA does not modify the animals’ performances nor the organoleptic traits of meat.HighlightsIntensive indoor fattening allows to obtain suckling, weaned and heavy lambs and mutton from culled ewes all-year-round beyond the seasonal supplies of traditional feeding system.The supplementation of indoor diets with rp-CLA does not modify the animals’ performances, slaughter traits or meat organoleptic traits. Intensive indoor fattening allows to obtain suckling, weaned and heavy lambs and mutton from culled ewes all-year-round beyond the seasonal supplies of traditional feeding system. The supplementation of indoor diets with rp-CLA does not modify the animals’ performances, slaughter traits or meat organoleptic traits.

In this study, we aimed to investigate differences in the genetics of fertility traits (heritability of traits and correlations between traits in divergent environments) in dairy cows of different production levels defined on the basis of the herd-average daily milk energy output (herd-dMEO). Data were obtained from Holstein-Friesian (n = 37,359 for fertility traits, 381,334 for dMEO), Brown Swiss (n = 79,638 for fertility traits, 665,697 for dMEO), and Simmental cows (n = 63,048 for fertility traits, 448,445 for dMEO) reared in northeastern Italy. Fertility traits under study were interval from calving to first service, interval from first service to conception, days open, calving interval, calving rate, and nonreturn rate at d 56. We classified herds into low and high productivity based on the herd-average dMEO (inferred using mixed effects models). We estimated genetic parameters using Bayesian bivariate animal models, where expressions of a phenotype in the low and high dMEO herds were taken as being different-albeit correlated-traits. Fertility traits were more favorable in Simmental than in Holstein-Friesian cows, whereas for all traits, Holstein-Friesian had the highest estimates of intraherd heritability [ranging from 0.021 (0.006-0.038) to 0.126 (0.10-0.15)] and Simmental the lowest [ranging from 0.008 (0.001-0.017) to 0.101 (0.08-0.12)]. The genetic correlations between fertility traits and dMEO were moderate and unfavorable, ranging, in absolute values, from 0.527 (0.37-0.68) to 0.619 (0.50-0.73) in Holstein-Friesian; from 0.339 (0.20-0.47) to 0.556 (0.45-0.66) in Brown Swiss; and from 0.340 (0.10-0.60) to 0.475 (0.33-0.61) in Simmental cattle. The only exception was the nonreturn rate at d 56, which had weak genetic correlations with dMEO in all 3 breeds. The herd correlations between fertility and dMEO tended to be modest and favorable and the residual correlations modest and variable. The heritability of fertility traits tended to be greater in the low dMEO than in the high dMEO herds in the case of the Holstein-Friesians, but not in the case of the Brown Swiss or Simmentals. The additive genetic correlations between fertility traits in the low and high dMEO herds were always lower than 1 [0.329 (-0.17 to 0.85) to 0.934 (0.86 to 0.99)] for all traits considered in all breeds. The correlation was particularly low for the threshold characters and the interval from first service to conception in Holstein-Friesian, suggesting that the relative performances of genotypes vary significantly between herds of different dMEO levels. Although there was large variability in the estimates, results might support making separate genetic evaluations of fertility in the different herd production groups. Our results also indicate that Simmental, a dual-purpose breed, has higher fertility and lower environmental sensitivity than Holstein-Friesian, with Brown Swiss being intermediate.

The yield, flavor, and texture of ripened cheese result from numerous interrelated microbiological, biochemical, and physical reactions that take place during ripening. The aims of the present study were to propose a 2-compartment first-order kinetic model of cheese weight loss over the ripening period; to test the variation in new informative phenotypes describing this process; and to assess the effects on these traits of dairy farming system, individual farms within dairy system, animal factors, and milk composition. A total of 1,211 model cheeses were produced in the laboratory using individual 1.5-L milk samples from Brown Swiss cows reared on 83 farms located in Trento Province. During ripening (60 d; temperature 15°C, relative humidity 85%), the weight of all model cheeses was measured, and cheese yield (cheese weight/processed milk weight, %CY) was calculated at 7 intervals from cheese-making (0, 1, 7, 14, 28, 42, and 60 d). Using these measures, a 2-compartment first-order kinetic model (3-parameter equation) was developed for modeling %CY during the ripening period, as follows: [Formula: see text] , where %CYt is the %CY at ripening time t; %CYi and %CYf are the modeled %CY traits at time 0 d (%CYi = initial %CY) and at the end of a ripening period sufficient to reach a constant wheel weight (%CYf = final %CY after 60 d ripening in the case of small model cheeses); kCY is the instant rate constant for cheese weight loss (%/d). Cheese weight and protein and fat losses were calculated as the % difference between the model cheeses at 0 and after 60 d of ripening. The variation in cheese pH was calculated as the % difference between pH at 0 and after 60 d. Dairy system, individual herd within dairy system, and the cow's parity and lactation stage (tested with a linear mixed model) strongly affected almost all the traits collected during model cheese ripening. Milk fat, protein, lactose, pH, and somatic cell score also greatly affected almost all the traits, although kCY was affected only by milk protein. After including milk composition in the linear mixed model, the importance of all the herd and animal sources of variation was greatly reduced for all traits. The proposed model and novel traits could be tested, first, with the aim of establishing new monitoring procedures enabling the dairy industry to improve milk quality-based payment systems at the herd level and, second, with a view to exploring possible genetic improvements to dairy cow populations.

The use of sexed semen to produce purebred replacement heifers allows a large proportion of dairy cows to be mated to double-muscled sires and for quantitative and qualitative improvements to be made to beef production from dairy herds. A survey of 1,285 crossbred calves showed that they are destined not only for veal production (like purebred dairy calves) but also for beef production. Crossbred calves obtained from Belgian Blue sires (sold on average at 34 d of age and 64 kg of weight) had a higher market value (€363/calf) than those obtained from double-muscled INRA 95 sires (€297/calf, used for veal production) and from Limousin sires (€216/calf for veal and €271/calf for beef production). As a sire breed, Simmental did not differ significantly from Belgian Blue, but as a dam breed, the crossbred calves fetched a higher price (€5.11/kg) than when Holstein was the dam breed (€4.50/kg). Compared with heifer calves, crossbred bull calves at sale were younger (34.1 vs. 37.2 d) and heavier (64.0 vs. 62.6 kg of live weight), fetched a higher price (€5.13 vs. €4.99/kg), and had a greater value (€328 vs. €312/calf). As the value of purebred dairy calves was about €80 to €100/calf, we are able to confirm that the combined use of sexed semen for pure-breeding and conventional beef semen for terminal crossbreeding can increase the income from dairy farms, especially when the sires are double-muscled beef bulls.

Authentication of dairy systems is of growing interest for the dairy industry and we investigated the potentiality of using volatile fingerprinting of ripened cheeses by proton transfer reaction time-of-flight mass spectrometry. A total of 1,075 individual model cheeses made from milk of individual Brown Swiss cows of 72 farms were analysed. Using a linear discriminant analysis, cows and herds were assigned to 3 or 5 dairy systems differing in management, available facilities, and diets. We obtained variable discrimination abilities (up to 77% of correct classification of cheeses and 70% of farms with cross-validation). We found m/z 61,028 (acetic acid), 109,070 (pyrazine), and m/z 137,132 (terpene) characterising model cheeses from traditional dairy systems and m/z 71,086 (3-methyl-butan-1-ol, 3-methyl-3-buten-1-ol, pentan-1-ol), m/z 101,097 (hexan-2-one, hexanal), m/z 123,117 (nonenal), m/z 129,127 (octan-1-one, octanal), and two unidentified peaks m/z 83,071 and m/z 93,090 characterising model cheeses from the modern farms. In conclusion, it seems possible to discriminate between a range of dairy systems using fast volatile fingerprinting of ripened cheeses but a proper validation of results obtained is needed.Highlights Mass spectrometry technique (PTR-ToF-MS) was able to discriminate between dairy systems. We found m/z 61,028 (acetic acid), 109,070 (pyrazine), and m/z 137,132 (terpene) characterising model cheeses from traditional dairy systems. We found m/z 71,086 (3-methyl-butan-1-ol, 3-methyl-3-buten-1-ol, pentan-1-ol), m/z 101,097 (hexan-2-one, hexanal), m/z 123,117 (nonenal), m/z 129,127 (octan-1-one, octanal), and two unidentified peaks m/z 83,071 and m/z 93,090 characterising model cheeses from the modern farms.

Dairy cows have high incidences of metabolic disturbances, which often lead to disease, having a subsequent significant impact on productivity and reproductive performance. As the milk fatty acid (FA) profile represents a fingerprint of the cow’s nutritional and metabolic status, it could be a suitable indicator of metabolic status at the cow level. In this study, we obtained milk FA profile and a set of metabolic indicators (body condition score, ultrasound liver measurements, and 29 hematochemical parameters) from 297 Holstein–Friesian cows. First, we applied a multivariate factor analysis to detect latent structure among the milk FAs. We then explored the associations between these new synthetic variables and the morphometric, ultrasonographic and hematic indicators of immune and metabolic status. Significant associations were exhibited by the odd-chain FAs, which were inversely associated with β-hydroxybutyrate and ceruloplasmin, and positively associated with glucose, albumin, and γ-glutamyl transferase. Short-chain FAs were inversely related to predicted triacylglycerol liver content. Rumen biohydrogenation intermediates were associated with glucose, cholesterol, and albumin. These results offer new insights into the potential use of milk FAs as indicators of variations in energy and nutritional metabolism in early lactating dairy cows.

In general, Fourier-transform infrared (FTIR) predictions are developed using a single-breed population split into a training and a validation set. However, using populations formed of different breeds is an attractive way to design cross-validation scenarios aimed at increasing prediction for difficult-to-measure traits in the dairy industry. This study aimed to evaluate the potential of FTIR prediction using training set combining specialized and dual-purpose dairy breeds to predict different phenotypes divergent in terms of biological meaning, variability, and heritability, such as body condition score (BCS), serum β-hydroxybutyrate (BHB), and kappa casein (k-CN) in the major cattle breed, i.e., Holstein-Friesian. Data were obtained from specialized dairy breeds: Holstein (468 cows) and Brown Swiss (657 cows), and dual-purpose breeds: Simmental (157 cows), Alpine Grey (75 cows), and Rendena (104 cows), giving a total of 1461 cows from 41 multi-breed dairy herds. The FTIR prediction model was developed using a gradient boosting machine (GBM), and predictive ability for the target phenotype in Holstein cows was assessed using different cross-validation (CV) strategies: a within-breed scenario using 10-fold cross-validation, for which the Holstein population was randomly split into 10 folds, one for validation and the remaining nine for training (10-fold_HO), an across-breed scenario (BS_HO) where the Brown Swiss cows were used as the training set and the Holstein cows as the validation set, a specialized multi-breed scenario (BS+HO_10-fold), where the entire Brown Swiss and Holstein populations were combined then split into 10 folds, and a multi-breed scenario (Multi-breed), where the training set comprised specialized (Holstein and Brown Swiss) and dual-purpose (Simmental, Alpine Grey, and Rendena) dairy cows, combined with nine folds of the Holstein cows. Lastly a Multi-breed CV2 scenario was implemented, assuming the same number of records as the reference scenario and using the same proportions as the multi-breed. Within-Holstein, FTIR predictions had a predictive ability of 0.63 for BCS, 0.81 for BHB, and 0.80 for k-CN. Using a specific breed (Brown Swiss) as the training set for prediction in the Holstein population reduced the prediction accuracy by 10% for BCS, 7% for BHB, and 11% for κ-CN. Notably, the combination of Holstein and Brown Swiss cows in the training set increased the predictive ability of the model by 6%, which was 0.66 for BCS, 0.85 for BHB, and 0.87 for k-CN. Using multiple specialized and dual-purpose animals in the training set outperforms the 10-fold_HO (standard) approach, with an increase in predictive ability of 8% for BCS, 7% for BHB, and 10% for k-CN. When the Multi-breed CV2 was implemented, no improvement was observed. Our findings suggest that FTIR prediction of different phenotypes in the Holstein breed can be improved by including different specialized and dual-purpose breeds in the training population. Our study also shows that predictive ability is enhanced when the size of the training population and the phenotypic variability are increased

The aim of this study was to quantify the major sources of variation in the levels of 15 minerals in individual milk samples collected from cows raised in multibreed dairy herds. The herds (n = 27) were classified into 2 categories, according to milk productivity. Milk productivity was based on the net energy of lactating cows' average daily milk yield. Milk samples were collected from 240 cows of 6 different breeds: 3 specialized dairy (Holstein-Friesian, Brown Swiss, and Jersey) and 3 dual-purpose (Simmental, Rendena, and Alpine Grey). The samples were analyzed for macro-elements (Na, Mg, P, S, K, and Ca), essential micro-elements (Mn, Fe, Cu, Zn, and Se), and environmental micro-elements (B, Si, Sr, and Sn), using inductively coupled plasma-optical emission spectrometry. Data were analyzed using a linear mixed model that included fixed effects of days in milk (DIM), parity, breed, and herd productivity, and a random effect of herd-date within productivity level. Results showed that the effect of herd-date varied across minerals. It was especially large for environmental minerals (ranging from 47 to 91% of total variance) and ranged from 11 to 61% for macrominerals and essential microminerals. Milk samples collected from farms with a high level of herd productivity had a richer mineral profile than samples from low-productivity herds. Parity only influenced macrominerals, with the exception of S and Ca, while DIM influenced almost all minerals, with a few exceptions among the environmental elements. Differences in mineral profile were small between specialized and dual-purpose breeds, but they were large within the group of the specialized cows. These breed differences were reduced after adjusting for milk quality and yield, particularly in the case of milk Mg, S, Ca, Mn, and Zn levels. Milk samples from the Jersey and Brown Swiss cows had higher mineral levels (Sn excluded) than milk from the Holstein-Friesian cows; the other breeds of Alpine origin produced milk of intermediate quality. Our findings suggest that breed has a stronger effect on macrominerals and some of the essential microminerals than herd productivity, parity, and DIM. The modification of the mineral profile in milk seems possible for many minerals, but it likely depends on genetics (e.g., breed, selection) and on environmental and management factors in variable proportions according to the mineral considered.

Fourier-transform infrared (FTIR) spectra collected during milk recording schemes at population level can be used for predicting novel traits of interest for farm management, cows' genetic improvement, and milk payment systems. The aims of this study were as follows. (1) To predict cheese yield traits using FTIR spectra from routine milk recordings exploiting previously developed calibration equations. (2) To compare the predicted cheese-making abilities of different dairy and dual-purpose breeds. (3) To analyze the effects of herds' level of intensiveness (HL) and of the cow's level of productivity (CL). (4) To compare the patterns of predicted cheese yields with the patterns of milk composition in different breeds to discern the drivers of cheese-making efficiency. The major sources of variation of FTIR predictions of cheese yield ability (fresh cheese or cheese solids produced per unit milk) of individual milk samples were studied on 115,819 cows of 4 breeds (2 specialized dairy breeds, Holstein and Brown Swiss, and 2 dual-purpose breeds, Simmental and Alpine Grey) from 6,430 herds and exploiting 1,759,706 FTIR test-day spectra collected over 7 yr of milk sampling. Calibration equations used were previously developed on 1,264 individual laboratory model cheese procedures (cross-validation R2 0.85 and 0.95 for fresh and solids cheese yields, respectively). The linear model used for statistical analysis included the effects of parity, lactation stage, year of calving, month of sampling, HL, CL, breed of cow, and the interactions breed × HL and breed × CL. The HL and CL stratifications (5 classes each) were based on average daily secretion of milk net energy per cow. All effects were highly significant. The major conclusions were as follows. (1) The FTIR-based prediction of cheese yield of milk goes beyond the knowledge of fat and protein content, partially explaining differences in cheese-making ability in different cows, breeds and herds. (2) Differences in cheese yields of different breeds are only partially explained by milk fat and protein composition, and less productive breeds are characterized by a higher milk nutrient content as well as a higher recovery of nutrients in the cheese. (3) High-intensive herds not only produce much more milk, but the milk has a higher nutrient content and a higher cheese yield, whereas within herds, compared with less productive cows, the more productive cows have a much greater milk yield, milk with a greater content of fat but not of protein, and a moderate improvement in cheese yield, differing little from expectations based on milk composition. Finally, (4) the effects of HL and CL on milk quality and cheese-making ability are similar but not identical in different breeds, the less productive ones having some advantage in terms of cheese-making ability. We can obtain FTIR-based prediction of cheese yield from individual milk samples retrospectively at population level, which seems to go beyond the simple knowledge of milk composition, incorporating information on nutrient retention ability in cheese, with possible advantages for management of farms, genetic improvement of dairy cows, and milk payment systems.

Background The possibility of assessing meat quality traits over the meat chain is strongly limited, especially in the context of selective breeding which requires a large number of phenotypes. The main objective of this study was to investigate the suitability of portable infrared spectrometers for phenotyping beef cattle aiming to genetically improving the quality of their meat. Meat quality traits (pH, color, water holding capacity, tenderness) were appraised on rib eye muscle samples of 1,327 Piemontese young bulls using traditional (i.e., reference/gold standard) laboratory analyses; the same traits were also predicted from spectra acquired at the abattoir on the intact muscle surface of the same animals 1 d after slaughtering. Genetic parameters were estimated for both laboratory measures of meat quality traits and their spectra-based predictions. Results The prediction performances of the calibration equations, assessed through external validation, were satisfactory for color traits (R2 from 0.52 to 0.80), low for pH and purge losses (R2 around 0.30), and very poor for cooking losses and tenderness (R2 below 0.20). Except for lightness and purge losses, the heritability estimates of most of the predicted traits were lower than those of the measured traits while the genetic correlations between measured and predicted traits were high (average value 0.81). Conclusions Results showed that NIRS predictions of color traits, pH, and purge losses could be used as indicator traits for the indirect genetic selection of the reference quality phenotypes. Results for cooking losses were less effective, while the NIR predictions of tenderness were affected by a relatively high uncertainty of estimate. Overall, genetic selection of some meat quality traits, whose direct phenotyping is difficult, can benefit of the application of infrared spectrometers technology.

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

Several local sheep breeds of the Alpine regions are in danger of extinction as a result of mountain farms being abandoned. Three fattening trials were carried out on 4 breeds of the Veneto region (Foza, Lamon, Alpagota, and Brogna) on a total of 115 ewes and suckling lambs, weaned lambs and heavy lambs slaughtered at 1356 ± 267, 104 ± 24, 219 ± 26, and 342 ± 14 d of age, respectively. The results confirmed that the adoption of total mixed diets resulted in good in vivo and post mortem performances of culled ewes (64 kg slaughter weight, 48% dressing percentage), suckling lambs (21 kg, 50%), weaned lambs (30 kg, 42%) and heavy lambs (61 kg, 48%). Carcass compactness and fatness increased with the age of the lambs, whereas the proportions of lean meat and bone decreased. Meat lightness and cooking losses decreased with age, whereas redness and chroma increased. The meat quality of culled ewes was like that of heavy lambs. Carcass conformation and composition (aside from size traits), and the quality of the meat differed little among the four Alpine breeds. These results show that local breeds are potentially valuable sources of lamb and mutton meat across different seasons, paving the way for strategies to develop and promote the meat with the involvement of high-quality restaurants and meat retailers.Highlights Compared to outdoor grazing, the adoption of indoor diets resulted in good in vivo and post-mortem performances of ewes, suckling lambs and heavy lambs. Autochthonous Alpine breeds are potentially valuable sources of high-quality lamb and mutton across different seasons of the year. This potential could be exploited in developing a prestigious product aimed at high quality restaurants and butchers. The subsequent demand for the product would, in turn, contribute to the survival of the local sheep breeds of Alpine regions.

The aims of this study were to investigate potential functional relationships among milk protein fractions in dairy cattle and to carry out a structural equation model (SEM) GWAS to provide a decomposition of total SNP effects into direct effects and effects mediated by traits that are upstream in a phenotypic network. To achieve these aims, we first fitted a mixed Bayesian multitrait genomic model to infer the genomic correlations among 6 milk nitrogen fractions [4 caseins (CN), namely κ-, β-, αS1-, and αS2-CN, and 2 whey proteins, namely β-lactoglobulin (β-LG) and α-lactalbumin (α-LA)], in a population of 989 Italian Brown Swiss cows. Animals were genotyped with the Illumina BovineSNP50 Bead Chip v.2 (Illumina Inc.). A Bayesian network approach using the max-min hill-climbing (MMHC) algorithm was implemented to model the dependencies or independence among traits. Strong and negative genomic correlations were found between β-CN and αS1-CN (-0.706) and between β-CN and κ-CN (-0.735). The application of the MMHC algorithm revealed that κ-CN and β-CN seemed to directly or indirectly influence all other milk protein fractions. By integrating multitrait model GWAS and SEM-GWAS, we identified a total of 127 significant SNP for κ-CN, 89 SNP for β-CN, 30 SNP for αS1-CN, and 14 SNP for αS2-CN (mostly shared among CN and located on Bos taurus autosome 6) and 15 SNP for β-LG (mostly located on Bos taurus autosome 11), whereas no SNP passed the significance threshold for α-LA. For the significant SNP, we assessed and quantified the contribution of direct and indirect paths to total marker effect. Pathway analyses confirmed that common regulatory mechanisms (e.g., energy metabolism and hormonal and neural signals) are involved in the control of milk protein synthesis and metabolism. The information acquired might be leveraged for setting up optimal management and selection strategies aimed at improving milk quality and technological characteristics in dairy cattle.

Mastitis is one of the most prevalent diseases in dairy cattle and is the cause of considerable economic losses. Alongside somatic cell count (SCC), differential somatic cell count (DSCC) has been recently introduced as a new indicator of intramammary infection. The DSCC is expressed as a count or a proportion (%) of polymorphonuclear neutrophils plus lymphocytes (PMN-LYM) in milk somatic cells. These numbers are complemented to total somatic cell count or to 100 by macrophages (MAC). The aim of this study was to investigate the genetic variation and heritability of DSCC, and its correlation with milk composition, udder health indicators, milk composition, and technological traits in Holstein cattle. Data used in the analysis consisted in single test-day records from 2,488 Holstein cows reared in 36 herds located in northern Italy. Fourier-transform infrared (FTIR) spectroscopy was used to predict missing information for some milk coagulation and cheese-making traits, to increase sample size and improve estimation of the genetic parameters. Bayesian animal models were implemented via Gibbs sampling. Marginal posterior means of the heritability estimates were 0.13 for somatic cell score (SCS); 0.11 for DSCC, MAC proportion, and MAC count; and 0.10 for PMN-LYM count. Posterior means of additive genetic correlations between SCS and milk composition and udder health were low to moderate and unfavorable. All the relevant genetic correlations between the SCC traits considered and the milk traits (composition, coagulation, cheese yield and nutrients recovery) were unfavorable. The SCS showed genetic correlations of -0.30 with the milk protein proportion, -0.56 with the lactose proportion and -0.52 with the casein index. In the case of milk technological traits, SCS showed genetic correlations of 0.38 with curd firming rate (k20), 0.45 with rennet coagulation time estimated using the curd firming over time equation (RCTeq), -0.39 with asymptotic potential curd firmness, -0.26 with maximum curd firmness (CFmax), and of -0.31 with protein recovery in the curd. Differential somatic cell count expressed as proportion was correlated with SCS (0.60) but had only 2 moderate genetic correlations with milk traits: with lactose (-0.32) and CFmax (-0.33). The SCS was highly correlated with the log PMN-LYM count (0.79) and with the log MAC count (0.69). The 2 latter traits were correlated with several milk traits: fat (-0.38 and -0.43 with PMN-LYM and MAC counts, respectively), lactose percentage (-0.40 and -0.46), RCTeq (0.53 and 0.41), tmax (0.38 and 0.48). Log MAC count was correlated with k20 (+0.34), and log PMN-LYM count was correlated with CFmax (-0.26) and weight of water curd as percentage of weight of milk processed (-0.26). The results obtained offer new insights into the relationships between the indicators of udder health and the milk technological traits in Holstein cows.

Our study investigated the inbreeding load for fertility traits in the Italian Brown Swiss dairy cattle breed. Fertility traits included continuous traits (i.e., interval from calving to first service, days open, and calving interval) and categorical traits (i.e., calving rate at first insemination and nonreturn date at d 56). We included only records of the first 3 parities of cows that calved between 2010 and 2018. We traced up the pedigree of the cows with records as far as possible, ending up with a total of 73,246 animals. The final data set consisted of 59,864 records from 34,921 cows. We analyzed all models using a Bayesian approach that included a covariate with total inbreeding in addition to systematic, permanent environment, additive genetic, and inbreeding load effects. We then evaluated the trends in heritabilities and ratios of the inbreeding load using a continuum of partial inbreeding coefficients from 0.001 to 0.100 as reference. Posterior estimates of heritabilities tended to decrease across the continuum, whereas ratios of the inbreeding load tended to increase, more noticeably in categorical traits (calving rate at first insemination and nonreturn date at d 56). From the results obtained, we confirmed the presence of heterogeneity in inbreeding depression. We then predicted the inbreeding load effects, which had a low reliability of prediction, explained by having only 513 ancestors generating inbreeding. However, reliability of prediction was high enough for some of the individuals, obtaining a favorable prediction of inbreeding load for a relevant percentage, which improved the phenotypic performance of their inbred descendants. These results make it feasible to implement breeding and management strategies that select ancestors with a favorable inbreeding load prediction. In addition, it opens the possibility to define a global index for the expected consequences of the inbreeding generated by each individual.

This study aimed to investigate the genetic and putative causal relationships between fertility traits [i.e., days open (DO) and calving rate (CR)] and milk quality, composition, and fatty acid contents (milk composition traits) in Holstein-Friesian, Brown Swiss, and Simmental cattle, using recursive models within a Bayesian framework. Trivariate animal models were run, each including one fertility trait, one milk composition trait, and, in all models, milk yield. The DO and CR data were merged with the test days closest to the insemination date for milk composition traits. After editing, 16,468 to 23,424 records for Holstein-Friesian, 23,424 to 46,660 for Brown Swiss, and 26,105 to 35,574 for Simmental were available for the analyses. Recursive animal models were applied to investigate the possible causal influences of milk composition traits on fertility and the genetic relationships among these traits. The results suggested a potential cause-and-effect relationship between milk composition traits and fertility traits, with the first trait influencing the latter. We also found greater recursive effects of milk composition traits on DO than on CR, the latter with some putative differences among breeds in terms of sensitivity. For instance, the putative causal effects of somatic cell score on CR (on the observed scale, %) varied from -0.96 to -1.39%, depending on the breed. Concerning fatty acids, we found relevant putative effects of C18:0 on CR, with estimates varying from -7.8 to -9.9%. Protein and casein percentages, and short-chain fatty acid showed larger recursive effects on CR, whereas fat, protein, and casein percentages, C16:0, C18:0, and long-chain fatty acid had larger effects on DO. The results obtained suggested that these milk traits could be considered as effective indicators of the effects of animal metabolic and physiological status on fertility.

This study aimed to infer the genetic parameters of five enteric methane emissions (EME) predicted from milk infrared spectra (13 models). The reference values were estimated from milk fatty acid profiles (chromatography), individual model-cheese, and daily milk yield of 1158 Brown Swiss cows (85 farms). Genetic parameters were estimated, under a Bayesian framework, for EME reference traits and their infrared predictions. Heritability of predicted EME traits were similar to EME reference values for methane yield (CH4/DM: 0.232&ndash, 0.317) and methane intensity per kg of corrected milk (CH4/CM: 0.177&ndash, 0.279), smaller per kg cheese solids (CH4/SO: 0.093&ndash, 0.165), but greater per kg fresh cheese (CH4/CU: 0.203&ndash, 0.267) and for methane production (dCH4: 0.195&ndash, 0.232). We found good additive genetic correlations between infrared-predicted methane intensities and the reference values (0.73 to 0.93), less favorable values for CH4/DM (0.45&ndash, 0.60), and very variable for dCH4 according to the prediction method (0.22 to 0.98). Easy-to-measure milk infrared-predicted EME traits, particularly CH4/CM, CH4/CU and dCH4, could be considered in breeding programs aimed at the improvement of milk ecological footprint.

Maedi–Visna virus (MVV) and Mycobacterium avium ssp. paratuberculosis (MAP) are two pathogens that cause chronic, production-limiting diseases in dairy sheep. Although they are present worldwide, there are no detailed reports on their actual effects on milk traits in the literature. This study was designed to investigate the effects of test positivity to MVV and MAP on ovine milk yield, composition and coagulation properties, and curd-firming over time (CFt) variables in clinically healthy animals at the field level. The additive genetic variation and heritabilities of MVV and MAP positivity were also estimated. Milk samples were collected from 1,079 Sarda sheep kept on 23 farms, and pedigree information was obtained from the flock book. Milk yield was also recorded on the sampling date. Positivity for MVV and MAP was determined from milk samples using indirect ELISA test kits. Milk composition traits were measured by spectroscopy, milk coagulation properties were measured with a Formagraph (Foss Italia, Padua, Italy), and CFt traits were calculated using the data from the Formagraph diagram. The effects of MVV and MAP positivity on milk traits were determined through a set of mixed linear models, which took into account various sources of variation, such as days in milk, parity, and flock effects, and included the effects (positive or negative) of the 2 pathogens. A Bayesian threshold sire model with sire relationship was used to estimate genetic variation and heritability. The overall animal prevalence of MVV-positive ewes was 43.6%; on only 1 farm of the 23 tested were all sampled ewes negative. An overall animal prevalence of 10.6% was recorded for MAP, with 4 farms at 0%. Positivity for MVV significantly affected the logarithmic score of the bacterial count, curd firmness after 30 min and 45 min, and the curd-firming instant rate constant. We found significant effects of MAP infection on milk composition, pH, and rennet coagulation time. The mean of the posterior distributions of heritability estimates on the liability scale was 0.15 for MAP and 0.07 for MVV. Our results demonstrate that only a few traits are negatively affected by MVV and MAP positivity, and that there is exploitable genetic variation in MVV and MAP susceptibility in dairy sheep.

Genomic selection (GS) reports on milk fatty acid (FA) profiles have been published quite recently and are still few despite this trait represents the most important aspect of milk nutritional and sensory quality. Reasons for this can be found in the high costs of phenotype recording but also in issues related to its nature of complex trait constituted by multiple genetically correlated variables with low heritabilities. One possible strategy to deal with such constraint is represented by the use of dimension reduction methods. We analysed 40 individual FAs from Italian Brown Swiss, Holstein and Simmental milk through multivariate factor analysis (MFA) to study the genetics of milk FA-related latent variables (factors) and assess their potential use in breeding. A total of nine factors were obtained, and their genetic parameters were inferred under a Bayesian framework using two statistical approaches: the classical pedigree best linear unbiased prediction (ABLUP) and the single-step genomic BLUP (ssGBLUP). The resulting factorial solutions were able to represent groups of FAs with common origin and function and can be considered concise pathway-level phenotypes. The heritability (h2 ) values showed relevant variations across different factors in each breed (0.03 ≤ h2 ≤ 0.38). The accuracies of breeding values predicted were low to high, ranging from 0.13 to 0.72 and from 0.18 to 0.74 considering the pedigree and the genomic model, respectively. The gain in accuracy in genetic prediction due to the addition of genomic information was ~30% and ~5% in validation and training groups respectively, confirming the contribution of genomic information in yielding more accurate predictions compared to the traditional ABLUP methodology. Our results suggest that MFA in combination with GS can be a valuable tool in dairy cattle breeding and deserves to be further investigated for use in future breeding programs to improve cow's milk FA-related traits.

Different fractions of milk nitrogenous compounds (not only caseins) have different effects on the nutritional value of milk, its coagulation and curd firming properties, and its cheese-making efficiency. To assess different sources of variation, especially the cows' breed and genetic variants of the main protein fractions, milk samples were collected from 1,504 cows belonging to 3 dairy breeds (Holstein-Friesian, Brown Swiss, and Jersey) and 3 dual-purpose breeds (Simmental, Rendena, and Alpine Grey) reared in 41 multibreed herds. Beyond crude protein, casein (CN), and urea, 7 protein fractions were analyzed using HPLC, and 5 other N fraction traits were calculated. All 15 traits were measured qualitatively (% of milk N) and quantitatively (g/L of milk). The HPLC technique allowed us to discriminate between the main genetic variants of β-CN, κ-CN, and β-lactoglobulin and thus to genotype the cows for the CSN2, CSN3, and BLG genes, respectively. Data were analyzed using 2 mixed models, both including the effects of herd-date, breed, parity, and lactation stage, and only one also including the effects of the genotypes of the milk proteins. Breed of cow explained 2 to 36% of phenotypic variability for all the N fractions, with the exception of the urea and total casein contents of milk and the urea and β-CN proportions of total milk N. Lactation stage had a considerable influence on the amount (g/L) of almost all the protein fractions in milk, but neither the nonprotein N fractions nor the percentage of milk N protein profile were affected. The inclusion of the CSN2, CSN3, and BLG genotypes in the model explained a large part of the total variability in all the milk protein and nonprotein fractions except urea. It also reduced the variance explained by breed and residual factors. An exception was shown by the proportion of αS1-CN variance explained by breed that moved from 13 to 28%. Similarly, for amount (g/L) of β-CN, the effect of breed became significant (12%), whereas it was almost null before inclusion of genotypes. In terms of percentage of milk N, the genotypes of CSN3 notably affected all the casein fractions, whereas the BLG genotypes had a much greater influence on most noncasein traits. The genotypes of the CSN2 gene exerted an appreciable effect on αS2-CN and not β-CN, as expected. Comparing the 2 models, we were also able to discriminate the effect of the breed on a milk N fraction, both quantitatively and qualitatively, in 2 quotas: the first due to the milk protein polymorphisms (major genes) and the second due to other genetic factors (polygene), after correcting for the effect of herd-date of sampling, parity, and lactation stage. The knowledge about the detailed milk protein profile of different cattle breeds provided by this study could be of great benefit for the dairy industry, providing new tools for the enhancement of milk payment systems and breeding program designs.

Depending on whether milk protein fractions are evaluated qualitatively or quantitatively, different genetic outcomes may emerge. In this study, we compared the genetic parameters for the major milk protein fractions&mdash, caseins (&alpha, S1-, &alpha, S2-, &beta, and к-CN), and whey proteins (&beta, lactoglobulin, &beta, LG, &alpha, lactalbumin, &alpha, LA)&mdash, estimated using the multi-trait genomic best linear unbiased prediction method and expressed variously as milk content (g/100g milk), percentage of milk nitrogen (%N) and daily yield per cow (g/d). The results showed that the genetic parameter estimates varied according to how the milk protein fractions were expressed. Heritability estimates for the caseins and whey protein fractions expressed as daily yields were lower than when they were expressed as proportions and contents, revealing important differences in genetic outcomes. The proportion and the content of &beta, CN were negatively correlated with the proportions and contents of &alpha, S1-CN, &alpha, S2-CN, and к-CN, while the daily yield of &beta, &ndash, CN was negatively correlated with the daily yields of &alpha, S1-CN and &alpha, S2-CN. The Spearman&rsquo, s rank correlations and the coincidence rates between the various predicted genomic breeding values (GEBV) for the milk protein fractions expressed in different ways indicated that these differences had a significant effect on the ranking of the animals. The results suggest that the way milk protein fractions are expressed has implications for breeding programs aimed at improving milk nutritional and technological characteristics.

The aim of this study was to assess the influence of the amounts of the αS1-, αS2-, β-, and κ-casein (CN) and the α-lactalbumin and β-lactoglobulin protein fractions on the efficiency of the cheese-making process independently of their genetic polymorphisms. The study was carried out on milk samples from 1,271 Brown Swiss cows from 85 herds classified into 4 categories according to management, feeding, and housing characteristics (traditional and modern systems). To assess the efficiency of the cheese-making process, we processed the milk samples according to a laboratory cheese-making procedure (1,500 mL/sample) and obtained the following measures: (1) 3 percentage cheese yields (%CYcurd, %CYsolids, %CYwater), (2) 2 daily cheese yields obtained by multiplying %CY (curd and total solids) by daily milk yields (dCYcurd, dCYsolids), (3) 4 measures of nutrient recovery in the curd (RECfat, RECprotein, RECsolids, RECenergy), and (4) 2 measures of cheese-making efficiency in terms of the ratio between the observed and theoretical %CY (Ef-%CYcurd, Ef-%CYsolids). All the aforementioned traits were analyzed by fitting 2 linear mixed models with protein fractions as fixed effects expressed as percentage in the milk (model M-%milk) and as percentage of the total casein content (model M-%cas) together with the effects of total casein content (only in model M-%cas), daily milk yield (only in model M-%milk; not for dCY traits), dairy system, herd (random effect), days in milk, parity, and vat. The efficiency of overall cheese yield (Ef-%CYcurd) was mostly positively associated with β-CN content in the milk, whereas Ef-%CYsolids was greater with higher amounts of κ-CN and αS1-CN (M-%milk) due to the strong influence of both fractions on the recovery rate of milk components in the curd (fat and total solids, protein with αS1-CN only) when expressed as percentage of milk and of total casein; only β-CN was more important for RECprotein. In contrast, we found β-lactoglobulin to be highly negatively related to all the traits related to the cheese-making process and to the daily cheese yield per cow, whereas α-lactalbumin was positively associated with the latter traits. Additional research on this topic is needed, with particular focus on the genetic and genomic aspects of the role of protein fractions in the cheese-making process and on the associations between genetic polymorphisms in milk protein and milk nutrient recovery in the curd.

The present study investigated the effect of somatic cell count, lactose, and pH on sheep milk composition, coagulation properties (MCP), and curd firming (CF) parameters. Individual milk samples were collected from 1,114 Sarda ewes reared in 23 farms. Milk composition, somatic cell count, single point MCP (rennet coagulation time, RCT; curd firming time, k 20 ; and curd firmness, a 30 , a 45 , and a 60 ), and CF model parameters were achieved. Phenotypic traits were statistically analyzed using a mixed model to estimate the effects of the different levels of milk somatic cell score (SCS), lactose, and pH, respectively. Additive genetic, herd, and residual correlations among these 3 traits, and with milk composition, MCP and CF parameters, were inferred using a Bayesian approach. From a phenotypic point of view, higher SCS levels caused a delayed gelification of milk. Lactose concentration and pH were significant for many milk quality traits, with a very intense effect on both coagulation times and curd firming. These traits (RCT, RCT estimated using the curd firming over time equation, and k 20 ) showed an unfavorable increase of about 20% from the highest to the lowest level of lactose. Milk samples with pH values lower than 6.56 versus higher than 6.78 were characterized by an increase of RCT (from 6.00 to 14.3 min) and k 20 (from 1.65 to 2.65 min) and a decrease of all the 3 curd firmness traits. From a genetic point of view, the marginal posterior distribution of heritability estimates evidenced a large and exploitable variability for all 3 phenotypes. The mean intra-farm heritability estimates were 0.173 for SCS, 0.418 for lactose content, and 0.206 for pH. Lactose (favorably), and SCS and pH (unfavorably), at phenotypic and genetic levels, were correlated mainly with RCT and RCT estimated using the curd firming over time equation and scarcely with the other curd firming traits. The SCS, lactose, and pH were significantly correlated with each other's. In conclusion, results reported in the present study suggest that SCS, pH, and lactose affect, contemporarily and independently, milk quality and MCP. These phenotypes, easily available during milk recording schemes measured by infrared spectra prediction, could be used as potential indicators traits for improving cheese-making ability of ovine milk.

Data on Holstein (16,890), Brown Swiss (31,441), Simmental (25,845), and Alpine Grey (12,535) cows reared in northeastern Italy were used to assess the ability of milk components (fat, protein, casein, and lactose) and Fourier transform infrared (FTIR) spectral data to diagnose pregnancy. Pregnancy status was defined as whether a pregnancy was confirmed by a subsequent calving and no other subsequent inseminations within 90 d of the breeding of specific interest. Milk samples were analyzed for components and FTIR full-spectrum data using a MilkoScan FT+ 6000 (Foss Electric, Hillerod, Denmark). The spectrum covered 1,060 wavenumbers (wn) from 5,010 to 925 cm −1 . Pregnancy status was predicted using generalized linear models with fat, protein, lactose, casein, and individual FTIR spectral bands or wavelengths as predictors. We also fitted a generalized linear model as a simultaneous function of all wavelengths (1,060 wn) with a Bayesian variable selection model using the BGLR R-package (https://r-forge.r-project.org/projects/bglr/). Prediction accuracy was determined using the area under a receiver operating characteristic curve based on a 10-fold cross-validation (CV-AUC) assessment based on sensitivities and specificities of phenotypic predictions. Overall, the best prediction accuracies were obtained for the model that included the complete FTIR spectral data. We observed similar patterns across breeds with small differences in prediction accuracy. The highest CV-AUC value was obtained for Alpine Grey cows (CV-AUC = 0.645), whereas Brown Swiss and Simmental cows had similar performance (CV-AUC = 0.630 and 0.628, respectively), followed by Holsteins (CV-AUC = 0.607). For single-wavelength analyses, important peaks were detected at wn 2,973 to 2,872 cm −1 where Fat-B (C-H stretch) is usually filtered, wn 1,773 cm −1 where Fat-A (C=O stretch) is filtered, wn 1,546 cm −1 where protein is filtered, wn 1,468 cm −1 associated with urea and fat, wn 1,399 and 1,245 cm −1 associated with acetone, and wn 1,025 to 1,013 cm −1 where lactose is filtered. In conclusion, this research provides new insight into alternative strategies for pregnancy screening of dairy cows.

Ruminants (and milk production) contribute to global climate change through enteric methane emissions (EME), and any attempt to reduce them is complicated by the fact that they are difficult and expensive to measure directly. In the case of dairy cows, a promising indirect method of estimating EME is to use the milk fatty acid profile as a proxy, as a relationship exists between microbial activity in the rumen and the molecules available for milk synthesis in the mammary gland. In the present study, we analyzed the detailed fatty acid profiles (through gas chromatography) of a large number of milk samples from 1,158 Brown Swiss cows reared on 85 farms with the aim of testing in the field 2 equations for estimating EME taken from a published meta-analysis. The average estimated methane yield (CH4 emission per kg of dry matter intake, 21.34 ± 1.60 g/kg) and methane intensity (per kg of corrected milk, 14.17 ± 1.78 g/kg), and the derived methane production (CH4 emissions per day per cow, 357 ± 109 g/d) were similar to those previously published. Using data from model cheese makings from individual cows, we also calculated estimated methane intensity per kilogram of fresh cheese (99.7 ± 16.4 g/kg) and cheese solids (207.5 ± 30.9 g/kg). Dairy system affected all EME estimates. Traditional dairy farms, and modern farms including corn silage in the TMR exhibited greater estimated methane intensities. We found very wide variability in estimated EME traits among different farms within dairy system (0.33 to 0.61 of total variance), suggesting the need to modify the farms' feeding regimens and management practices to mitigate emissions. Among the individual factors, parity order affected all estimated EME traits excepted methane yield, with an increase from first lactation to the following ones. Lactation stage exhibited more favorable estimated EME traits during early lactation, concomitant with the availability of nutrients from body tissue mobilization for mammary synthesis of milk. Our results showed a coherence between the EME traits estimated from the analysis of milk fatty acids and the expectations according to current knowledge. Further research is needed to validate the results obtained in this study in other breeds and populations, to assess the magnitude of the genetic variation and the potential of these phenotypes to be exploited in breeding programs with the aim to mitigate emissions.

The aim of this study was to analyze milk protein composition in purebred and crossbred dairy cattle and estimate the effects of individual sources of variation on the investigated traits. Milk samples were collected from 505 cows from three commercial farms located in Northern Italy, some of which had originated from crossbreeding programs, although most were purebred Holsteins (HO). The basic crossbreeding scheme was a three-breed rotational system using Swedish Red (SR) semen on HO cows (SR×HO), Montbeliarde (MO) semen on SR×HO cows (MO×(SR×HO)) and HO semen again on MO×(SR×HO) cows. A smaller number of purebred HO from each of the herds were mated inverting the breed order (MO×HO and SR×(MO×HO)) or using Brown Swiss (BS) bulls (BS×HO) then MO bulls (MO×(BS×HO)). Milk samples were analyzed by reverse-phase HPLC to obtain protein fraction amounts (g/l) and proportions (% of total true protein). Traits were analyzed using a linear model, which included the fixed effects of herd-test-day (HTD), parity, days in milk and breed combination. Results showed that milk protein fractions were influenced by HTD, stage of lactation, parity and breed combination. The increase in protein concentration during lactation was due in particular to β-casein (β-CN), α S1-CN and β-lactoglobulin (β-LG). The higher protein content of primiparous milk was mainly due to higher concentrations of all casein fractions. The milk from crossbred cows had higher contents and proportions of κ-CN and α-lactalbumin (α-LA), lower proportions of β-LG and greater proportion of caseins/smaller in whey proteins on milk true protein than purebred HO. The three-way crossbreds differed from two-way crossbreds only in having greater proportions of α-LA in their milk. Of the three-way crossbreds, the SR sired cows yielded milk with a smaller content and proportion of β-LG than the MO sired cows, and, consequently, a higher proportion of caseins than whey proteins. Results from this study support the feasibility of using crossbreeding programs to alter milk protein profiles with the aim of improving milk quality and cheese-making properties.

The aim of this study was to apply Bayesian models to the Fourier-transform infrared spectroscopy spectra of individual sheep milk samples to derive calibration equations to predict traditional and modeled milk coagulation properties (MCP), and to assess the repeatability of MCP measures and their predictions. Data consisted of 1,002 individual milk samples collected from Sarda ewes reared in 22 farms in the region of Sardinia (Italy) for which MCP and modeled curd-firming parameters were available. Two milk samples were taken from 87 ewes and analyzed with the aim of estimating repeatability, whereas a single sample was taken from the other 915 ewes. Therefore, a total of 1,089 analyses were performed. For each sample, 2 spectra in the infrared region 5,011 to 925 cm −1 were available and averaged before data analysis. BayesB models were used to calibrate equations for each of the traits. Prediction accuracy was estimated for each trait and model using 20 replicates of a training-testing validation procedure. The repeatability of MCP measures and their predictions were also compared. The correlations between measured and predicted traits, in the external validation, were always higher than 0.5 (0.88 for rennet coagulation time). We confirmed that the most important element for finding the prediction accuracy is the repeatability of the gold standard analyses used for building calibration equations. Repeatability measures of the predicted traits were generally high (≥95%), even for those traits with moderate analytical repeatability. Our results show that Bayesian models applied to Fourier-transform infrared spectra are powerful tools for cheap and rapid prediction of important traits in ovine milk and, compared with other methods, could help in the interpretation of results.

This study aimed at investigating the effects of growth rate (GR) on live performance, carcass and dressed green ham (DRHAM) traits of finishing heavy pigs managed in accordance with the regulations governing the production of Protected Denomination of Origin (PDO) dry-cured hams. Data were obtained from 233 crossbred pigs restrictively fed from 90 to 165 kg BW. Individual BW, daily feed intake, and weights of carcass and main cuts were recorded. Thighs were measured for subcutaneous fat thickness, analysed for fatty acid composition of subcutaneous fat, and dry-cured hams were weighed. Records were classified according to GR quartile, calculated using the residuals of a statistical model that included the effects of pen, sex and diet. Data were analysed using the same model with the addition of the effect of GR quartile. Average GR ranged from 578 in the first quartile to 748 g/d in the fourth. Final BW and weights of carcass, DRHAM and dry-cured ham increased as the GR quartile increased (p

Using sexed semen to produce purebred replacement heifers makes it possible to mate a large proportion of dairy cows to double-muscled sires and to quantitatively and qualitatively improve beef production and increase the income from dairy herds. Net profit first depends on changes in the farm's overall fertility rate. The objective of this study was to analyze the conception rate in herds using a combination of conventional dairy semen (for pure- and crossbreeding), X-sorted dairy semen (to produce purebred replacement heifers), and conventional beef semen (for terminal crossbreeding). Data were obtained from 50,785 inseminations of 15,580 dairy cows (78% Holstein-Friesian, 15% Brown Swiss, 2% Simmental, and 5% crossbreds) from 106 dairy farms (average milk yield 35.1 ± 9.4 kg/d, with 3.76 ± 0.83% fat and 3.32 ± 0.39% protein contents). To account for the main potential confounders, we used separate generalized linear mixed-effects models for cows and virgin heifers. The results showed that the odds ratio of conception improved (1.00 to 1.34) with an increase in the average milk yield of the herd but worsened (1.12 to 0.70) with an increase in the milk yield of individual cows within herd. The summer months showed a strong reduction in the odds ratio of conception in cows (0.56 in July and August) but not in virgin heifers. Multiparous cows had a lower odds ratio of conception (0.85) than primiparous cows (1.00). The order of insemination did not affect the fertility of the cows or heifers, whereas the odds ratio of conception improved with advancing lactation (1.00 to 2.12). The Simmental cows were more fertile than Holstein-Friesians (1.37 vs. 1.00), whereas the fertility of the heifers was not affected by breed. Taking all these possible confounders into account simultaneously, in pure-breeding the odds ratio of conception using sexed semen did not differ from that using conventional dairy semen in cows (0.90 vs. 1.00) or in virgin heifers (0.95 vs. 1.00). However, crossbreeding using conventional beef and dairy semen improved the odds ratio of conception (1.10 and 1.17, respectively) in cows (1.37 using beef semen) and heifers (1.25 using dairy semen). The proportion of newborn heifer calves was ≥90% using sexed dairy semen. The combined use of sexed semen, especially on heifers, to produce purebred replacement females and beef semen to produce terminal crossbred calves was shown to have the potential to increase overall herd fertility, which could be further improved using sexed dairy semen to produce dairy crossbreds instead of purebred replacement heifers.

A key concern in beef production is how to improve carcass and meat quality traits. Identifying the genomic regions and biological pathways that contribute to explaining variability in these traits is of great importance for selection purposes. In this study, genome wide-association (GWAS) and pathway-based analyses of carcass traits (age at slaughter (AS), carcass weight (CW), carcass daily gain (CDG), conformation score and rib-eye muscle area) and meat quality traits (pH, Warner-Bratzler shear force, purge loss, cooking loss and colour parameters (lightness, redness, yellowness, chroma, hue)) were conducted using genotype data from the ‘GeneSeek Genomic Profiler Bovine LD’ array in a cohort of 1166 double-muscled Piemontese beef cattle. The genome wide-association analysis was based on the GRAMMAR-GC approach and identified 37 significant single nucleotide polymorphisms (SNPs), which were associated with 12 traits (P

Fourier-transform near- and mid-infrared (FTIR) milk spectral data are routinely collected in many countries worldwide. Establishing an optimal strategy to use spectral data in genetic evaluations requires knowledge of the heritabilities of individual FTIR wavelength absorbances. Previous FTIR heritability estimates have been based on relatively small sample sizes and have not considered the possibility that heritability may vary across parities and stages of the lactation. We used data from ∼370,000 test-day records of Canadian Holstein cows to produce a landscape of the heritability of FTIR spectra, 1,060 wavelengths in the near- and mid-infrared spectrum (5,011-925 cm-1), by parity and month of the lactation (mo 1 to 3 and mo 1 to 6, respectively). The 2 regions of the spectrum associated with absorption of electromagnetic energy by water molecules were estimated to have very high phenotypic variances, very low heritabilities, and very low proportion of variance explained by herd-year-season (HYS) subclasses. The near- or short-wavelength infrared (SWIR: 5,066-3,672 cm-1) region was also characterized by low heritability estimates, whereas the estimated proportion of the variance explained by HYS was high. The mid-wavelength infrared region (MWIR: 3,000-2,500 cm-1) and the transition between mid and long-wavelength infrared region (MWIR-LWIR: 1,500-925 cm-1) harbor several waves characterized by moderately high (≥0.4) heritabilities. Most of the high-heritability regions contained wavelengths that are reported to be associated with important milk metabolites and components. Interestingly, these 2 same regions tended to show more variability in heritabilities between parity and lactation stage. Second parity showed heritability patterns that were distinctly different from those of the first and third parities, whereas the first 2 mo of the lactation had clearly distinct heritability patterns compared with mo 3 to 6.

In the present study we estimated the genetic parameters of enteric methane emissions (EME) traits predicted from milk fatty acid profile (FA) and those of their predictors in 1,091 Brown Swiss cows reared on 85 farms in order to assess the potential of using EME-related phenotypes in selective breeding. Univariate and bivariate genetic models were fitted in a Bayesian framework. The means of the marginal posterior distribution of intra-herd heritability ranged from 0.12 for estimated methane production (g/d/cow) to 0.24 for estimated methane yield (g/kg dry matter intake [DMI]), with intermediate values for estimated methane intensity, increasingly higher when expressed per kg of corrected milk (0.13), fresh cheese (0.16), or cheese solids (0.20). Regarding the correlations, the milk quality traits and percentage cheese yields were generally moderately correlated with the estimated EME traits, and were variable in terms of sign. Daily milk and cheese yield traits were, as expected, all highly positively correlated with estimated daily methane production. In contrast, they were negatively correlated with estimated methane yield and intensity, the estimates being large in the case of phenotypic and herd correlations, and low in the case of additive genetic and residual correlations. With the exception of the negative correlations with daily methane production, EME traits exhibited trivial correlations with body size and BCS of cows, which, in turn, were negatively correlated with milk yield. Although the results should be validated on a larger population and different breeds, our study demonstrate the presence of additive genetic variation of EME traits, which could be exploited in breeding programmes for the improvement in both milk production and the ecological footprint of dairy farming.HighlightsEnteric methane emissions (EME) of dairy cows can be estimated on the basis of milk fatty acid profile.EME exhibited exploitable genetic variation.Genetic selection could be preferentially based on predicted methane intensity per kg of milk, or per kg of cheese in countries where milk production is used mainly for cheese-making. Enteric methane emissions (EME) of dairy cows can be estimated on the basis of milk fatty acid profile. EME exhibited exploitable genetic variation. Genetic selection could be preferentially based on predicted methane intensity per kg of milk, or per kg of cheese in countries where milk production is used mainly for cheese-making.

The aim of this study was to perform genetic, genome-wide association (GWAS), and gene-set enrichment analyses with latent variables related to milk fatty acid profile (i.e., fatty acids factor scores; FAF), milk composition, and udder health in a cohort of 1,158 Italian Brown Swiss cows. The phenotypes under study were 12 FAF previously identified through factor analysis and classified as follows: de novo FA (F1), branched-chain FA-milk yield (F2), biohydrogenation (F3), long-chain fatty acids (F4), desaturation (F5), short-chain fatty acids (F6), milk protein and fat contents (F7), odd fatty acids (F8), conjugated linoleic acids (F9), linoleic acid (F10), udder health (F11) and vaccelenic acid (F12). (Co)variance components were estimated for factor scores using a Bayesian linear animal model via Gibbs sampling. The animals were genotyped with the Illumina BovineSNP50 BeadChip v.2 (Illumina Inc., San Diego, CA). A single marker regression model was fitted for GWAS analysis. The gene-set enrichment analysis was run on the GWAS results using the Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes pathway databases to identify the ontologies and pathways associated with the FAF. Marginal posterior means of the heritabilities of the aforementioned FAF ranged from 0.048 for F12 to 0.310 for F5. Factors F1 and F6 had the highest number of relevant genetic correlations with the other traits. The genomic analysis detected a total of 39 significant SNP located on 17 Bos taurus autosomes. All latent variables produced signals except for F2 and F10. The traits with the highest number of significant associations were F11 (17) and F12 (7). Gene-set enrichment analyses identified significant pathways (false discovery rate 5%) for F3 and F7. In particular, systemic lupus erythematosus was enriched for F3, whereas the MAPK (mitogen-activated protein kinase) signaling pathway was overrepresented for F7. The results support the existence of important and exploitable genetic and genomic variation in these latent explanatory phenotypes. Information acquired might be exploited in selection programs and when designing further studies on the role of the putative candidate genes identified in the regulation of milk composition and udder health.

We investigated the influence of rumen-protected conjugated linoleic acid (rpCLA) on growth performances, and carcass and meat quality traits in beef. Twenty-four young bulls and 30 heifers obtained from double-muscled beef sires and dairy cows were fed a low-protein ration (110 g/kg DM of crude protein) supplemented with 0, 8 or 80 g/d of a commercial rpCLA product. The animals were monthly weighed and scored for body muscularity and fatness. Blood samples were collected after 140 days on feed. Animals were slaughtered when they reached average in vivo fatness scores of around 2.5 (heifers) and 2.0 (bulls) points respectively. At slaughter, carcasses, various organs and parts of the gastrointestinal tract were weighed; the 5th rib was dissected and its tissue and muscle chemical composition was determined. The rpCLA had little influence on growth performance but decreased the blood urea content by 28% (p < 0.01). The rpCLA × sex interactions for daily gain (p < 0.05), conformation scores (p < 0.01), and blood creatinine content (p < 0.05) suggest that males were more responsive to rpCLA than females when fed a low-protein ration, probably because of the metabolic protein-sparing effect of CLA. Only slight differences were observed in carcass weight and quality at slaughter. The results indicate that the response of beef cattle to rpCLA is dependent on sex or on their propensity for lean and fat accretion. It is also possible that counteracting feedback mechanisms compensate for the influence of rpCLA administration over the course of growth.

Breast cancer is the most common cancer in women worldwide. Cancer metastases are responsible for the high mortality rate. A small but distinct subset of cells, cancer stem cells (CSCs), have the capacity to self-renew, initiate tumour formation, and develop metastases. The CSC content in human breast cancer correlates with the Hippo tumour suppressor signalling pathway. Specifically, the activity of YAP/TAZ, transcription co-activators of the Hippo pathway, sustains the self-renewal and tumour-initiation capacities of CSCs. Little is known about YAP/TAZ in canine and feline mammary tumours, which are very common tumours. The preliminary aim of the study was to investigate the expression of YAP/TAZ in canine and feline mammary tumours by Western blot and immunohistochemistry. Increased cytoplasmic and nuclear expression of YAP/TAZ was observed in all carcinomas compared to normal tissues, indicating neoplastic deregulation of the Hippo pathway. Nuclear expression significantly increased in grade III (high grade carcinomas) compared to grade I (low grade carcinomas) tumours, suggesting that YAP/TAZ play a role in the increased aggressiveness of these tumours. Moreover, different scoring systems for immunohistochemical analyses were compared and the H index and the Allred scores were the most significant. In conclusion, YAP/TAZ are expressed in aggressive canine and feline mammary tumours as reported in some human cancers. Further studies might better elucidate the role of the Hippo pathway in prognosis and as a target for new therapies. In addition, tumours in dogs and cats may be a useful model to study this pathway.

There is an ever-growing interest in research oriented towards the improvement of quality of animal products. In this context, one major operational bottleneck is the possibility to collect quality indicators over the meat and dairy chains and for selective breeding purposes. The use of near-infrared (NIR) and the Fourier-transformed infrared (FTIR) spectroscopy techniques have been proven to be powerful precision phenotyping tools for high-throughput meat and milk quality assessment. Such technologies allow scoring large number of animals and/or derived-products for novel (predicted) phenotypes and indicator traits to set-up potential new payment systems and boost the genetic improvement. One important step in the use of NIR and FTIR tools is the definition of the “gold standard” as the infrared-based predictions could act only as indicators traits. Indeed, the definition of a robust calibration set, the assessment of repeatability and reproducibility of the reference (i.e., gold standard) as well as the detection of random and systematic errors are crucial steps. Once the reference phenotype has been defined, different statistical methodologies could be applied to infrared spectra data. For instance, the partial least squares regression (PLS) is a multivariate regression method commonly used to build up prediction models using NIR and FTIR spectra data. However, the implementation of advanced statistical approaches, such as Bayesian approaches and machine learning methods, might allow us to achieve more robust and accurate predictions. In this talk, we will describe and discuss some of the challenges and potentials of NIR and FTIR tools for large-scale precision phenotyping. Some examples include the use of NIR and Visible-NIR (Vis-NIR) for assessing meat quality parameters (also using portable instruments able to collect spectra directly from the muscle surface at the slaughterhouse) and the use of FTIR for predicting several traits related to fine milk composition and technological traits in dairy cattle.

The aim of this study was to investigate in Holstein cows the genetic basis of blood serum metabolites [i.e., total protein, albumin, globulin, albumin:globulin ratio (A:G), and blood β-hydroxybutyrate (BHB)], a set of milk phenotypes related to udder health, milk quality technological characteristics, and genetic relationships among them. Samples of milk were collected from 498 Holstein cows belonging to 28 herds. All animal welfare and milk phenotypes were assessed using standard analytical methodology. A set of Bayesian univariate and bivariate animal models was implemented via Gibbs sampling, and statistical inference was based on the marginal posterior distributions of parameters of concern. We observed a small additive genetic influence for serum albumin concentrations, moderate heritability (≥0.20) for total proteins, globulins, and A:G, and high heritability (0.37) for blood BHB. Udder health traits (somatic cell score, milk lactose, and milk pH) showed low or moderate heritabilities (0.15–0.20), whereas variations in milk protein fraction concentrations were confirmed as mostly under genetic control (heritability: 0.21–0.71). The moderate and high heritabilities observed for milk coagulation properties and curd firming modeling parameters provided confirmation that genetic background exerts a strong influence on the cheese-making ability of milk, largely due to genetic polymorphisms in the major milk protein genes. Blood BHB showed strong negative genetic correlations with globulins (−0.619) but positive correlations with serum albumin (0.629) and A:G (0.717), which suggests that alterations in the serum protein pattern and BHB blood levels are likely to be genetically related. Strong relationships were found between albumin and fat percentages (−0.894), between globulin and αS2-CN (−0.610), and, to a lesser extent, between serum protein pattern and milk technological characteristics. Genetic relationships between blood BHB and traits related to udder health and milk quality and technological characteristics were mostly weak. This study provides evidence that there is exploitable additive genetic variation for traits related to animal health and welfare and throws light on the shared genetic basis of these traits and the phenotypes related to the quality and cheese-making ability of milk.

The relationship of the estrous cycle to milk composition and milk physical properties was assessed on Holstein (n = 10,696), Brown Swiss (n = 20,501), Simmental (n = 17,837), and Alpine Grey (n = 8,595) cows reared in northeastern Italy. The first insemination after calving for each cow was chosen to be the day of estrus and insemination. Test days surrounding the insemination date (from 10 d before to 10 d after the day of the estrus) were selected and categorized in phases relative to estrus as diestrus high-progesterone, proestrus, estrus, metestrus, and diestrus increasing-progesterone phases. Milk components and physical properties were predicted on the basis of Fourier-transform infrared spectra of milk samples and were analyzed using a linear mixed model, which included the random effects of herd, the fixed classification effects of year-month, parity number, breed, estrous cycle phase, day nested within the estrous cycle phase, conception, partial regressions on linear and quadratic effects of days in milk nested within parity number, as well as the interactions between conception outcome with estrous cycle phase and breed with estrous cycle phase. Milk composition, particularly fat, protein, and lactose, showed clear differences among the estrous cycle phases. Fat increased by 0.14% from diestrus high-progesterone to estrous phase, whereas protein concomitantly decreased by 0.03%. Lactose appeared to remain relatively constant over diestrus high-progesterone, rising 1 d before the day of estrus followed by a gradual reduction over the subsequent phases. Specific fatty acids were also affected across the estrous cycle phases: C14:0 and C16:0 decreased (-0.34 and -0.48%) from proestrus to estrus with a concomitant increase in C18:0 and C18:1 cis-9 (0.40 and 0.73%). More general categories of fatty acids showed a similar behavior; that is, unsaturated fatty acids, monounsaturated fatty acids, polyunsaturated fatty acids, trans fatty acids, and long-chain fatty acids increased, whereas the saturated fatty acids, medium-chain fatty acids, and short-chain fatty acids decreased during the estrous phase. Finally, urea, somatic cell score, freezing point, pH, and homogenization index were also affected indicating variation associated with the hormonal and behavioral changes of cows in standing estrus. Hence, the variation in milk profiles of cows showing estrus should potentially be taken into account for precision dairy farming management.