Mid-infrared prediction of lactoferrin content in bovine milk: potential indicator of mastitis

Genetic parameters for mid-infrared-spectroscopy-predicted mastitis phenotypes and related traits

Genetic improvement of udder health in dairy cows is of high relevance as mastitis is one of the most prevalent diseases. Since it is known that the heritability of mastitis is low and direct data on mastitis cases are often not available in large numbers, auxiliary traits, such as somatic cell count (SCC) are used for the genetic evaluation of udder health. In previous studies, models to predict clinical mastitis based on mid-infrared (MIR) spectral data and a somatic cell count-derived score (SCS) were developed. Those models can provide a probability of mastitis for each cow at every test-day, which is potentially useful as an additional auxiliary trait for the genetic evaluation of udder health. Furthermore, MIR spectral data were used to estimate contents of lactoferrin, a glycoprotein positively associated with immune response. The present study aimed to estimate heritabilities (h2) and genetic correlations (ra) for clinical mastitis diagnosis (CM), SCS, MIR-predicted mastitis probability (MIRprob), MIR + SCS-predicted mastitis probability (MIRSCSprob) and lactoferrin estimates (LF). Data for this study were collected within the routine milk recording and health monitoring system of Austria from 2014 to 2021 and included records of approximately 54,000 Fleckvieh cows. Analyses were performed in two datasets, including test-day records from 5 to 150 or 5 to 305 days in milk. Prediction models were applied to obtain MIR- and SCS-based phenotypes (MIRprob, MIRSCSprob, LF). To estimate heritabilities and genetic correlations bivariate linear animal models were applied for all traits. A lactation model was used for CM, defined as a binary trait, and a test-day model for all other continuous traits. In addition to the random animal genetic effect, the fixed effects year-season of calving and parity-age at calving and the random permanent environmental effect were considered in all models. For CM the random herd-year effect, for continuous traits the random herd-test day effect and the covariate days in milk (linear and quadratic) were additionally fitted. The obtained genetic parameters were similar in both datasets. The heritability found for CM was expectedly low (h2 = 0.02). For SCS and MIRSCSprob, heritability estimates ranged from 0.23 to 0.25, and for MIRprob and LF from 0.15 to 0.17. CM was highly correlated with SCS and MIRSCSprob (ra = 0.85 to 0.88). Genetic correlations of CM were moderate with MIRprob (ra = 0.26 and 0.37) during 150 and 305 days in milk, respectively and low with LF (h2 = 0.10 and 0.11). However, basic selection index calculations indicate that the added value of the new MIR-predicted phenotypes is limited for genetic evaluation of udder health.

The Genetic Characteristics of FT-MIRS-Predicted Milk Fatty Acids in Chinese Holstein Cows

Fourier Transform Mid-Infrared Spectroscopy (FT-MIRS) can be used for quantitative detection of milk components. Here, milk samples of 458 Chinese Holstein cows from 11 provinces in China were collected and we established a total of 22 quantitative prediction models in milk fatty acids by FT-MIRS. The coefficient of determination of the validation set ranged from 0.59 (C18:0) to 0.76 (C4:0). The models were adopted to predict the milk fatty acids from 2138 cows and a new high-throughput computing software HiBLUP was employed to construct a multi-trait model to estimate and analyze genetic parameters in dairy cows. Finally, genome-wide association analysis was performed and seven novel SNPs significantly associated with fatty acid content were selected, investigated, and verified with the FarmCPU method, which stands for "Fixed and random model Circulating Probability Unification". The findings of this study lay a foundation and offer technical support for the study of fatty acid trait breeding and the screening and grouping of characteristic dairy cows in China with rich, high-quality fatty acids. It is hoped that in the future, the method established in this study will be able to screen milk sources rich in high-quality fatty acids.

Genetic analysis of milk minerals in dairy cattle: a review

Mineral composition in milk can affect its nutritional as well as physio-chemical properties of milk and is considered an important trait in the evaluation of milk quality. The composition and concentration of milk minerals could be altered with natural source of variation including nutrition and genetics. The effect of diet on milk minerals is well studied. However, genetic effects on the milk minerals have recently gained the attention. This review provides an overview of the genetic variation of milk minerals, and the genomic regions associated with mineral concentration in the milk are also discussed. The difference of milk minerals between breeds and the genetic parameters including heritability estimates and correlation among minerals indicates that milk minerals are under strong genetic control. Recently, the genome-wide association study (GWAS) has explored several regions associated with milk minerals and thus provides a new genetic source for improving the milk quality through genomics-assisted breeding. Hence, a combination of the qualitative and molecular approaches can be exploited to improving the nutritional quality of cattle milk in terms of its mineral composition.

Genetic aspects of immunoglobulins and cyclophilin A in milk as potential indicators of mastitis resistance in Holstein cows

Mastitis is one of the most frequent and costly diseases affecting dairy cattle. Natural antibodies (immunoglobulins) and cyclophilin A (CyPA), the most abundant member of the family of peptidyl prolyl cis/trans isomerases, in milk may serve as indicators of mastitis resistance in dairy cattle. However, genetic information for CyPA is not available, and knowledge on the genetic and nongenetic relationships between these immune-related traits and somatic cell score (SCS) and milk yield in dairy cattle is sparse. Therefore, we aimed to comprehensively evaluate whether immune-related traits consisting of 5 Ig classes (IgG, IgG1, IgG2, IgA, and IgM) and CyPA in the test-day milk of Holstein cows can be used as genetic indicators of mastitis resistance by evaluating the genetic and nongenetic relationships with SCS in milk. The nongenetic factors affecting immune-related traits and the effects of these traits on SCS were evaluated. Furthermore, the genetic parameters of immune-related traits according to health status and genetic relationships under different SCS environments were estimated. All immune-related traits were significantly associated with SCS and directly proportional. Additionally, evaluation using a classification tree revealed that IgA, IgG2, and IgG were associated with SCS levels. Genetic factor analyses indicated that heritability estimates were low for CyPA (0.08) but moderate for IgG (0.37), IgA (0.44), and IgM (0.44), with positive genetic correlations among Ig (0.25-0.96). We also evaluated the differences in milk yield and SCS of cows between the low and high groups according to their sires' estimated breeding value for immune-related traits. In the high group, IgA had a significantly lower SCS in milk at 7 to 30 d compared with that in the low group. Furthermore, the Ig in milk had high positive genetic correlations between healthy and infected conditions (0.82-0.99), suggesting that Ig in milk under healthy conditions could interact with those under infected conditions, owing to the genetic ability based on the level of Ig in milk. Thus, Ig in milk are potential indicators for the genetic selection of mastitis resistance. However, because only the relationship between immune-related traits and SCS was investigated in this study, further study on the relationship between clinical mastitis and Ig in milk is needed before Ig can be used as an indicator of mastitis resistance.

Minerals and essential amino acids of bovine colostrum: Phenotypic variability and predictive ability of mid- and near-infrared spectroscopy

Colostrum quality and volume are fundamental for calves because it is the primary supplier of antibodies and the first source of energy, carbohydrates, lipids, proteins, minerals, and vitamins for the newborn. Assessing the detailed composition (i.e., AA and mineral content) of bovine colostrum (BC) on-line and at a reasonable cost would help dairy stakeholders such as farmers or veterinarians for precision feeding purposes and industries producing preparations containing BC such as foodstuff, supplements, and medicaments. In the present study we evaluated mid- (MIRS) and near-infrared spectroscopy (NIRS) prediction ability for AA and mineral composition of individual BC. Second, we the investigated the major factors affecting the phenotypic variability of such traits also evaluating the correlations with the Ig concentration. Results demonstrated that MIRS and NIRS were able to provide sufficiently accurate predictions for all the AA. The coefficient of determination in external validation (R2V) fell, in fact, within the range of 0.70 to 0.86, with the exception of Ile, His, and Met. Only some minerals reached a sufficient accuracy (i.e., Ca, P, S, and Mg; R2V ≥ 0.66) using MIRS, and also S (R2V = 0.87) using NIRS. Phenotypically, both parity and calving season affected the variability of these BC composition traits. Heifers' colostrum was the one with the greatest concentration of Ca and P, the 2 most abundant minerals. These minerals were however very low in cows calving in summer compared with the rest of the year. The pattern of AA across parities and calving season was not linear, likely because their variability was scarcely (or not) affected by these effects. Finally, samples characterized by high IgG concentration were those presenting on average greater concentration of AA. Findings suggest that infrared spectroscopy has the potential to be used to predict certain AA and minerals, outlining the possibility of implementing on-site analyses for the evaluation of the broad-sense BC quality.

Mid-infrared spectroscopy for large-scale phenotyping of bovine colostrum gross composition and immunoglobulin concentration

Immunoglobulin G is the fundamental antibody for acquisition of passive transfer of immunity in ruminant newborns. Colostrum, in fact, must be administered as soon as possible after birth to ensure a successful transfer of IgG from the dam to the calf. Assessment of colostrum Ig concentration and gross composition via gold standards is expensive, time consuming, and hardly implementable for large-scale investigations. Therefore, in the present study we evaluated the predictive ability of mid-infrared spectroscopy (MIRS) as an indirect determination method. A total of 714 colostrum samples collected within 6 h from parturition from Italian Holstein cows, 30% primiparous and 70% pluriparous, were scanned using a benchtop spectrometer after dilution in pure water. The prediction models were developed by correlating spectral information with the reference measurements: IgG concentration (93.54 ± 33.87 g/L), total Ig concentrations (102.82 ± 35.04 g/L), and content of protein (14.71 ± 3.51%), fat (4.61 ± 3.04%), and lactose (2.36 ± 0.51 mg/100 mg). We found a good to excellent performance in prediction of colostrum IgG concentration and traditional composition traits in cross-validation (R2CV ≥ 0.92) and a promising and good predictive ability in external validation with R2V equal to 0.84, 0.89, and 0.74 for IgG, protein, and fat, respectively. In the case of IgG and protein content, for example, the coefficient of determination in external validation was greater than 0.84. The other Ig fractions, A and M, presented insufficient prediction accuracy likely due to their extremely low concentration compared with IgG (4.56 and 5.06 g/L vs. 93.54 g/L). The discriminant ability of MIRS-predicted IgG and protein content was outstanding when trying to classify samples according to the quality level (i.e., low vs. high concentration of IgG). In particular, the cut-off that better discriminate low- from high-quality colostrum was 75.40 g/L in the case of the MIRS-predicted IgG and 13.32% for the MIRS-predicted protein content. Therefore, MIRS is proposed as a rapid and cheap tool for large-scale punctual IgG, protein, and lactose quantification and for the screening of low-quality samples. From a practical perspective, there is the possibility to install colostrum models in the MIRS benchtop machineries already present in laboratories in charge of official milk testing. Colostrum phenotypes collected on an individual basis will be useful to breeders for the definition of specific selection strategies and to farmers for management scopes. Finally, our findings may be relevant for other stakeholders, given the fact that colostrum is an emerging ingredient for the animal and human food and pharmaceutical industry.

Recent Advances in the Determination of Milk Adulterants and Contaminants by Mid-Infrared Spectroscopy

The presence of chemical contaminants, toxins, or veterinary drugs in milk, as well as the adulteration of milk from different species, has driven the development of new tools to ensure safety and quality. Several analytical procedures have been proposed for the rapid screening of hazardous substances or the selective confirmation of the authenticity of milk. Mid-infrared spectroscopy and Fourier-transform infrared have been two of the most relevant technologies conventionally employed in the dairy industry. These fingerprint methodologies can be very powerful in determining the trait of raw material without knowing the identity of each constituent, and several aspects suggest their potential as a screening method to detect adulteration. This paper reviews the latest advances in applying mid-infrared spectroscopy for the detection and quantification of adulterants, milk dilution, the presence of pathogenic bacteria, veterinary drugs, and hazardous substances in milk.

Effectiveness of mid-infrared spectroscopy for the prediction of cow milk metabolites

Proton nuclear magnetic resonance (¹H NMR) spectroscopy is acknowledged as one of the most powerful analytical methods with cross-cutting applications in dairy foods. To date, the use of ¹H NMR spectroscopy for the collection of milk metabolic profile is hindered by costly and time-consuming sample preparation and analysis. The present study aimed at evaluating the accuracy of mid-infrared spectroscopy (MIRS) as a rapid method for the prediction of cow milk metabolites determined through ¹H NMR spectroscopy. Bulk milk (n = 72) and individual milk samples (n = 482) were analyzed through one-dimensional ¹H NMR spectroscopy and MIRS. Nuclear magnetic resonance spectroscopy identified 35 milk metabolites, which were quantified in terms of relative abundance, and MIRS prediction models were developed on the same 35 milk metabolites, using partial least squares regression analysis. The best MIRS prediction models were developed for galactose-1-phosphate, glycerophosphocholine, orotate, choline, galactose, lecithin, glutamate, and lactose, with coefficient of determination in external validation from 0.58 to 0.85, and ratio of performance to deviation in external validation from 1.50 to 2.64. The remaining 27 metabolites were poorly predicted. This study represents a first attempt to predict milk metabolome. Further research is needed to specifically address whether developed prediction models may find practical application in the dairy sector, with particular regard to the screening of dairy cows' metabolic status, the quality control of dairy foods, and the identification of processed milk or incorrectly stored milk.

Candidate genes for mastitis resistance in dairy cattle: a data integration approach

BACKGROUND

Inflammation of the mammary tissue (mastitis) is one of the most detrimental health conditions in dairy ruminants and is considered the most economically important infectious disease of the dairy sector. Improving mastitis resistance is becoming an important goal in dairy ruminant breeding programmes. However, mastitis resistance is a complex trait and identification of mastitis-associated alleles in livestock is difficult. Currently, the only applicable approach to identify candidate loci for complex traits in large farm animals is to combine different information that supports the functionality of the identified genomic regions with respect to a complex trait.

METHODS

To identify the most promising candidate loci for mastitis resistance we integrated heterogeneous data from multiple sources and compiled the information into a comprehensive database of mastitis-associated candidate loci. Mastitis-associated candidate genes reported in association, expression, and mouse model studies were collected by searching the relevant literature and databases. The collected data were integrated into a single database, screened for overlaps, and used for gene set enrichment analysis.

RESULTS

The database contains candidate genes from association and expression studies and relevant transgenic mouse models. The 2448 collected candidate loci are evenly distributed across bovine chromosomes. Data integration and analysis revealed overlaps between different studies and/or with mastitis-associated QTL, revealing promising candidate genes for mastitis resistance.

CONCLUSION

Mastitis resistance is a complex trait influenced by numerous alleles. Based on the number of independent studies, we were able to prioritise candidate genes and propose a list of the 22 most promising. To our knowledge this is the most comprehensive database of mastitis associated candidate genes and could be helpful in selecting genes for functional validation studies.

Providing New Insights on the Molecular Properties and Thermal Stability of Ovotransferrin and Lactoferrin

Ovotransferrin (OVT) is a multi-functional protein showing over 50% homology with Bovine lactoferrin (BLF) and human lactoferrin (HLF), which have the potential to be a substitute for lactoferrin (LF) due to the limited production of LF. To explore the substitutability of OVT, the molecular properties and thermal stability of OVT, BLF and HLF were characterized because these properties will affect the processing quality and biological activities of protein products when exposed to different processing conditions (e.g., temperature, pH, ion strength). The results showed that although obviously different isoelectric point (5.31, 9.12 and 8.75 for OVT, BLF and HLF, respectively), particle size distribution and hydrophobicity were found, they exhibited good dispersity because of high potential value. They showed an endothermic peak at 80.64 °C, 65.71 °C and 90.01 °C, respectively, and the denaturation temperature varied at different pH and ionic strength. OVT and BLF were more susceptible to heating at pH 5.0 as reflected by the decline of denaturation temperature (21.78 °C shift for OVT and 5.81 °C shift for BLF), while HLF could remain stable. Compared with BLF, OVT showed higher secondary structure stability at pH 7.0 and 9.0 with heating. For example, the α-helix content of OVT changed from 20.35% to 15.4% at pH 7.0 after heating, while that of BLF changed from 20.05% to 6.65%. The increase on fluorescence intensity and redshifts on the maximum wavelength after heating indicated the changes of tertiary structure of them. The turbidity measurements showed that the thermal aggregation degree of OVT was lower than BLF and HLF at pH 7.0 (30.98%, 59.53% and 35.66%, respectively) and pH 9.0 (4.83%, 12.80% and 39.87%, respectively). This work demonstrated the similar molecular properties and comparable thermal stability of OVT to BLF and HLF, which can offer a useful reference for the substitute of LF by OVT.

Predictions of Milk Fatty Acid Contents by Mid-Infrared Spectroscopy in Chinese Holstein Cows

Genetic improvement of milk fatty acid content traits in dairy cattle is of great significance. However, chromatography-based methods to measure milk fatty acid content have several disadvantages. Thus, quick and accurate predictions of various milk fatty acid contents based on the mid-infrared spectrum (MIRS) from dairy herd improvement (DHI) data are essential and meaningful to expand the amount of phenotypic data available. In this study, 24 kinds of milk fatty acid concentrations were measured from the milk samples of 336 Holstein cows in Shandong Province, China, using the gas chromatography (GC) technique, which simultaneously produced MIRS values for the prediction of fatty acids. After quantification by the GC technique, milk fatty acid contents expressed as g/100 g of milk (milk-basis) and g/100 g of fat (fat-basis) were processed by five spectral pre-processing algorithms: first-order derivative (DER1), second-order derivative (DER2), multiple scattering correction (MSC), standard normal transform (SNV), and Savitzky-Golsy convolution smoothing (SG), and four regression models: random forest regression (RFR), partial least square regression (PLSR), least absolute shrinkage and selection operator regression (LassoR), and ridge regression (RidgeR). Two ranges of wavebands (4000~400 cm^-1 and 3017~2823 cm^-1/1805~1734 cm^-1) were also used in the above analysis. The prediction accuracy was evaluated using a 10-fold cross validation procedure, with the ratio of the training set and the test set as 3:1, where the determination coefficient (R²) and residual predictive deviation (RPD) were used for evaluations. The results showed that 17 out of 31 milk fatty acids were accurately predicted using MIRS, with RPD values higher than 2 and R² values higher than 0.75. In addition, 16 out of 31 fatty acids were accurately predicted by RFR, indicating that the ensemble learning model potentially resulted in a higher prediction accuracy. Meanwhile, DER1, DER2 and SG pre-processing algorithms led to high prediction accuracy for most fatty acids. In summary, these results imply that the application of MIRS to predict the fatty acid contents of milk is feasible.

Immune mechanisms, resistance genes, and their roles in the prevention of mastitis in dairy cows

Mastitis is one of the most important diseases of the mammary gland. The increased incidence of this disease in cows is due to the breeding of dairy cattle for higher yields, which is accompanied by an increased susceptibility to mastitis. Therefore, the difficulty involved with preventing this disease has increased. An integral part of current research is the elimination of mastitis in order to reduce the consumption of antibiotic drugs, thereby reducing the resistance of microorganisms and decreasing companies' economic losses due to mastitis (i.e. decreased milk yield, increased drug costs, and reduced milk supply). Susceptibility to mastitis is based on dairy cows' immunity, health, nutrition, and welfare. Thus, it is important to understand the immune processes in the body in order to increase the resistance of animals. Recently, various studies have focused on the selection of mastitis resistance genes. An important point is also the prevention of mastitis. This publication aims to describe the physiology of the mammary gland along with its immune mechanisms and to approximate their connection with potential mastitis resistance genes. This work describes various options for mastitis elimination and focuses on genetic selection and a closer specification of resistance genes to mastitis. Among the most promising resistance genes for mastitis, we consider CD14, CXCR1, lactoferrin, and lactoglobulin.

Prediction of Acute and Chronic Mastitis in Dairy Cows Based on Somatic Cell Score and Mid-Infrared Spectroscopy of Milk

Simple Summary

Mid-infrared (MIR) spectroscopy is the method of choice to determine milk components like fat, protein and urea. We examined the potential of MIR spectra analyses for the prediction of clinical mastitis events of dairy cows additionally, or alternatively, to somatic cell count, which is routinely used as an indicator for mastitis monitoring. Prediction models based on MIR spectra and a somatic cell count-derived score (SCS) were developed and compared. A model based on MIR spectra and SCS proved more accurate at predicting mastitis than models based on either indicator alone. Consequently, MIR spectra analyses add extra value in the prediction of clinical mastitis, making them potentially useful for dairy farm management and as an auxiliary trait for the genetic evaluation of udder health.

Abstract

Monitoring for mastitis on dairy farms is of particular importance, as it is one of the most prevalent bovine diseases. A commonly used indicator for mastitis monitoring is somatic cell count. A supplementary tool to predict mastitis risk may be mid-infrared (MIR) spectroscopy of milk. Because bovine health status can affect milk composition, this technique is already routinely used to determine standard milk components. The aim of the present study was to compare the performance of models to predict clinical mastitis based on MIR spectral data and/or somatic cell count score (SCS), and to explore differences of prediction accuracies for acute and chronic clinical mastitis diagnoses. Test-day data of the routine Austrian milk recording system and diagnosis data of its health monitoring, from 59,002 cows of the breeds Fleckvieh (dual purpose Simmental), Holstein Friesian and Brown Swiss, were used. Test-day records within 21 days before and 21 days after a mastitis diagnosis were defined as mastitis cases. Three different models (MIR, SCS, MIR + SCS) were compared, applying Partial Least Squares Discriminant Analysis. Results of external validation in the overall time window (−/+21 days) showed area under receiver operating characteristic curves (AUC) of 0.70 when based only on MIR, 0.72 when based only on SCS, and 0.76 when based on both. Considering as mastitis cases only the test-day records within 7 days after mastitis diagnosis, the corresponding areas under the curve were 0.77, 0.83 and 0.85. Hence, the model combining MIR spectral data and SCS was performing best. Mastitis probabilities derived from the prediction models are potentially valuable for routine mastitis monitoring for farmers, as well as for the genetic evaluation of the trait udder health.

Multiple Breeds and Countries' Predictions of Mineral Contents in Milk from Milk Mid-Infrared Spectrometry

Measuring the mineral composition of milk is of major interest in the dairy sector. This study aims to develop and validate robust multi-breed and multi-country models predicting the major minerals through milk mid-infrared spectrometry using partial least square regressions. A total of 1281 samples coming from five countries were analyzed to obtain spectra and in ICP-AES to measure the mineral reference contents. Models were built from records coming from four countries (n = 1181) and validated using records from the fifth country, Austria (n = 100). The importance of including local samples was tested by integrating 30 Austrian samples in the model while validating with the remaining 70 samples. The best performances were achieved using this second set of models, confirming the need to cover the spectral variability of a country before making a prediction. Validation root mean square errors were 54.56, 63.60, 7.30, 59.87, and 152.89 mg/kg for Na, Ca, Mg, P, and K, respectively. The built models were applied on the Walloon milk recording large-scale spectral database, including 3,510,077. The large-scale predictions on this dairy herd improvement database provide new insight regarding the minerals' variability in the population, as well as the effect of parity, stage of lactation, breeds, and seasons.

Invited review: Toward a common language in data-driven mastitis detection research

Sensor technologies for mastitis detection have resulted in the collection and availability of a large amount of data. As a result, scientific publications reporting mastitis detection research have become less driven by approaches based on biological assumptions and more by data-driven modeling. Most of these approaches try to predict mastitis events from (combinations of) raw sensor data to which a wide variety of methods are applied originating from machine learning and classical statistical approaches. However, an even wider variety in terminologies is used by researchers for methods that are similar in nature. This makes it difficult for readers from other disciplines to understand the specific methods that are used and how these differ from each other. The aim of this paper was to provide a framework (filtering, transformation, and classification) for describing the different methods applied in sensor data-based clinical mastitis detection research and use this framework to review and categorize the approaches and underlying methods described in the scientific literature on mastitis detection. We identified 40 scientific publications between 1992 and 2020 that applied methods to detect clinical mastitis from sensor data. Based on these publications, we developed and used the framework and categorized these scientific publications into the 2 data processing techniques of filtering and transformation. These data processing techniques make raw data more amendable to be used for the third step in our framework, that of classification, which is used to distinguish between healthy and nonhealthy (mastitis) cows. Most publications (n = 34) used filtering or transformation, or a combination of these 2, for data processing before classification, whereas the remaining publications (n = 6) classified the observations directly from raw data. Concerning classification, applying a simple threshold was the most used method (n = 19 publications). Our work identified that within approaches several different methods and terminologies for similar methods were used. Not all publications provided a clear description of the method used, and therefore it seemed that different methods were used between publications, whereas in fact just a different terminology was used, or the other way around. This paper is intended to serve as a reference for people from various research disciplines who need to collaborate and communicate efficiently about the topic of sensor-based mastitis detection and the methods used in this context. The framework used in this paper can support future research to correctly classify approaches and methods, which can improve the understanding of scientific publication. We encourage future research on sensor-based animal disease detection, including that of mastitis detection, to use a more coherent terminology for methods, and clearly state which technique (e.g., filtering) and approach (e.g., moving average) are used. This paper, therefore, can serve as a starting point and further stimulates the interdisciplinary cooperation in sensor-based mastitis research.

Improving robustness and accuracy of predicted daily methane emissions of dairy cows using milk mid-infrared spectra

BACKGROUND

A robust proxy for estimating methane (CH₄ ) emissions of individual dairy cows would be valuable especially for selective breeding. This study aimed to improve the robustness and accuracy of prediction models that estimate daily CH₄ emissions from milk Fourier transform mid-infrared (FT-MIR) spectra by (i) increasing the reference dataset and (ii) adjusting for routinely recorded phenotypic information. Prediction equations for CH₄ were developed using a combined dataset including daily CH₄ measurements (n = 1089; g d^-1 ) collected using the SF₆ tracer technique (n = 513) and measurements using respiration chambers (RC, n = 576). Furthermore, in addition to the milk FT-MIR spectra, the variables of milk yield (MY) on the test day, parity (P) and breed (B) of cows were included in the regression analysis as explanatory variables.

RESULTS

Models developed based on a combined RC and SF₆ dataset predicted the expected pattern in CH₄ values (in g d^-1 ) during a lactation cycle, namely an increase during the first weeks after calving followed by a gradual decrease until the end of lactation. The model including MY, P and B information provided the best prediction results (cross-validation statistics: R² = 0.68 and standard error = 57 g CH₄ d^-1 ).

CONCLUSIONS

The models developed accounted for more of the observed variability in CH₄ emissions than previously developed models and thus were considered more robust. This approach is suitable for large-scale studies (e.g. animal genetic evaluation) where robustness is paramount for accurate predictions across a range of animal conditions. © 2020 Society of Chemical Industry.

Influence of Estrus on the Milk Characteristics and Mid-Infrared Spectra of Dairy Cows

Simple Summary

Some studies have confirmed the variation in milk profiles when dairy cows show estrus. However, only a few milk components, such as fat, protein, and lactose, have been investigated so far, and thus any changes in the many other parts of milk’s composition due to estrus are unknown. Milk mid-infrared (MIR) spectra consist of wavenumbers, which provide insight into the chemical composition of milk. The MIR spectrum reflects the global composition of milk, but this information is currently underused. In this study, we considered MIR wavenumbers as traits, and directly studied the spectral information as a way to study the estrus of dairy cows linked to milk composition. This research provides a deeper understanding of the milk MIR spectrum and may lead to new approaches for estrus detection in dairy cows from routine milk analysis, thereby guiding an opportune insemination time.

Abstract

Milk produced by dairy cows is a complex combination of many components. However, at present, changes in only a few milk components (e.g., fat, protein, and lactose) during the estrus cycle in dairy cows have been documented. Mid-infrared (MIR) spectroscopy is a worldwide method routinely used for milk analysis, as MIR spectra reflect the global composition of milk. Therefore, this study aimed to investigate the changes in milk MIR spectra and milk production traits (fat, protein, lactose, urea, total solids (TS), and solid not fat (SnF)) due to estrus. Cows that were successfully inseminated, leading to conception, were included. Cows confirmed to be pregnant were considered to be in estrus at the day of insemination (day 0). A general linear mixed model, which included the random effect of cows, the fixed classification effects of parity number, days in relation to estrus, as well as the interaction between parity number and days in relation to estrus, was applied to investigate the changes in milk production traits and 1060 milk infrared wavenumbers, ranging from 925 to 5011 cm⁻¹, of 371 records from 162 Holstein cows on the days before (day −3, day −2, and day −1) and on the day of estrus (day 0). The days in relation to estrus had a significant effect on fat, protein, urea, TS, and SnF, whose contents increased from day −3 to day 0. Lactose did not seem to be significantly influenced by the occurrence of estrus. The days in relation to estrus had significant effects on the majority of the wavenumbers. Besides, we found that some of the wavenumbers in the water absorption regions were significantly changed on the days before and on the day of estrus. This suggests that these wavenumbers may contain useful information. In conclusion, the changes in the milk composition due to estrus can be observed through the analysis of the milk MIR spectrum. Further analyses are warranted to more deeply explore the potential use of milk MIR spectra in the detection of estrus.

Fourier Transform Infrared Spectroscopy as a Tool to Study Milk Composition Changes in Dairy Cows Attributed to Housing Modifications to Improve Animal Welfare

Animal welfare status is assessed today through visual evaluations requiring an on-farm visit. A convenient alternative would be to detect cow welfare status directly in milk samples, already routinely collected for milk recording. The objective of this study was to propose a novel approach to demonstrate that Fourier transform infrared (FTIR) spectroscopy can detect changes in milk composition related to cows subjected to movement restriction at the tie stall with four tie-rail configurations varying in height and position (TR1, TR2, TR3 and TR4). Milk mid-infrared spectra were collected on weekly basis. Long-term average spectra were calculated for each cow using spectra collected in weeks 8-10 of treatment. Principal component analysis was applied to spectral averages and the scores of principal components (PCs) were tested for treatment effect by mixed modelling. PC7 revealed a significant treatment effect (p = 0.01), particularly for TR3 (configuration with restricted movement) vs. TR1 (recommended configuration) (p = 0.03). The loading spectrum of PC7 revealed high loadings at wavenumbers that could be assigned to biomarkers related to negative energy balance, such as β-hydroxybutyrate, citrate and acetone. This observation suggests that TR3 might have been restrictive for cows to access feed. Milk FTIR spectroscopy showed promising results in detecting welfare status and housing conditions in dairy cows.

Development of a subacute ruminal acidosis risk score and its prediction using milk mid-infrared spectra in early-lactation cows

A routine monitoring for subacute ruminal acidosis (SARA) on the individual level could support the minimization of economic losses and the ensuring of animal welfare in dairy cows. The objectives of this study were (1) to develop a SARA risk score (SRS) by combining information from different data acquisition systems to generate an integrative indicator trait, (2) the investigation of associations of the SRS with feed analysis data, blood characteristics, performance data, and milk composition, including the fatty acid (FA) profile, (3) the development of a milk mid-infrared (MIR) spectra-based prediction equation for this novel reference trait SRS, and (4) its application to an external data set consisting of MIR data of test day records to investigate the association between the MIR-based predictions of the SRS and the milk FA profile. The primary data set, which was used for the objectives (1) to (3), consisted of data collected from 10 commercial farms with a total of 100 Holstein cows in early lactation. The data comprised barn climate parameters, pH and temperature logging from intrareticular measurement boluses, as well as jaw movement and locomotion behavior recordings of noseband-sensor halters and pedometers. Further sampling and data collection included feed samples, blood samples, milk performance, and milk samples, whereof the latter were used to get the milk MIR spectra and to estimate the main milk components, the milk FA profile, and the lactoferrin content. Because all measurements were characterized by different temporal resolutions, the data preparation consisted of an aggregation into values on a daily basis and merging it into one data set. For the development of the SRS, a total of 7 traits were selected, which were derived from measurements of pH and temperature in the reticulum, chewing behavior, and milk yield. After adjustment for fixed effects and standardization, these 7 traits were combined into the SRS using a linear combination and directional weights based on current knowledge derived from literature studies. The secondary data set was used for objective (4) and consisted of test day records of the entire herds, including performance data, milk MIR spectra and MIR-predicted FA. At farm level, it could be shown that diets with higher proportions of concentrated feed resulted in both lower daily mean pH and higher SRS values. On the individual level, an increased SRS could be associated with a modified FA profile (e.g., lower levels of short- and medium-chain FA, higher levels of C17:0, odd- and branched-chain FA). Furthermore, a milk MIR-based partial least squares regression model with a moderate predictability was established for the SRS. This work provides the basis for the development of routine SARA monitoring and demonstrates the high potential of milk composition-based assessment of the health status of lactating cows.

Predicting pregnancy status from mid-infrared spectroscopy in dairy cow milk using deep learning

Accurately identifying pregnancy status is imperative for a profitable dairy enterprise. Mid-infrared (MIR) spectroscopy is routinely used to determine fat and protein concentrations in milk samples. Mid-infrared spectra have successfully been used to predict other economically important traits, including fatty acid content, mineral content, body energy status, lactoferrin, feed intake, and methane emissions. Machine learning has been used in a variety of fields to find patterns in vast quantities of data. This study aims to use deep learning, a sub-branch of machine learning, to establish pregnancy status from routinely collected milk MIR spectral data. Milk spectral data were obtained from National Milk Records (Chippenham, UK), who collect large volumes of data continuously on a monthly basis. Two approaches were followed: using genetic algorithms for feature selection and network design (model 1), and transfer learning with a pretrained DenseNet model (model 2). Feature selection in model 1 showed that the number of wave points in MIR data could be reduced from 1,060 to 196 wave points. The trained model converged after 162 epochs with validation accuracy and loss of 0.89 and 0.18, respectively. Although the accuracy was sufficiently high, the loss (in terms of predicting only 2 labels) was considered too high and suggested that the model would not be robust enough to apply to industry. Model 2 was trained in 2 stages of 100 epochs each with spectral data converted to gray-scale images and resulted in accuracy and loss of 0.97 and 0.08, respectively. Inspection on inference data showed prediction sensitivity of 0.89, specificity of 0.86, and prediction accuracy of 0.88. Results indicate that milk MIR data contains features relating to pregnancy status and the underlying metabolic changes in dairy cows, and such features can be identified by means of deep learning. Prediction equations from trained models can be used to alert farmers of nonviable pregnancies as well as to verify conception dates.

A comparison of 4 different machine learning algorithms to predict lactoferrin content in bovine milk from mid-infrared spectra

Lactoferrin (LF) is a glycoprotein naturally present in milk. Its content varies throughout lactation, but also with mastitis; therefore it is a potential additional indicator of udder health beyond somatic cell count. Condequently, there is an interest in quantifying this biomolecule routinely. First prediction equations proposed in the literature to predict the content in milk using milk mid-infrared spectrometry were built using partial least square regression (PLSR) due to the limited size of the data set. Thanks to a large data set, the current study aimed to test 4 different machine learning algorithms using a large data set comprising 6,619 records collected across different herds, breeds, and countries. The first algorithm was a PLSR, as used in past investigations. The second and third algorithms used partial least square (PLS) factors combined with a linear and polynomial support vector regression (PLS + SVR). The fourth algorithm also used PLS factors, but included in an artificial neural network with 1 hidden layer (PLS + ANN). The training and validation sets comprised 5,541 and 836 records, respectively. Even if the calibration prediction performances were the best for PLS + polynomial SVR, their validation prediction performances were the worst. The 3 other algorithms had similar validation performances. Indeed, the validation root mean squared error (RMSE) ranged between 162.17 and 166.75 mg/L of milk. However, the lower standard deviation of cross-validation RMSE and the better normality of the residual distribution observed for PLS + ANN suggest that this modeling was more suitable to predict the LF content in milk from milk mid-infrared spectra (R²v = 0.60 and validation RMSE = 162.17 mg/L of milk). This PLS +ANN model was then applied to almost 6 million spectral records. The predicted LF showed the expected relationships with milk yield, somatic cell score, somatic cell count, and stage of lactation. The model tended to underestimate high LF values (higher than 600 mg/L of milk). However, if the prediction threshold was set to 500 mg/L, 82% of samples from the validation having a content of LF higher than 600 mg/L were detected. Future research should aim to increase the number of those extremely high LF records in the calibration set.

Predicting bovine tuberculosis status of dairy cows from mid-infrared spectral data of milk using deep learning

Bovine tuberculosis (bTB) is a zoonotic disease in cattle that is transmissible to humans, distributed worldwide, and considered endemic throughout much of England and Wales. Mid-infrared (MIR) analysis of milk is used routinely to predict fat and protein concentration, and is also a robust predictor of several other economically important traits including individual fatty acids and body energy. This study predicted bTB status of UK dairy cows using their MIR spectral profiles collected as part of routine milk recording. Bovine tuberculosis data were collected as part of the national bTB testing program for Scotland, England, and Wales; these data provided information from over 40,500 bTB herd breakdowns. Corresponding individual cow life-history data were also available and provided information on births, movements, and deaths of all cows in the study. Data relating to single intradermal comparative cervical tuberculin (SICCT) skin-test results, culture, slaughter status, and presence of lesions were combined to create a binary bTB phenotype labeled 0 to represent nonresponders (i.e., healthy cows) and 1 to represent responders (i.e., bTB-affected cows). Contemporaneous individual milk MIR spectral data were collected as part of monthly routine milk recording and matched to bTB status of individual animals on the single intradermal comparative cervical tuberculin test date (±15 d). Deep learning, a sub-branch of machine learning, was used to train artificial neural networks and develop a prediction pipeline for subsequent use in national herds as part of routine milk recording. Spectra were first converted to 53 × 20-pixel PNG images, then used to train a deep convolutional neural network. Deep convolutional neural networks resulted in a bTB prediction accuracy (i.e., the number of correct predictions divided by the total number of predictions) of 71% after training for 278 epochs. This was accompanied by both a low validation loss (0.71) and moderate sensitivity and specificity (0.79 and 0.65, respectively). To balance data in each class, additional training data were synthesized using the synthetic minority over sampling technique. Accuracy was further increased to 95% (after 295 epochs), with corresponding validation loss minimized (0.26), when synthesized data were included during training of the network. Sensitivity and specificity also saw a 1.22- and 1.45-fold increase to 0.96 and 0.94, respectively, when synthesized data were included during training. We believe this study to be the first of its kind to predict bTB status from milk MIR spectral data. We also believe it to be the first study to use milk MIR spectral data to predict a disease phenotype, and posit that the automated prediction of bTB status at routine milk recording could provide farmers with a robust tool that enables them to make early management decisions on potential reactor cows, and thus help slow the spread of bTB.

Large-scale phenotyping in dairy sector using milk MIR spectra: Key factors affecting the quality of predictions

Methods and technologies enabling the estimation at large scale of important traits for the dairy sector are of great interest. Those phenotypes are necessary to improve herd management, animal genetic evaluation, and milk quality control. In the recent years, the research was very active to predict new phenotypes from the mid-infrared (MIR) analysis of milk. Models were developed to predict phenotypes such as fine milk composition, milk technological properties or traits related to cow health, fertility and environmental impact. Most of models were developed within research contexts and often not designed for routine use. The implementation of models at a large scale to predict new traits of interest brings new challenges as the factors influencing the robustness of models are poorly documented. The first objective of this work is to highlight the impact on prediction accuracy of factors such as the variability of the spectral and reference data, the spectral regions used and the complexity of models. The second objective is to emphasize methods and indicators to evaluate the quality of models and the quality of predictions generated under routine conditions. The last objective is to outline the issues and the solutions linked with the use and transfer of models on large number of instruments. Based on partial least square regression and 10 datasets including milk MIR spectra and reference quantitative values for 57 traits of interest, the impact of the different factors is illustrated by evaluating the influence on the validation root mean square error of prediction (RMSEP). In the displayed examples, all factors, when well set up, increase the quality of predictions, with an improvement of the RMSEP ranging from 12% to 43%. This work also aims to underline the need for and the complementarity between different validation procedures, statistical parameters and quality assurance methods. Finally, when using and transferring models, the impact of the spectral standardization on the prediction reproducibility is highlighted with an improvement up to 86% with the tested models, and the monitoring of individual spectrometer stability over time appears essential. This list inspired from our experience is of course not exhaustive. The displayed results are only examples and not general rules and other aspects play a role in the quality of final predictions. However, this work highlights good practices, methods and indicators to increase and evaluate quality of phenotypes predicted at a large scale. The results obtained argue for the development of guidelines at international levels, as well as international collaborations in order to constitute large and robust datasets and enable the use of models in routine conditions.

Phenotypic investigation of fine milk components in bovine milk and their prediction using mid-infrared spectroscopy

This study aimed to examine the phenotypic variation observed in fine milk components, the use of mid-infrared (MIR) spectroscopy to predict these components, and the correlations with other milk production traits. A total of 4116 milk samples were collected from 414 dairy cows from 44 Canadian herds. Overall 3117 samples were analyzed for casein micelle size, 986 and 937 for total and soluble calcium, respectively, and 2054 for lactoferrin. Individual milk MIR spectra were obtained. Linear mixed models with repeated records were used to test the effects of herd nested within breed, breed, days in milk class, parity, and season on milk components. Lactoferrin content in milk increased with stage of lactation and parity number. Prediction equations were developed by partial least squares regression to predict each milk component using MIR spectroscopy, with the greatest coefficient of determination achieved for lactoferrin (0.55). Pearson correlation coefficients for lactoferrin with somatic cell score and protein percentage were 0.49% and 0.40%, respectively. The estimated correlation increased to 0.53% for predicted lactoferrin and protein percentage. The examined fine milk components exhibited variation, but low accuracies for their MIR prediction may limit the utility of this technology for their phenotyping.

Proteomic Analyses of Mammary Glands Provide Insight into the Immunity and Metabolism Pathways Associated with Clinical Mastitis in Meat Sheep

Simple Summary

Clinical mastitis is one of the most common diseases in sheep and is of major economic concern due to treatment costs, inadequate lamb growth and premature eliminate of ewes. To preliminarily explore possible regulatory roles of proteins involved in the host-pathogen interactions during intramammary infection triggered by this disease in meat sheep, mammary tissues were harvested from sheep with healthy and clinical mastitis caused by natural infection, and the differentially expressed proteins were identified in an infected group when compared to a healthy group, using comparative proteomics based on two-dimensional electrophoresis. Further enrichment analyses indicated that most of the differentially expressed proteins mainly engaged in regulating immune responses and metabolisms. These findings offer candidate proteins for further studies on molecular mechanisms of host defense response and metabolism in sheep cases.

Abstract

Clinical mastitis is still an intractable problem for sheep breeding. The natural immunologic mechanisms of the mammary gland against infections are not yet understood. For a better understanding of the disease-associated proteins during clinical mastitis in meat sheep, we performed two-dimensional electrophoresis (2-DE)-based comparative proteomic analyses of mammary tissues, including from healthy mammary tissues (HMTs) and from mammary tissues with clinical mastitis (CMMTs). The 2-DE results showed that a total of 10 up-regulated and 16 down-regulated proteins were identified in CMMTs when compared to HMTs. Of these, Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) enrichment analyses revealed that most proteins were associated with immune responses or metabolisms. The results of qRT-PCR and Western blot for randomly selected four differentially expressed proteins (DEPs) including superoxide dismutase [Mn] (SOD2), annexin A2 (ANAX2), keratin 10 (KRT10) and endoplasmic reticulum resident protein 29 (ERP29) showed that their expression trends were consistent with 2-DE results except ANXA2 mRNA levels. This is an initial report describing the 2-DE-based proteomics study of the meat sheep mammary gland with clinical mastitis caused by natural infection, which provides additional insight into the immune and metabolic mechanisms during sheep mastitis.

A landscape of the heritability of Fourier-transform infrared spectral wavelengths of milk samples by parity and lactation stage in Holstein cows

Fourier-transform near- and mid-infrared (FTIR) milk spectral data are routinely collected in many countries worldwide. Establishing an optimal strategy to use spectral data in genetic evaluations requires knowledge of the heritabilities of individual FTIR wavelength absorbances. Previous FTIR heritability estimates have been based on relatively small sample sizes and have not considered the possibility that heritability may vary across parities and stages of the lactation. We used data from ∼370,000 test-day records of Canadian Holstein cows to produce a landscape of the heritability of FTIR spectra, 1,060 wavelengths in the near- and mid-infrared spectrum (5,011-925 cm^-1), by parity and month of the lactation (mo 1 to 3 and mo 1 to 6, respectively). The 2 regions of the spectrum associated with absorption of electromagnetic energy by water molecules were estimated to have very high phenotypic variances, very low heritabilities, and very low proportion of variance explained by herd-year-season (HYS) subclasses. The near- or short-wavelength infrared (SWIR: 5,066-3,672 cm^-1) region was also characterized by low heritability estimates, whereas the estimated proportion of the variance explained by HYS was high. The mid-wavelength infrared region (MWIR: 3,000-2,500 cm^-1) and the transition between mid and long-wavelength infrared region (MWIR-LWIR: 1,500-925 cm^-1) harbor several waves characterized by moderately high (≥0.4) heritabilities. Most of the high-heritability regions contained wavelengths that are reported to be associated with important milk metabolites and components. Interestingly, these 2 same regions tended to show more variability in heritabilities between parity and lactation stage. Second parity showed heritability patterns that were distinctly different from those of the first and third parities, whereas the first 2 mo of the lactation had clearly distinct heritability patterns compared with mo 3 to 6.

Short communication: Development of an equation for estimating methane emissions of dairy cows from milk Fourier transform mid-infrared spectra by using reference data obtained exclusively from respiration chambers

Evaluation and mitigation of enteric methane (CH₄) emissions from ruminant livestock, in particular from dairy cows, have acquired global importance for sustainable, climate-smart cattle production. Based on CH₄ reference measurements obtained with the SF₆ tracer technique to determine ruminal CH₄ production, a current equation permits evaluation of individual daily CH₄ emissions of dairy cows based on milk Fourier transform mid-infrared (FT-MIR) spectra. However, the respiration chamber (RC) technique is considered to be more accurate than SF₆ to measure CH₄ production from cattle. This study aimed to develop an equation that allows estimating CH₄ emissions of lactating cows recorded in an RC from corresponding milk FT-MIR spectra and to challenge its robustness and relevance through validation processes and its application on a milk spectral database. This would permit confirming the conclusions drawn with the existing equation based on SF₆ reference measurements regarding the potential to estimate daily CH₄ emissions of dairy cows from milk FT-MIR spectra. A total of 584 RC reference CH₄ measurements (mean ± standard deviation of 400 ± 72 g of CH₄/d) and corresponding standardized milk mid-infrared spectra were obtained from 148 individual lactating cows between 7 and 321 d in milk in 5 European countries (Germany, Switzerland, Denmark, France, and Northern Ireland). The developed equation based on RC measurements showed calibration and cross-validation coefficients of determination of 0.65 and 0.57, respectively, which is lower than those obtained earlier by the equation based on 532 SF₆ measurements (0.74 and 0.70, respectively). This means that the RC-based model is unable to explain the variability observed in the corresponding reference data as well as the SF₆-based model. The standard errors of calibration and cross-validation were lower for the RC model (43 and 47 g/d vs. 66 and 70 g/d for the SF₆ version, respectively), indicating that the model based on RC data was closer to actual values. The root mean squared error (RMSE) of calibration of 42 g/d represents only 10% of the overall daily CH₄ production, which is 23 g/d lower than the RMSE for the SF₆-based equation. During the external validation step an RMSE of 62 g/d was observed. When the RC equation was applied to a standardized spectral database of milk recordings collected in the Walloon region of Belgium between January 2012 and December 2017 (1,515,137 spectra from 132,658 lactating cows in 1,176 different herds), an average ± standard deviation of 446 ± 51 g of CH₄/d was estimated, which is consistent with the range of the values measured using both RC and SF₆ techniques. This study confirmed that milk FT-MIR spectra could be used as a potential proxy to estimate daily CH₄ emissions from dairy cows provided that the variability to predict is covered by the model.

The potential of Fourier transform infrared spectroscopy of milk samples to predict energy intake and efficiency in dairy cows

Knowledge of animal-level and herd-level energy intake, energy balance, and feed efficiency affect day-to-day herd management strategies; information on these traits at an individual animal level is also useful in animal breeding programs. A paucity of data (especially at the individual cow level), of feed intake in particular, hinders the inclusion of such attributes in herd management decision-support tools and breeding programs. Dairy producers have access to an individual cow milk sample at least once daily during lactation, and consequently any low-cost phenotyping strategy should consider exploiting measureable properties in this biological sample, reflecting the physiological status and performance of the cow. Infrared spectroscopy is the study of the interaction of an electromagnetic wave with matter and it is used globally to predict milk quality parameters on routinely acquired individual cow milk samples and bulk tank samples. Thus, exploiting infrared spectroscopy in next-generation phenotyping will ensure potentially rapid application globally with a negligible additional implementation cost as the infrastructure already exists. Fourier-transform infrared spectroscopy (FTIRS) analysis is already used to predict milk fat and protein concentrations, the ratio of which has been proposed as an indicator of energy balance. Milk FTIRS is also able to predict the concentration of various fatty acids in milk, the composition of which is known to change when body tissue is mobilized; that is, when the cow is in negative energy balance. Energy balance is mathematically very similar to residual energy intake (REI), a suggested measure of feed efficiency. Therefore, the prediction of energy intake, energy balance, and feed efficiency (i.e., REI) from milk FTIRS seems logical. In fact, the accuracy of predicting (i.e., correlation between predicted and actual values; root mean square error in parentheses) energy intake, energy balance, and REI from milk FTIRS in dairy cows was 0.88 (20.0MJ), 0.78 (18.6MJ), and 0.63 (22.0MJ), respectively, based on cross-validation. These studies, however, are limited to results from one research group based on data from 2 contrasting production systems in the United Kingdom and Ireland and would need to be replicated, especially in a range of production systems because the prediction equations are not accurate when the variability used in validation is not represented in the calibration data set. Heritable genetic variation exists for all predicted traits. Phenotypic differences in energy intake also exists among animals stratified based on genetic merit for energy intake predicted from milk FTIRS, substantiating the usefulness of such FTIR-predicted phenotypes not only for day-to-day herd management, but also as part of a breeding strategy to improve cow performance.

Correlations between the content of selected whey proteins in cow milk

The aim of the study was to determine the correlations between the concentrations of major whey proteins in cow milk. A total of 2,278 milk samples from Polish Holstein-Friesian (Black-and-White and Red-and-White varieties), Simmental and Jersey cows were analysed. In each sample the content of major whey proteins, i.e. α-lactalbumin, β-lactoglobulin, bovine serum albumin, lactoferrin and lysozyme were determined by the RP-HPLC method. Matrix scatter plots were prepared to determine the correlations between the concentrations of individual whey proteins. In the vast majority of cases a significant relationship was found between the content of individual whey proteins. Taking into account the production season and breed of cow, highly significant (p=0.001) negative correlation coefficients were obtained for the content of α-lactalbumin and bovine serum albumin, for α-lactalbumin and lysozyme, for β-lactoglobulin and bovine serum albumin, and for β-lactoglobulin and lysozyme. Positive correlations were observed for the concentrations of α-lactalbumin with β-lactoglobulin, lactoferrin and lysozyme, as well as for bovine serum albumin with lysozyme.

Short communication: Mid-infrared spectroscopy prediction of fine milk composition and technological properties in Italian Simmental

The objective of this study was to evaluate the ability of mid-infrared predictions of fine milk composition and technological traits to serve as a tool for large-scale phenotyping of the Italian Simmental population. Calibration equations accurately predicted the fatty acid profile of the milk, but we obtained moderate or poor accuracy for detailed protein composition, coagulation properties, curd yield and composition, lactoferrin, and concentration of major minerals. To evaluate the role of infrared predictions as indicator traits of fine milk composition in indirect selective breeding programs, the genetic parameters of the traits predicted using mid-infrared spectra need to be estimated.

A proteomics-based identification of putative biomarkers for disease in bovine milk

The objective of this study was to identify and characterize potential biomarkers for disease resistance in bovine milk that can be used to indicate dairy cows at risk to develop future health problems. We selected high- and low-resistant cows i.e. cows that were less or more prone to develop diseases according to farmers' experience and notifications in the disease registration data. The protein composition of milk serum samples of these high- and low-resistant cows were compared using NanoLC-MS/MS. In total 78 proteins were identified and quantified of which 13 were significantly more abundant in low-resistant cows than high-resistant cows. Quantification of one of these proteins, lactoferrin (LF), by ELISA in a new and much larger set of full fat milk samples confirmed higher LF levels in low- versus high-resistant cows. These high- and low-resistant cows were selected based on comprehensive disease registration and milk recording data, and absence of disease for at least 4 weeks. Relating the experienced diseases to LF levels in milk showed that lameness was associated with higher LF levels in milk. Analysis of the prognostic value of LF showed that low-resistant cows with higher LF levels in milk had a higher risk of being culled within one year after testing than high-resistant cows. In conclusion, LF in milk are higher in low-resistant cows, are associated with lameness and may be a prognostic marker for risk of premature culling.

Milk mid-infrared spectra enable prediction of lactation-stage-dependent methane emissions of dairy cattle within routine population-scale milk recording schemes

Mitigating the proportion of energy intake lost as methane could improve the sustainability and profitability of dairy production. As widespread measurement of methane emissions is precluded by current in vivo methods, the development of an easily measured proxy is desirable. An equation has been developed to predict methane from the mid-infrared (MIR) spectra of milk within routine milk-recording programs. The main goals of this study were to improve the prediction equation for methane emissions from milk MIR spectra and to illustrate its already available usefulness as a high throughput phenotypic screening tool. A total of 532 methane measurements considered as reference data (430 ± 129 g of methane/day) linked with milk MIR spectra were obtained from 165 cows using the SF6 technique. A first derivative was applied to the MIR spectra. Constant (P0), linear (P1) and quadratic (P2) modified Legendre polynomials were computed from each cows stage of lactation (days in milk), at the day of SF6 methane measurement. The calibration model was developed using a modified partial least-squares regression on first derivative MIR data points × P0, first derivative MIR data points × P1, and first derivative MIR data points × P2 as variables. The MIR-predicted methane emissions (g/day) showed a calibration coefficient of determination of 0.74, a cross-validation coefficient of determination of 0.70 and a standard error of calibration of 66 g/day. When applied to milk MIR spectra recorded in the Walloon Region of Belgium (≈2 000 000 records), this equation was useful to study lactational, annual, seasonal, and regional methane emissions. We conclude that milk MIR spectra has potential to be used to conduct high throughput screening of lactating dairy cattle for methane emissions. The data generated enable monitoring of methane emissions and production characteristics across and within herds. Milk MIR spectra could now be used for widespread screening of dairy herds in order to develop management and genetic selection tools to reduce methane emissions.

Invited review: overview of new traits and phenotyping strategies in dairy cattle with a focus on functional traits

For several decades, breeding goals in dairy cattle focussed on increased milk production. However, many functional traits have negative genetic correlations with milk yield, and reductions in genetic merit for health and fitness have been observed. Herd management has been challenged to compensate for these effects and to balance fertility, udder health and metabolic diseases against increased production to maximize profit without compromising welfare. Functional traits, such as direct information on cow health, have also become more important because of growing concern about animal well-being and consumer demands for healthy and natural products. There are major concerns about the impact of drugs used in veterinary medicine on the spread of antibiotic-resistant strains of bacteria that can negatively impact human health. Sustainability and efficiency are also increasingly important because of the growing competition for high-quality, plant-based sources of energy and protein. Disruptions to global environments because of climate change may encourage yet more emphasis on these traits. To be successful, it is vital that there be a balance between the effort required for data recording and subsequent benefits. The motivation of farmers and other stakeholders involved in documentation and recording is essential to ensure good data quality. To keep labour costs reasonable, existing data sources should be used as much as possible. Examples include the use of milk composition data to provide additional information about the metabolic status or energy balance of the animals. Recent advances in the use of mid-infrared spectroscopy to measure milk have shown considerable promise, and may provide cost-effective alternative phenotypes for difficult or expensive-to-measure traits, such as feed efficiency. There are other valuable data sources in countries that have compulsory documentation of veterinary treatments and drug use. Additional sources of data outside of the farm include, for example, slaughter houses (meat composition and quality) and veterinary labs (specific pathogens, viral loads). At the farm level, many data are available from automated and semi-automated milking and management systems. Electronic devices measuring physiological status or activity parameters can be used to predict events such as oestrus, and also behavioural traits. Challenges concerning the predictive biology of indicator traits or standardization need to be solved. To develop effective selection programmes for new traits, the development of large databases is necessary so that high-reliability breeding values can be estimated. For expensive-to-record traits, extensive phenotyping in combination with genotyping of females is a possibility.

Breeding the dairy cow of the future: what do we need?

Genetics is responsible for approximately half the observed changes in animal performance in well structured breeding programs. Key characteristics of the dairy cow of the future include (1) production of a large quantity of high-value output (i.e. milk and meat), (2) good reproductive performance, (3) good health status, (4) good longevity, (5) no requirement for a large quantity of feed, yet being able to eat sufficient feed to meet its requirements, (6) easy to manage (i.e. easy calving, docile), (7) good conformation (over and above reflective of health, reproductive performance and longevity), (8) low environmental footprint, and (9) resilience to external perturbations. Pertinent and balanced breeding goals must be developed and implemented to achieve this type of animal; excluding any characteristic from the breeding goal could be detrimental for genetic gain in this characteristic. Attributes currently not explicitly considered in most dairy-cow breeding objectives include product quality, feed intake and efficiency, and environmental footprint; animal health is poorly represented in most breeding objectives. Lessons from the past deterioration in reproductive performance in the global Holstein population remind us of the consequences of ignoring or failing to monitor certain animal characteristics. More importantly, however, current knowledge clearly demonstrates that once unfavourable trends have been identified and the appropriate breeding strategy implemented, the reversal of genetic trends is achievable, even for low-heritability traits such as reproductive performance. Genetic variation exists in all the characteristics described. In the genomics era, the relevance of heritability statistics for most traits is less; the exception is traits not amenable to routine measurement in large populations. Phenotyping strategies (e.g. more detailed phenotypes, larger population) will remain a key component of an animal breeding strategy to achieve the cow of the future as well as providing the necessary tools and information to monitor performance. The inclusion of genomic information in genetic evaluations is, and will continue, to improve the accuracy of genetic evaluations, which, in turn, will augment genetic gain; genomics, however, can also contribute to gains in performance over and above support of increased genetic gain. Nonetheless, the faster genetic gain and thus reduced ability to purge out unfavourable alleles necessitates the appropriate breeding goal and breeding scheme and very close monitoring of performance, in particular for traits not included in the breeding goals. Developments in other disciplines (e.g. reproductive technologies), coupled with commercial struggle for increased market share of the breeding industry, imply a possible change in the landscape of dairy-cow breeding in the future.

Evaluating somatic cell scores with a Bayesian Gaussian linear state-space model

Because accurate characterization of health state is important for managing dairy herds, we propose to validate the use of a linear state-space model (LSSM) for evaluating monthly somatic cell scores (SCSs). To do so, we retrieved SCS from a dairy database and collected reports on clinical mastitis collected in 20 farms, during the period from January 2008 to December 2011 in the Walloon region of Belgium. The dependent variable was the SCS, and the independent variables were the number of days from calving, year of calving and parity. The LSSM also incorporated an error-free underlying variable that described the trend across time as a function of previous clinical and subclinical status. We computed the mean sum of squared differences between observed SCS and median values of the posterior SCS distribution and constructed the receiver operating characteristic (ROC) curve for SCS thresholds going from 0 to 6. Our results show SCS estimates are close to observed SCS and area under the ROC curve is higher than 90%. We discuss the meaning of the parameters in light of our current knowledge of the disease and propose methods to incorporate, in LSSM, this knowledge often expressed in the form of ordinary differential equations.