Search

Abstract Background Various blood metabolites are known to be useful indicators of health status in dairy cattle, but their routine assessment is time-consuming, expensive, and stressful for the cows at the herd level. Thus, we evaluated the effectiveness of combining in-line near infrared (NIR) milk spectra with on-farm (days in milk [DIM] and parity) and genetic markers for predicting blood metabolites in Holstein cattle. Data were obtained from 388 Holstein cows from a farm with an AfiLab system. NIR spectra, on-farm information, and single nucleotide polymorphisms (SNP) markers were blended to develop calibration equations for blood metabolites using the elastic net (ENet) approach, considering 3 models: (1) Model 1 (M1) including only NIR information, (2) Model 2 (M2) with both NIR and on-farm information, and (3) Model 3 (M3) combining NIR, on-farm and genomic information. Dimension reduction was considered for M3 by preselecting SNP markers from genome-wide association study (GWAS) results. Results Results indicate that M2 improved the predictive ability by an average of 19% for energy-related metabolites (glucose, cholesterol, NEFA, BHB, urea, and creatinine), 20% for liver function/hepatic damage, 7% for inflammation/innate immunity, 24% for oxidative stress metabolites, and 23% for minerals compared to M1. Meanwhile, M3 further enhanced the predictive ability by 34% for energy-related metabolites, 32% for liver function/hepatic damage, 22% for inflammation/innate immunity, 42.1% for oxidative stress metabolites, and 41% for minerals, compared to M1. We found improved predictive ability of M3 using selected SNP markers from GWAS results using a threshold of > 2.0 by 5% for energy-related metabolites, 9% for liver function/hepatic damage, 8% for inflammation/innate immunity, 22% for oxidative stress metabolites, and 9% for minerals. Slight reductions were observed for phosphorus (2%), ferric-reducing antioxidant power (1%), and glucose (3%). Furthermore, it was found that prediction accuracies are influenced by using more restrictive thresholds (−log10(P-value) > 2.5 and 3.0), with a lower increase in the predictive ability. Conclusion Our results highlighted the potential of combining several sources of information, such as genetic markers, on-farm information, and in-line NIR infrared data improves the predictive ability of blood metabolites in dairy cattle, representing an effective strategy for large-scale in-line health monitoring in commercial herds.

Abstract Background Metabolic disturbances adversely impact productive and reproductive performance of dairy cattle due to changes in endocrine status and immune function, which increase the risk of disease. This may occur in the post-partum phase, but also throughout lactation, with sub-clinical symptoms. Recently, increased attention has been directed towards improved health and resilience in dairy cattle, and genomic selection (GS) could be a helpful tool for selecting animals that are more resilient to metabolic disturbances throughout lactation. Hence, we evaluated the genomic prediction of serum biomarkers levels for metabolic distress in 1353 Holsteins genotyped with the 100K single nucleotide polymorphism (SNP) chip assay. The GS was evaluated using parametric models best linear unbiased prediction (GBLUP), Bayesian B (BayesB), elastic net (ENET), and nonparametric models, gradient boosting machine (GBM) and stacking ensemble (Stack), which combines ENET and GBM approaches. Results The results show that the Stack approach outperformed other methods with a relative difference (RD), calculated as an increment in prediction accuracy, of approximately 18.0% compared to GBLUP, 12.6% compared to BayesB, 8.7% compared to ENET, and 4.4% compared to GBM. The highest RD in prediction accuracy between other models with respect to GBLUP was observed for haptoglobin (hapto) from 17.7% for BayesB to 41.2% for Stack; for Zn from 9.8% (BayesB) to 29.3% (Stack); for ceruloplasmin (CuCp) from 9.3% (BayesB) to 27.9% (Stack); for ferric reducing antioxidant power (FRAP) from 8.0% (BayesB) to 40.0% (Stack); and for total protein (PROTt) from 5.7% (BayesB) to 22.9% (Stack). Using a subset of top SNPs (1.5k) selected from the GBM approach improved the accuracy for GBLUP from 1.8 to 76.5%. However, for the other models reductions in prediction accuracy of 4.8% for ENET (average of 10 traits), 5.9% for GBM (average of 21 traits), and 6.6% for Stack (average of 16 traits) were observed. Conclusions Our results indicate that the Stack approach was more accurate in predicting metabolic disturbances than GBLUP, BayesB, ENET, and GBM and seemed to be competitive for predicting complex phenotypes with various degrees of mode of inheritance, i.e. additive and non-additive effects. Selecting markers based on GBM improved accuracy of GBLUP.

ABSTRACT: Udder health has a crucial role in sustainable milk production, and various reports have pointed out that changes in udder condition seem to affect milk mineral content. The somatic cell count (SCC) is the most recognized indicator for the determination of udder health status. Recently, a new parameter, the differential somatic cell count (DSCC), has been proposed for a more detailed evaluation of intramammary infection patterns. Specifically, the DSCC is the combined proportions of polymorphonuclear neutrophils and lymphocytes (PMN-LYM) on the total SCC, with macrophages (MAC) representing the remainder proportion. In this study, we evaluated the association between DSCC in combination with SCC on a detailed milk mineral profile in 1,013 Holstein-Friesian cows reared in 5 herds. An inductively coupled plasma–optical emission spectrometry was used to quantify 32 milk mineral elements. Two different linear mixed models were fitted to explore the associations between the milk mineral elements and first, the DSCC combined with SCC, and second, DSCC expressed as the PMN-LYM and MAC counts, obtained by multiplying the proportion of PMN-LYM and MAC by SCC. We observed a significant positive association between SCC and milk Na, S, and Fe levels. Differential somatic cell count showed an opposite behavior to the one displayed by SCC, with a negative association with Na and positive association with K milk concentrations. When considering DSCC as count, Na and K showed contrasting behavior when associated with PMN-LYM or MAC counts, with decreasing of Na content and increasing K when associated with increasing PMN-LYM counts, and increasing Na and decreasing K when associated with increasing MAC count. These findings confirmed that an increase in SCC is associated with altered milk Na and K amounts. Moreover, MAC count seemed to mirror SCC patterns, with the worsening of inflammation. Differently, PMN-LYM count exhibited patterns of associations with milk Na and K contents attributable more to LYM than PMN, given the nonpathological condition of the majority of the investigated population. An interesting association was observed for milk S content, which increased with increasing of inflammatory conditions (i.e., increased SCC and MAC count) probably attributable to its relationship with milk proteins, especially whey proteins. Moreover, milk Fe content showed positive associations with the PMN-LYM population, highlighting its role in immune regulation during inflammation. Further studies including individuals with clinical condition are needed to achieve a comprehensive view of milk mineral behavior during udder health impairment.

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.HIGHLIGHTS 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.

The use of sexed semen to produce replacement heifers in dairy farms allows many cows to be inseminated with bulls of beef breeds with the aim of improving meat production. The major dairy and beef breeds are undergoing rapid genetic improvement as a result of more efficient selection methods, necessitating evaluation of the quality of meat produced by beef × dairy crosses obtained using current genetics. As part of a larger project, we carried out quantitative and qualitative evaluations of the veal and beef produced from crossbred calves produced by dairy farms that use sexed semen to obtain purebred replacement heifers. A total of 231 bull-calves were reared in the dairy farms or in specialised fattening farms, slaughtered, and sampled for meat quality analysis. These comprised: 104 crossbred calves, mainly from Belgian Blue sires and Holstein dams, which were intensively fattened for veal production; 84 beef × dairy crossbred calves, destined for beef production on the dairy farm of origin (12 calves) or in a specialised intensive fattening centre (72 calves); 26 purebred Charolais and 17 beef crosses from suckler cows, included for comparison, which were fattened in the same fattening centre and fed the same diet as the beef × dairy crosses. Veal and beef meat differed in all quality traits; the beef produced by dairy farmers was better than the beef produced by specialised fatteners in terms of cooking losses and tenderness, even though these animals were much older and weighed less at slaughter. There were minor differences in the quality traits of veal between the various beef × dairy breed combinations. In the case of intensive beef production, minor differences were observed between Belgian Blue-sired crosses from Brown Swiss cows (lipid content and shear force of meat) and Belgian Blue-sired crosses from crossbred cows (haem iron, lightness, yellowness, and hue) in comparison with crosses from Holstein cows. The beef of young bulls from suckler cows (Charolais and crosses) had a lower pH, dry matter and lipid contents, slightly higher cholesterol and haem iron contents, and higher cooking losses than meat from beef × dairy crosses. Taking into account the results from the entire project, we conclude that the combined use of sexed semen for pure-breeding and conventional beef semen for terminal crossbreeding improves meat production from dairy herds, both quantitatively and qualitatively.HIGHLIGHTS The meat from beef × dairy crossbred calves was found to have valuable qualitative characteristics. With regard to veal production, the meat quality characteristics of crossbred calves did not differ appreciably between breed combinations. With regard to intensive beef production, the meat from beef × dairy young bulls had greater quality than that from suckler cows. Crosses fattened by dairy farmers grew slower than those by specialised fatteners, but their beef was better in terms of cooking losses and tenderness.

The milk Fourier-transform infra-red spectrometry (FTIRS) can be used to predict the enteric methane emissions (EME) at population level. In this study, the variability in FTIRS predicted EME traits due to the breed of cow, farm, the production level of individual cows within herds, and their interactions were evaluated. A dataset obtained from milk recordings, which covered four breeds (Holstein, Brown Swiss, Simmental and Alpine Grey), 6,430 herds, 115,819 cows, and 1,759,706 test-day milk/spectra records was used. The herds were stratified into 5 production levels considering their average daily milk energy production; individual cows within herd were similarly stratified considering their individual production levels. The EME traits were predicted directly from milk FTIR spectra or indirectly from six informative milk FA predicted from milk spectra. The statistical model included, separately for each trait and method, breed, herd intensiveness level, cow production level, and their interactions, year, month, parity, and lactation stage. The direct and indirect methods yielded similar results in predicting CH4 yield per kg of DMI and CH4 intensity per kg of fat-protein corrected milk. The indirect method was reliable in predicting daily EME production per cow, the indirect one did not. EME was affected by the breed x herd production intensiveness interaction, and to a lesser degree, by the breed x cow production level interaction. A better understanding of the complex interactions influencing EME in dairy herds was achieved. This would be useful for the genetic improvement, the environmental certification of farms, and for setting prices in milk payment schemes.HIGHLIGHTS Causes of variation of the FTIR predicted EME were studied at population level using data from 6,430 herds and 1,759,706 milk spectra records Reliable daily EME estimates were achieved from 6 six informative fatty acids predicted from milk spectra EME was influenced by the interaction breed × level of herd intensiveness and, to a lesser extent, by the individual production level

ABSTRACT: The causes of variation in the milk mineral profile of dairy cattle during the first phase of lactation were studied under the hypothesis that the milk mineral profile partially reflects the animals' metabolic status. Correlations between the minerals and the main milk constituents (i.e., protein, fat, and lactose percentages), and their associations with the cows' metabolic status indicators were explored. The metabolic status indicators (MET) that we used were blood energy-protein metabolites [nonesterified fatty acids, β-hydroxybutyrate (BHB), glucose, cholesterol, creatinine, and urea], and liver ultrasound measurements (predicted triacylglycerol liver content, portal vein area, portal vein diameter and liver depth). Milk and blood samples, and ultrasound measurements were taken from 295 Holstein cows belonging to 2 herds and in the first 120 d in milk (DIM). Milk mineral contents were determined by ICP-OES; these were considered the response variable and analyzed through a mixed model which included DIM, parity, milk yield, and MET as fixed effects, and the herd/date as a random effect. The MET traits were divided in tertiles. The results showed that milk protein was positively associated with body condition score (BCS) and glucose, and negatively associated with BHB blood content; milk fat was positively associated with BHB content; milk lactose was positively associated with BCS; and Ca, P, K and S were the minerals with the greatest number of associations with the cows' energy indicators, particularly BCS, predicted triacylglycerol liver content, glucose, BHB and urea. We conclude that the protein, fat, lactose, and mineral contents of milk partially reflect the metabolic adaptation of cows during lactation and within 120 DIM. Variations in the milk mineral profile were consistent with changes in the major milk constituents and the metabolic status of cows.

Abstract Background Subclinical intramammary infection (IMI) represents a significant problem in maintaining dairy cows’ health. Disease severity and extent depend on the interaction between the causative agent, environment, and host. To investigate the molecular mechanisms behind the host immune response, we used RNA-Seq for the milk somatic cells (SC) transcriptome profiling in healthy cows (n = 9), and cows naturally affected by subclinical IMI from Prototheca spp. (n = 11) and Streptococcus agalactiae (S. agalactiae; n = 11). Data Integration Analysis for Biomarker discovery using Latent Components (DIABLO) was used to integrate transcriptomic data and host phenotypic traits related to milk composition, SC composition, and udder health to identify hub variables for subclinical IMI detection. Results A total of 1,682 and 2,427 differentially expressed genes (DEGs) were identified when comparing Prototheca spp. and S. agalactiae to healthy animals, respectively. Pathogen-specific pathway analyses evidenced that Prototheca’s infection upregulated antigen processing and lymphocyte proliferation pathways while S. agalactiae induced a reduction of energy-related pathways like the tricarboxylic acid cycle, and carbohydrate and lipid metabolism. The integrative analysis of commonly shared DEGs between the two pathogens (n = 681) referred to the core-mastitis response genes, and phenotypic data evidenced a strong covariation between those genes and the flow cytometry immune cells (r 2 = 0.72), followed by the udder health (r 2 = 0.64) and milk quality parameters (r 2 = 0.64). Variables with r ≥ 0.90 were used to build a network in which the top 20 hub variables were identified with the Cytoscape cytohubba plug-in. The genes in common between DIABLO and cytohubba (n = 10) were submitted to a ROC analysis which showed they had excellent predictive performances in terms of discriminating healthy and mastitis-affected animals (sensitivity > 0.89, specificity > 0.81, accuracy > 0.87, and precision > 0.69). Among these genes, CIITA could play a key role in regulating the animals’ response to subclinical IMI. Conclusions Despite some differences in the enriched pathways, the two mastitis-causing pathogens seemed to induce a shared host immune-transcriptomic response. The hub variables identified with the integrative approach might be included in screening and diagnostic tools for subclinical IMI detection.

Abstract 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.

The transformation of milk into cheese largely depends on the technological properties of the raw material, with breed being a crucial factor that influences both the composition and coagulation properties of the milk used for cheesemaking. This study uses canonical correlation analysis to explore the relationships between physicochemical traits and coagulation properties in milk from various Spanish breeds, aiming to identify both common and breed-specific patterns that impact milk technological aptitude. A total of 832 milk samples from Manchega, Assaf, Merino de Grazalema, and Merino de Los Pedroches breeds were analyzed. The milk characteristics investigated included pH, composition (fat, protein, lactose, total solids), and coagulation properties (curd firmness—A60, rennet clotting time—RCT, curd firming time—k20, and individual laboratory curd yield—ILCY). The results reveal a shared correlation structure across breeds and unique covariation patterns in some breeds that deviate from the general trend. While Assaf and Merino de Los Pedroches follow the common correlation pattern, Manchega and Merino de Grazalema exhibit distinct patterns. This research underscores the need for in-depth study and suggests that the dairy industry could benefit from shifting from the traditional focus on maximizing fat and protein for higher curd yields to considering technological traits for selective breeding.

Sheep milk from local breeds is important for the production of high-quality cheeses throughout the Mediterranean region, such as Manchego cheese in Spain. To maintain sustainable and efficient production, it is necessary to reach a better understanding of how the composition and hygiene of the milk affect the coagulation process, with the aim of optimizing production yield. This study implemented a stochastic production frontier function to estimate the potential production of curd and efficiency using data from the four seasons of a study of 77 Manchega sheep farms. The Cobb–Douglas production frontier model was estimated using the maximum likelihood estimation method. The results showed that the content of protein, lactose, and fat exhibited increasing returns to scale, with protein content being the most significant factor for curd production. Approximately half of the inefficiency was due to factors related to the technological properties and the hygiene of the milk. The pH, curd firmness, and concentration of lactic acid bacteria improved the efficiency of coagulation, while the concentration of spores of lactate-fermenting Clostridium spp., Pseudomonas spp., staphylococci, and catalase-negative gram-positive cocci favored the inefficiency of the coagulation process. To date, this is the first study to evaluate the effect of different factors, such as microbial groups, milk composition, and technological properties, on the efficiency of the coagulation process in dairy sheep.

This study conducted a seasonal analysis of bulk tank milk from 77 sheep farms to establish relationships between the concentration of major microbial groups and milk coagulation properties. The investigated milk traits included composition (pH, fat, casein, lactose), coagulation properties (curd firmness: A60-, rennet clotting time: RCT-, curd firming time: k20-, curd yield: CY-), and somatic cell score (SCS). The main microbial groups analyzed were total mesophilic bacteria (SPC), thermodurics (THERMO), psychrotrophs (PSYCHRO), Pseudomonas spp. (PSEUDO), lactic acid bacteria (LAB), catalase-negative gram-positive cocci (GPCNC), Escherichia coli (ECOLI), coliforms other than Escherichia coli (COLI), coagulase-positive staphylococci (CPS), coagulase-negative staphylococci (CNS), and spores of lactate-fermenting Clostridium (BAB). Mixed linear models were used to explore associations between coagulation properties and the aforementioned variables. Results demonstrated that incorporating microbial loads into the models improves their fit and the relative quality of the outcomes. An important seasonality is demonstrated by an increase in CY and A60, along with a decrease in RCT and k20 during autumn and winter, contrasting with spring and summer. BAB concentration resulted in a reduction of A60 and an increase in RCT, whereas SPC concentration led to an enhancement of A60 and a reduction in RCT. An increase in GPCNC concentration was associated with an increase in k20 and a decrease in CY.

Intensive selection for improved productivity has been accompanied by an increase in inbreeding rates and a reduction in genetic diversity. The increase in inbreeding tends to impact performance, especially fitness-related traits such as male fertility. Inbreeding can be monitored using runs of homozygosity (ROH), defined as contiguous lengths of homozygous genotypes observed in an individual’s chromosome. The goal of this study was to evaluate the presence of ROH in Italian Brown Swiss cattle and assess its association with bull fertility. First, we evaluated the association between ROH and male fertility using 1,102 Italian Brown Swiss bulls with sire conception rate records and 572 K SNPs spanning the entire genome. Second, we split the entire population into 100 high-fertility and 100 low-fertility bulls to investigate the potential enrichment of ROH segments in the low-fertility group. Finally, we mapped the significant ROH regions to the bovine genome to identify candidate genes associated with sperm biology and male fertility. Notably, there was a negative association between bull fertility and the amount of homozygosity. Four different ROH regions located in chromosomes 6, 10, 11, and 24 were significantly overrepresented in low-fertility bulls (Fisher’s exact test, p-value

ABSTRACT: 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 (αS1-CN, αS2-CN, β-CN, κ-CN, β-LG, α-LA) expressed qualitatively as percentages of total nitrogen content (% N), quantitatively as contents in milk (g/L), and as daily production levels (g/d). Individual milk samples were collected from 1,264 Brown Swiss cows reared in 85 commercial herds in Trento Province (northeast Italy). A total of 989 cows were successfully genotyped using the Illumina Bovine SNP50 v.2 BeadChip (Illumina Inc.), and a genomic relationship matrix was constructed using the 37,519 SNP markers obtained. Milk protein fractions were quantified and the β-CN, κ-CN, and β-LG genetic variants were identified by reversed-phase HPLC (RP-HPLC). All protein fractions were analyzed through a Bayesian multitrait animal model implemented via Gibbs sampling. The effects of days in milk, parity order, and the CSN2, CSN3, and BLG genotypes were assigned flat priors in this model, whereas the effects of herd and animal additive genetic were assigned Gaussian prior distributions, and inverse Wishart distributions were assumed for the respective co-variance matrices. Marginal posterior distributions of the parameters of interest were compared before and after the inclusion of the effects of the 3 major genes in the model. The results showed that a high portion of the genetic variance was controlled by the major genes. This was particularly apparent in the qualitative protein profile, which was found to have a higher heritability than the protein fraction contents in milk and their daily yields. When the genes were included individually in the model, CSN2 was the major gene controlling all the casein fractions except for κ-CN, which was controlled directly by the CSN3 gene. The BLG gene had the most influence on the 2 whey proteins. The genetic correlations showed the major genes had only a small effect on the relationships between the protein fractions, but through comparison of the correlation coefficients of the proteins expressed in different ways they revealed potential mechanisms of regulation and competitive synthesis in the mammary gland. The estimates for the effects of the CSN2 and CSN3 genes on protein profiles showed overexpression of protein synthesis in the presence of the B allele in the genotype. Conversely, the β-LG B variant was associated with a lower concentration of β-LG compared with the β-LG A variant, independently of how the protein fractions were expressed, and it was followed by downregulation (or upregulation in the case of the β-LG B) of all other protein fractions. These results should be borne in mind when seeking to design more efficient selection programs aimed at improving milk quality for the efficiency of dairy industry and the effect of dairy products on human health.

Abstract Improving reproductive performance remains a major goal in dairy cattle worldwide. Service sire has been recognized as an important factor affecting herd fertility. The main objective of this study was to reveal the genetic basis of male fertility in Italian Brown Swiss dairy cattle. Dataset included 1102 Italian Brown Swiss bulls with sire conception rate records genotyped with 454k single nucleotide polymorphisms. The analysis included whole-genome scans and gene-set analyses to identify genomic regions, individual genes and genetic mechanisms affecting Brown Swiss bull fertility. One genomic region on BTA1 showed significant additive effects. This region harbors gene RABL3 which is implicated cell proliferation and motility. Two genomic regions, located on BTA6 and BTA26, showed marked non-additive effects. These regions harbor genes, such as WDR19 and ADGRA1, that are directly involved in male fertility, including sperm motility, acrosome reaction, and embryonic development. The gene-set analysis revealed functional terms related to cell adhesion, cellular signaling, cellular transport, immune system, and embryonic development. Remarkably, a gene-set analysis also including Holstein and Jersey data, revealed significant processes that are common to the three dairy breeds, including cell migration, cell–cell interaction, GTPase activity, and the immune function. Overall, this comprehensive study contributes to a better understanding of the genetic basis of male fertility in cattle. In addition, our findings may guide the development of novel genomic strategies for improving service sire fertility in Brown Swiss cattle.

ABSTRACT: Cheese-making traits in dairy cattle are important to the dairy industry but are difficult to measure at the individual level because there are limitations on collecting phenotypic information. Mid-infrared spectroscopy has its advantages, but it can only be used during monthly milk recordings. Recently, in-line devices for real-time analysis of milk quality have been developed. The AfiLab recording system (Afimilk) offers significant benefits as phenotypes can be collected from each cow at each milking session. The objective of this study was to assess the potential of integrating AfiLab real-time milk analyzer measures with the stacking ensemble learning technique using heterogeneous base learners for the in-line daily monitoring of cheese-making traits in Holstein cattle with a view to developing a precision livestock farming system for monitoring the technological quality of milk. Data and samples for wet-laboratory analyses were collected from 499 Holstein cows belonging to 2 farms where the AfiLab system was installed. The traits of concern were 9 milk coagulation traits [3 milk coagulation properties (MCP), and 6 curd firming traits (CFt)], and 7 cheese-making traits [3 cheese yield (CY) traits, and 4 milk nutrient recovery in the curd (REC) traits]. The near-infrared AfiLab spectral data and on-farm information (days in milk and parity) were used to assess the predictive ability of different statistical methods [elastic net (EN), gradient boosting machine (GBM), extreme gradient boosting (XGBoost), and artificial neural network (ANN)] across different cross-validation scenarios. These statistical methods were considered the base learners, which were then combined in a stacking ensemble learning. Results indicate that including information on the cows (days in milk and parity) in the AfiLab infrared prediction increased its accuracy by 10.3% for traditional MCP, 13.8% for curd firming, 9.8% for CY, and 11.2% for REC traits compared with those obtained from near-infrared AfiLab alone. The statistical approaches exhibited high prediction accuracies (R2) averaged across the cross-validation scenarios for traditional MCP (0.58 for ANN, 0.55 for EN and GBM, 0.52 for XGBoost, and 0.62 for stacking ensemble), CFt (0.55 for ANN, 0.54 for EN and GBM, 0.53 for XGBoost, and 0.61 for stacking ensemble), and similar R2 averages for CY and REC (0.55 for ANN, 0.54 for EN and GBM, 0.53 for XGBoost, and 0.61 for stacking ensemble). The ANN approach was more accurate than the other base learners (EN, GBM, and XGBoost) and improved accuracy across cross-validation scenarios on average by 7% for traditional MCP, 5% for CFt, 8% for CY, and 7% for REC. The stacking ensemble method improved prediction accuracy by 3% to 31% for traditional MCP, 2% to 26% for CFt, 1% to 38% for CY traits, and 2% to 27% for REC traits compared with the base learners. The prediction accuracies of the different approaches evaluated tended to decrease from the 10-fold cross-validation to the independent validation scenario, although there was a smaller reduction in prediction accuracy with the stacking ensemble learning technique across all the cross-validation scenarios. Our results show that combining in-line on-farm information with stacking ensemble machine learning represents an effective alternative for obtaining robust daily predictions of milk cheese-making traits.

Abstract Precision livestock farming technologies are used to monitor animal health and welfare parameters continuously and in real time in order to optimize nutrition and productivity and to detect health issues at an early stage. The possibility of predicting blood metabolites from milk samples obtained during routine milking by means of infrared spectroscopy has become increasingly attractive. We developed, for the first time, prediction equations for a set of blood metabolites using diverse machine learning methods and milk near-infrared spectra collected by the AfiLab instrument. Our dataset was obtained from 385 Holstein Friesian dairy cows. Stacking ensemble and multi-layer feedforward artificial neural network outperformed the other machine learning methods tested, with a reduction in the root mean square error of between 3 and 6% in most blood parameters. We obtained moderate correlations (r) between the observed and predicted phenotypes for γ-glutamyl transferase (r = 0.58), alkaline phosphatase (0.54), haptoglobin (0.66), globulins (0.61), total reactive oxygen metabolites (0.60) and thiol groups (0.57). The AfiLab instrument has strong potential but may not yet be ready to predict the metabolic stress of dairy cows in practice. Further research is needed to find out methods that allow an improvement in accuracy of prediction equations.

ABSTRACT: Substantial research has been carried out on rapid, nondestructive, and inexpensive techniques for predicting cheese composition using spectroscopy in the visible and near-infrared radiation range. Moreover, in recent years, new portable and handheld spectrometers have been used to predict chemical composition from spectra captured directly on the cheese surface in dairies, storage facilities, and food plants, removing the need to collect, transport, and process cheese samples. For this review, we selected 71 papers (mainly dealing with prediction of the chemical composition of cheese) and summarized their results, focusing our attention on the major sources of variation in prediction accuracy related to cheese variability, spectrometer and spectra characteristics, and chemometrics techniques. The average coefficient of determination obtained from the validation samples ranged from 86 to 90% for predicting the moisture, fat, and protein contents of cheese, but was lower for predicting NaCl content and cheese pH (79 and 56%, respectively). There was wide variability with respect to all traits in the results of the various studies (standard deviation: 9–30%). This review draws attention to the need for more robust equations for predicting cheese composition in different situations; the calibration data set should consist of representative cheese samples to avoid bias due to an overly specific field of application and ensure the results are not biased for a particular category of cheese. Different spectrometers have different accuracies, which do not seem to depend on the spectrum extension. Furthermore, specific areas of the spectrum—the visible, infrared-A, or infrared-B range—may yield similar results to broad-range spectra; this is because several signals related to cheese composition are distributed along the spectrum. Small, portable instruments have been shown to be viable alternatives to large bench-top instruments. Last, chemometrics (spectra pre-treatment and prediction models) play an important role, especially with regard to difficult-to-predict traits. A proper, fully independent, validation strategy is essential to avoid overoptimistic results.

Fusarium mycotoxins are inactivated by rumen flora; however, a certain amount can pass the rumen and reticulum or be converted into biological active metabolites. Limited scientific evidence is available on the impact and mitigation of Fusarium mycotoxins on dairy cows’ performance and health, particularly when cows are exposed for an extended period (more than 2 months). The available information related to these mycotoxin effects on milk cheese-making parameters is also very poor. The objective of this study was to evaluate a commercially available mycotoxin mitigation product (MMP, i.e., TOXO® HP-R, Selko, Tilburg, The Netherlands) in lactating dairy cows fed a Fusarium mycotoxin-contaminated diet, and the repercussions on the dry matter intake, milk yield, milk quality, cheese-making traits and health status of cows. The MMP contains smectite clays, yeast cell walls and antioxidants. In the study, 36 lactating Holstein cows were grouped based on the number of days of producing milk, milk yield, body condition score and those randomly assigned to specific treatments. The study ran over 2 periods (March/May–May/July 2022). In each period, six animals/treatment were considered. The experimental periods consisted of 9 days of adaptation and 54 days of exposure. The physical activity, rumination time, daily milk production and milk quality were measured. The cows were fed once daily with the same total mixed ration (TMR) composition. The experimental groups consisted of a control (CTR) diet, with a TMR with low contamination, high moisture corn (HMC), and beet pulp; a mycotoxins (MTX) diet, with a TMR with highly contaminated HMC, and beet pulp; and an MTX diet supplemented with 100 g/cow/day of the mycotoxin mitigation product (MMP). The trial has shown that the use of MMP reduced the mycotoxin’s negative effects on the milk yield and quality (protein, casein and lactose). The MTX diet had a lower milk yield and feed efficiency than the CTR and MMP HP-R diets. The MMP limited the negative effect of mycotoxin contamination on clotting parameters, mitigating the variations on some coagulation properties; however, the MMP inclusion tended to decrease the protein and apparent starch digestibility of the diet. These results provide a better understanding of mycotoxin risk on dairy cows’ performances and milk quality. The inclusion of an MMP product mitigated some negative effects of the Fusarium mycotoxin contamination during this trial. The major effects were on the milk yield and quality in both the experimental periods. These results provide better insight on the effects of mycotoxins on the performance and quality of milk, as well as the cheese-making traits. Further analyses should be carried out to evaluate MMP’s outcome on immune–metabolic responses and diet digestibility.

The composition of raw milk is of major importance for dairy products, especially fat, protein, and casein (CN) contents, which are used worldwide in breeding programs for dairy species because of their role in human nutrition and in determining cheese yield (%CY). The aim of the study was to develop formulas based on detailed milk composition to disentangle the role of each milk component on %CY traits. To this end, 1,271 individual milk samples (1.5 L/cow) from Brown Swiss cows were processed according to a laboratory model cheese-making procedure. Fresh %CY (%CYCURD), total solids and water retained in the fresh cheese (%CYSOLIDS and %CYWATER), and 60-days ripened cheese (%CYRIPENED) were the reference traits and were used as response variables. Training-testing linear regression modeling was performed: 80% of observations were randomly assigned to the training set, 20% to the validation set, and the procedure was repeated 10 times. Four groups of predictive equations were identified, in which different combinations of predictors were tested separately to predict %CY traits: (i) basic composition, i.e., fat, protein, and CN, tested individually and in combination; (ii) udder health indicators (UHI), i.e., fat + protein or CN + lactose and/or somatic cell score (SCS); (iii) detailed protein profile, i.e., fat + protein fractions [CN fractions, whey proteins, and nonprotein nitrogen (NPN) compounds]; (iv) detailed protein profile + UHI, i.e., fat + protein fractions + NPN compounds and/or UHI. Aside from the positive effect of fat, protein, and total casein on %CY, our results allowed us to disentangle the role of each casein fraction and whey protein, confirming the central role of β-CN and κ-CN, but also showing α-lactalbumin (α-LA) to have a favorable effect, and β-lactoglobulin (β-LG) a negative effect. Replacing protein or casein with individual milk protein and NPN fractions in the statistical models appreciably increased the validation accuracy of the equations. The cheese industry would benefit from an improvement, through genetic selection, of traits related to cheese yield and this study offers new insights into the quantification of the influence of milk components in composite selection indices with the aim of directly enhancing cheese production.

The infrared spectrum of bovine milk is used to predict many interesting traits, whereas there have been few studies on goat milk in this regard. The objective of this study was to characterize the major sources of variation in the absorbance of the infrared spectrum in caprine milk samples. A total of 657 goats belonging to 6 breeds and reared on 20 farms under traditional and modern dairy systems were milk-sampled once. Fourier-transform infrared (FTIR) spectra were taken (2 replicates per sample, 1314 spectra), and each spectrum contained absorbance values at 1060 different wavenumbers (5000 to 930 × cm−1), which were treated as a response variable and analyzed one at a time (i.e., 1060 runs). A mixed model, including the random effects of sample/goat, breed, flock, parity, stage of lactation, and the residual, was used. The pattern and variability of the FTIR spectrum of caprine milk was similar to those of bovine milk. The major sources of variation in the entire spectrum were as follows: sample/goat (33% of the total variance); flock (21%); breed (15%); lactation stage (11%); parity (9%); and the residual unexplained variation (10%). The entire spectrum was segmented into five relatively homogeneous regions. Two of them exhibited very large variations, especially the residual variation. These regions are known to be affected by the absorbance of water, although they also exhibited wide variations in the other sources of variation. The average repeatability of these two regions were 45% and 75%, whereas for the other three regions it was about 99%. The FTIR spectrum of caprine milk could probably be used to predict several traits and to authenticate the origin of goat milk.

In its 800-year history, the University of Padova (Padua, Italy) has come to play an important role in the development of animal science in Italy and Europe. Having founded the oldest university botanical garden (1545; UNESCO World Heritage Site) and anatomical theatre (1595), and awarded the first university degree to a woman (Elena-Lucrezia Cornaro-Piscopia, 1678), the University instituted the Public School of Agriculture (Orto Agrario) in 1765 and the Collegium Zooiatricum in 1773. Between 1986 and 1995, under the leadership of Professor Mario Bonsembiante, animal scientist and rector, the scientific-technological center of Agripolis was established, including the Faculties of Agriculture and Veterinary Medicine (with the experimental farm and veterinary hospital), the Istituto Zooprofilattico delle Venezie, and Veneto Agricoltura. Agripolis rapidly became the most innovative center in Italy for teaching, research, and extension services in the fields of agricultural, animal, veterinary, food and environmental sciences. In the last 35 years, the number of Italian animal scientists has remained almost constant, but they have moved on to more innovative research topics closer to society’s demand and concerns, and significantly enhanced their international standing, reaching the first position in Europe and the fourth in the world in terms of citations.

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.