The effects of dietary iodine content, milking system, and farming practices on milk iodine concentration and quality traits

Various management practices can influence milk quality traits in dairy cattle. As an example, an increasing investment in automatic milking system to substitute milking parlors has been observed in the last 2 decades in dairy farms which could have affected certain bulk milk quality traits. What is more, milking practices can also affect certain milk parameters; as an example, teat disinfectants containing I are used in commercial farms where pre- or postdipping is performed, leading to presence of some I in the bulk milk. However, this trace mineral is also supplied in cows' diet to fulfill their nutritional requirements, partly contributing to the milk I final concentration. Therefore, the aim of this study was to evaluate the sources of variation of milk I along with other traditional milk quality traits. A total of 91 dairy farms in northeastern Italy were enrolled in the study. In each farm, diet and bulk milk samples were collected on the same day for chemical analysis. Concentration of I, in particular, was determined in both milk and feed with gold standard. Pearson correlations were calculated among the traits available for milk and diet, and a general linear model was used to test significance of fixed effects (feeding system, milking system, farming system, herd size, herd stage of lactation, and sampling month) on milk quality traits including the I concentration. In the case of milk I, diet I and presence of I-based predipping and postdipping teat disinfect application were also tested as fixed effects. Results showed a positive linear correlation between milk and diet I content (correlation coefficient [r] = 0.78). Although milk I was also positively correlated with lactose content (r = 0.25), dietary I was not correlated with other milk traits. Milk I content was significantly affected by dietary I, I-based predipping teat disinfectant application, and herd composition. Compared with conventional farms, organic farms showed lower protein content and greater somatic cell score (SCS) but similar milk I. Milking system significantly affected only lactose content and SCS of milk. Sampling month was only significant for milk urea nitrogen and herd composition, feeding system, herd size, and herd average days in milk did not modify milk gross composition and SCS. In conclusion, dietary supply of I is the main factor affecting milk I concentration and findings suggest that I level in milk can be naturally improved in dairy cows by modulating the I content in the diet administered. However, further research is needed to evaluate the effect of I-based sanitizers on milk I.

Modeling and Optimization of Herb-Fortified Fresh Kombucha Cheese: An Artificial Neural Network Approach for Enhancing Quality Characteristics

In this study, an Artificial Neural Network (ANN) model is used to solve the complex task of producing fresh cheese with the desired quality parameters. The study focuses on kombucha fresh cheese samples fortified with ground wild thyme, supercritical fluid extract of wild thyme, ground sage and supercritical fluid extract of sage and optimizes the parameters of chemical composition, antioxidant potential and microbiological profile. The ANN models demonstrate robust generalization capabilities and accurately predict the observed results based on the input parameters. The optimal neural network model (MLP 6-10-16) with 10 neurons provides high r2 values (0.993 for training, 0.992 for testing, and 0.992 for validation cycles). The ANN model identified the optimal sample, a supercritical fluid extract of sage, on the 20th day of storage, showcasing specific favorable process parameters. These parameters encompass dry matter, fat, ash, proteins, water activity, pH, antioxidant potential (TP, DPPH, ABTS, FRAP), and microbiological profile. These findings offer valuable insights into producing fresh cheese efficiently with the desired quality attributes. Moreover, they highlight the effectiveness of the ANN model in optimizing diverse parameters for enhanced product development in the dairy industry.

A Highly Selective Analytical Method Based on Salt-Assisted Liquid-Liquid Extraction for Trace-Level Enrichment of Multiclass Pesticide Residues in Cow Milk for Quantitative Liquid Chromatographic Analysis

In this study, a simple, inexpensive, selective, and fast salting-out assisted liquid-liquid extraction (SALLE) technique coupled with high-pressure liquid chromatography-diode array detection (HPLC-DAD) was developed for the extraction, preconcentration, and analysis of trace level seven multiclass pesticide residues in pasteurized and raw cow milk samples. The significant factors that affect the extent to which the target analytes are extracted, such as the type of extraction solvent and its volume, the type and concentration of salting-out salts, the pH of the solution, and the extraction time, have been investigated. Under optimum conditions, the correlation coefficient (r2) was obtained within a range of 0.9982-0.9997 for a broad linear range concentration of 2-1500 ng·mL-1. Reliable sensitivity was achieved with limits of detection (LODs) and limits of quantification (LOQs) ranging from 0.58-2.56 ng·mL-1 and 1.95-8.51 ng·mL-1, respectively. While precision with interday and intraday in terms of relative standard deviations (RSDs) was observed in the range of 1.97 - 7.88% and 4.52 - 8.04%, respectively. The results of the precision studies reveal that good repeatability and reproducibility (RSDs <9) were achieved, thus showing a low variability extraction of the developed method. Finally, the proposed and validated approach was effectively used to extract and determine pesticide residues in real milk matrices; however, the target analytes were not detected in all samples.

β-Casein A1 and A2: Effects of polymorphism on the cheese-making process

Of late, "A2 milk" has gained prominence in the dairy sector due to its potential implications in human health. Consequently, the frequency of A2 homozygous animals has considerably increased in many countries. To elucidate the potential implications that beta casein (β-CN) A1 and A2 may have on cheese-making traits, it is fundamental to investigate the relationships between the genetic polymorphisms and cheese-making traits at the dairy plant level. Thus, the aim of the present study was to evaluate the relevance of the β-CN A1/A2 polymorphism on detailed protein profile and cheese-making process in bulk milk. Based on the β-CN genotype of individual cows, 5 milk pools diverging for presence of the 2 β-CN variants were obtained: (1) 100% A1; (2) 75% A1 and 25% A2; (3) 50% A1 and 50% A2; (4) 25% A1 and 75% A2; and (5) 100% A2. For each cheese-making day (n = 6), 25 L of milk (divided into 5 pools, 5 L each) were processed, for a total of 30 cheese-making processes. Cheese yield, curd nutrient recovery, whey composition, and cheese composition were assessed. For every cheese-making process, detailed milk protein fractions were determined through reversed-phase HPLC. Data were analyzed by fitting a mixed model, which included the fixed effects of the 5 different pools, the protein and fat content as a covariate, and the random effect of the cheese-making sessions. Results showed that the percentage of κ-CN significantly decreased up to 2% when the proportion of β-CN A2 in the pool was ≥25%. An increase in the relative content of β-CN A2 (≥50% of total milk processed) was also associated with a significantly lower cheese yield both 1 and 48 h after cheese production, whereas no effects were observed after 7 d of ripening. Concordantly, recovery of nutrients reflected a more efficient process when the inclusion of β-CN A2 was ≤75%. Finally, no differences in the final cheese composition obtained by the different β-CN pools were observed.

Seasonal Variation in Raw Milk VOC Profile within Intensive Feeding Systems

The study aimed to assess the seasonal variation in raw milk volatile organic compounds (VOCs) from three indoor feeding systems based on maize silage (n = 31), silages/hay (n = 19) or hay (n = 16). After headspace solid-phase microextraction (HS-SPME), VOC profiles were determined by gas chromatography (GC). Chemical and VOC (log₁₀ transformations of the peak areas) data were submitted to a two-way ANOVA to assess the feeding system (FS) and season (S) effects; an interactive principal component analysis (iPCA) was also performed. The interaction FS × S was never significant. The FS showed the highest (p < 0.05) protein and casein content for hay-milk samples, while it did not affect any VOCs. Winter milk had higher (p < 0.05) proportions of protein, casein, fat and some carboxylic acids, while summer milk was higher (p < 0.05) in urea and 2-pentanol and methyl aldehydes. The iPCA confirmed a seasonal spatial separation. Carboxylic acids might generate from incomplete esterification in the mammary gland and/or milk lipolytic activity, while aldehydes seemed to be correlated with endogenous lipid or amino acid oxidation and/or feed transfer. The outcomes suggested that VOCs could be an operative support to trace raw milk for further mild processing.

Screening of antibiotic residues in raw bovine milk in Lombardy, Italy: Microbial growth inhibition assay and LC-HRMS technique integration for an accurate monitoring

Antibiotic residues in food of animal origin is a great concern for public health worldwide in terms of antibiotic resistance development, potential allergic reactions and disruption of intestinal flora equilibrium. In this study the presence of antibiotic residues in raw bovine milk samples collected from farms located in Lombardy region in Italy from 2018 to 2022 was assessed in the context of the national milk quality payment system. Samples were screened with microbiological growth inhibition test Delvotest ® SP NT and a very low positivity rate ranging from 0.1% to 0.07% over the four years was determined. A total of 79 positive samples were further analysed by LC-HRMS screening technique to confirm positivity and detect the specific antibiotic compound contaminating the sample. The β-lactam antibiotics resulted to be the most frequently detected, with the penicillin G being the most abundant compound. The data suggested that low levels of antibiotic contamination are consistently maintained over the last four years and the integration of the techniques used in this study is a valuable tool for a deep and precise monitoring of antibiotic residues in milk.

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.

Genomic Evaluation of Primiparous High-Producing Dairy Cows: Inbreeding Effects on Genotypic and Phenotypic Production-Reproductive Traits

Simple Summary

Improving the genomic prediction methodologies in high-producing dairy cattle is a key factor for the selection of suitable individuals to ensure better productivity. However, the most advanced prediction tools based on genotyping show ~75% reliability. Nowadays, the incorporation of new indices to genomic prediction methods, such as the Inbreeding Index (II), can significantly facilitate the selection of reliable production and reproductive traits for progeny selection. Thus, the objective of this study was to determine the impact of II (low: LI and high: HI), based on genomic analysis, and its effect on production and reproductive phenotypic traits in high-producing primiparous dairy cows. Individuals with II between ≥2.5 and ≤5.0 have shown up to a two-fold increase in negative correlations comparing LI versus HI genomic production and reproductive parameters, severely affecting important traits such as Milk Production at 305 d, Protein Production at 305 d, Fertility Index, and Daughter Pregnancy Rate. Therefore, high-producing dairy cows face an increased risk of negative II-derived effects in their selection programs, particularly at II ≥ 2.5.

Abstract

The main objective of this study was to analyze the effects of the inbreeding degree in high-producing primiparous dairy cows genotypically and phenotypically evaluated and its impacts on production and reproductive parameters. Eighty Holstein–Friesian primiparous cows (age: ~26 months; ~450 kg body weight) were previously genomically analyzed to determine the Inbreeding Index (II) and were divided into two groups: low inbreeding group (LI: <2.5; n = 40) and high inbreeding group (HI: ≥2.5 and ≤5.0; n = 40). Genomic determinations of production and reproductive parameters (14 in total), together with analyses of production (12) and reproductive (11) phenotypic parameters (23 in total) were carried out. Statistically significant differences were obtained between groups concerning the genomic parameters of Milk Production at 305 d and Protein Production at 305 d and the reproductive parameter Daughter Calving Ease, the first two being higher in cows of the HI group and the third lower in the LI group (p < 0.05). For the production phenotypic parameters, statistically significant differences were observed between both groups in the Total Fat, Total Protein, and Urea parameters, the first two being higher in the LI group (p < 0.05). Also, significant differences were observed in several reproductive phenotypic parameters, such as Number of Services per Conception, Calving to Conception Interval, Days Open Post Service, and Current Inter-Partum Period, all of which negatively influenced the HI group (p < 0.05). In addition, correlation analyses were performed between production and reproductive genomic parameters separately and in each consanguinity group. The results showed multiple positive and negative correlations between the production and reproductive parameters independently of the group analyzed, being these correlations more remarkable for the reproductive parameters in the LI group and the production parameters in the HI group (p < 0.05). In conclusion, the degree of inbreeding significantly influenced the results, affecting different genomic and phenotypic production and reproductive parameters in high-producing primiparous cows. The determination of the II in first-calf heifers is crucial to evaluate the negative effects associated with homozygosity avoiding an increase in inbreeding depression on production and reproductive traits.

Mid-infrared spectrometry prediction of the cheese-making properties of raw Montbéliarde milks from herds and cheese dairy vats used for the production of Protected Designation of Origin and Protected Geographical Indication cheeses in Franche-Comté

Franche-Comté is the primary producing region of Protected Designation of Origin cheeses in France. Normally, mid-infrared (MIR) prediction models for cheese-making property (CMP) traits are developed using individual bovine milks. However, considering the requests of all actors in the dairy sector, the present study aimed to assess the feasibility of MIR spectroscopy to develop CMP equations of Montbéliarde herd and dairy vat milks. For this purpose, 22 CMP traits were analyzed on samples collected in 2016 (half in February-March and half in May-June) from 100 commercial herds and 70 dairy vats (55 cheese dairies) located in Franche-Comté. These characteristics included 11 rennet coagulation traits and 8 lactic acidification traits measured in either soft cheese or pressed cooked cheese conditions and 3 laboratory curd yields. Models of MIR prediction for each of the 22 CMP traits were built using partial least squares regression with external validation by dividing the data set into calibration (70%) and validation (30%) sets. We confirmed that the variability of milk traits depends largely on the production scale and is higher for individual milk than for herd milk and even higher for vat milk. The best prediction models were obtained in herd milk samples for curd yields expressed in dry matter or fresh, with a coefficient of determination (R²) in external validation of 0.78 and 0.77, respectively. As with individual milk, these traits are closely related to the gross composition of the milk and therefore easier to predict by MIR spectroscopy. However, these curd yield traits were poorly predicted (R² = 0.58) in vat milk samples due to their lower variability. In herd milk samples, prediction models of other CMP traits were poorly accurate except for the ratio of the time to obtain a standard firmness to the rennet coagulation time in soft cheese or pressed cooked cheese conditions, which showed R² > 0.66 in external validation. Such trait is important in qualifying the behavior of milk during cheese production. Prediction models of other CMP traits for either herd or vat milk samples had poor accuracy, and further work is needed to improve their performance.

RESEARCH OF ORGANOLEPTIC PARAMETERS OF DUTCH CHEESE, PRODUCED FROM MILK OF COWS OF DIFFERENT BREEDS

The aim of this research is to study an influence of types of different breeds of milk cows on organoleptic properties of Dutch cheese at different rations of their food. It allows to receive cheese with prognosticated quality characteristics and to manage marketing strategies at cheese realization. The work presents the results of studying organoleptic properties of Dutch cheese, produced of milk of cows of different breeds, traditionally bred in the Northern-Eastern region of Ukraine. For studying an influence of cows feeding on exclusive properties of hard cheese, traditional food rations of cows were added with forages of a silage-hay type with Lucerne as a main source of proteins (44 % of the daily norm). Rations were practically identical with traditional ones by the content of energy and main food value factors and in general corresponded to norms of cows feeding. It has been demonstrated, that the high mark of cheese, up to 45 points for the taste and smell and total point mark 99 of 100 possible was obtained as a result of changing food rations at the expanse of introducing Lucerne silage in them. The data on optimization of food rations of certain cow breeds allow managing quality characteristics of milk and products of and are expedient for making cheese with unique regional characteristics. At conducting the comparative organoleptic assessment, there has been revealed a distinct dependence between a point mark of quality parameters of cheese on a breed of milk cows, and also on food rations of them. At the silage-hay food ration of animals, the received Dutch cheese had higher quality characteristics after 60 days of storage than cheese samples, obtained of cow milk at the traditional food ration.

Development of Fourier-transformed mid-infrared spectroscopy prediction models for major constituents of fractions of delactosated, defatted milk obtained through ultra- and nanofiltration

Milk filtration procedures are gaining relevance in the dairy industry because milk ultra- and nanofiltrates are used to increase milk processing efficiency, and as additives for products with improved nutraceutical properties. This study aimed to develop Fourier-transformed mid-infrared spectroscopy calibrations for ultra- and nanopermeate and retentate fractions of defatted and delactosated milk. A total of 154 samples from different milk fractions were collected and analyzed using reference methods to determine protein, solids-not-fat, glucose, and galactose content. The obtained values were matched with their respective Fourier-transformed mid-infrared spectroscopy spectra to develop new prediction models. Calibrations for each trait were built following 3 different approaches to get the best prediction models: (1) using the entire data set, (2) using 3 subsets based on component concentrations (level approach), and (3) using hierarchical clusters calculated with pairwise Mahalanobis distance among spectra (cluster approach). Calibrations were developed using partial least squares regression, after removing low signal-to-noise ratio wavelengths, and validated through a leave-one-out cross-validation procedure. In addition, the accuracy of the predicted values within each fraction was checked for each approach. Dividing the data set into subsets improved prediction models for each trait and for the samples in each milk fraction. Without considering milk fraction, the best improvement was observed for glucose and galactose. Glucose ratio performance deviation in cross-validation (RPD) increased from 7.42 to 11.31 and 11.06, for cluster and level approaches, respectively, whereas galactose RPD increased from 8.86 to 11.69 and 11.27 for cluster and level approaches, respectively. Considering milk fractions, the best improvement was observed for protein content, where RPD ranged from 0.08 to 6.06 for the whole data set calibration, whereas it ranged from 0.43 to 40.34 for the subset calibration approaches. Cluster and level approaches to build calibration models were comparable for samples from different fractions, suggesting that the 2 subsetting protocols should be both investigated to get the best prediction performances.