G. Visentin

Prediction of bovine milk technological traits from mid-infrared spectroscopy analysis in dairy cows

Rapid, cost-effective monitoring of milk technological traits is a significant challenge for dairy industries specialized in cheese manufacturing. The objective of the present study was to investigate the ability of mid-infrared spectroscopy to predict rennet coagulation time, curd-firming time, curd firmness at 30 and 60min after rennet addition, heat coagulation time, casein micelle size, and pH in cow milk samples, and to quantify associations between these milk technological traits and conventional milk quality traits. Samples (n=713) were collected from 605 cows from multiple herds; the samples represented multiple breeds, stages of lactation, parities, and milking times. Reference analyses were undertaken in accordance with standardized methods, and mid-infrared spectra in the range of 900 to 5,000cm(-1) were available for all samples. Prediction models were developed using partial least squares regression, and prediction accuracy was based on both cross and external validation. The proportion of variance explained by the prediction models in external validation was greatest for pH (71%), followed by rennet coagulation time (55%) and milk heat coagulation time (46%). Models to predict curd firmness 60min from rennet addition and casein micelle size, however, were poor, explaining only 25 and 13%, respectively, of the total variance in each trait within external validation. On average, all prediction models tended to be unbiased. The linear regression coefficient of the reference value on the predicted value varied from 0.17 (casein micelle size regression model) to 0.83 (pH regression model) but all differed from 1. The ratio performance deviation of 1.07 (casein micelle size prediction model) to 1.79 (pH prediction model) for all prediction models in the external validation was <2, suggesting that none of the prediction models could be used for analytical purposes. With the exception of casein micelle size and curd firmness at 60min after rennet addition, the developed prediction models may be useful as a screening method, because the concordance correlation coefficient ranged from 0.63 (heat coagulation time prediction model) to 0.84 (pH prediction model) in the external validation.

Prediction of individual milk proteins including free amino acids in bovine milk using mid-infrared spectroscopy and their correlations with milk processing characteristics

The aim of this study was to evaluate the effectiveness of mid-infrared spectroscopy in predicting milk protein and free amino acid (FAA) composition in bovine milk. Milk samples were collected from 7 Irish research herds and represented cows from a range of breeds, parities, and stages of lactation. Mid-infrared spectral data in the range of 900 to 5,000 cm(-1) were available for 730 milk samples; gold standard methods were used to quantify individual protein fractions and FAA of these samples with a view to predicting these gold standard protein fractions and FAA levels with available mid-infrared spectroscopy data. Separate prediction equations were developed for each trait using partial least squares regression; accuracy of prediction was assessed using both cross validation on a calibration data set (n=400 to 591 samples) and external validation on an independent data set (n=143 to 294 samples). The accuracy of prediction in external validation was the same irrespective of whether undertaken on the entire external validation data set or just within the Holstein-Friesian breed. The strongest coefficient of correlation obtained for protein fractions in external validation was 0.74, 0.69, and 0.67 for total casein, total β-lactoglobulin, and β-casein, respectively. Total proteins (i.e., total casein, total whey, and total lactoglobulin) were predicted with greater accuracy then their respective component traits; prediction accuracy using the infrared spectrum was superior to prediction using just milk protein concentration. Weak to moderate prediction accuracies were observed for FAA. The greatest coefficient of correlation in both cross validation and external validation was for Gly (0.75), indicating a moderate accuracy of prediction. Overall, the FAA prediction models overpredicted the gold standard values. Near-unity correlations existed between total casein and β-casein irrespective of whether the traits were based on the gold standard (0.92) or mid-infrared spectroscopy predictions (0.95). Weaker correlations among FAA were observed than the correlations among the protein fractions. Pearson correlations between gold standard protein fractions and the milk processing characteristics of rennet coagulation time, curd firming time, curd firmness, heat coagulating time, pH, and casein micelle size were weak to moderate and ranged from -0.48 (protein and pH) to 0.50 (total casein and a30). Pearson correlations between gold standard FAA and these milk processing characteristics were also weak to moderate and ranged from -0.60 (Val and pH) to 0.49 (Val and K20). Results from this study indicate that mid-infrared spectroscopy has the potential to predict protein fractions and some FAA in milk at a population level.

Predictive ability of mid-infrared spectroscopy for major mineral composition and coagulation traits of bovine milk by using the uninformative variable selection algorithm

Milk minerals and coagulation properties are important for both consumers and processors, and they can aid in increasing milk added value. However, large-scale monitoring of these traits is hampered by expensive and time-consuming reference analyses. The objective of the present study was to develop prediction models for major mineral contents (Ca, K, Mg, Na, and P) and milk coagulation properties (MCP: rennet coagulation time, curd-firming time, and curd firmness) using mid-infrared spectroscopy. Individual milk samples (n=923) of Holstein-Friesian, Brown Swiss, Alpine Grey, and Simmental cows were collected from single-breed herds between January and December 2014. Reference analysis for the determination of both mineral contents and MCP was undertaken with standardized methods. For each milk sample, the mid-infrared spectrum in the range from 900 to 5,000cm(-1) was stored. Prediction models were calibrated using partial least squares regression coupled with a wavenumber selection technique called uninformative variable elimination, to improve model accuracy, and validated both internally and externally. The average reduction of wavenumbers used in partial least squares regression was 80%, which was accompanied by an average increment of 20% of the explained variance in external validation. The proportion of explained variance in external validation was about 70% for P, K, Ca, and Mg, and it was lower (40%) for Na. Milk coagulation properties prediction models explained between 54% (rennet coagulation time) and 56% (curd-firming time) of the total variance in external validation. The ratio of standard deviation of each trait to the respective root mean square error of prediction, which is an indicator of the predictive ability of an equation, suggested that the developed models might be effective for screening and collection of milk minerals and coagulation properties at the population level. Although prediction equations were not accurate enough to be proposed for analytic purposes, mid-infrared spectroscopy predictions could be evaluated as phenotypic information to genetically improve milk minerals and MCP on a large scale.

Effectiveness of mid-infrared spectroscopy to predict the color of bovine milk and the relationship between milk color and traditional milk quality traits

The color of milk affects the subsequent color features of the resulting dairy products; milk color is also related to milk fat concentration. The objective of the present study was to quantify the ability of mid-infrared spectroscopy (MIRS) to predict color-related traits in milk samples and to estimate the correlations between these color-related characteristics and traditional milk quality traits. Mid-infrared spectral data were available on 601 milk samples from 529 cows, all of which had corresponding gold standard milk color measures determined using a Chroma Meter (Konica Minolta Sensing Europe, Nieuwegein, the Netherlands); milk color was expressed using the CIELAB uniform color space. Separate prediction equations were developed for each of the 3 color parameters (L*=lightness, a*=greenness, b*=yellowness) using partial least squares regression. Accuracy of prediction was determined using both cross validation on a calibration data set (n=422 to 457 samples) and external validation on a data set of 144 to 152 samples. Moderate accuracy of prediction was achieved for the b* index (coefficient of correlation for external validation=0.72), although poor predictive ability was obtained for both a* and L* indices (coefficient of correlation for external validation of 0.30 and 0.55, respectively). The linear regression coefficient of the gold standard values on the respective MIRS-predicted values of a*, L*, and b* was 0.81, 0.88, and 0.96, respectively; only the regression coefficient on L* was different from 1. The mean bias of prediction (i.e., the average difference between the MIRS-predicted values and gold standard values in external validation) was not different from zero for any of 3 parameters evaluated. A moderate correlation (0.56) existed between the MIRS-predicted L* and b* indices, both of which were weakly correlated with the a* index. Milk fat, protein, and casein were moderately correlated with both the gold standard and MIRS-predicted values for b*. Results from the present study indicate that MIRS data provides an efficient, low-cost screening method to determine the b* color of milk at a population level.

Phenotypic characterisation of milk technological traits,protein fractions,and major mineral and fatty acid composition of Burlina cattle breed

Abstract The aim of the present study was to characterise milk of Burlina local cattle breed for traits of technological and nutritional relevance, such as milk coagulation properties (MCP), and protein, major mineral and fatty acid (FA) composition. Burlina is mainly reared in mountain areas of Veneto Region (Italy) and it has been inserted in conservation plans aiming to avoid biodiversity loss and marginal pasture areas abandonment. Eighty-one individual milk samples were collected in four farms. Milk coagulation properties were determined using Formagraph, and protein, mineral and FA composition were analysed in high performance liquid chromatography, inductively coupled plasma optical emission spectrometry and gas chromatography, respectively. Results evidenced good protein percentage (3.38%) and considerable casein content (28.89 mg/mL), as well as a desirable FA profile, with ω-6 to ω-3 ratio of 4.04. Somatic cell score, averaging 3.13, is a trait that should be enhanced through the improvement of farm management. This would have positive effects on MCP. Among milk minerals, the most and less abundant were K (1493.53 mg/kg) and Mg (110.07 mg/kg), respectively. Overall, herd, parity and lactation stage explained moderate to low variation of the studied traits. Results of the present study could be useful to valorise Burlina local breed and preserve biodiversity in marginal areas.

Processing characteristics of dairy cow milk are moderately heritable

Milk processing attributes represent a group of milk quality traits that are important to the dairy industry to inform product portfolio. However, because of the resources required to routinely measure such quality traits, precise genetic parameter estimates from a large population of animals are lacking for these traits. Milk processing characteristics considered in the present study-rennet coagulation time, curd-firming time, curd firmness at 30 and 60 min after rennet addition, heat coagulation time, casein micelle size, and milk pH-were all estimated using mid-infrared spectroscopy prediction equations. Variance components for these traits were estimated using 136,807 test-day records from 5 to 305 d in milk (DIM) from 9,824 cows using random regressions to model the additive genetic and within-lactation permanent environmental variances. Heritability estimates ranged from 0.18 ± 0.01 (26 DIM) to 0.38 ± 0.02 (180 DIM) for rennet coagulation time; from 0.26 ± 0.02 (5 DIM) to 0.57 ± 0.02 (174 DIM) for curd-firming time; from 0.16 ± 0.01 (30 DIM) to 0.56 ± 0.02 (271 DIM) for curd firmness at 30 min; from 0.13 ± 0.01 (30 DIM) to 0.48 ± 0.02 (271 DIM) for curd firmness at 60 min; from 0.08 ± 0.01 (17 DIM) to 0.24 ± 0.01 (180 DIM) for heat coagulation time; from 0.23 ± 0.02 (30 DIM) to 0.43 ± 0.02 (261 DIM) for casein micelle size; and from 0.20 ± 0.01 (30 DIM) to 0.36 ± 0.02 (151 DIM) for milk pH. Within-trait genetic correlations across DIM weakened as the number of days between compared intervals increased but were mostly >0.4 except between the peripheries of the lactation. Eigenvalues and associated eigenfunctions of the additive genetic covariance matrix for all traits revealed that at least the 80% of the genetic variation among animals in lactation profiles was associated with the height of the lactation profile. Curd-firming time and curd firmness at 30 min were weakly to moderately genetically correlated with milk yield (from 0.33 ± 0.05 to 0.59 ± 0.05 for curd-firming time, and from -0.62 ± 0.03 to -0.21 ± 0.06 for curd firmness at 30 min). Milk protein concentration was strongly genetically correlated with curd firmness at 30 min (0.84 ± 0.02 to 0.94 ± 0.01) but only weakly genetically correlated with milk heat coagulation time (-0.27 ± 0.07 to 0.19 ± 0.06). Results from the present study indicate the existence of exploitable genetic variation for milk processing characteristics. Because of possible indirect deterioration in milk processing characteristics due to selection for greater milk yield, emphasis on milk processing characteristics is advised.

Cow and environmental factors associated with protein fractions and free amino acids predicted using mid-infrared spectroscopy in bovine milk

The objective of the present study was to identify the factors associated with both the protein composition and free amino acid (FAA) composition of bovine milk predicted using mid-infrared spectroscopy. Milk samples were available from 7 research herds and 69 commercial herds. The spectral data from the research herds comprised 94,286 separate morning and evening milk samples; the spectral data from the commercial herds comprised 40,260 milk samples representing a composite sample of both the morning and evening milkings. Mid-infrared spectroscopy prediction models developed in a previous study were applied to all spectra. Factors associated with the predicted protein and FAA composition were quantified using linear mixed models. Factors considered in the model included the fixed effects of calendar month of the test, milking time (i.e., morning, evening, or both combined), parity (1, 2, 3, 4, 5, and ≥6), stage of lactation, the interaction between parity and stage of lactation, breed proportion of the cow (Friesian, Jersey, Norwegian Red, Montbéliarde, and other), and both the general heterosis and recombination coefficients of the cow. Contemporary group as well as both within- and across-lactation permanent environmental effects were included in all models as random effects. Total proteins (i.e., total casein, CN; total whey; and total β-lactoglobulin) and protein fractions (with the exception of α-lactalbumin) decreased postcalving until 36 to 65 days in milk and increased thereafter. After adjusting the statistical model for differences in crude protein content and milk yield separately, irrespective of stage of lactation, younger animals produced more total proteins (i.e., total CN, total whey, and total β-lactoglobulin) as well as more total FAA, Glu, and Asp than their older contemporaries. The concentration of all protein fractions (except β-CN) in milk was greatest in the evening milk, even after adjusting for differences in the crude protein content of the milk. Relative to a purebred Holstein cow, Jersey cows, on average, produced a greater concentration of all CN fractions but less total FAA, Glu, Gly, Asp, and Val in milk. Relative to their respective purebred parental average, first-cross cows produced more total CN and more β-CN. Results from the present study indicate that many cow-level factors, as well as other factors, are associated with protein composition and FAA composition of bovine milk.

Phenotypic characterisation of major mineral composition predicted by mid-infrared spectroscopy in cow milk

Abstract Population-level phenotyping is, for some traits, hardly achievable due to limitation of reference analyses. Mid-infrared spectroscopy (MIRS) is a quick and cost-effective laboratory technique allowing the prediction of innovative milk quality features on a large scale. The present study aimed at investigating factors associated with milk Ca, K, Mg, Na and P contents predicted by MIRS models on a large multi-breed spectral data set of 123,240 test-day records. Two dairy (Holstein–Friesian and Brown Swiss) and two dual-purpose (Simmental and Alpine Grey) cattle breeds were considered. Sources of variation of predicted milk mineral contents were studied using linear mixed models, including the fixed effects of breed, month of sampling, stage of lactation, parity and interactions between the main effects. Random factors were herd nested within breed, cow nested within breed and the residual. Calcium content was greater in milk of dual-purpose than dairy cattle breeds. Simmental cows produced milk with the greatest content of Ca and Na, and milk of Holstein-Friesian had the lowest P content. Variation of content of Ca, Mg, Na and P throughout lactations exhibited an opposite trend to that of milk yield, with the lowest values around the peak of lactation. On the other hand, K content mirrored the trend of milk yield. Multiparous cows had lower content of milk minerals, with the exception of Na, compared with first parity animals. Environmental factors identified in the present study can be considered for within-breed genetic evaluation to adjust records of milk mineral contents for these effects.

Genetic and nongenetic factors associated with milk color in dairy cows

Milk color is one of the sensory properties that can influence consumer choice of one product over another and it influences the quality of processed dairy products. This study aims to quantify the cow-level genetic and nongenetic factors associated with bovine milk color traits. A total of 136,807 spectra from Irish commercial and research herds (with multiple breeds and crosses) were used. Milk lightness (Lˆ*), red-green index (aˆ*) and yellow-blue index (bˆ*) were predicted for individual milk samples using only the mid-infrared spectrum of the milk sample. Factors associated with milk color were breed, stage of lactation, parity, milking-time, udder health status, pasture grazing, and seasonal calving. (Co)variance components for Lˆ*,aˆ*, and bˆ* were estimated using random regressions on the additive genetic and within-lactation permanent environmental effects. Greater bˆ* value (i.e., more yellow color) was evident in milk from Jersey cows. Milk Lˆ* increased consistently with stage of lactation, whereas aˆ* increased until mid lactation to subsequently plateau. Milk bˆ* deteriorated until 31 to 60 DIM, but then improved thereafter until the end of lactation. Relative to multiparous cows, milk yielded by primiparae was, on average, lighter (i.e., greater Lˆ*), more red (i.e., greater aˆ*), and less yellow (i.e., lower bˆ*). Milk from the morning milk session had lower Lˆ*,aˆ*, and bˆ* Heritability estimates (±SE) for milk color varied between 0.15 ± 0.02 (30 DIM) and 0.46 ± 0.02 (210 DIM) for Lˆ*, between 0.09 ± 0.01 (30 DIM) and 0.15 ± 0.02 (305 DIM) for aˆ*, and between 0.18 ± 0.02 (21 DIM) and 0.56 ± 0.03 (305 DIM) for bˆ* For all the 3 milk color features, the within-trait genetic correlations approached unity as the time intervals compared shortened and were generally <0.40 between the peripheries of the lactation. Strong positive genetic correlations existed between bˆ* value and milk fat concentration, ranging from 0.82 ± 0.19 at 5 DIM to 0.96 ± 0.01 at 305 DIM and confirming the observed phenotypic correlation (0.64, SE = 0.01). Results of the present study suggest that breeding strategies for the enhancement of milk color traits could be implemented for dairy cattle populations. Such strategies, coupled with the knowledge of milk color traits variation due to nongenetic factors, may represent a tool for the dairy processors to reduce, if not eliminate, the use of artificial pigments during milk manufacturing.

Short communication: Prediction of milk coagulation and acidity traits in Mediterranean buffalo milk using Fourier-transform mid-infrared spectroscopy

Milk coagulation and acidity traits are important factors to inform the cheesemaking process. Those traits have been deeply studied in bovine milk, whereas scarce information is available for buffalo milk. However, the dairy industry is interested in a method to determine milk coagulation and acidity features quickly and in a cost-effective manner, which could be provided by Fourier-transform mid-infrared (FT-MIR) spectroscopy. The aim of this study was to evaluate the potential of FT-MIR to predict coagulation and acidity traits of Mediterranean buffalo milk. A total of 654 records from 36 herds located in central Italy with information on milk yield, somatic cell score, milk chemical composition, milk acidity [pH, titratable acidity (TA)], and milk coagulation properties (rennet coagulation time, curd firming time, and curd firmness) were available for statistical analysis. Reference measures of milk acidity and coagulation properties were matched with milk spectral information, and FT-MIR prediction models were built using partial least squares regression. The data set was divided into a calibration set (75%) and a validation set (25%). The capacity of FT-MIR spectroscopy to correctly classify milk samples based on their renneting ability was evaluated by a canonical discriminant analysis. Average values for milk coagulation traits were 13.32 min, 3.24 min, and 39.27 mm for rennet coagulation time, curd firming time, and curd firmness, respectively. Milk acidity traits averaged 6.66 (pH) and 7.22 Soxhlet-Henkel degrees/100 mL (TA). All milk coagulation and acidity traits, except for pH, had high variability (17 to 46%). Prediction models of coagulation traits were moderately to scarcely accurate, whereas the coefficients of determination of external validation were 0.76 and 0.66 for pH and TA, respectively. Canonical discriminant analysis indicated that information on milk coagulating ability is present in the MIR spectra, and the model correctly classified as noncoagulating the 91.57 and 67.86% of milk samples in the calibration and validation sets, respectively. In conclusion, our results can be relevant to the dairy industry to classify buffalo milk samples before processing.