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.

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.

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.