Variation in caprine milk composition and coagulation as affected by udder health indicators

Understanding the phenotypic and genetic background of the lactose content in Sarda dairy sheep

Lactose, the principal carbohydrate found in milk, plays an important role in the physiological processes of milk production because it is related to milk volume, and it is responsible for the osmotic equilibrium between blood and milk in the mammary gland. In this study, factors affecting lactose content (LC) in sheep milk are investigated. For this purpose, 2,358 test-day records were sampled from 509 ewes (3-7 records per animal). The LC and other main milk traits were analyzed using a mixed linear model that included days in milk (DIM) class, parity, lambing month, and type of lambing as fixed effects and animal, permanent environment, and flock test day as random effects. The pedigree-based approach was used to estimate the heritability and repeatability of LC. Moreover, the genomic background of LC was investigated through a GWAS. The LC was affected by all tested factors (i.e., DIM class, parity, lambing month, and type of lambing). Low heritability (0.10 ± 0.05) and moderate repeatability (0.42 ± 0.02) were estimated for LC. High negative genetic correlations were estimated between LC and NaCl (-0.99 ± 0.01) and between LC and somatic cell count (-0.94 ± 0.05). Only 2 markers passed the chromosome-wide Bonferroni threshold. Results of the present study, although obtained on a relatively small sample, suggest the possibility to include LC in the breeding programs, particularly because of its strong relationship with NaCl and somatic cell count.

Predictive formulas for different measures of cheese yield using milk composition from individual goat samples

The objective of this study was to develop formulas based on milk composition of individual goat samples for predicting cheese yield (%CY) traits (fresh curd, milk solids, and water retained in the curd). The specific aims were to assess and quantify (1) the contribution of major milk components (fat, protein, and casein) and udder health indicators (lactose, somatic cell count, pH, and bacterial count) on %CY traits (fresh curd, milk solids, and water retained in the curd); (2) the cheese-making method; and (3) goat breed effects on prediction accuracy of the %CY formulas. The %CY traits were analyzed in duplicate from 600 goats, using an individual laboratory cheese-making procedure (9-MilCA method; 9 mL of milk per observation) for a total of 1,200 observations. Goats were reared in 36 herds and belonged to 6 breeds (Saanen, Murciano-Granadina, Camosciata delle Alpi, Maltese, Sarda, and Sarda Primitiva). Fresh %CY (%CY_CURD), total solids (%CY_SOLIDS), and water retained (%CY_WATER) in the curd were used as response variables. Single and multiple linear regression models were tested via different combinations of standard milk components (fat, protein, casein) and indirect udder health indicators (UHI; lactose, somatic cell count, pH, and bacterial count). The 2 %CY observations within animal were averaged, and a cross-validation (CrV) scheme was adopted, in which 80% of observations were randomly assigned to the calibration (CAL) set and 20% to the validation (VAL) set. The procedure was repeated 10 times to account for sampling variability. Further, the model presenting the best prediction accuracy in CrV (i.e., comprehensive formula) was used in a secondary analysis to assess the accuracy of the %CY predictive formulas as part of the laboratory cheese-making procedure (within-animal validation, WAV), in which the first %CY observation within animal was assigned to CAL, and the second to the VAL set. Finally, a stratified CrV (SCrV) was adopted to assess the %CY traits prediction accuracy across goat breeds, again using the best model, in which 5 breeds were included in CAL and the remaining one in the VAL set. Fitting statistics of the formulas were assessed by coefficient of determination of validation (R²_VAL) and the root mean square error of validation (RMSE_VAL). In CrV, the formula with the best prediction accuracy for all %CY traits included fat, casein, and UHI (R²_VAL = 0.65, 0.96, and 0.23 for %CY_CURD, %CY_SOLIDS, and %CY_WATER, respectively). The WAV procedure showed R²_VAL higher than those obtained in CrV, evidencing a low effect of the 9-MilCA method and, indirectly, its high repeatability. In the SCrV, large differences for %CY_CURD and %CY_WATER among breeds evidenced that the breed is a fundamental factor to consider in %CY predictive formulas. These results may be useful to monitor milk composition and quantify the influence of milk traits in the composite selection indices of specific breeds, and for the direct genetic improvement of cheese production.

Pointing Out Some Issues Regarding Reproduction Management in Murciano-Granadina Goats

Simple Summary

The hypothesis of this experiment proposes that it could be possible to identify pregnant goats through maximum progesterone milk levels at any time in the pregnancy, and that there is an optimal moment to apply a lactation inhibitor to dry off lactating goats. The maximum progesterone concentration in milk varied depending on the season of the year, and those concentrations were similar for pregnant and non-pregnant goats, but significantly higher in the case of gestating goats with four foetuses, for which it would be possible to distinguish the pregnancy. The milk yield of goats at mating does not affect fertility until a value of at least 3250 mL/day. If using lactation inhibitors, their application up to the 10th week post-mating would be optimal for drying off lactating goats.

Abstract

Two of the most important problems in high-yielding dairy goat farms are early and accurate pregnancy diagnosis and the appropriate dry off of lactating does before the next kidding. The hypothesis posits that it could be possible to identify pregnant does through maximum progesterone milk levels at any time during the pregnancy, and that there is an optimal time to apply a lactation inhibitor to help dry off lactating does. Therefore, 114 Murciano-Granadina breed goats were used, from which 74 goats were inseminated at week 20 of lactation and samples of milk from pregnant and non-pregnant goats were taken at two-week intervals. The average maximum progesterone milk levels were higher outside the natural breeding season (40° latitude) than in the breeding season (11.6 ± 1.13 vs. 8.6 ± 1.02 ng/mL), although the levels from pregnant and non-pregnant goats were similar (10.85 ± 1.3 vs. 9.74 ± 1.6 ng/mL), except in the case of pregnancy with four foetuses (12.5 ± 1.3 ng/mL). Milk yield at mating does not affect fertility until a value of at least 3250 mL/day. Pregnancy started to affect milk yield up to the +7th week and was 59.9% lower in the +10th week after mating, so the use of lactation inhibitors could be more effective from this latter week. In conclusion, the results show that it is not possible to detect gestation in goats reliably through the maximum concentration of progesterone in milk at any time during lactation, except in the case of goats gestating four foetuses, that the milk yield of goats at mating does not affect fertility until a value of at least 3250 mL/day, and that from the 10th week post-mating, the application of lactation inhibitors would be optimal.

Rule Discovery in Milk Content towards Mastitis Diagnosis: Dealing with Farm Heterogeneity over Multiple Years through Classification Based on Associations

Simple Summary

Invisible (subclinical) mastitis decreases milk quality and production. Invisible mastitis is linked to an increased use of antimicrobials. The risk of the emergence of antimicrobial-resistant bacteria is a major public health concern worldwide. Therefore, early detection of infected cows is of great importance. Machine learning has opened a new avenue for early mastitis prediction based on simple and accessible milking parameters, such as milk volume, fat, protein, lactose, electrical conductivity (EC), milking time, and milking peak flow. However, farm heterogeneity is a major challenge where multiple patterns can predict mastitis. Here, we employed a classification based on associations and scaling approach for multiple pattern discovery over multiple years. The approach we have developed helps to address farm heterogeneity and generalise machine learning-based diagnosis of mastitis worldwide.

Abstract

Subclinical mastitis, an economically challenging disease of dairy cattle, is associated with an increased use of antimicrobials which reduces milk quantity and quality. It is more common than clinical mastitis and far more difficult to detect. Recently, much attention has been paid to the development of machine-learning expert systems for early detection of subclinical mastitis from milking features. However, differences between animals within a farm as well as between farms, particularly across multiple years, are major obstacles to the generalisation of machine learning models. Here, for the first time, we integrated scaling by quartiling with classification based on associations in a multi-year study to deal with farm heterogeneity by discovery of multiple patterns towards mastitis. The data were obtained from one farm comprising Holstein Friesian cows in Ongaonga, New Zealand, using an electronic automated monitoring system. The data collection was repeated annually over 3 consecutive years. Some discovered rules, such as when the milking peak flow is low, electrical conductivity (EC) of milk is low, milk lactose is low, milk fat is high, and milk volume is low, the cow has subclinical mastitis, reached high confidence (>70%) in multiple years. On averages, over 3 years, low level of milk lactose and high value of milk EC were part of 93% and 83.8% of all subclinical mastitis detecting rules, offering a reproducible pattern of subclinical mastitis detection. The scaled year-independent combinational rules provide an easy-to-apply and cost-effective machine-learning expert system for early detection of hidden mastitis using milking parameters.

The mineral profile affects the coagulation pattern and cheese-making efficiency of bovine milk

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

Somatic Cell Count in Goat Milk: An Indirect Quality Indicator

A high somatic cell count (SCC) impacts dairy quality to a large extent. The goal of this work was to investigate differences in goat milk composition and technological parameters according to SCC cut-off (600, 700, 800, and 1000.10³/mL). Thirty-four individual milk samples of White Shorthair goats in a similar stage of lactation were investigated. The first differences in milk quality appeared already at SCC cut-off of 600.10³/mL (5.58 LSCS-linear somatic cell score), yet the most striking differences were found for SCC over 1000.10³/mL (6.32 LSCS), which was expressed by lowering heat stability (126 vs. 217 s, p = 0.034), increasing protein (3.41 vs. 3.04%, p = 0.009), casein (2.80 vs. 2.44%, p = 0.034) and chloride (164 vs. 147 mg/100 mL, p = 0.004) levels, as well as non-fat dry matter (8.79 vs. 8.45%, p = 0.045). It has been shown that low levels of Staphylococcus spp. bacteria (120-1600 CFU/mL) in the mammary gland correlated with decreased lactose content (4.60 vs. 4.47 g/100 g, p = 0.022). Since our results indicate that even low SCC values may significantly affect the technological properties of goat milk, SCC should therefore be routinely screened and reported to dairy manufacturers to assure the consumer of high end-product quality.

Breed of goat affects the prediction accuracy of milk coagulation properties using Fourier-transform infrared spectroscopy

The prediction of traditional goat milk coagulation properties (MCP) and curd firmness over time (CF_t) parameters via Fourier-transform infrared (FTIR) spectroscopy can be of significant economic interest to the dairy industry and can contribute to the breeding objectives for the genetic improvement of dairy goat breeds. Therefore, the aims of this study were to (1) explore the variability of milk FTIR spectra from 4 goat breeds (Camosciata delle Alpi, Murciano-Granadina, Maltese, and Sarda), and to assess the possible discriminant power of milk FTIR spectra among breeds, (2) assess the viability to predict coagulation traits by using milk FTIR spectra, and (3) quantify the effect of the breed on the prediction accuracy of MCP and CF_t parameters. In total, 611 individual goat milk samples were used. Analysis of variance of measured MCP and CF_t parameters was carried out using a mixed model including the farm and pendulum as random factors, and breed, parity, and days in milk as fixed factors. Milk spectra for each goat were collected over the spectral range from wavenumber 5,011 to 925 × cm^-1. Discriminant analysis of principal components was used to assess the ability of FTIR spectra to identify breed of origin. A Bayesian model was used to calibrate equations for each coagulation trait. The accuracy of the model and the prediction equation was assessed by cross-validation (CRV; 80% training and 20% testing set) and stratified CRV (SCV; 3 breeds in the training set, one breed in the testing set) procedures. Prediction accuracy was assessed by using coefficient of determination of validation (R²_VAL), the root mean square error of validation (RMSE_VAL), and the ratio performance deviation. Moreover, measured and FTIR predicted traits were compared in the SCV procedure by assessing their least squares means for the breed effect, Pearson correlations, and variance heteroscedasticity. Results showed the feasibility of using FTIR spectra and multivariate analyses to correctly assign milk samples to their breeds of origin. The R²_VAL values obtained with the CRV procedure were moderate to high for the majority of coagulation traits, with RMSE_VAL and ratio performance deviation values increasing as the coagulation process progresses from rennet addition. Prediction accuracy obtained with the SCV were strongly influenced by the breed, presenting general low values restricting a practical application. In addition, the low Pearson correlation coefficients of Sarda breed for all the traits analyzed, and the heteroscedastic variances of Camosciata delle Alpi, Murciano-Granadina, and Maltese breeds, further indicated that it is fundamental to consider the differences existing among breeds for the prediction of milk coagulation traits.

Goat farm variability affects milk Fourier-transform infrared spectra used for predicting coagulation properties

Driven by the large amount of goat milk destined for cheese production, and to pioneer the goat cheese industry, the objective of this study was to assess the effect of farm in predicting goat milk-coagulation and curd-firmness traits via Fourier-transform infrared spectroscopy. Spectra from 452 Sarda goats belonging to 14 farms in central and southeast Sardinia (Italy) were collected. A Bayesian linear regression model was used, estimating all spectral wavelengths' effects simultaneously. Three traditional milk-coagulation properties [rennet coagulation time (min), time to curd firmness of 20 mm (min), and curd firmness 30 min after rennet addition (mm)] and 3 curd-firmness measures modeled over time [rennet coagulation time estimated according to curd firmness change over time (RCT_eq), instant curd-firming rate constant, and asymptotical curd firmness] were considered. A stratified cross validation (SCV) was assigned, evaluating each farm separately (validation set; VAL) and keeping the remaining farms to train (calibration set) the statistical model. Moreover, a SCV, where 20% of the goats randomly taken (10 replicates per farm) from the VAL farm entered the calibration set, was also considered (SCV₈₀). To assess model performance, coefficient of determination (R²_VAL) and the root mean squared error of validation were recorded. The R²_VAL varied between 0.14 and 0.45 (instant curd-firming rate constant and RCT_eq, respectively), albeit the standard deviation was approximating half of the mean for all the traits. Although average results of the 2 SCV procedures were similar, in SCV₈₀, the maximum R²_VAL increased at about 15% across traits, with the highest observed for time to curd firmness of 20 mm (20%) and the lowest for RCT_eq (6%). Further investigation evidenced important variability among farms, with R²_VAL for some of them being close to 0. Our work outlined the importance of considering the effect of farm when developing Fourier-transform infrared spectroscopy prediction equations for coagulation and curd-firmness traits in goats.

Characterization of milk composition, coagulation properties, and cheese-making ability of goats reared in extensive farms

The aims of this study were to explore the variability of milk composition, coagulation properties, and cheese-making traits of the Sarda goat breed, and to investigate the effects of animal and farm factors, and the geographic area (Central-East vs. South-West) of an insular region of Italy, Sardinia. A total of 570 Sarda goats reared in 21 farms were milk-sampled during morning milking. Individual milk samples were analyzed for composition, traditional milk coagulation properties (MCP), modeled curd-firming over time parameters (CF_t), and cheese-making traits (cheese yield, %CY; recovery of nutrients, %REC; daily cheese yield, dCY). Farms were classified into 2 categories based on milk energy level (MEL; high or low), defined according to the average net energy of milk daily produced by the lactating goats. Milk yield and composition were analyzed using a mixed model including the fixed effects of MEL, geographic area, days in milk, and parity, and the random effect of farm within MEL and geographic area. Data about MCP, CF_t, and the cheese-making process were analyzed using the same model, with the inclusion of the effects of animal and pendulum of the lactodynamograph instrument, allowing the measure of repeatability of these traits. Results showed that animal had greater influence on coagulation and cheese-making traits compared with farm effect. Days in milk influenced milk composition, whose changes partly reflected the modifications of %CY traits. Moreover, large differences were observed between primiparous and multiparous goats: primiparous goats produced less milk of better quality (higher fat, lower somatic cell and bacterial counts) and less cheese, but with higher recovery of fat and protein in the curd, compared with multiparous goats. The repeatability was very high, for both coagulation (84.0 to 98.8%) and cheese-making traits (89.7 to 99.9%). The effect of MEL was significant for daily productions of milk and cheese, coagulation time, and recovery of protein in the curd, which were better in high-MEL farms. As regards geographic area, milk composition and percentage cheese yield were superior in the Central-East area, whereas daily milk and cheese production and MCP were better in the South-West. This result was explainable by the phenomenon of crossbreeding Sarda goats with Maltese bucks, which occurred with greater intensity in the South-West than in the Central-East area of the island. The results provided by this study could be of great interest for the goat dairy sector. Indeed, the methods described in the present study could be applicable for other farming methods, goat breeds, and geographic areas. The collection of a wide range of phenotypes at individual animal level is fundamental for the characterization of local populations and can be used to guarantee breed conservation and the persistence of traditional farming systems, and to increase farmers' profit.

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.