Presence of pathogen DNA in milk harvested from quarters is associated to changes in cows' milk yield and composition

BACKGROUND

Intramammary infection is the result of invasion and multiplication of microorganisms in the mammary gland and commonly leads to mastitis in dairy animals. Although much has been done to improve cows' udder health, mastitis remains a significant and costly health issue for dairy farmers, especially if subclinical. In this study, quarter milk samples from clinically healthy cows were harvested to detect pathogens via quantitative PCR (qPCR) and evaluate changes in individual milk traits according to the number of quarters infected and the type of microorganism(s). A commercial qPCR kit was used for detection of Mycoplasma bovis, Mycoplasma spp., Staphylococcus aureus, coagulase-negative staphylococci (CNS), Streptococcus agalactiae, Streptococcus dysgalactiae, Streptococcus uberis, Prototheca spp., Escherichia coli, Klebsiella spp., Enterococcus spp. and Lactococcus lactis ssp. lactis. Quarter and pooled milk information of 383 Holstein, 132 Simmental, 129 Rendena, and 112 Jersey cows in 9 Italian single-breed herds was available.

RESULTS

Among the cows with pathogen(s) present in at least 1 quarter, CNS was the most commonly detected DNA, followed by Streptococcus uberis, Mycoplasma bovis, and Streptococcus agalactiae. Cows negative to qPCR were 206 and had the lowest milk somatic cell count. Viceversa, cows with DNA isolated in ≥ 3 quarters were those with the highest somatic cell count. Moreover, when major pathogens were isolated in ≥ 3 quarters, milk had the lowest casein index and lactose content. In animals with pathogen(s) DNA isolated, the extent with whom milk yield and major solids were impaired did not significantly differ between major and minor pathogens.

CONCLUSIONS

The effect of the number of affected quarters on the pool milk quality traits was investigated in clinically healthy cows using a commercial kit. Results remark the important negative effect of subclinical udder inflammations on milk yield and quality, but more efforts should be made to investigate the presence of untargeted microorganisms, as they may be potentially dangerous for cows. For a smarter use of antimicrobials, analysis of milk via qPCR is advisable - especially in cows at dry off - to identify quarters at high risk of inflammation and thus apply a targeted/tailored treatment.

Genetic parameters for novel mastitis traits defined by combining test-day somatic cell score and differential somatic cell count in the first lactation of Japanese Holsteins

In this study, we aimed to improve current udder health genetic evaluations by addressing the limitations of monthly sampled somatic cell score (SCS) for distinguishing cows with robust innate immunity from those susceptible to chronic infections. The objectives were to (1) establish novel somatic cell traits by integrating SCS and the differential somatic cell count (DSCC), which represents the combined proportion of polymorphonuclear leukocytes and lymphocytes in somatic cells and (2) estimate genetic parameters for the new traits, including their daily heritability and genetic correlations with milk production traits and SCS, using a random regression test-day model (RRTDM). We derived 3 traits, termed ML_SCS_DSCC, SCS_4_DSCC_65_binary, and ML_SCS_DSCC_binary, by using milk loss (ML) estimates at corresponding SCS and DSCC levels, thresholds established in previous studies, and a threshold established from milk loss estimates, respectively. Data consisted of test-day records collected during January 2021 through March 2022 from 265 herds in Hokkaido, Japan. From these records, we extracted records between 7 to 305 d in milk (DIM) in the first lactation to fit the RRTDM. The model included the random effect of herd-test-day, the fixed effect of year-month, fixed lactation curves nested with calving age groups, and random regressions with Legendre polynomials of order 3 for additive genetic and permanent environmental effects. The analysis was performed using Gibbs sampling with Gibbsf90+ software. The averages (ranges) of the daily heritability estimates over lactation were 0.086 (0.075-0.095) for SCS, 0.104 (0.073-0.127) for ML_SCS_DSCC, 0.137 (0.014-0.297) for SCS_4_DSCC_65_binary, and 0.138 (0.115-0.185) for ML_SCS_DSCC_binary; the heritability curve for SCS_4_DSCC_65_binary was erratic. Genetic correlations within the trait decreased as the DIM interval widened, especially for those integrating DSCC, indicating that these traits should be analyzed using RRTDM rather than repeatability models. The averages (ranges) of genetic correlations with milk yield over lactation were 0.01 (-0.22 to 0.28) for SCS, -0.05 (-0.40 to 0.13) for ML_SCS_DSCC, -0.08 (-0.17 to 0.09) for SCS_4_DSCC_65_binary, and -0.08 (-0.22 to 0.27) for ML_SCS_DSCC_binary. Compared with SCS, the newly defined traits exhibited slightly stronger negative genetic correlations with milk yield. Especially in late lactation stages, the genetic correlation between ML_SCS_DSCC and milk yield was significantly below zero, with a posterior median of -0.40. Furthermore, the new traits showed positive correlations with SCS, having estimates varying from 0.68 to 0.85 for ML_SCS_DSCC, 0.14 to 0.47 for SCS_4_DSCC_65_binary, and 0.61 to 0.66 for ML_SCS_DSCC_binary, depending on DIM. Considering that ML_SCS_DSCC and ML_SCS_DSCC_binary have relatively high heritability (compared with SCS) and favorable genetic correlations with milk production traits and SCS, their incorporation into breeding programs appears promising. Nevertheless, their genetic relationships with (sub)clinical mastitis require further investigation.

The use of milk Fourier-transform infrared spectra for predicting cheesemaking traits in Grana Padano Protected Designation of Origin cheese

The prediction of the cheese yield (%CY) traits for curd, solids, and retained water and the amount of fat, protein, solids, and energy recovered from the milk into the curd (%REC) by Bayesian models, using Fourier-transform infrared spectroscopy (FTIR), can be of significant economic interest to the dairy industry and can contribute to the improvement of the cheese process efficiency. The yields give a quantitative measure of the ratio between weights of the input and output of the process, whereas the nutrient recovery allows to assess the quantitative transfer of a component from milk to cheese (expressed in % of the initial weight). The aims of this study were: (1) to investigate the feasibility of using bulk milk spectra to predict %CY and %REC traits, and (2) to quantify the effect of the dairy industry and the contribution of single-spectrum wavelengths on the prediction accuracy of these traits using vat milk samples destined to the production of Grana Padano Protected Designation of Origin cheese. Information from 72 cheesemaking days (in total, 216 vats) from 3 dairy industries were collected. For each vat, the milk was weighed and analyzed for composition (total solids [TS], lactose, protein, and fat). After 48 h from cheesemaking, each cheese was weighed, and the resulting whey was sampled for composition as well (TS, lactose, protein, and fat). Two spectra from each milk sample were collected in the range between 5,011 and 925 cm-1 and averaged before the data analysis. The calibration models were developed via a Bayesian approach by using the BGLR (Bayesian Generalized Linear Regression) package of R software. The performance of the models was assessed by the coefficient of determination (R2VAL) and the root mean squared error (RMSEVAL) of validation. Random cross-validation (CVL) was applied [80% calibration and 20% validation set] with 10 replicates. Then, a stratified cross-validation (SCV) was performed to assess the effect of the dairy industry on prediction accuracy. The study was repeated using a selection of informative wavelengths to assess the necessity of using whole spectra to optimize prediction accuracy. Results showed the feasibility of using FTIR spectra and Bayesian models to predict cheesemaking traits. The R2VAL values obtained with the CVL procedure were promising in particular for the %CY and %REC for protein, ranging from 0.44 to 0.66 with very low RMSEVAL (from 0.16 to 0.53). Prediction accuracy obtained with the SCV was strongly influenced by the dairy factory industry. The general low values gained with the SCV do not permit a practical application of this approach, but they highlight the importance of building calibration models with a dataset covering the largest possible sample variability. This study also demonstrated that the use of the full FTIR spectra may be redundant for the prediction of the cheesemaking traits and that a specific selection of the most informative wavelengths led to improved prediction accuracy. This could lead to the development of dedicated spectrometers using selected wavelengths with built-in calibrations for the online prediction of these innovative traits.

Exploring the sources of variation of electrical conductivity and total and differential somatic cell count in Italian Mediterranean buffaloes

In the buffalo dairy sector, a huge effort is still needed to improve mastitis prevention, detection, and management. Electrical conductivity (EC) and total somatic cell count (SCC) are well-known indirect indicators of mastitis. Differential somatic cell count (DSCC), which represents the proportion of neutrophils and lymphocytes on the total SCC, is instead a novel phenotype collected in the dairy cattle sector in the last lustrum. As little is known about this novel trait in dairy buffalo, in the present study we explored the nongenetic factors affecting DSCC, as well as EC and total somatic cell score (SCS), in the Italian Mediterranean buffalo. The data set used for the analysis included 14,571 test-day (TD) records of 1,501 animals from 6 herds, and climatic information of the sampling locations. The original data were filtered to exclude animals with less than 3 TD per lactation and, for the investigated traits, outliers beyond 4 standard deviations. In the statistical model we included the fixed effects of herd (6 classes), days in milk (DIM; 10 classes of 30 d, with the last being an open class until 360 d), parity (6 classes, from 1 to 6+), year-season of calving (11 classes, from summer 2019 to winter 2021/2022), year-season of sampling (9 classes, from spring 2020 to spring 2022), production level (4 classes based on quartiles of average milk yield by herd), and temperature-humidity index (THI; 4 classes based on quartiles, calculated using the average temperature and relative humidity of the 5 d before sampling). Average EC, SCS, and DSCC vary across herds. Considering DIM, greater EC values were observed at the beginning and the end of lactation; SCS was slightly lower, but DSCC was greater around the lactation peak. Increased EC, SCS, and DSCC levels with increasing parity were reported. Year-season calving and year-season sampling only slightly affected the variation of the investigated traits. Milk of high-producing buffaloes was characterized by lower EC and SCS mean values, nevertheless it had slightly greater DSCC percentages. Buffaloes grouped in the highest THI classes (classes 3 and 4) showed, on average, greater EC, SCS, and DSCC in comparison to the lower classes, especially to class 2. Results of the present study represent a preliminary as well as necessary step for the possible future inclusion of EC, SCS, or DSCC in breeding programs aimed to improve mastitis resistance in dairy buffaloes.

A First Investigation into the Use of Differential Somatic Cell Count as a Predictor of Udder Health in Sheep

Differential somatic cell count (DSCC), the percentage of somatic cell count (SCC) due to polymorphonuclear leukocytes (PMNs) and lymphocytes (LYMs), is a promising effective diagnostic marker for dairy animals with infected mammary glands. Well-explored in dairy cows, DSCC is also potentially valid in sheep, where clinical and subclinical mastitis outbreaks are among the principal causes of culling. We pioneered the application of DSCC in dairy ewes by applying receiver-operating characteristic (ROC) curve analysis to define the most accurate thresholds to facilitate early discrimination of sheep with potential intramammary infection (IMI) from healthy animals. We tested four predefined SCC cut-offs established in previous research. Specifically, we applied SCC cut-offs of 265 × 103 cells/mL, 500 × 103 cells/mL, 645 × 103 cells/mL, and 1000 × 103 cells/mL. The performance of DSCC as a diagnostic test was assessed by examining sensitivity (Se), specificity (Sp), positive predictive value (PPV), negative predictive value (NPV), and area under curve (AUC) analyses. The designated threshold value for DSCC in the detection of subclinical mastitis is established at 79.8%. This threshold exhibits Se and Sp of 0.84 and 0.81, accompanied by an AUC of 0.88. This study represents the inaugural exploration of the potential use of DSCC in sheep's milk as an early indicator of udder inflammation.

Heritability and genetic correlations of total and differential somatic cell count with milk yield and composition traits in Italian Simmental cows

Costs of production have deeply increased each year in the last decades, breeders are continuously looking for more cost effective and more efficient ways to produce milk. Despite the major signs of progress in productivity, it is fundamental to optimize rather than maximize the performances of the dairy cows. Mastitis is still a highly prevalent disease in the dairy sector which causes several economic losses and environmental effect. Its accurate and early diagnosis is crucial to improve profitability of dairy cows and contribute to a more sustainable dairy industry. Among mastitis reduction strategies, there is the urgent need to implement breeding objectives to select cows displaying mastitis resistance by investigating the genetic mechanisms at the base of the inflammatory response. Therefore, in this study we aimed to further understand the genetic background of the differential somatic cell count (DSCC), which provides thorough insights on the actual inflammatory status of the mammary glands. The objectives of this study were to estimate on a cohort of 20,215 Italian Simmental cows over a 3-yr period: (1) the heritability and repeatability values of somatic cell score (SCS) and DSCC, (2) the genetic and phenotypic correlations between these 2 traits and milk production and milk composition traits, (3) the heritability and repeatability values of SCS and DSCC within class of udder health status. Heritability was low both for SCS (0.06) and DSCC (0.08), whereas the repeatability values for these traits were 0.43 and 0.36, suggesting that the magnitude of cow permanent environmental effect for these traits is remarkable. The genetic and phenotypic correlation of SCS with DSCC was 0.612 and 0.605, respectively. Because both significantly differed from the unit, we must consider those traits as different ones. This latter aspect corroborates the need to consider the DSCC as a further indicator of inflammatory status which might be implemented in the Simmental breed genetic evaluation. It is worthy to mention that heritability estimates for SCS and DSCC were the highest in healthy cows compared with the other udder health classes. This implies that when the udder health status changes, it is most likely due to environmental factors rather than aspects related to the animal's genetics. In contrast, the highest additive genetic variance and heritability found for SCS and DSCC in the healthy group might reveal the potential to further implement breeding strategies to select for healthier animals.

Udder health-related traits in cow milk: phenotypic variability and effect on milk yield and composition

The milk differential somatic cell count (DSCC) has been proposed in recent years as a mean by which to better monitor the udder health status (UHS) in dairy cows. Milk DSCC is the amount of polymorphonuclear neutrophils and lymphocytes contributing to the total somatic cell count (SCC) and can be determined on a routine basis in individual milk samples subjected to official analysis. In the present study, 522 865 milk test-day records of 77 143 cows were scrutinised to identify factors affecting the variability of both DSCC and SCC in Holstein Friesian, Jersey, Simmental and Rendena cows through linear mixed models. The fixed effects were breed, parity, lactation stage, sampling season, and all the first-order interactions of breed. Cow and herd-test-date were considered as random. Subsequently, four UHS groups were created (1: SCC ≤ 200 000 cells/mL and DSCC ≤ 65%; 2: SCC ≤ 200 000 cells/mL and DSCC > 65%; 3: SCC > 200 000 cells/mL and DSCC > 65%; 4: SCC > 200 000 cells/mL and DSCC ≤ 65%) to compare milk yield and quality. Milk SCS and DSCC differed across lactation, parity, sampling season and breed. In particular, Simmental cows had the lowest SCC and Jersey the lowest DSCC. Depending on the breed, UHS affected daily milk yield and composition to a different extent. The UHS group 4, i.e. the one grouping test-day records with high SCC and low DSCC, presented the lowest estimate of milk yield and lactose content no matter the breeds. Our findings support that udder health-related traits (SCS and DSCC) are useful information to improve udder health at individual cow and herd levels. Moreover, the combination of SCS and DSCC is useful to monitor milk yield and composition.

Exploiting machine learning methods with monthly routine milk recording data and climatic information to predict subclinical mastitis in Italian Mediterranean buffaloes

Mastitis has detrimental effects on the world's dairy industry, reducing animal health, milk production and quality, as well as income for farmers. In addition, consumers' growing interest in food safety and rational usage of antibiotics highlights the need to develop novel strategies to improve mastitis detection, prevention, and management. In the present study we applied machine learning (ML) analyses to predict presence or absence of subclinical mastitis in Italian Mediterranean buffaloes, exploiting information collected the previous month during routine milk recording procedures, as well as climatic data. The data set included 3,891 records of 1,038 buffaloes from 6 herds located in Basilicata Region (South Italy). Prediction models were developed using 4 different ML algorithms (Generalized Linear Model, Support Vector Machines, Random Forest, and Neural Network) and 2 data set splitting approaches for the creation of the training and test sets (by record or by animal ID number, always with 80% of the data used for model training and the remaining 20% for model testing). Support Vector Machine was the best method to predict high or low somatic cell count at the subsequent test-day record in the validation set, and therefore it was used to estimate the contribution of each feature to the best model. Independently from the data set splitting approach, the most important features were somatic cell score, differential somatic cell count, electrical conductivity, and milk production. Among climatic data, the most informative were temperature and relative humidity. When the data were split by animal ID, an improvement in models' predictive performance on the test set was observed, suggesting this as the most appropriate data splitting approach in data sets with repeated measures to avoid data leakage. According to different metrics, Neural Network was the best method for making predictions on the test set. Our findings confirmed the promising role of ML methods to improve prevention and surveillance of subclinical mastitis, exploiting the large amount of data currently available to identify animals that would possibly have high somatic cell count the subsequent month.

Genomic-based genetic parameters for milkability traits derived from automatic milking systems in North American Holstein cattle

The number of dairy farms adopting automatic milking systems (AMS) has considerably increased around the world aiming to reduce labor costs, improve cow welfare, increase overall performance, and generate a large amount of daily data, including production, behavior, health, and milk quality records. In this context, this study aimed to (1) estimate genomic-based variance components for milkability traits derived from AMS in North American Holstein cattle based on random regression models; and (2) derive and estimate genetic parameters for novel behavioral indicators based on AMS-derived data. A total of 1,752,713 daily records collected using 36 milking robot stations and 70,958 test-day records from 4,118 genotyped Holstein cows were used in this study. A total of 57,600 SNP remained after quality control. The daily-measured traits evaluated were milk yield (MY, kg), somatic cell score (SCS, score unit), milk electrical conductivity (EC, mS), milking efficiency (ME, kg/min), average milk flow rate (FR, kg/min), maximum milk flow rate (FRM, kg/min), milking time (MT, min), milking failures (MFAIL), and milking refusals (MREF). Variance components and genetic parameters for MY, SCS, ME, FR, FRM, MT, and EC were estimated using the AIREMLF90 software under a random regression model fitting a third-order Legendre orthogonal polynomial. A threshold Bayesian model using the THRGIBBS1F90 software was used for genetically evaluating MFAIL and MREF. The daily heritability estimates across days in milk (DIM) ranged from 0.07 to 0.28 for MY, 0.02 to 0.08 for SCS, 0.38 to 0.49 for EC, 0.45 to 0.56 for ME, 0.43 to 0.52 for FR, 0.47 to 0.58 for FRM, and 0.22 to 0.28 for MT. The estimates of heritability (± SD) for MFAIL and MREF were 0.02 ± 0.01 and 0.09 ± 0.01, respectively. Slight differences in the genetic correlations were observed across DIM for each trait. Strong and positive genetic correlations were observed among ME, FR, and FRM, with estimates ranging from 0.94 to 0.99. Also, moderate to high and negative genetic correlations (ranging from -0.48 to -0.86) were observed between MT and other traits such as SCS, ME, FR, and FRM. The genetic correlation (± SD) between MFAIL and MREF was 0.25 ± 0.02, indicating that both traits are influenced by different sets of genes. High and negative genetic correlations were observed between MFAIL and FR (-0.58 ± 0.02) and MFAIL and FRM (-0.56 ± 0.02), indicating that cows with more MFAIL are those with lower FR. The use of random regression models is a useful alternative for genetically evaluating AMS-derived traits measured throughout the lactation. All the milkability traits evaluated in this study are heritable and have demonstrated selective potential, suggesting that their use in dairy cattle breeding programs can improve dairy production efficiency in AMS.

Herd and animal factors affect the variability of total and differential somatic cell count in bovine milk

The aim of this study was to quantify some environmental (individual herds, herd productivity, milking system, and season) and animal factors [individual animals, breed, days in milk (DIM) and parity] on the variability of the log-10 transformation of somatic cell count (LSCC) and differential somatic cell count (DSCC) on individual bovine milk. A total of 159,360 test-day records related to milk production and composition were extracted from 12,849 Holstein-Friesian and 9,275 Simmental cows distributed across 223 herds. Herds were classified into high and low productivity, defined according to the average daily milk net energy output (DMEO) yielded by the cows. Data included daily milk yield (DYM; kg/d), milk fat, protein, lactose, SCC, and DSCC, and information on herds (i.e., productivity, milking system). The daily production of total and differential somatic cells in milk was calculated and then log-10 transformed, obtaining DLSCC and DLDSCC, respectively. Data were analyzed using a mixed model including the effects of individual herd, animal, repeated measurements intra animal as random, and herd productivity, milking system, season, breed, DIM, parity, DIM × parity, breed × season, DIM × milking system and parity × milking system as fixed factors. Herds with a high DMEO were characterized by a lower content of LSCC and DSCC, and higher DLSCC and DLDSCC, compared to the low DMEO herds. The association between milking system and somatic cell traits suggested that the use of the automatic milking systems would not allow for a rapid intervention on the cow, as evidenced by the higher content of all somatic cell traits compared to the other milking systems. Season was an important source of variation, as evidenced by high LSCC and DSCC content in milk during summer. Breed of cow had a large influence, with Holstein-Friesian having greater LSCC, DSCC, DLSCC, and DLDSCC compared to Simmental. With regard to DIM, the variability of LSCC was mostly related to that of DSCC, showing an increase from calving to the end of lactation, and suggesting the higher occurrence of chronic mastitis in cows toward the end of lactation. All the somatic cell traits increased across number of parities, possibly because older cows may have increased susceptibility to intramammary infections.

Immune mechanisms, resistance genes, and their roles in the prevention of mastitis in dairy cows

Mastitis is one of the most important diseases of the mammary gland. The increased incidence of this disease in cows is due to the breeding of dairy cattle for higher yields, which is accompanied by an increased susceptibility to mastitis. Therefore, the difficulty involved with preventing this disease has increased. An integral part of current research is the elimination of mastitis in order to reduce the consumption of antibiotic drugs, thereby reducing the resistance of microorganisms and decreasing companies' economic losses due to mastitis (i.e. decreased milk yield, increased drug costs, and reduced milk supply). Susceptibility to mastitis is based on dairy cows' immunity, health, nutrition, and welfare. Thus, it is important to understand the immune processes in the body in order to increase the resistance of animals. Recently, various studies have focused on the selection of mastitis resistance genes. An important point is also the prevention of mastitis. This publication aims to describe the physiology of the mammary gland along with its immune mechanisms and to approximate their connection with potential mastitis resistance genes. This work describes various options for mastitis elimination and focuses on genetic selection and a closer specification of resistance genes to mastitis. Among the most promising resistance genes for mastitis, we consider CD14, CXCR1, lactoferrin, and lactoglobulin.

Formation of Blood Neutrophil Extracellular Traps Increases the Mastitis Risk of Dairy Cows During the Transition Period

The current study was conducted to analyze the functions of blood neutrophils in transition cows and their association with postpartum mastitis risk as indicated by somatic cell counts (SCCs) in milk. Seventy-six healthy Holstein dairy cows were monitored from Week 4 prepartum to Week 4 postpartum. Five dairy cows with low SCCs (38 ± 6.0 × 10³/mL) and five with high SCCs (3,753 ± 570.0 × 10³/mL) were selected based on milk SCCs during the first three weeks of lactation. At Week 1 pre- and postpartum, serum samples were obtained from each cow to measure neutrophil extracellular trap (NET)-related variables, and blood neutrophils were collected for transcriptome analysis by RNA sequencing. The serum concentration of NETs was significantly higher (P < 0.05) in cows with high SCCs than in cows with low SCCs (36.5 ± 2.92 vs. 18.4 ± 1.73 ng/mL). The transcriptomic analysis revealed that the transcriptome differences in neutrophils between high- and low-SCC cows were mainly in cell cycle-related pathways (42.6%), including the cell cycle, DNA damage, and chromosomal conformation, at Week 1 prepartum. The hub genes of these pathways were mainly involved in both the cell cycle and NETosis. These results indicated that the formation of NETs in the blood of transition dairy cows was different between cows with low and high SCCs, which may be used as a potential indicator for the prognosis of postpartum mastitis risk and management strategies of perinatal dairy cows.

Genetic parameters of differential somatic cell count, milk composition, and cheese-making traits measured and predicted using spectral data in Holstein cows

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 (k₂₀), 0.45 with rennet coagulation time estimated using the curd firming over time equation (RCT_eq), -0.39 with asymptotic potential curd firmness, -0.26 with maximum curd firmness (CF_max), 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 CF_max (-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), RCT_eq (0.53 and 0.41), t_max (0.38 and 0.48). Log MAC count was correlated with k₂₀ (+0.34), and log PMN-LYM count was correlated with CF_max (-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.

Comparison of machine learning methods to predict udder health status based on somatic cell counts in dairy cows

Bovine mastitis is one of the most important economic and health issues in dairy farms. Data collection during routine recording procedures and access to large datasets have shed the light on the possibility to use trained machine learning algorithms to predict the udder health status of cows. In this study, we compared eight different machine learning methods (Linear Discriminant Analysis, Generalized Linear Model with logit link function, Naïve Bayes, Classification and Regression Trees, k-Nearest Neighbors, Support Vector Machines, Random Forest and Neural Network) to predict udder health status of cows based on somatic cell counts. Prediction accuracies of all methods were above 75%. According to different metrics, Neural Network, Random Forest and linear methods had the best performance in predicting udder health classes at a given test-day (healthy or mastitic according to somatic cell count below or above a predefined threshold of 200,000 cells/mL) based on the cow’s milk traits recorded at previous test-day. Our findings suggest machine learning algorithms as a promising tool to improve decision making for farmers. Machine learning analysis would improve the surveillance methods and help farmers to identify in advance those cows that would possibly have high somatic cell count in the subsequent test-day.

Genome-Wide Associations for Microscopic Differential Somatic Cell Count and Specific Mastitis Pathogens in Holstein Cows in Compost-Bedded Pack and Cubicle Farming Systems

Simple Summary

New free walk housing systems such as compost-bedded pack barns might positively influence animal welfare. However, udder health can be negatively affected due to the microbial environment in the pack. Udder health depends on many factors, such as the environment, the feed, the pathogen species, and the genetic mechanisms of the cow’s immune system. For a more precise evaluation of udder health, we examined novel traits including specific mastitis pathogens and differential somatic cell fractions in milk. In order to identify possible candidate genes for udder health, a genome-wide association study, including single-nucleotide polymorphisms (SNP) by housing system interactions (compost-bedded pack barn and conventional cubicle barn), was performed. We identified two potential candidate genes for the interaction effect in relation to udder health. The identified potential candidate gene HEMK1 (HemK methyltransferase family member 1) is involved in immune system development, and CHL1 (cell adhesion molecule L1 like) has an immunosuppressive effect during stress conditions. The results suggest housing system-specific breeding strategies in order to improve udder health in compost-bedded pack and conventional cubicle barns.

Abstract

The aim of the present study was to detect significant SNP (single-nucleotide polymorphism) effects and to annotate potential candidate genes for novel udder health traits in two different farming systems. We focused on specific mastitis pathogens and differential somatic cell fractions from 2198 udder quarters of 537 genotyped Holstein Friesian cows. The farming systems comprised compost-bedded pack and conventional cubicle barns. We developed a computer algorithm for genome-wide association studies allowing the estimation of main SNP effects plus consideration of SNPs by farming system interactions. With regard to the main effect, 35 significant SNPs were detected on 14 different chromosomes for the cell fractions and the pathogens. Six SNPs were significant for the interaction effect with the farming system for most of the udder health traits. We inferred two possible candidate genes based on significant SNP interactions. HEMK1 plays a role in the development of the immune system, depending on environmental stressors. CHL1 is regulated in relation to stress level and influences immune system mechanisms. The significant interactions indicate that gene activity can fluctuate depending on environmental stressors. Phenotypically, the prevalence of mastitis indicators differed between systems, with a notably lower prevalence of minor bacterial indicators in compost systems.

Dairy Cows' Health during Alpine Summer Grazing as Assessed by Milk Traits, Including Differential Somatic Cell Count: A Case Study from Italy

Simple Summary

Dairy herds in alpine areas usually adopt summer grazing, mainly to reduce feeding costs. This practice is related to the maintenance of local traditions and to the manufacturing of niche dairy products. However, it is important to assess the impact of this practice on cattle health. This case study investigated how milk-related health traits vary across extensive grazing during the summer period, using data collected in a dairy herd whose cows were repeatedly controlled for individual milk samples. Although the transition from barn farming to pasture led to a reduction in milk production, proper grazing management can make dairy cows more resilient in terms of udder health and metabolic efficiency. Findings of the present research report suggested that pasture can be adopted to maintain dairy herd sustainability without impairing animal health.

Abstract

Extensive summer grazing is a dairy herd management practice frequently adopted in mountainous areas. Nowadays, this activity is threatened by its high labour demand, but it is fundamental for environmental, touristic and economic implications, as well as for the preservation of social and cultural traditions. Scarce information on the effects of such low-input farming systems on cattle health is available. Therefore, the present case study aimed at investigating how grazing may affect the health status of dairy cows by using milk traits routinely available from the national milk recording scheme. The research involved a dairy herd of 52 Simmental and 19 Holstein × Simmental crossbred cows. The herd had access to the pasture according to a rotational grazing scheme from late spring up to the end of summer. A total of 616 test day records collected immediately before and during the grazing season were used. Individual milk yield was registered during the milking procedure. Milk samples were analysed for composition (fat, protein, casein and lactose contents) and health-related milk indicators (electrical conductivity, urea and β-hydroxybutyrate) using mid-infrared spectroscopy. Somatic cell count (SCC) and differential SCC were also determined. Data were analysed with a linear mixed model, which included the fixed effects of the period of sampling, cow breed, stage of lactation and parity, and the random effects of cow nested within breed and the residual. The transition from barn farming to pasture had a negative effect on milk yield, together with a small deterioration of fat and protein percentages. Health-related milk indicators showed a minor deterioration of the fat to protein ratio, differential SCC and electrical conductivity, particularly towards the end of the grazing season, whereas the somatic cell score and β-hydroxybutyrate were relatively constant. Overall, the study showed that, when properly managed, pasture grazing does not have detrimental effects on dairy cows in terms of udder health and efficiency. Therefore, the proper management of cows on pasture can be a valuable solution to preserve the economic, social and environmental sustainability of small dairy farms in the alpine regions, without impairing cows’ health.

Total and Differential Cell Counts as a Tool to Identify Intramammary Infections in Cows after Calving

Simple Summary

Mastitis is a costly disease and needs to be identified as soon as possible to reduce the negative effect on milk quality and quantity and to maximize the chance of cure when an antimicrobial therapy is applied. Bacteriological diagnosis is expensive and not easily available in some areas, therefore approaches to reduce the number of samples to be taken, focusing the interest on cows with higher chances to have an intramammary infections are desirable. The results of our study based on a large database of quarter milk samples analyses including bacteriological analysis, total (SCC) and differential (DSCC) cell count in the first 5–30 days after calving suggest a new and sustainable approach. Indeed, a marker (PLCC) calculated by multiplying SCC and DSCC showed to have the lowest cost when applied to identify udder quarters at risk to have an intramammary infection due to major pathogens. Moreover, this approach as well as the one based on SCC became a benefit when the prevalence of these infections exceeds 10%, and it be of high interest, when selective dry cow therapy is applied, to improve animal health at the herd level.

Abstract

Milk differential somatic cells count (DSCC), made possible under field conditions by the recent availability of a high-throughput milk analyzer may represent an improvement in mastitis diagnosis. While an increasing number of studies reports data on DSCC on individual cow samples, very few concerns DSCC from quarter milk samples. This paper reports for the first time the results of a retrospective study aiming to assess the performance of total (SCC), DSCC, and a novel calculated marker (PLCC) measured on quarter milk samples as a method to identify cows at risk for intramammary infection (IMI) in the first 30 days after calving. Overall, 14,586 valid quarter milk samples (3658 cows) taken in the first 30 days of lactation were considered. Quarters with major pathogens (MP) IMI, as expected, showed significantly higher means for SCC, DSCC, and PLCC. The accuracy, sensitivity, and specificity of the diagnosis based on different cut-offs calculated by ROC analysis are relatively close among DSCC, PLCC, and SCC (up to cut-off of 200,000 cells/mL). However, decision-tree analysis which includes the costs of analysis, but also the costs of the actions taken after test results showed as PLCC has the lowest cost among the three markers, and PLCC and SCC are cost effective when MP prevalence is higher than 6–10%. This diagnostic approach is of high interest particularly when selective dry cow therapy is applied to improve animal health at the herd level.

Associations between differential somatic cell count and milk yield, quality, and technological characteristics in Holstein cows

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.

Differential Somatic Cell Count: Value for Udder Health Management

Intramammary infection (IMI) can cause mastitis, which is one of the costliest and most prevalent diseases in dairy cattle herds. Somatic cell count (SCC) is a well-established parameter to indicate IMI, and it represents the total count of immune cells in the milk. The differential somatic cell count (DSCC) has also long been suggested to indicate IMI, but no machine was available until recently to provide this parameter automatically. Two new machines have recently been introduced to measure the milk DSCC as an additional indicator of IMI. Here we provide insights about the DSCC measured by these two machines and the value it may provide for udder health management, based on the available literature. We also provide perspectives for future research to investigate potential value in using the DSCC to improve udder health.

Genetic Parameters of Bovine Milk Fatty Acid Profile, Yield, Composition, Total and Differential Somatic Cell Count

The growing interest of consumers for milk and dairy products of high nutritional value has pushed researchers to evaluate the feasibility of including fatty acids (FA) in selection programs to modify milk fat profile and improve its nutritional quality. Therefore, the aim of this study was to estimate genetic parameters of FA profile predicted by mid-infrared spectroscopy, milk yield, composition, and total and differential somatic cell count. Edited data included 35,331 test-day records of 25,407 Italian Holstein cows from 652 herds. Variance components and heritability were estimated using single-trait repeatability animal models, whereas bivariate repeatability animal models were used to estimate genetic and phenotypic correlations between traits, including the fixed effects of stage of lactation, parity, and herd-test-date, and the random effects of additive genetic animal, cow permanent environment and the residual. Heritabilities and genetic correlations obtained in the present study reflected both the origins of FA (extracted from the blood or synthesized de novo by the mammary gland) and their grouping according to saturation or chain length. In addition, correlations among FA groups were in line with correlation among individual FA. Moderate negative genetic correlations between FA and milk yield and moderate to strong positive correlations with fat, protein, and casein percentages suggest that actual selection programs are currently affecting all FA groups, not only the desired ones (e.g., polyunsaturated FA). The absence of association with differential somatic cell count and the weak association with somatic cell score indicate that selection on FA profile would not affect selection on resistance to mastitis and vice versa. In conclusion, our findings suggest that genetic selection on FA content is feasible, as FA are variable and moderately heritable. Nevertheless, in the light of correlations with other milk traits estimated in this study, a clear breeding goal should first be established.

Basal Levels of CD18 Antigen Presenting Cells in Cow Milk Associate with Copy Number Variation of Fc Gamma Receptors

Differentiation of cells by flow cytometry provides informative somatic cell counts (SCCs) that allow analyzing leukocyte population patterns in udder infections of different etiologies. Postulating that this approach also enhances the statistical power to detect genetic variants linked to cell levels in milk of healthy mammary glands, we used monoclonal antibodies anti-CD18, anti-CD4, anti--CD14, and anti-PMN to count cells presenting these surface antigens, and performed a genome-wide association study of these counts in 125 Israeli Holsteins genotyped using SNP BeadChips. We identified an informative haplotype of 15 SNPs in the centromeric end of BTA3 that was strongly associated with CD18 cells (p < 2.3 × 10-9). Within this region, examination of the network of genes interacting with ITGB2 (CD18) indicated an Fc-γ-receptor gene cluster, including FCGR2A (CD32). Sanger-sequence analysis of FCGR2s-linked exon 3 variation to CD18 counts. Meta-analysis of RNA-Seq data revealed a significant negative correlation (R = -0.51) between expression of CD32 and CD18 in milk. Assembly of DNA-Seq reads uncovered FCGR copy-number variation and a variant, designated V7, was abundant in dairy cattle, probably reflecting adaptation to selection pressure for low SCC in Holstein milk.

Differential Somatic Cell Count as a Novel Indicator of Milk Quality in Dairy Cows

Recent available instruments allow to record the number of differential somatic cell count (DSCC), representing the combined proportion of polymorphonuclear leukocytes and lymphocytes, on a large number of milk samples. Milk DSCC provides indirect information on the udder health status of dairy cows. However, literature is limited regarding the effect of DSCC on milk composition at the individual cow level, as well as its relation to the somatic cell score (SCS). Hence, the aims of this study were to (i) investigate the effect of different levels of DSCC on milk composition (fat, protein, casein, casein index, and lactose) and (ii) explore the combined effect of DSCC and SCS on these traits. Statistical models included the fixed effects of days in milk, parity, SCS, DSCC and the interaction between SCS × DSCC, and the random effects of herd, animal within parity, and repeated measurements within cow. Results evidenced a decrease of milk fat and an increase in milk fatty acids at increasing DSCC levels, while protein, casein and their proportion showed their lowest values at the highest DSCC. A positive association was found between DSCC and lactose. The interaction between SCS and DSCC was important for lactose and casein index, as they varied differently upon high and low SCS and according to DSCC levels.