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.
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.
The objectives of this study were to explore the use of Fourier-transform infrared (FTIR) spectroscopy on individual sheep milk samples for predicting cheese-making traits, and to test the effect of ...the farm variability on their prediction accuracy. For each of 121 ewes from 4 farms, a laboratory model cheese was produced, and 3 actual cheese yield traits (fresh cheese, cheese solids, and cheese water) and 4 milk nutrient recovery traits (fat, protein, total solids, and energy) in the curd were measured. Calibration equations were developed using a Bayesian approach with 2 different scenarios: (1) a random cross-validation (80% calibration; 20% validation set), and (2) a leave-one-out validation (3 farms used as calibration, and the remaining one as validation set) to assess the accuracy of prediction of samples from external farms, not included in calibration set. The best performance was obtained for predicting the yield and recovery of total solids, justifying for the practical application of the method at sheep population and dairy industry levels. Performances for the remaining traits were lower, but still useful for the monitoring of the milk processing in the case of fresh curd and recovery of energy. Insufficient accuracies were found for the recovery of protein and fat, highlighting the complex nature of the relationships among the milk nutrients and their recovery in the curd. The leave-one-out validation procedure, as expected, showed lower prediction accuracies, as a result of the characteristics of the farming systems, which were different between calibration and validation sets. In this regard, the inclusion of information related to the farm could help to improve the prediction accuracy of these traits. Overall, a large contribution to the prediction of the cheese-making traits came from the areas known as “water” and “fingerprint” regions. These findings suggest that, according to the traits studied, the inclusion of water regions for the development of the prediction equation models is fundamental to maintain a high prediction accuracy. However, further studies are necessary to better understand the role of specific absorbance peaks and their contribution to the prediction of cheese-making traits, to offer reliable tools applicable along the dairy ovine chain.
The cheese-derived strains
84C isolated from Nostrano cheese, and
DSM 32386 isolated from Traditional Mountain Malga cheese have been previously reported as γ-aminobutyric acid (GABA)-producers
. In ...the present study, the ability of these strains to produce GABA was studied in experimental raw milk cheeses, with the aim to investigate the effect of the culture and the ripening time on the GABA concentration. The cultures used consisted on
84C alone (84C) or in combination with
DSM 32386 (84C-DSM). The control culture was a commercial
strain, which was tested alone (CTRL) or in combination with the
DSM 32386 (CTRL-DSM). The pH evolution, microbiological counts, MiSeq Illumina and UHPLC-HQOMS analysis on milk and cheese samples were performed after 2, 9, and 20 days ripening. During the whole ripening, the pH was always under 5.5 in all batches. The concentration of GABA increased during ripening, with the highest content in 84C after 9 days ripening (84 ± 37 mg/kg), in 84C-DSM and CTRL-DSM after 20 days ripening (91 ± 28 and 88 ± 24 mg/kg, respectively). The data obtained support the hypothesis that
84C and
DSM 32386 could be exploited as functional cultures, improving the
bio-synthesis of GABA during cheese ripening.
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 (%CYCURD), total solids (%CYSOLIDS), and water retained (%CYWATER) 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 (R2VAL) and the root mean square error of validation (RMSEVAL). In CrV, the formula with the best prediction accuracy for all %CY traits included fat, casein, and UHI (R2VAL = 0.65, 0.96, and 0.23 for %CYCURD, %CYSOLIDS, and %CYWATER, respectively). The WAV procedure showed R2VAL 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 %CYCURD and %CYWATER 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.
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.
Cheeses are produced by many different procedures, giving rise to many types differing in ripening time, size, shape, chemical composition, color, texture, and sensory properties. As the first step ...in a large project, our aim was to characterize and quantify the major sources of variation in cheese characteristics by sampling 1050 different cheeses manufactured by over 100 producers and grouped into 37 categories (16 with protected designation of origin, 4 traditional cheese categories, 3 pasta filata cheese categories, 5 flavored cheese categories, 2 goat milk categories, and 7 other categories ranging from very fresh to very hard cheeses). We obtained 17 traits from each cheese (shape, height, diameter, weight, moisture, fat, protein, water soluble nitrogen, ash, pH, 5 color traits, firmness, and adhesiveness). The main groups of cheese categories were characterized and are discussed in terms of the effects of the prevalent area of origin/feeding system, species of lactating females, main cheese-making technologies, and additives used. The results will allow us to proceed with the further steps, which will address the interrelationships among the different traits characterizing cheeses, detailed analyses of the nutrients affecting human health and sensorial fingerprinting.
Genetic diversity has become an urgent matter not only in small local breeds but also in more specialized ones. While the use of genomic data in livestock breeding programs increased genetic gain, ...there is increasing evidence that this benefit may be counterbalanced by the potential loss of genetic variability. Thus, in this study, we aimed to investigate the genetic diversity in the Italian Holstein dairy cattle using pedigree and genomic data from cows born between 2002 and 2020. We estimated variation in inbreeding, effective population size, and generation interval and compared those aspects prior to and after the introduction of genomic selection in the breed. The dataset contained 84,443 single-nucleotide polymorphisms (SNPs), and 74,485 cows were analyzed. Pedigree depth based on complete generation equivalent was equal to 10.67. A run of homozygosity (ROH) analysis was adopted to estimate SNP-based inbreeding (F
). The average pedigree inbreeding was 0.07, while the average F
was more than double, being equal to 0.17. The pattern of the effective population size based on pedigree and SNP data was similar although different in scale, with a constant decrease within the last five generations. The overall inbreeding rate (ΔF) per year was equal to +0.27% and +0.44% for F
and F
throughout the studied period, which corresponded to about +1.35% and +2.2% per generation, respectively. A significant increase in the ΔF was found since the introduction of genomic selection in the breed. This study in the Italian Holstein dairy cattle showed the importance of controlling the loss of genetic diversity to ensure the long-term sustainability of this breed, as well as to guarantee future market demands.
The assessment of the efficiency of the cheese-making process (ECMP) is crucial for the profitability of cheese-factories. A simple way to estimate the ECMP is the measure of the estimated ...cheese-making losses (ECL), expressed by the ratio between the concentration of each constituent in the residual whey and in the processed milk. The aim of this research was to evaluate the influence of the season and cheese factory on the efficiency of the cheese-making process in Parmigiano Reggiano cheese manufacture. The study followed the production of 288 Parmigiano Reggiano cheese on 12 batches in three commercial cheese factories. For each batch, samples of the processed milk and whey were collected. Protein, casein, and fat ECL resulted in an average of 27.01%, 0.72%, and 16.93%, respectively. Both milk crude protein and casein contents were negatively correlated with protein ECL, r = -0.141 (
≤ 0.05), and r = -0.223 (
≤ 0.001), respectively. The same parameters resulted in a negative correlation with casein ECL (
≤ 0.001) (r = -0.227 and -0.212, respectively). Moreover, fat ECL was correlated with worse milk coagulation properties and negatively correlated with casein content (r = -0.120;
≤ 0.05). In conclusion, ECLs depend on both milk characteristics and season.
The aims of this proof of principle study were to compare two different chemometric approaches using a Bayesian method, Partial Least Square (PLS) and PLS-discriminant analysis (DA), for the ...prediction of the chemical composition and texture properties of the Grana Padano (GP) and Parmigiano Reggiano (PR) PDO cheeses by using NIR and Raman spectra and quantify their ability to distinguish between the two PDO and among their ripening periods. For each dairy chain consortium, 9 cheese samples from 3 dairy industries were collected for a total of 18 cheese samples. Three seasoning times were chosen for each dairy industry: 12, 20, and 36 months for GP and 12, 24, and 36 months for PR. A portable NIR instrument (spectral range: 950-1,650 nm) was used on 3 selected spots on the paste of each cheese sample, for a total of 54 spectra collected. An Alpha300 R confocal Raman microscope was used to collect 10 individual spectra for each cheese sample in each spot for a total of 540 Raman spectra collected. After the detection of eventual outliers, the spectra were also concatenated together (NIR + Raman). All the cheese samples were assessed in terms of chemical composition and texture properties following the official reference methods. A Bayesian approach and PLS-DA were applied to the NIR, Raman, and fused spectra to predict the PDO type and seasoning time. The PLS-DA reached the best performances, with 100% correctly identified PDO type using Raman only. The fusion of the data improved the results in 60% of the cases with the Bayesian and of 40% with the PLS-DA approach. A Bayesian approach and a PLS procedure were applied to the NIR, Raman, and fused spectra to predict the chemical composition of the cheese samples and their texture properties. In this case, the best performance in validation was reached with the Bayesian method on Raman spectra for fat (R2VAL = 0.74). The fusion of the data was not always helpful in improving the prediction accuracy. Given the limitations associated with our sample set, future studies will expand the sample size and incorporate diverse PDO cheeses.