Prediction and repeatability of milk coagulation properties and curd-firming modeling parameters of ovine milk using Fourier-transform infrared spectroscopy and Bayesian models.

The aim of this study was to apply Bayesian models to the Fourier-transform infrared spectroscopy spectra of individual sheep milk samples to derive calibration equations to predict traditional and modeled milk coagulation properties (MCP), and to assess the repeatability of MCP measures and their predictions. Data consisted of 1,002 individual milk samples collected from Sarda ewes reared in 22 farms in the region of Sardinia (Italy) for which MCP and modeled curd-firming parameters were available. Two milk samples were taken from 87 ewes and analyzed with the aim of estimating repeatability, whereas a single sample was taken from the other 915 ewes. Therefore, a total of 1,089 analyses were performed. For each sample, 2 spectra in the infrared region 5,011 to 925 cm-1 were available and averaged before data analysis. BayesB models were used to calibrate equations for each of the traits. Prediction accuracy was estimated for each trait and model using 20 replicates of a training-testing validation procedure. The repeatability of MCP measures and their predictions were also compared. The correlations between measured and predicted traits, in the external validation, were always higher than 0.5 (0.88 for rennet coagulation time). We confirmed that the most important element for finding the prediction accuracy is the repeatability of the gold standard analyses used for building calibration equations. Repeatability measures of the predicted traits were generally high (≥95%), even for those traits with moderate analytical repeatability. Our results show that Bayesian models applied to Fourier-transform infrared spectra are powerful tools for cheap and rapid prediction of important traits in ovine milk and, compared with other methods, could help in the interpretation of results.

Genetic parameters of cheese yield and curd nutrient recovery or whey loss traits predicted using Fourier-transform infrared spectroscopy of samples collected during milk recording on Holstein, Brown Swiss, and Simmental dairy cows.

Cheese yield is the most important technological parameter in the dairy industry in many countries. The aim of this study was to infer (co)variance components for cheese yields (CY) and nutrient recoveries in curd (REC) predicted using Fourier-transform infrared (FTIR) spectroscopy of samples collected during milk recording on Holstein, Brown Swiss, and Simmental dairy cows. A total of 311,354 FTIR spectra representing the test-day records of 29,208 dairy cows (Holstein, Brown Swiss, and Simmental) from 654 herds, collected over a 3-yr period, were available for the study. The traits of interest for each cow consisted of 3 cheese yield traits (%CY: fresh curd, curd total solids, and curd water as a percent of the weight of the processed milk), 4 curd nutrient recovery traits (REC: fat, protein, total solids, and the energy of the curd as a percent of the same nutrient in the processed milk), and 3 daily cheese production traits (daily fresh curd, total solids, and the water of the curd per cow). Calibration equations (freely available upon request to the corresponding author) were used to predict individual test-day observations for these traits. The (co)variance components were estimated for the CY, REC, milk production, and milk composition traits via a set of 4-trait analyses within each breed. All analyses were performed using REML and linear animal models. The heritabilities of the %CY were always higher for Holstein and Brown Swiss cows (0.22 to 0.33) compared with Simmental cows (0.14 to 0.18). In general, the fresh cheese yield (%CYCURD) showed genetic variation and heritability estimates that were slightly higher than those of its components, %CYSOLIDS and %CYWATER. The parameter RECPROTEIN was the most heritable trait in all the 3 breeds, with values ranging from 0.32 to 0.41. Our estimation of the genetic relationships of the CY and REC with milk production and composition revealed that the current selection strategies used in dairy cattle are expected to exert only limited effects on the REC traits. Instead, breeders may be able to exploit genetic variations in the %CY, particularly RECFAT and RECPROTEIN. This last component is not explained by the milk protein content, suggesting that its direct selection could be beneficial for cheese production aptitude. Collectively, our findings indicate that breeding strategies aimed at enhancing CY and REC could be easily and rapidly implemented for dairy cattle populations in which FTIR spectra are routinely acquired from individual milk samples.

Bayesian regression models outperform partial least squares methods for predicting milk components and technological properties using infrared spectral data.

The aim of this study was to assess the performance of Bayesian models commonly used for genomic selection to predict "difficult-to-predict" dairy traits, such as milk fatty acid (FA) expressed as percentage of total fatty acids, and technological properties, such as fresh cheese yield and protein recovery, using Fourier-transform infrared (FTIR) spectral data. Our main hypothesis was that Bayesian models that can estimate shrinkage and perform variable selection may improve our ability to predict FA traits and technological traits above and beyond what can be achieved using the current calibration models (e.g., partial least squares, PLS). To this end, we assessed a series of Bayesian methods and compared their prediction performance with that of PLS. The comparison between models was done using the same sets of data (i.e., same samples, same variability, same spectral treatment) for each trait. Data consisted of 1,264 individual milk samples collected from Brown Swiss cows for which gas chromatographic FA composition, milk coagulation properties, and cheese-yield traits were available. For each sample, 2 spectra in the infrared region from 5,011 to 925 cm(-1) were available and averaged before data analysis. Three Bayesian models: Bayesian ridge regression (Bayes RR), Bayes A, and Bayes B, and 2 reference models: PLS and modified PLS (MPLS) procedures, were used to calibrate equations for each of the traits. The Bayesian models used were implemented in the R package BGLR (http://cran.r-project.org/web/packages/BGLR/index.html), whereas the PLS and MPLS were those implemented in the WinISI II software (Infrasoft International LLC, State College, PA). Prediction accuracy was estimated for each trait and model using 25 replicates of a training-testing validation procedure. Compared with PLS, which is currently the most widely used calibration method, MPLS and the 3 Bayesian methods showed significantly greater prediction accuracy. Accuracy increased in moving from calibration to external validation methods, and in moving from PLS and MPLS to Bayesian methods, particularly Bayes A and Bayes B. The maximum R(2) value of validation was obtained with Bayes B and Bayes A. For the FA, C10:0 (% of each FA on total FA basis) had the highest R(2) (0.75, achieved with Bayes A and Bayes B), and among the technological traits, fresh cheese yield R(2) of 0.82 (achieved with Bayes B). These 2 methods have proven to be useful instruments in shrinking and selecting very informative wavelengths and inferring the structure and functions of the analyzed traits. We conclude that Bayesian models are powerful tools for deriving calibration equations, and, importantly, these equations can be easily developed using existing open-source software. As part of our study, we provide scripts based on the open source R software BGLR, which can be used to train customized prediction equations for other traits or populations.

Stress related changes during a half marathon in master endurance athletes.

AIM: The aim of the present study was to investigate heart rate (HR), salivary cortisol (sC) alpha-amylase (sAA) and rating of perceived exertion (RPE) in relation to competition outcome during a half marathon. METHODS: HR was monitored and salivary samples were collected during an official half marathon in five Master endurance runners (age 47 ± 7 years). RPE was collected using a 100-mm Visual Analogue Scale (VAS) 30 minutes after the end of competition. RESULTS: Performance corresponded to 94% of their personal best (PB). Athletes spent 53.7% of total race time at intensities above 95% HRmax. RPE showed values of 68 ± 8 mm. With respect to pre-competition values (25.54 ± 6.39 nmol/L), sC concentrations significantly increased (P=0.043) by 59% immediately after the race (40.54 ± 3.95 nmol/L) and remained elevated until 1 h post exercise. Pre-competition sAA concentrations (90.59 ± 42.86 U/mL) were 118% higher (P=0.043) with respect to time-matched baseline values (197.92 ± 132 U/mL). sAA increased (192%; P=0.043) immediately after the race and was higher than time-matched resting samples. The better each athlete performed the greater cortisol increase during exercise (P<0.001). Performance was not correlated to the anticipatory sAA (the percent difference between pre-competition values and time-matched baseline ones) or to the sAA increase during exercise. CONCLUSION: This is the first attempt to study the stress-related responses during official endurance competitions in master runners. Although the strict criteria of inclusion might have limited the statistical significance, the present findings indicate that endurance competition is a remarkable stressor for psycho-physiological aspects of master athletes.

Comparison between genetic parameters of cheese yield and nutrient recovery or whey loss traits measured from individual model cheese-making methods or predicted from unprocessed bovine milk samples using Fourier-transform infrared spectroscopy.

Cheese yield is an important technological trait in the dairy industry. The aim of this study was to infer the genetic parameters of some cheese yield-related traits predicted using Fourier-transform infrared (FTIR) spectral analysis and compare the results with those obtained using an individual model cheese-producing procedure. A total of 1,264 model cheeses were produced using 1,500-mL milk samples collected from individual Brown Swiss cows, and individual measurements were taken for 10 traits: 3 cheese yield traits (fresh curd, curd total solids, and curd water as a percent of the weight of the processed milk), 4 milk nutrient recovery traits (fat, protein, total solids, and energy of the curd as a percent of the same nutrient in the processed milk), and 3 daily cheese production traits per cow (fresh curd, total solids, and water weight of the curd). Each unprocessed milk sample was analyzed using a MilkoScan FT6000 (Foss, Hillerød, Denmark) over the spectral range, from 5,000 to 900 wavenumber × cm(-1). The FTIR spectrum-based prediction models for the previously mentioned traits were developed using modified partial least-square regression. Cross-validation of the whole data set yielded coefficients of determination between the predicted and measured values in cross-validation of 0.65 to 0.95 for all traits, except for the recovery of fat (0.41). A 3-fold external validation was also used, in which the available data were partitioned into 2 subsets: a training set (one-third of the herds) and a testing set (two-thirds). The training set was used to develop calibration equations, whereas the testing subsets were used for external validation of the calibration equations and to estimate the heritabilities and genetic correlations of the measured and FTIR-predicted phenotypes. The coefficients of determination between the predicted and measured values in cross-validation results obtained from the training sets were very similar to those obtained from the whole data set, but the coefficient of determination of validation values for the external validation sets were much lower for all traits (0.30 to 0.73), and particularly for fat recovery (0.05 to 0.18), for the training sets compared with the full data set. For each testing subset, the (co)variance components for the measured and FTIR-predicted phenotypes were estimated using bivariate Bayesian analyses and linear models. The intraherd heritabilities for the predicted traits obtained from our internal cross-validation using the whole data set ranged from 0.085 for daily yield of curd solids to 0.576 for protein recovery, and were similar to those obtained from the measured traits (0.079 to 0.586, respectively). The heritabilities estimated from the testing data set used for external validation were more variable but similar (on average) to the corresponding values obtained from the whole data set. Moreover, the genetic correlations between the predicted and measured traits were high in general (0.791 to 0.996), and they were always higher than the corresponding phenotypic correlations (0.383 to 0.995), especially for the external validation subset. In conclusion, we herein report that application of the cross-validation technique to the whole data set tended to overestimate the predictive ability of FTIR spectra, give more precise phenotypic predictions than the calibrations obtained using smaller data sets, and yield genetic correlations similar to those obtained from the measured traits. Collectively, our findings indicate that FTIR predictions have the potential to be used as indicator traits for the rapid and inexpensive selection of dairy populations for improvement of cheese yield, milk nutrient recovery in curd, and daily cheese production per cow.

The use of Fourier-transform infrared spectroscopy to predict cheese yield and nutrient recovery or whey loss traits from unprocessed bovine milk samples.

Cheese yield is an important technological trait in the dairy industry in many countries. The aim of this study was to evaluate the effectiveness of Fourier-transform infrared (FTIR) spectral analysis of fresh unprocessed milk samples for predicting cheese yield and nutrient recovery traits. A total of 1,264 model cheeses were obtained from 1,500-mL milk samples collected from individual Brown Swiss cows. Individual measurements of 7 new cheese yield-related traits were obtained from the laboratory cheese-making procedure, including the fresh cheese yield, total solid cheese yield, and the water retained in curd, all as a percentage of the processed milk, and nutrient recovery (fat, protein, total solids, and energy) in the curd as a percentage of the same nutrient contained in the milk. All individual milk samples were analyzed using a MilkoScan FT6000 over the spectral range from 5,000 to 900 wavenumber × cm(-1). Two spectral acquisitions were carried out for each sample and the results were averaged before data analysis. Different chemometric models were fitted and compared with the aim of improving the accuracy of the calibration equations for predicting these traits. The most accurate predictions were obtained for total solid cheese yield and fresh cheese yield, which exhibited coefficients of determination between the predicted and measured values in cross-validation (1-VR) of 0.95 and 0.83, respectively. A less favorable result was obtained for water retained in curd (1-VR=0.65). Promising results were obtained for recovered protein (1-VR=0.81), total solids (1-VR=0.86), and energy (1-VR=0.76), whereas recovered fat exhibited a low accuracy (1-VR=0.41). As FTIR spectroscopy is a rapid, cheap, high-throughput technique that is already used to collect standard milk recording data, these FTIR calibrations for cheese yield and nutrient recovery highlight additional potential applications of the technique in the dairy industry, especially for monitoring cheese-making processes and milk payment systems. In addition, the prediction models can be used to provide breeding organizations with information on new phenotypes for cheese yield and milk nutrient recovery, potentially allowing these traits to be enhanced through selection.

The use of session RPE for interval training in master endurance athletes: should rest be included?

AIM: The aim of this study was to evaluate the effectiveness of the session rate of perceived exertion (RPE) method as a tool to quantify internal training load during interval training in master athletes. In addition, we investigated whether it is appropriate to take into account rest periods when calculating the session-RPE. METHODS: Eight male master endurance athletes (age: 45.3±7.3 years; stature: 1.74±0.06 m; body mass: 64.9±9.1 kg) were monitored during an interval training session consisting of 5 x 1000 m performed at 95% of vVO2max with 5 min rest between bouts. Edwards' summated heart rate zone method was used as a reference measure and the session RPE rating was obtained using the CR10 Borg's scale modified by Foster. RESULTS: High (r: 0.82; R2: 0.67) and significant (P=0.013) correlation was observed between the Edwards' heart rate (HR) and the session-RPE method when rest periods are taken into account; meanwhile a higher significant correlation (r: 0.86; R2: 0.74; P=0.003) was found between Edwards' HR and the session-RPE methods when rest periods were eliminated for the session-RPE computation. CONCLUSION: Despite the rest period exclusion from the computation of session RPE seems more appropriate, the statistical analysis indicates that there is no significant difference between the two correlation coefficients. These findings suggest that the session-RPE can be a useful tool to monitor internal training load during interval training and that the inclusion/exclusion of rest periods in its computation needs further investigation.