Insights on early response to acute heat shock of bovine mammary epithelial cells through a multimethod approach.

Heat stress is a significant challenge in dairy cattle herds, affecting milk production and quality, and generating important changes at the cellular level. Most in vitro research on heat shock (HS) effects on dairy cow mammary cells was focused on medium-long-term effects. In recent years, Fourier transform-infrared (FT-IR) micro-spectroscopy has been increasingly used to study the effects of several external stresses on different cell lines, down to the level of single cellular components, such as DNA/RNA, lipids, and proteins. In this study, the possible changes at the biochemical and molecular level induced by acute (30 min-2 h) HS in bovine mammary epithelial (BME-UV1) cells were investigated. The cells were exposed to different temperatures, thermoneutral (TN, 37 °C) and HS (42 °C), and FT-IR spectra were acquired to analyse the effects of HS on biochemical characteristics of BME-UV1 cellular components (proteins, lipids, and DNA/RNA). Moreover, cell viability assay, reactive oxygen species production, and mRNA expression of heat shock proteins (HSPA1A, HSP90AA1, GRP78, GRP94) and antioxidant genes (SOD1, SOD2) by RT-qPCR were also analysed. The FT-IR results showed a change already at 30 min of HS exposure, in the content of long-chain fatty acids, which probably acted as a response to a modification of membrane fluidity in HS cells compared with TN cells. After 2 h of HS exposure, modification of DNA/RNA activity and accumulation of aggregated proteins was highlighted in HS cells. The gene expression analyses showed the overexpression of HSPA1A and HSP90AA1 starting from 30 min up to 2 h in HS cells compared with TN cells. At 2 h of HS exposure, also the overexpression of GRP94 was observed in HS cells. Acute HS did not affect cell viability, reactive oxygen species level, and SOD1 and SOD2 gene expression of BME-UV1 cells. According to the results obtained, cells initiate early defence mechanisms in case of acute HS and probably this efficient response capacity may be decisive for tolerance to heat stress of dairy cattle.

The use of milk Fourier-transform infrared spectra for predicting cheesemaking traits in Grana Padano Protected Designation of Origin cheese.

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.

Investigation of phenotypic, genetic and genomic background of Milk spectra in Sarda dairy sheep.

Aim of this study was to analyse the genetic background of milk Fourier transform infrared (FTIR) spectra in dairy sheep. Individual milk FTIR spectra, with 1060 wavenumbers each, were available for 793 adult Sarda breed ewes genotyped at 45,813 SNP. The absorbance values of each wavenumber was analysed using a linear mixed model that included dim class, parity and lambing month as fixed effects and flock-test date and animal as random effects. The model was applied to estimate variance components and heritability and to perform a genome-wide association study for each wavenumber. Average h2 of wavenumbers absorbance was 0.13 ± 0.08, with the largest values observed in the regions associated with the characteristic bonds of carbonylic and methylenic groups of milk fat (h2 = 0.57 at 1724-1728 cm-1; and h2 = 0.34 at 2811-2834 cm-1, respectively). The absorbance values of wavenumbers were moderately correlated with the estimated heritabilities. After the Bonferroni correction, a total of nine markers were found to be significantly associated with 32 different wavenumbers. Of particular interest was the SNP s63269.1, mapped on chromosome 2, that was found to be associated with 27 wavenumbers. Genes previously found to be related to traits of interest (e.g. disease resistance, milk yield and quality, cheese firmness) are located close to the significant SNP. As expected, the heritability estimated for the absorbance of each wavenumbers seems to be associated with the related milk components.

Prediction of detailed blood metabolic profile using milk infrared spectra and machine learning methods in dairy cattle.

The adoption of preventive management decisions is crucial to dealing with metabolic impairments in dairy cattle. Various serum metabolites are known to be useful indicators of the health status of cows. In this study, we used milk Fourier-transform mid-infrared (FTIR) spectra and various machine learning (ML) algorithms to develop prediction equations for a panel of 29 blood metabolites, including those related to energy metabolism, liver function/hepatic damage, oxidative stress, inflammation/innate immunity, and minerals. For most traits, the data set comprised observations from 1,204 Holstein-Friesian dairy cows belonging to 5 herds. An exception was represented by ß-hydroxybutyrate prediction, which contained observations from 2,701 multibreed cows pertaining to 33 herds. The best predictive model was developed using an automatic ML algorithm that tested various methods, including elastic net, distributed random forest, gradient boosting machine, artificial neural network, and stacking ensemble. These ML predictions were compared with partial least squares regression, the most commonly used method for FTIR prediction of blood traits. Performance of each model was evaluated using 2 cross-validation (CV) scenarios: 5-fold random (CVr) and herd-out (CVh). We also tested the best model's ability to classify values precisely in the 2 extreme tails, namely, the 25th (Q25) and 75th (Q75) percentiles (true-positive prediction scenario). Compared with partial least squares regression, ML algorithms achieved more accurate performance. Specifically, elastic net increased the R2 value from 5% to 75% for CVr and 2% to 139% for CVh, whereas the stacking ensemble increased the R2 value from 4% to 70% for CVr and 4% to 150% for CVh. Considering the best model, with the CVr scenario, good prediction accuracies were obtained for glucose (R2 = 0.81), urea (R2 = 0.73), albumin (R2 = 0.75), total reactive oxygen metabolites (R2 = 0.79), total thiol groups (R2 = 0.76), ceruloplasmin (R2 = 0.74), total proteins (R2 = 0.81), globulins (R2 = 0.87), and Na (R2 = 0.72). Good prediction accuracy in classifying extreme values was achieved for glucose (Q25 = 70.8%, Q75 = 69.9%), albumin (Q25 = 72.3%), total reactive oxygen metabolites (Q25 = 75.1%, Q75 = 74%), thiol groups (Q75 = 70.4%), total proteins (Q25 = 72.4%, Q75 = 77.2.%), globulins (Q25 = 74.8%, Q75 = 81.5%), and haptoglobin (Q75 = 74.4%). In conclusion, our study shows that FTIR spectra can be used to predict blood metabolites with relatively good accuracy, depending on trait, and are a promising tool for large-scale monitoring.

Novel prediction models for hyperketonemia using bovine milk Fourier-transform infrared spectroscopy.

Metabolic diseases driven by negative energy balance in dairy cattle contribute to reduced milk production, increased disease incidence, culling, and death. Cow side tests for negative energy balance markers are available but are labor-intensive. Milk sample analysis using Fourier transform infrared spectroscopy (FTIR) allows for sampling numerous cows simultaneously. FTIR prediction models have moderate accuracy for hyperketonemia diagnosis (beta-hydroxybutyrate (BHB) ≥ 1.2 mmol/L). Most research using FTIR has focused on homogenous datasets and conventional prediction models, including partial least squares, linear discriminant analysis, and ElasticNet. Our objective was to evaluate more diverse modeling options, such as deep learning, gradient boosting machine models, and model ensembles for hyperketonemia classification. We compiled a sizable, heterogeneous dataset including milk FTIR and concurrent blood samples. Blood samples were tested for blood BHB, and wavenumber data was obtained from milk FTIR analysis. Using this dataset, we trained conventional prediction models and other options listed above. We demonstrate prediction model performance is similar for convolutional neural networks and ensemble models to simpler algorithm options. Results obtained from this study indicate that deep learning and model ensembles are potential algorithm options for predicting hyperketonemia in dairy cattle. Additionally, our results indicate hyperketonemia prediction models can be developed using heterogeneous datasets.

Predicting milk protein fractions using infrared spectroscopy and a gradient boosting machine for breeding purposes in Holstein cattle.

In recent years, increasing attention has been focused on the genetic evaluation of protein fractions in cow milk with the aim of improving milk quality and technological characteristics. In this context, advances in high-throughput phenotyping by Fourier transform infrared (FTIR) spectroscopy offer the opportunity for large-scale, efficient measurement of novel traits that can be exploited in breeding programs as indicator traits. We took milk samples from 2,558 Holstein cows belonging to 38 herds in northern Italy, operating under different production systems. Fourier transform infrared spectra were collected on the same day as milk sampling and stored for subsequent analysis. Two sets of data (i.e., phenotypes and FTIR spectra) collected in 2 different years (2013 and 2019-2020) were compiled. The following traits were assessed using HPLC: true protein, major casein fractions [αS1-casein (CN), αS2-CN, ß-CN, κ-CN, and glycosylated-κ-CN], and major whey proteins (ß-lactoglobulin and α-lactalbumin), all of which were measured both in grams per liter (g/L) and proportion of total nitrogen (% N). The FTIR predictions were calculated using the gradient boosting machine technique and tested by 3 different cross-validation (CRV) methods. We used the following CRV scenarios: (1) random 10-fold, which randomly split the whole into 10-folds of equal size (9-folds for training and 1-fold for validation); (2) herd/date-out CRV, which assigned 80% of herd/date as the training set with independence of 20% of herd/date assigned as the validation set; (3) forward/backward CRV, which split the data set in training and validation set according with the year of milk sampling (FTIR and gold standard data assessed in 2013 or 2019-2020) using the "old" and "new" databases for training and validation, and vice-versa with independence among them; (4) the CRV for genetic parameters (CRV-gen), where animals without pedigree as assigned as a fixed training population and animals with pedigree information was split in 5-folds, in which 1-fold was assigned to the fixed training population, and 4-folds were assigned to the validation set (independent from the training set). The results (i.e., measures and predictions) of CRV-gen were used to infer the genetic parameters for gold standard laboratory measurements (i.e., proteins assessed with HPLC) and FTIR-based predictions considering the CRV-gen scenario from a bi-trait animal model using single-step genomic BLUP. We found that the prediction accuracies of the gradient boosting machine equations differed according to the way in which the proteins were expressed, achieving higher accuracy when expressed in g/L than when expressed as % N in all CRV scenarios. Concerning the reproducibility of the equations over the different years, the results showed no relevant differences in predictive ability between using "old" data as the training set and "new" data as the validation set and vice-versa. Comparing the additive genetic variance estimates for milk protein fractions between the FTIR predicted and HPLC measures, we found reductions of -19.7% for milk protein fractions expressed in g/L, and -21.19% expressed as % N. Although we found reductions in the heritability estimates, they were small, with values ranging from -1.9 to -7.25% for g/L, and -1.6 to -7.9% for % N. The posterior distributions of the additive genetic correlations (ra) between the FTIR predictions and the laboratory measurements were generally high (>0.8), even when the milk protein fractions were expressed as % N. Our results show the potential of using FTIR predictions in breeding programs as indicator traits for the selection of animals to enhance milk protein fraction contents. We expect acceptable responses to selection due to the high genetic correlations between HPLC measurements and FTIR predictions.

Using vibrational molecular spectroscopy to reveal carbohydrate molecular structure properties of faba bean partitions and faba bean silage before and after rumen incubation in relation to nutrient availability and supply to dairy cattle.

To our knowledge, there is limited study on the relationship between the molecular structure of feed and nutrient availability in the ruminant system. The objective of this study is to use advanced vibrational molecular spectroscopy (attenuated total reflection [ATR]-Fourier transform infrared [FT/IR]) to reveal carbohydrate molecular structure properties of faba bean partitions (stem, leaf, whole pods [WP], and whole plant) and faba bean silage before and after rumen incubation in relation to nutrient availability and supply to dairy cattle. The study included the correlation between carbohydrate-related spectral profiles and chemical profiles, feed energy values, Cornell Net Carbohydrate and Protein System carbohydrate fractions, and rumen degradation parameters of faba bean samples (whole crop, stem, leaf, WP, and silage) before and after rumen incubation. FTIR spectra of faba bean sample before and after 12 and 24 h rumen incubations were collected with JASCO FT/IR-4200 with ATR at mid-IR range (ca. 4000-700 cm-1 ) with 128 scans and at 4 cm-1 resolution. The univariate molecular spectral analysis was carried out using OMNIC software. The results show that ATR-FT/IR spectroscopic technique could detect the change of microbial digestion to carbohydrate-related molecular structure. The spectral parameters of feed rumen incubation residues had a stronger correlation with less degradable carbohydrate fractions (neutral detergent fiber, acid detergent fiber, acid detergent lignin, hemicellulose, and cellulose) while spectral profiles of original faba samples had a stronger correlation with easily degradable carbohydrate fractions (starch). In conclusion, rumen degradation of carbohydrate contents can be reflected in the change of its molecular spectral profiles. The study shows that vibrational molecular spectroscopy (ATR-FT/IR) shows high potential as a fast analytical tool to evaluate and predict nutrient supply in the ruminant system.

Effect of CaCl2 on 2 heat-induced whey protein concentrate fibrillation pathways: Spontaneous and nuclear induction.

Amyloid fibrils have many excellent functional properties that facilitate their applications in the food industry. There are 2 pathways for whey protein concentrate (WPC) to form amyloid fibril aggregates: spontaneous pathway and nuclear induction pathway. Low ionic strength is a necessary condition for the spontaneous pathway to proceed successfully. In this paper, the effect of salt ions on 2 WPC fibrillation pathways was investigated by adding CaCl2. The results demonstrated WPC fibrils were unable to form normally through spontaneous pathway as adding CaCl2; but still could form through nuclear induction pathway with 20 to 30 mM CaCl2, the nuclei accelerated the fibrillation process led to the resistance to the disordered aggregation brought by CaCl2. Moreover, divalent cations (Ca2+, Mg2+) had much stronger effects than monovalent cations (Na+) on fibril formation, and the results of X-ray photoelectron spectrum together with Fourier-transform infrared spectroscopy suggested that Ca2+ had a greater effect on the fibril formation than Cl-.

Association of housing and management practices with milk yield, milk composition, and fatty acid profile, predicted using Fourier transform mid-infrared spectroscopy, in farms with automated milking systems.

Milk fatty acid (FA) profile can be divided into (1) de novo (C4-C14) that are synthesized in the mammary gland; (2) preformed (≥C18) that are absorbed from blood and originate from mobilized adipose tissues or dietary fat; and (3) mixed (C16), which have both origins. Our objectives were to describe the FA profile, as predicted using Fourier transform mid-infrared spectroscopy, of bulk tank milk from automated milking system (AMS) farms and to assess the association of management and housing factors with the bulk tank milk composition and FA profile of those AMS farms. The data used were collected from 124 commercial Canadian Holstein dairy farms with AMS, located in the provinces of Ontario (n = 68) and Quebec (n = 56). The farms were visited once from April to September 2019, and information were collected on barn design and herd management practices. Information regarding individual cow milk yield (kg/d), days in milk, parity, and the number of milking cows were automatically collected by the AMS units on each farm. These data were extracted for the entire period that the bulk tank milk samples were monitored, from April 2019 to April 2020 in Quebec and from August 2019 to May 2020 in Ontario. Across herds, milk yield averaged (mean ± standard error) 35.9 ± 0.4 kg/d, with 3.97 ± 0.01% fat and 3.09 ± 0.01% protein, whereas FA profile averaged 26.2 ± 0.1, 33.1 ± 0.1, and 40.7 ± 0.2 g/100 g of FA for de novo, mixed, and preformed, respectively. The FA yield averaged 0.34 ± 0.01, 0.44 ± 0.01, and 0.54 ± 0.01 kg/d for de novo, mixed, and preformed, respectively. Multivariable regression models were used to associate herd-level housing factors and management practices with milk production, composition, and FA profile. Milk yield was positively associated with using a robot feed pusher (+2.1 kg/d) and the use of deep bedding (+2.6 kg/d). The use of a robot feed pusher, deep bedding, and greater stall raking frequency were positively associated with greater yield (kg/d) of de novo, mixed, preformed, and de novo + mixed FA. Use of deep bedding was negatively associated with concentration of fat, de novo FA, mixed FA, and de novo + mixed FA, expressed in grams per 100 g (%) of milk. A wider lying alley width (≥305 cm) was associated with a greater concentration (g/100 g of milk) of de novo and de novo + mixed FA. Greater frequency of partial mixed ration delivery (>2×/d vs. 1 and 2×/d) was positively associated with a greater proportion (g/100 g of FA) of de novo, mixed, and de novo + mixed FA and negatively associated with the proportion of preformed FA. Overall, these associations indicated that bulk tank FA profile can be used as a tool to monitor and adjust management and housing in AMS farms.

Synthesis and physicochemical characterization of zinc-lactoferrin complexes.

One trend of the modern world is the search for new biologically active substances based on renewable resources. Milk proteins can be a solution for such purposes as they have been known for a long time as compounds that can be used for the manufacturing of multiple food and non-food products. Thus, the goal of the work was to investigate the parameters of Zn-bovine lactoferrin (bLTF) interactions, which enables the synthesis of Zn-rich protein complexes. Zinc-bLTF complexes can be used as food additives or wound-healing agents. Methodology of the study included bLTF characterization by sodium dodecyl sulfate-PAGE, MALDI-TOF, and MALDI-TOF/TOF mass spectrometry as well Zn-bLTF interactions by attenuated total reflection-Fourier-transform infrared, Raman spectroscopy, scanning and transmission microscopy, and zeta potential measurements. The obtained results revealed that the factors that affect Zn-bLTF interactions most significantly were found to be pH and ionic strength of the solution and, in particular, the concentration of Zn2+. These findings imply that these factors should be considered when aiming at the synthesis of Zn-bLTF metallocomplexes.

Inhibition of porcine epidemic diarrhea virus by Alpiniae oxyphyllae fructus polysaccharide 3.

Porcine epidemic diarrhea virus (PEDV) is a deadly pathogen that still plagues suckling piglets. However, there is still no anti-PEDV drug available in clinics. To develop potential anti-PEDV drugs, the antiviral activity of Alpiniae oxyphyllae fructus polysaccharide 3 (AOFP3) against PEDV infection in IPEC-J2 cells were assessed in our present study. The structural characterization of AOFP3 was studied by using HPAEC, GC-MS, FT-IR and NMR techniques. At the same time, the anti-PEDV activity of AOFP3 was investigated by performing RT-qPCR, Western blot and immunofluorescence assays. The results showed that AOFP3 (44.4 kDa) was composed of glucose and galacturonic acid at a molar ratio of 77.54:22.46 and consisted of →4)-α-D-Glcp-(1→, →4,6)-α-D-Glcp-(1→, T-α-D-Glcp-(1→ and →4)-α-D-GalAp-(1→. AOFP3 significantly decreased PEDV titer in IPEC-J2 cells and prevented cellular damage of IPEC-J2 cells caused by PEDV infection. Furthermore, AOFP3 showed an antioxidative activity in inhibiting PEDV reproduction. Therefore, AOFP3 was expected to be a material of anti-PEDV drug.

Developing a predictive model for beta-hydroxybutyrate and non-esterified fatty acids using milk fourier-transform infrared spectroscopy in dairy cows.

Negative energy balance following parturition predisposes dairy cattle to numerous metabolic disorders. Current diagnostics are limited by their labor requirements and inability to measure multiple metabolic markers simultaneously. Fourier-transform Infrared spectroscopy (FTIR) data, measured from milk samples, could improve the detection of metabolic disorders using routine milk samples from dairy farms. The objective of this study was to develop a predictive model for numeric values of blood beta-hydroxybutyrate (BHB) and blood non-esterified fatty acids (NEFA). The study utilized a dataset comprised of 622 observations with known blood BHB and blood NEFA values measured concurrently with the milk FTIR data. Using ElasticNet regression on milk FTIR data and production information, we built regression models for numeric blood BHB and blood NEFA prediction and classification models for blood BHB values greater than 1.2 mmol/L and blood NEFA values greater than 0.7 mmol/L. The R2 of the best fitting model was 0.56 and 0.51 for log-transformed BHB and log-transformed NEFA, respectively. The BHB classification model had a 90 % sensitivity and 83 % specificity and the NEFA classification model achieved a sensitivity of 73 % and specificity of 74 %. We applied our numeric prediction models to an external dataset (n = 9660) with known blood metabolites to validate the prediction accuracy of log-transformed blood BHB and log-transformed blood NEFA models. Log-transformed BHB root mean square error (RMSE) was 0.4018 and log-transformed NEFA RMSE was 0.4043. The second objective of this study was to develop a categorization for cows as either metabolically disordered or healthy. Responses to negative energy balance in transition cows are related to blood levels of BHB and NEFA. Cows suffering from metabolic disorders without elevated blood BHB values remain unidentified when detection is focused on blood BHB alone. To account for this differentiated metabolic response, we categorized cows as either 'metabolically healthy' or suffering a 'metabolic disorder' by using a combination of blood BHB, blood NEFA, and milk fat to protein quotient. We obtained a balanced accuracy of 94 % for the prediction of cow metabolic status. Direct prediction of metabolic status can be used to identify hyperketonemic cows in addition to cows exhibiting metabolic response patterns not detected by elevated blood BHB alone.

Evaluating the performance of machine learning methods and variable selection methods for predicting difficult-to-measure traits in Holstein dairy cattle using milk infrared spectral data.

Fourier-transform infrared (FTIR) spectroscopy is a powerful high-throughput phenotyping tool for predicting traits that are expensive and difficult to measure in dairy cattle. Calibration equations are often developed using standard methods, such as partial least squares (PLS) regression. Methods that employ penalization, rank-reduction, and variable selection, as well as being able to model the nonlinear relations between phenotype and FTIR, might offer improvements in predictive ability and model robustness. This study aimed to compare the predictive ability of 2 machine learning methods, namely random forest (RF) and gradient boosting machine (GBM), and penalized regression against PLS regression for predicting 3 phenotypes differing in terms of biological meaning and relationships with milk composition (i.e., phenotypes measurable directly and not directly in milk, reflecting different biological processes which can be captured using milk spectra) in Holstein-Friesian cattle under 2 cross-validation scenarios. The data set comprised phenotypic information from 471 Holstein-Friesian cows, and 3 target phenotypes were evaluated: (1) body condition score (BCS), (2) blood ß-hydroxybutyrate (BHB, mmol/L), and (3) κ-casein expressed as a percentage of nitrogen (κ-CN, % N). The data set was split considering 2 cross-validation scenarios: samples-out random in which the population was randomly split into 10-folds (8-folds for training and 1-fold for validation and testing); and herd/date-out in which the population was randomly assigned to training (70% herd), validation (10%), and testing (20% herd) based on the herd and date in which the samples were collected. The random grid search was performed using the training subset for the hyperparameter optimization and the validation set was used for the generalization of prediction error. The trained model was then used to assess the final prediction in the testing subset. The grid search for penalized regression evidenced that the elastic net (EN) was the best regularization with increase in predictive ability of 5%. The performance of PLS (standard model) was compared against 2 machine learning techniques and penalized regression using 2 cross-validation scenarios. Machine learning methods showed a greater predictive ability for BCS (0.63 for GBM and 0.61 for RF), BHB (0.80 for GBM and 0.79 for RF), and κ-CN (0.81 for GBM and 0.80 for RF) in samples-out cross-validation. Considering a herd/date-out cross-validation these values were 0.58 (GBM and RF) for BCS, 0.73 (GBM and RF) for BHB, and 0.77 (GBM and RF) for κ-CN. The GBM model tended to outperform other methods in predictive ability around 4%, 1%, and 7% for EN, RF, and PLS, respectively. The prediction accuracies of the GBM and RF models were similar, and differed statistically from the PLS model in samples-out random cross-validation. Although, machine learning techniques outperformed PLS in herd/date-out cross-validation, no significant differences were observed in terms of predictive ability due to the large standard deviation observed for predictions. Overall, GBM achieved the highest accuracy of FTIR-based prediction of the different phenotypic traits across the cross-validation scenarios. These results indicate that GBM is a promising method for obtaining more accurate FTIR-based predictions for different phenotypes in dairy cattle.

Breed of goat affects the prediction accuracy of milk coagulation properties using Fourier-transform infrared spectroscopy.

The prediction of traditional goat milk coagulation properties (MCP) and curd firmness over time (CFt) parameters via Fourier-transform infrared (FTIR) spectroscopy can be of significant economic interest to the dairy industry and can contribute to the breeding objectives for the genetic improvement of dairy goat breeds. Therefore, the aims of this study were to (1) explore the variability of milk FTIR spectra from 4 goat breeds (Camosciata delle Alpi, Murciano-Granadina, Maltese, and Sarda), and to assess the possible discriminant power of milk FTIR spectra among breeds, (2) assess the viability to predict coagulation traits by using milk FTIR spectra, and (3) quantify the effect of the breed on the prediction accuracy of MCP and CFt parameters. In total, 611 individual goat milk samples were used. Analysis of variance of measured MCP and CFt parameters was carried out using a mixed model including the farm and pendulum as random factors, and breed, parity, and days in milk as fixed factors. Milk spectra for each goat were collected over the spectral range from wavenumber 5,011 to 925 × cm-1. Discriminant analysis of principal components was used to assess the ability of FTIR spectra to identify breed of origin. A Bayesian model was used to calibrate equations for each coagulation trait. The accuracy of the model and the prediction equation was assessed by cross-validation (CRV; 80% training and 20% testing set) and stratified CRV (SCV; 3 breeds in the training set, one breed in the testing set) procedures. Prediction accuracy was assessed by using coefficient of determination of validation (R2VAL), the root mean square error of validation (RMSEVAL), and the ratio performance deviation. Moreover, measured and FTIR predicted traits were compared in the SCV procedure by assessing their least squares means for the breed effect, Pearson correlations, and variance heteroscedasticity. Results showed the feasibility of using FTIR spectra and multivariate analyses to correctly assign milk samples to their breeds of origin. The R2VAL values obtained with the CRV procedure were moderate to high for the majority of coagulation traits, with RMSEVAL and ratio performance deviation values increasing as the coagulation process progresses from rennet addition. Prediction accuracy obtained with the SCV were strongly influenced by the breed, presenting general low values restricting a practical application. In addition, the low Pearson correlation coefficients of Sarda breed for all the traits analyzed, and the heteroscedastic variances of Camosciata delle Alpi, Murciano-Granadina, and Maltese breeds, further indicated that it is fundamental to consider the differences existing among breeds for the prediction of milk coagulation traits.

Economic value of information from an alert system on physiological imbalance in fresh cows.

Physiological imbalance is an abnormal physiological condition that cannot be directly observed but is assumed to precede subclinical and clinical diseases in the beginning of lactation. Alert systems to detect the physiological imbalance in a cow using Fourier transform mid-infrared spectroscopy in milk have been developed. The objective of this study was to estimate the value of information provided from such system with different indicator accuracies, herd prevalence and prices. A decision tree was created to model the probabilities of detection and associated costs of test outcome, intervention and occurrence of disease. We assumed that the negative effect of physiological imbalance was the development of subclinical ketosis and that this negative effect was prevented by drenching the cows with propylene glycol for 5 days. We simulated the economic impact of subclinical ketosis mediated through physiological imbalance to be $194 per case. The results showed that if the alert system was highly accurate (Se = 0.99/Sp = 0.99), and the prevalence of physiological imbalance was 30 %, the value of information provided from the system is $19 per cow-year. In case the prevalence is 5 % or 50 %, the value of information is $3 and $13, respectively. These estimates for the value do not cover the capital costs and operational costs of the alert system. This study furthermore clearly demonstrated that in order to estimate the value of information correctly, it is important to consider that drenching all cows and not drenching any of the cows are the two relevant alternative options in the absence of the alert system. In conclusion, the decision tree and sensitivity analysis developed in this study show that final economic results are highly variable to the prevalence of physiological imbalance and highest at an intermediate prevalence. Other relevant factors are the costs associated with drenching and the cost associated with treating false positives and not treating false negatives. In addition, this study highlights the benefits of simulation to pinpoint where additional information is needed to further quantify the economic value and required accuracy of an indication-based intervention system.

Opportunities and limitations of milk mid-infrared spectra-based estimation of acetone and ß-hydroxybutyrate for the prediction of metabolic stress and ketosis in dairy cows.

Subclinical (SCK) and clinical (CK) ketosis are metabolic disorders responsible for big losses in dairy production. Although Fourier-transform mid-infrared spectrometry (FTIR) to predict ketosis in cows exposed to great metabolic stress was studied extensively, little is known about its suitability in predicting hyperketonemia using individual samples, e.g. in small dairy herds or when only few animals are at risk of ketosis. The objective of the present research was to determine the applicability of milk metabolites predicted by FTIR spectrometry in the individual screening for ketosis. In experiment 1, blood and milk samples were taken every two weeks after calving from Holstein (n = 80), Brown Swiss (n = 72) and Swiss Fleckvieh (n = 58) cows. In experiment 2, cows diagnosed with CK (n = 474) and 420 samples with blood ß-hydroxybutyrate [BHB] <1.0 mmol/l were used to investigate if CK could be detected by FTIR-predicted BHB and acetone from a preceding milk control. In experiment 3, correlations between data from an in farm automatic milk analyser and FTIR-predicted BHB and acetone from the monthly milk controls were evaluated. Hyperketonemia occurred in majority during the first eight weeks of lactation. Correlations between blood BHB and FTIR-predicted BHB and acetone were low (r = 0.37 and 0.12, respectively, P < 0.0001), as well as the percentage of true positive values (11.9 and 16.6%, respectively). No association of FTIR predicted ketone bodies with the interval of milk sampling relative to CK diagnosis was found. Data obtained from the automatic milk analyser were moderately correlated with the same day FTIR-predicted BHB analysis (r = 0.61). In conclusion, the low correlations with blood BHB and the small number of true positive samples discourage the use of milk mid-infrared spectrometry analyses as the only method to predict hyperketonemia at the individual cow level.

Pathophysiology and pathological remodelling associated with dilated cardiomyopathy in broiler chickens predisposed to heart pump failure.

Broiler chickens selected for rapid growth are highly susceptible to dilated cardiomyopathy (DCM). In order to elucidate the pathophysiology of DCM, the present study examines the fundamental features of pathological remodelling associated with DCM in broiler chickens using light microscopy, transmission electron microscopy (TEM), and synchrotron Fourier Transform Infrared (FTIR) micro-spectroscopy. The morphological features and FTIR spectra of the left ventricular myocardium were compared among broiler chickens affected by DCM with clinical signs of heart pump failure, apparently normal fast-growing broiler chickens showing signs of subclinical DCM (high risk of heart failure), slow-growing broiler chickens (low risk of heart failure) and Leghorn chickens (resistant to heart failure, used here as physiological reference). The findings indicate that DCM and heart pump failure in fast-growing broiler chickens are a result of a complex metabolic syndrome involving multiple catabolic pathways. Our data indicate that a good deal of DCM pathophysiology in chickens selected for rapid growth is associated with conformational changes of cardiac proteins, and pathological changes indicative of accumulation of misfolded and aggregated proteins in the affected cardiomyocytes. From TEM image analysis it is evident that the affected cardiomyocytes demonstrate significant difficulty in the disposal of damaged proteins and maintenance of proteostasis, which leads to pathological remodelling of the heart and contractile dysfunction. It appears that the underlying causes of accumulation of damaged proteins are associated with dysregulated auto phagosome and proteasome systems, which, in susceptible individuals, create a milieu conducive for the development of DCM and heart failure. RESEARCH HIGHLIGHTS The light and electron microscopy image analyses revealed degenerative changes and protein aggregates in the cardiomyocytes of chickens affected by DCM. The analyses of FTIR spectra of the myocardium revealed that DCM and heart pump failure in broiler chickens are associated with conformational changes of myocardial proteins. The morphological changes in cardiomyocytes and conformational changes in myocardial proteins architecture are integral constituents of pathophysiology of DCM in fast-growing broiler chickens.

Potential of milk mid-infrared spectra to predict nitrogen use efficiency of individual dairy cows in early lactation.

Improving nitrogen use efficiency (NUE) at both the individual cow and the herd level has become a key target in dairy production systems, for both environmental and economic reasons. Cost-effective and large-scale phenotyping methods are required to improve NUE through genetic selection and by feeding and management strategies. The aim of this study was to evaluate the possibility of using mid-infrared (MIR) spectra of milk to predict individual dairy cow NUE during early lactation. Data were collected from 129 Holstein cows, from calving until 50 d in milk, in 3 research herds (Denmark, Ireland, and the UK). In 2 of the herds, diets were designed to challenge cows metabolically, whereas a diet reflecting local management practices was offered in the third herd. Nitrogen intake (kg/d) and nitrogen excreted in milk (kg/d) were calculated daily. Nitrogen use efficiency was calculated as the ratio between nitrogen in milk and nitrogen intake, and expressed as a percentage. Individual daily values for NUE ranged from 9.7 to 81.7%, with an average of 36.9% and standard deviation of 10.4%. Milk MIR spectra were recorded twice weekly and were standardized into a common format to avoid bias between apparatus or sampling periods. Regression models predicting NUE using milk MIR spectra were developed on 1,034 observations using partial least squares or support vector machines regression methods. The models were then evaluated through (1) a cross-validation using 10 subsets, (2) a cow validation excluding 25% of the cows to be used as a validation set, and (3) a diet validation excluding each of the diets one by one to be used as validation sets. The best statistical performances were obtained when using the support vector machines method. Inclusion of milk yield and lactation number as predictors, in combination with the spectra, also improved the calibration. In cross-validation, the best model predicted NUE with a coefficient of determination of cross-validation of 0.74 and a relative error of 14%, which is suitable to discriminate between low- and high-NUE cows. When performing the cow validation, the relative error remained at 14%, and during the diet validation the relative error ranged from 12 to 34%. In the diet validation, the models showed a lack of robustness, demonstrating difficulties in predicting NUE for diets and for samples that were not represented in the calibration data set. Hence, a need exists to integrate more data in the models to cover a maximum of variability regarding breeds, diets, lactation stages, management practices, seasons, MIR instruments, and geographic regions. Although the model needs to be validated and improved for use in routine conditions, these preliminary results showed that it was possible to obtain information on NUE through milk MIR spectra. This could potentially allow large-scale predictions to aid both further genetic and genomic studies, and the development of farm management tools.

Genetic analysis of orotic acid predicted with Fourier transform infrared milk spectra.

Fourier transform infrared spectral analysis is a cheap and fast method to predict milk composition. A not very well studied milk component is orotic acid. Orotic acid is an intermediate in the biosynthesis pathway of pyrimidine nucleotides and is an indicator for the metabolic cattle disorder deficiency of uridine monophosphate synthase. The function of orotic acid in milk and its effect on calf health, health of humans consuming milk or milk products, manufacturing properties of milk, and its potential as an indicator trait are largely unknown. The aims of this study were to determine if milk orotic acid can be predicted from infrared milk spectra and to perform a large-scale phenotypic and genetic analysis of infrared-predicted milk orotic acid. An infrared prediction model for orotic acid was built using a training population of 292 Danish Holstein and 299 Danish Jersey cows, and a validation population of 381 Danish Holstein cows. Milk orotic acid concentration was determined with nuclear magnetic resonance spectroscopy. For genetic analysis of infrared orotic acid, 3 study populations were used: 3,210 Danish Holstein cows, 3,360 Danish Jersey cows, and 1,349 Dutch Holstein Friesian cows. Using partial least square regression, a prediction model for orotic acid was built with 18 latent variables. The error of the prediction for the infrared model varied from 1.0 to 3.2 mg/L, and the accuracy varied from 0.68 to 0.86. Heritability of infrared orotic acid predicted with the standardized prediction model was 0.18 for Danish Holstein, 0.09 for Danish Jersey, and 0.37 for Dutch Holstein Friesian. We conclude that milk orotic acid can be predicted with moderate to good accuracy based on infrared milk spectra and that infrared-predicted orotic acid is heritable. The availability of a cheap and fast method to predict milk orotic acid opens up possibilities to study the largely unknown functions of milk orotic acid.

Symposium review: Big data, big predictions: Utilizing milk Fourier-transform infrared and genomics to improve hyperketonemia management.

Negative animal health and performance outcomes are associated with disease incidences that can be labor-intensive, costly, and cumbersome for many farms. Amelioration of unfavorable outcomes through early detection and treatment of disease has emphasized the value of improving health monitoring. Although the value is recognized, detecting hyperketonemia (HYK) is still difficult for many farms to do practically and efficiently. Increasing data streams available to farms presents opportunities to use data to better monitor cow and herd health; however, challenges remain with regard to validating, integrating, and interpreting data. During the transition to lactation period, useful data are presented in the form of milk production and composition, milk Fourier-transform infrared (FTIR) wavelength absorbance, cow management records, and genomics, which have been employed to monitor postpartum onset of HYK. Attempts to predict postpartum HYK from test-day milk and performance variables incorporated into multiple linear regression models have demonstrated sufficient accuracy to monitor monthly herd prevalence; however, they lacked the sensitivity and specificity for individual cow diagnostics. Subsequent artificial neural network prediction models employing FTIR data and milk composition variables achieved 83 and 81% sensitivity and specificity for individual cow diagnostics. Although these results fail to reach the diagnostic goals of 90%, they are achieved without individual cow blood samples, which may justify acceptance of lower performance. The caveat is that these models depend on milk analysis, which is traditionally done every 4 weeks. This infrequent sampling allows for a single diagnostic sample for about half of the fresh cows. Benefits to farms are greatly improved if postpartum cows can be milk-tested weekly. Additionally, this allows for close monitoring of somatic cell count and may open the door for use of other herd health monitoring tools. Future improvements in these models may be achievable by maximizing sensitivity at the expense of specificity and may be most economical in disorders for which the cost of treatment is less than that of mistreating (e.g., HYK). Genomic predictions for HYK may be improved by incorporating genome-wide associated SNP and further utilized for precision management of HYK risk groups. Development and validation of HYK prediction models may provide producers with individual cow and herd-level management tools.