Microwave Dielectric Response of Bovine Milk as Pregnancy Detection Tool in Dairy Cows

The most reliable methods for pregnancy diagnosis in dairy herds include rectal palpation, ultrasound examination, and evaluation of plasma progesterone concentrations. However, these methods are expensive, labor-intensive, and invasive. Thus, there is a need to develop a practical, non-invasive, cost-effective method that can be implemented on the farm to detect pregnancy. This study suggests employing microwave dielectric spectroscopy (MDS, 0.5-40 GHz) as a method to evaluate reproduction events in dairy cows. The approach involves the integration of MDS data with information on milk solids to detect pregnancy and identify early embryonic loss in dairy cows. To test the ability to predict pregnancy according to these measurements, milk samples were collected from (i) pregnant and non-pregnant randomly selected cows, (ii) weekly from selected cows (n = 12) before insemination until a positive pregnancy test, and (iii) daily from selected cows (n = 10) prior to insemination until a positive pregnancy test. The results indicated that the dielectric strength of Δε and the relaxation time, τ, exhibited reduced variability in the case of a positive pregnancy diagnosis. Using principal component analysis (PCA), a clear distinction between pregnancy and nonpregnancy status was observed, with improved differentiation upon a higher sampling frequency. Additionally, a neural network machine learning technique was employed to develop a prediction algorithm with an accuracy of 73%. These findings demonstrate that MDS can be used to detect changes in milk upon pregnancy. The developed machine learning provides a broad classification that could be further enhanced with additional data.

Assessing the relationship between somatic cell count and the milk mid-infrared spectrum in Chinese Holstein cows

BACKGROUND

Milk produced by dairy cows is a complex combination of many components, but the effect of mastitis has only been investigated for a few of these components. Milk mid-infrared (MIR) spectra can reflect the global composition of milk, and this study aimed to detect the relationships between milk MIR spectral wavenumbers and milk somatic cell count (SCC)-a sensitive biomarker for mastitis.

METHODS

Pearson correlation analysis was used to calculate the correlation coefficient between somatic count score (SCS) and spectral wavenumbers. A general linear mixed model was applied to investigate the effect of three different classes of SCC (low, middle and high) on spectral wavenumbers.

RESULTS

The mean correlation coefficient between the 'fingerprint region' (wavenumbers 925-1582 cm-1 ) and the SCS was higher than that for other regions of the MIR spectrum, and the specific wavenumber with the strongest correlation with the SCS was within the 'fingerprint region'. SCC class had a significant (p < 0.05) effect on 639 spectral wavenumbers. In particular, some spectral wavenumbers within the 'fingerprint region' were highly affected by the SCC class.

LIMITATION

The data were collected from only one province in China, so the generalisability of the findings may be limited.

CONCLUSION

SCC had close relationships with milk spectral wavenumbers related to important milk components or chemical bonds.

A machine learning approach using partitioning around medoids clustering and random forest classification to model groups of farms in regard to production parameters and bulk tank milk antibody status of two major internal parasites in dairy cows

Fasciola hepatica and Ostertagia ostertagi are internal parasites of cattle compromising physiology, productivity, and well-being. Parasites are complex in their effect on hosts, sometimes making it difficult to identify clear directions of associations between infection and production parameters. Therefore, unsupervised approaches not assuming a structure reduce the risk of introducing bias to the analysis. They may provide insights which cannot be obtained with conventional, supervised methodology. An unsupervised, exploratory cluster analysis approach using the k–mode algorithm and partitioning around medoids detected two distinct clusters in a cross-sectional data set of milk yield, milk fat content, milk protein content as well as F. hepatica or O. ostertagi bulk tank milk antibody status from 606 dairy farms in three structurally different dairying regions in Germany. Parasite–positive farms grouped together with their respective production parameters to form separate clusters. A random forests algorithm characterised clusters with regard to external variables. Across all study regions, co–infections with F. hepatica or O. ostertagi, respectively, farming type, and pasture access appeared to be the most important factors discriminating clusters (i.e. farms). Furthermore, farm level lameness prevalence, herd size, BCS, stage of lactation, and somatic cell count were relevant criteria distinguishing clusters. This study is among the first to apply a cluster analysis approach in this context and potentially the first to implement a k–medoids algorithm and partitioning around medoids in the veterinary field. The results demonstrated that biologically relevant patterns of parasite status and milk parameters exist between farms positive for F. hepatica or O. ostertagi, respectively, and negative farms. Moreover, the machine learning approach confirmed results of previous work and shed further light on the complex setting of associations a between parasitic diseases, milk yield and milk constituents, and management practices.

Pregnancy status predicted using milk mid-infrared spectra from dairy cattle

Accurate and timely pregnancy diagnosis is an important component of effective herd management in dairy cattle. Predicting pregnancy from Fourier-transform mid-infrared (FT-MIR) spectroscopy data is of particular interest because the data are often already available from routine milk testing. The purpose of this study was to evaluate how well pregnancy status could be predicted in a large data set of 1,161,436 FT-MIR milk spectra records from 863,982 mixed-breed pasture-based New Zealand dairy cattle managed within seasonal calving systems. Three strategies were assessed for defining the nonpregnant cows when partitioning the records according to pregnancy status in the training population. Two of these used records for cows with a subsequent calving only, whereas the third also included records for cows without a subsequent calving. For each partitioning strategy, partial least squares discriminant analysis models were developed, whereby spectra from all the cows in 80% of herds were used to train the models, and predictions on cows in the remaining herds were used for validation. A separate data set was also used as a secondary validation, whereby pregnancy diagnosis had been assigned according to the presence of pregnancy-associated glycoproteins (PAG) in the milk samples. We examined different ways of accounting for stage of lactation in the prediction models, either by including it as an effect in the prediction model, or by pre-adjusting spectra before fitting the model. For a subset of strategies, we also assessed prediction accuracies from deep learning approaches, utilizing either the raw spectra or images of spectra. Across all strategies, prediction accuracies were highest for models using the unadjusted spectra as model predictors. Strategies for cows with a subsequent calving performed well in herd-independent validation with sensitivities above 0.79, specificities above 0.91 and area under the receiver operating characteristic curve (AUC) values over 0.91. However, for these strategies, the specificity to predict nonpregnant cows in the external PAG data set was poor (0.002-0.04). The best performing models were those that included records for cows without a subsequent calving, and used unadjusted spectra and days in milk as predictors, with consistent results observed across the training, herd-independent validation and PAG data sets. For the partial least squares discriminant analysis model, sensitivity was 0.71, specificity was 0.54 and AUC values were 0.68 in the PAG data set; and for an image-based deep learning model, the sensitivity was 0.74, specificity was 0.52 and the AUC value was 0.69. Our results demonstrate that in pasture-based seasonal calving herds, confounding between pregnancy status and spectral changes associated with stage of lactation can inflate prediction accuracies. When the effect of this confounding was reduced, prediction accuracies were not sufficiently high enough to use as a sole indicator of pregnancy status.

Sequence-based genome-wide association study of individual milk mid-infrared wavenumbers in mixed-breed dairy cattle

BACKGROUND

Fourier-transform mid-infrared (FT-MIR) spectroscopy provides a high-throughput and inexpensive method for predicting milk composition and other novel traits from milk samples. While there have been many genome-wide association studies (GWAS) conducted on FT-MIR predicted traits, there have been few GWAS for individual FT-MIR wavenumbers. Using imputed whole-genome sequence for 38,085 mixed-breed New Zealand dairy cattle, we conducted GWAS on 895 individual FT-MIR wavenumber phenotypes, and assessed the value of these direct phenotypes for identifying candidate causal genes and variants, and improving our understanding of the physico-chemical properties of milk.

RESULTS

Separate GWAS conducted for each of 895 individual FT-MIR wavenumber phenotypes, identified 450 1-Mbp genomic regions with significant FT-MIR wavenumber QTL, compared to 246 1-Mbp genomic regions with QTL identified for FT-MIR predicted milk composition traits. Use of mammary RNA-seq data and gene annotation information identified 38 co-localized and co-segregating expression QTL (eQTL), and 31 protein-sequence mutations for FT-MIR wavenumber phenotypes, the latter including a null mutation in the ABO gene that has a potential role in changing milk oligosaccharide profiles. For the candidate causative genes implicated in these analyses, we examined the strength of association between relevant loci and each wavenumber across the mid-infrared spectrum. This revealed shared association patterns for groups of genomically-distant loci, highlighting clusters of loci linked through their biological roles in lactation and their presumed impacts on the chemical composition of milk.

CONCLUSIONS

This study demonstrates the utility of FT-MIR wavenumber phenotypes for improving our understanding of milk composition, presenting a larger number of QTL and putative causative genes and variants than found from FT-MIR predicted composition traits. Examining patterns of significance across the mid-infrared spectrum for loci of interest further highlighted commonalities of association, which likely reflects the physico-chemical properties of milk constituents.

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.

Integrating genomic and infrared spectral data improves the prediction of milk protein composition in dairy cattle

Background

Over the past decade, Fourier transform infrared (FTIR) spectroscopy has been used to predict novel milk protein phenotypes. Genomic data might help predict these phenotypes when integrated with milk FTIR spectra. The objective of this study was to investigate prediction accuracy for milk protein phenotypes when heterogeneous on-farm, genomic, and pedigree data were integrated with the spectra. To this end, we used the records of 966 Italian Brown Swiss cows with milk FTIR spectra, on-farm information, medium-density genetic markers, and pedigree data. True and total whey protein, and five casein, and two whey protein traits were analyzed. Multiple kernel learning constructed from spectral and genomic (pedigree) relationship matrices and multilayer BayesB assigning separate priors for FTIR and markers were benchmarked against a baseline partial least squares (PLS) regression. Seven combinations of covariates were considered, and their predictive abilities were evaluated by repeated random sub-sampling and herd cross-validations (CV).

Results

Addition of the on-farm effects such as herd, days in milk, and parity to spectral data improved predictions as compared to those obtained using the spectra alone. Integrating genomics and/or the top three markers with a large effect further enhanced the predictions. Pedigree data also improved prediction, but to a lesser extent than genomic data. Multiple kernel learning and multilayer BayesB increased predictive performance, whereas PLS did not. Overall, multilayer BayesB provided better predictions than multiple kernel learning, and lower prediction performance was observed in herd CV compared to repeated random sub-sampling CV.

Conclusions

Integration of genomic information with milk FTIR spectral can enhance milk protein trait predictions by 25% and 7% on average for repeated random sub-sampling and herd CV, respectively. Multiple kernel learning and multilayer BayesB outperformed PLS when used to integrate heterogeneous data for phenotypic predictions.

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.

Association between days open and milk spectral data in dairy cows

Data on 19,489 Brown Swiss cows reared in northeastern Italy were used to associate absorbances of individual wavenumbers within the mid-infrared range with days open (DO). Different postcalving days in milk (DIM) intervals were studied to determine the most informative milk sampling periods for predicting DO. Milk samples were analyzed using a MilkoScan (Foss Electric, Hillerød, Denmark) Fourier-transform infrared (FTIR) spectrometer for 1,060 wavenumbers (wn) ranging from 5,011 to 925 cm^-1. To determine DO, we considered an insemination to lead to conception when there was no return of heat (i.e., no successive insemination) and the cow had a subsequent calving date whereby gestation length was required to be within ±30 d of 290 d. Only milk records within the first 90 DIM were considered. Associations were inferred by (1) fitting linear regression models between the DO and each individual wavenumber or milk component, and (2) fitting a Bayesian regression model that included the complete FTIR spectral data. The effects of including systematic effects (parity number, year-season, herd) in the model on these associations were also studied. These analyses were performed for the complete data (5-90 DIM) and for data stratified by DIM period (5 to 30, 31 to 60, and 61 to 90 DIM). Overall, regions of wavenumbers of the milk FTIR spectra that were associated with DO included wn 2,973 to 2,830 cm^-1 [related to fat-B (C-H stretch)], wn 2,217 to 1,769 cm^-1 [related to fat-A (C = O stretch)], wn 1,546 cm^-1 (related to protein), wn 1,465 cm^-1 (related to urea and fat), wn 1,399 to 1,245 cm^-1 (related to acetone), and wn 1,110 cm^-1 (related to lactose). Estimated effects depended on the DIM period, with milk samples drawn during DIM intervals 31 to 60 d and 61 to 90 d being most strongly associated with DO. These DIM intervals are also typically most associated with negative energy balance and peak lactation.

Phenotypic and genetic effects of pregnancy on milk production traits in Holstein-Friesian cattle

Pregnancy is a prerequisite for the initiation of lactation and for maintaining the milk production cycle. Pregnancy affects milk production and therefore should be accounted for in the genetic evaluation. Furthermore, there might be genetic differences in pregnancy effects on milk composition. The objective of this study was to estimate phenotypic and genetic effects of pregnancy on milk production traits. For this purpose, test-day records and conception dates of 1,359 first-parity Holstein-Friesian cows were analyzed. Significant effects of pregnancy on all milk production traits were detected except somatic cell score (e.g., the cumulative effects of pregnancy on milk yield were -247 kg). The pregnancy effects on milk yield, lactose yield, protein yield, fat yield, and fat content were small during early gestation (<150 d) and substantially increased in late gestation. The effects of pregnancy on milk protein yield were relatively stronger than those on fat yield. The effects of pregnancy on milk production traits differed for DGAT1 genotypes. Milk yield, lactose yield, protein yield, and fat yield of DGAT1 AA cows were more affected by pregnancy than that of DGAT1 KK cows (e.g., the cumulative effects of pregnancy on milk yield were negligible for DGAT1 KK cows and were -443 kg for DGAT1 AA cows). These results suggest that DGAT1 KK cows may be more suitable for shortening or omitting the dry period than DGAT1 AA cows.

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.

Predicting milk mid-infrared spectra from first-parity Holstein cows using a test-day mixed model with the perspective of herd management

The use of test-day models to model milk mid-infrared (MIR) spectra for genetic purposes has already been explored; however, little attention has been given to their use to predict milk MIR spectra for management purposes. The aim of this paper was to study the ability of a test-day mixed model to predict milk MIR spectra for management purposes. A data set containing 467,496 test-day observations from 53,781 Holstein dairy cows in first lactation was used for model building. Principal component analysis was implemented on the selected 311 MIR spectral wavenumbers to reduce the number of traits for modeling; 12 principal components (PC) were retained, explaining approximately 96% of the total spectral variation. Each of the retained PC was modeled using a single trait test-day mixed model. The model solutions were used to compute the predicted scores of each PC, followed by a back-transformation to obtain the 311 predicted MIR spectral wavenumbers. Four new data sets, containing altogether 122,032 records, were used to test the ability of the model to predict milk MIR spectra in 4 distinct scenarios with different levels of information about the cows. The average correlation between observed and predicted values of each spectral wavenumber was 0.85 for the modeling data set and ranged from 0.36 to 0.62 for the scenarios. Correlations between milk fat, protein, and lactose contents predicted from the observed spectra and from the modeled spectra ranged from 0.83 to 0.89 for the modeling set and from 0.32 to 0.73 for the scenarios. Our results demonstrated a moderate but promising ability to predict milk MIR spectra using a test-day mixed model. Current and future MIR traits prediction equations could be applied on the modeled spectra to predict all MIR traits in different situations instead of developing one test-day model separately for each trait. Modeling MIR spectra would benefit farmers for cow and herd management, for instance through prediction of future records or comparison between observed and expected wavenumbers or MIR traits for the detection of health and management problems. Potential resulting tools could be incorporated into milk recording systems.

Use of Large and Diverse Datasets for 1H NMR Serum Metabolic Profiling of Early Lactation Dairy Cows

Most livestock metabolomic studies involve relatively small, homogenous populations of animals. However, livestock farming systems are non-homogenous, and large and more diverse datasets are required to ensure that biomarkers are robust. The aims of this study were therefore to (1) investigate the feasibility of using a large and diverse dataset for untargeted proton nuclear magnetic resonance (¹H NMR) serum metabolomic profiling, and (2) investigate the impact of fixed effects (farm of origin, parity and stage of lactation) on the serum metabolome of early-lactation dairy cows. First, we used multiple linear regression to correct a large spectral dataset (707 cows from 13 farms) for fixed effects prior to multivariate statistical analysis with principal component analysis (PCA). Results showed that farm of origin accounted for up to 57% of overall spectral variation, and nearly 80% of variation for some individual metabolite concentrations. Parity and week of lactation had much smaller effects on both the spectra as a whole and individual metabolites (< 3% and < 20%, respectively). In order to assess the effect of fixed effects on prediction accuracy and biomarker discovery, we used orthogonal partial least squares (OPLS) regression to quantify the relationship between NMR spectra and concentrations of the current gold standard serum biomarker of energy balance, β-hydroxybutyrate (BHBA). Models constructed using data from multiple farms provided reasonably robust predictions of serum BHBA concentration (0.05 ≤ RMSE ≤ 0.18). Fixed effects influenced the results biomarker discovery; however, these impacts could be controlled using the proposed method of linear regression spectral correction.

The evolving role of Fourier-transform mid-infrared spectroscopy in genetic improvement of dairy cattle

Over the last 100 years, significant advances have been made in the characterisation of milk composition for dairy cattle improvement programs. Technological progress has enabled a shift from labour intensive, on-farm collection and processing of samples that assess yield and fat levels in milk, to large-scale processing of samples through centralised laboratories, with the scope extended to include quantification of other traits. Fourier-transform mid-infrared (FT-MIR) spectroscopy has had a significant role in the transformation of milk composition phenotyping, with spectral-based predictions of major milk components already being widely used in milk payment and animal evaluation systems globally. Increasingly, there is interest in analysing the individual FT-MIR wavenumbers, and in utilising the FT-MIR data to predict other novel traits of importance to breeding programs. This includes traits related to the nutritional value of milk, the processability of milk into products such as cheese, and traits relevant to animal health and the environment. The ability to successfully incorporate these traits into breeding programs is dependent on the heritability of the FT-MIR predicted traits, and the genetic correlations between the FT-MIR predicted and actual trait values. Linking FT-MIR predicted traits to the underlying mutations responsible for their variation can be difficult because the phenotypic expression of these traits are a function of a diverse range of molecular and biological mechanisms that can obscure their genetic basis. The individual FT-MIR wavenumbers give insights into the chemical composition of milk and provide an additional layer of granularity that may assist with establishing causal links between the genome and observed phenotypes. Additionally, there are other molecular phenotypes such as those related to the metabolome, chromatin accessibility, and RNA editing that could improve our understanding of the underlying biological systems controlling traits of interest. Here we review topics of importance to phenotyping and genetic applications of FT-MIR spectra datasets, and discuss opportunities for consolidating FT-MIR datasets with other genomic and molecular data sources to improve future dairy cattle breeding programs.

Genome-wide association study for genotype by lactation stage interaction of milk production traits in dairy cattle

Substantial evidence demonstrates that the genetic background of milk production traits changes during lactation. However, most GWAS for milk production traits assume that genetic effects are constant during lactation and therefore might miss those quantitative trait loci (QTL) whose effects change during lactation. The GWAS for genotype by lactation stage interaction are aimed at explicitly detecting the QTL whose effects change during lactation. The purpose of this study was to perform GWAS for genotype by lactation stage interaction for milk yield, lactose yield, lactose content, fat yield, fat content, protein yield, and somatic cell score to detect QTL with changing effects during lactation. For this study, 19,286 test-day records of 1,800 first-parity Dutch Holstein cows were available and cows were genotyped using a 50K SNP panel. A total of 7 genomic regions with effects that change during lactation were detected in the GWAS for genotype by lactation stage interaction. Two regions on Bos taurus autosome (BTA)14 and BTA19 were also significant based on a GWAS that assumed constant genetic effects during lactation. Five regions on BTA4, BTA10, BTA11, BTA16, and BTA23 were only significant in the GWAS for genotype by lactation stage interaction. The biological mechanisms that cause these changes in genetic effects are still unknown, but negative energy balance and effects of pregnancy may play a role. These findings increase our understanding of the genetic background of lactation and may contribute to the development of better management indicators based on milk composition.

Diagnosing the pregnancy status of dairy cows: How useful is milk mid-infrared spectroscopy?

Pregnancy diagnosis is an essential part of successful breeding programs on dairy farms. Milk composition alters with pregnancy, and this is well documented. Fourier-transform mid-infrared (MIR) spectroscopy is a rapid and cost-effective method for providing milk spectra that reflect the detailed composition of milk samples. Therefore, the aim of this study was to assess the ability of MIR spectroscopy to predict the pregnancy status of dairy cows. The MIR spectra and insemination records were available from 8,064 Holstein cows of 19 commercial dairy farms in Australia. Three strategies were studied to classify cows as open or pregnant using partial least squares discriminant analysis models with random cow-independent 10-fold cross-validation and external validation on a cow-independent test set. The first strategy considered 6,754 MIR spectra after insemination used as independent variables in the model. The results showed little ability to detect the pregnancy status as the area under the receiver operating characteristic curve was 0.63 and 0.65 for cross-validation and testing, respectively. The second strategy, involving 1,664 records, aimed to reduce noise in the MIR spectra used as predictors by subtracting a spectrum before insemination (i.e., open spectrum) from the spectrum after insemination. The accuracy was comparable with the first approach, showing no superiority of the method. Given the limited results for these models when using combined data from all stages after insemination, the third strategy explored separate models at 7 stages after insemination comprising 348 to 1,566 records each (i.e., progressively greater gestation) with single MIR spectra after insemination as predictors. The models developed using data recorded after 150 d of pregnancy showed promising prediction accuracy with the average value of area under the receiver operating characteristic curve of 0.78 and 0.76 obtained through cross-validation and testing, respectively. If this can be confirmed on a larger data set and extended to somewhat earlier stages after insemination, the model could be used as a complementary tool to detect fetal abortion.

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.

Genetic Analysis of Milk Production Traits and Mid-Infrared Spectra in Chinese Holstein Population

Simple Summary

Usually, spectral data are used as predictors to predict milk components, animal characteristics, and even reproductive status. Another innovative way to use spectral data involves considering spectral wavenumbers as traits and then analyzing from the genetic perspective. In this study, we considered milk spectral data directly as traits, then detected the influence of some non-genetic factors on spectral wavenumbers and estimated the genetic parameters of spectral points. The result of the present study could be used as a management tool for dairy farm and also provides a further understanding of genetic background of milk mid-infrared (MIR) spectra. In future, milk spectral data could be applied more effective. For example, some sub-clinical diseases might be detected based on the difference between the expected and observed values of the spectral traits. In addition, we could also use genetic correlation between wavenumbers and a trait of interest, which are difficult and expensive to measure, to apply for the genetic improvement of dairy species.

Abstract

Milk composition always serves as an indicator for the cow’s health status and body condition. Some non-genetic factors such as parity, days in milk (DIM), and calving season, which obviously affect milk performance, therefore, need to be considered in dairy farm management. However, only a few milk compositions are used in the current animal selection programs. The mid-infrared (MIR) spectroscopy can reflect the global composition of milk, but this information is currently underused. The objectives of this study were to detect the effect of some non-genetic factors on milk production traits as well as 1060 individual spectral points covering from 925.92 cm⁻¹ to 5011.54 cm⁻¹, estimate heritabilities of milk production traits and MIR spectral wavenumbers, and explore the genetic correlations between milk production traits and 1060 individual spectral points in a Chinese Holstein population. The mixed models procedure of SAS software was used to test the non-genetic factors. Single-trait animal models were used to estimate heritabilities and bivariate animal models were used to estimate genetic correlations using the package of ASReml in R software. The results showed that herd, parity, calving season, and lactation stage had significant effects on the percentages of protein and lactose, whereas herd and lactation stage had significant effects on fat percentage. Moreover, the herd showed a significant effect on all of the 1060 individual wavenumbers, whereas lactation stage, parity, and calving season had significant effect on most of the wavenumbers of the lactose-region (925 cm⁻¹ to 1200 cm⁻¹), protein-region (1240 cm⁻¹ to 1600 cm⁻¹), and fat-regions (1680 cm⁻¹ to 1770 cm⁻¹ and 2800 cm⁻¹ to 3015 cm⁻¹). The estimated heritabilities for protein percentage (PP), fat percentage (FP), and lactose percentage (LP) were 0.08, 0.05, and 0.09, respectively. Further, the milk spectrum was heritable but low for most individual points. Heritabilities of 1060 individual spectral points were 0.04 on average, ranging from 0 to 0.11. In particular, heritabilities for wavenumbers of spectral regions related to water absorption were very low and even null, and heritabilities for wavenumbers of specific MIR regions associated with fat-I, fat-II, protein, and lactose were 0.04, 0.06, 0.05, and 0.06 on average, respectively. The genetic correlations between PP and FP, PP and LP, FP, and LP were 0.78, −0.29, and −0.14, respectively. In addition, PP, FP, and LP shared the similar patterns of genetic correlations with the spectral wavenumbers. The genetic correlations between milk production traits and spectral regions related to important milk components varied from weak to very strong (0.01 to 0.94, and −0.01 to −0.96). The current study could be used as a management tool for dairy farms and also provides a further understanding of the genetic background of milk MIR spectra.

Predicting milk fatty acids and energy balance of dairy cows in Australia using milk mid-infrared spectroscopy

Mid-infrared spectroscopy (MIRS) is traditionally used for analysing milk fat, protein and lactose concentrations in dairy production, but there is growing interest in using it to predict difficult, or expensive-to-measure, phenotypes on a large scale. The resulting prediction equations can be applied to MIRS data from commercial herd-testing, to facilitate management and feeding decisions, or for genomic selection purposes. We investigated the ability of MIRS of milk samples to predict milk fatty acids (FAs) and energy balance (EB) of dairy cows in Australia. Data from 240 Holstein lactating cows that were part of two 32-day experiments, were used. Milk FAs were measured twice during the experimental period. Prediction models were developed using partial least-square regression with a 10-fold cross-validation. Measures of prediction accuracy included the coefficient of determination (R2cv) and root mean-square error. Milk FAs with a chain length of ≤16 were accurately predicted (0.89 ≤ R2cv ≤ 0.95), while prediction accuracy for FAs with a chain length of ≥17 was slightly lower (0.72 ≤ R2cv ≤ 0.82). The accuracy of the model prediction was moderate for EB, with the value of R2cv of 0.48. In conclusion, the ability of MIRS to predict milk FAs was high, while EB was moderately predicted. A larger dataset is needed to improve the accuracy and the robustness of the prediction models.

Energy profiling of dairy cows from routine milk mid-infrared analysis

The balance of body energy within and across lactations can have health and fertility consequences for the dairy cow. This study aimed to create a large calibration data set of dairy cow body energy traits across the cow's productive life, with concurrent milk mid-infrared (MIR) spectral data, to generate a prediction tool for use in commercial dairy herds. Detailed phenotypic data from 1,101 Holstein Friesian cows from the Langhill research herd (SRUC, Scotland) were used to generate energy balance (EB) and effective energy intake (EI), both in megajoules per day. Pretreatment of spectral data involved standardization to account for drift over time and machine. Body energy estimates were aligned with their spectral data to generate a prediction of these traits based on milk MIR spectroscopy. After data edits, partial least squares analysis generated prediction equations with a coefficient of determination from split sample 10-fold cross validation of 0.77 and 0.75 for EB and EI, respectively. These prediction equations were applied to national milk MIR spectra on over 11 million animal test dates (January 2013 to December 2016) from 4,453 farms. The predictions generated from these were subject to phenotypic analyses with a fixed regression model highlighting differences between the main dairy breeds in terms of energy traits. Genetic analyses generated heritability estimates for EB and EI ranging from 0.12 to 0.17 and 0.13 to 0.15, respectively. This study shows that MIR-based predictions from routinely collected national data can be used to generate predictions of dairy cow energy turnover profiles for both animal management and genetic improvement of such difficult and expensive-to-record traits.

Strategies for noise reduction and standardization of milk mid-infrared spectra from dairy cattle

The use of Fourier-transform mid-infrared (FTIR) spectroscopy is of interest to the dairy industry worldwide for predicting milk composition and other novel traits that are difficult or expensive to measure directly. Although there are many valuable applications for FTIR spectra, noise from differences in spectral responses between instruments is problematic because it reduces prediction accuracy if ignored. The purpose of this study was to develop strategies to reduce the impact of noise and to compare methods for standardizing FTIR spectra in order to reduce between-instrument variability in multiple-instrument networks. Noise levels in bands of the infrared spectrum caused by the water content of milk were characterized, and a method for identifying and removing outliers was developed. Two standardization methods were assessed and compared: piecewise direct standardization (PDS), which related spectra on a primary instrument to spectra on 5 other (secondary) instruments using identical milk-based reference samples (n = 918) analyzed across the 6 instruments; and retroactive percentile standardization (RPS), whereby percentiles of observed spectra from routine milk test samples (n = 2,044,094) were used to map and exploit primary- and secondary-instrument relationships. Different applications of each method were studied to determine the optimal way to implement each method across time. Industry-standard predictions of milk components from 2,044,094 spectra records were regressed against predictions from spectra before and after standardization using PDS or RPS. The PDS approach resulted in an overall decrease in root mean square error between industry-standard predictions and predictions from spectra from 0.190 to 0.071 g/100 mL for fat, from 0.129 to 0.055 g/100 mL for protein, and from 0.143 to 0.088 g/100 mL for lactose. Reductions in prediction error for RPS were similar but less consistent than those for PDS across time, but similar reductions were achieved when PDS coefficients were updated monthly and separate primary instruments were assigned for the North and South Islands of New Zealand. We demonstrated that the PDS approach is the most consistent method to reduce prediction errors across time. We also showed that the RPS approach is sensitive to shifts in milk composition but can be used to reduce prediction errors, provided that secondary-instrument spectra are standardized to a primary instrument with samples of broadly equivalent milk composition. Appropriate implementation of either of these approaches will improve the quality of predictions based on FTIR spectra for various downstream applications.

Using geometric morphometrics for the genetics analysis of shape and size of lactation curves in Israeli first-parity Holstein cattle

Our objective was to combine the methods of geometric morphometrics and multivariate quantitative genetics to genetic evaluation of the size and shape of lactation curves of milk of 3,492 Israeli first-parity Holstein cattle. Lactation records were treated as morphological data, for which 2 different lactation shape functions were evaluated, one depicted by a line graph and the other by an orbital graph. The lactation curves from both shape functions were represented by 2-dimensional Cartesian landmark coordinates. The 2 sets of landmarks were then analyzed individually for each shape function with geometric morphometrics to separate variation into components of size and shape. The analysis yielded 2 size measures and 2 sets of shape variables, and they were the inputs to estimate variance components using the MTC REML individual animal model program. Variance components were also estimated for the 305-d lactation production as a reference. Shape variables showed negligible correlation with 305-d production, providing evidence of size and shape of lactation curve as separate characters. The size measure derived from the orbital-depicted lactation curve had equal heritability (0.39 ± 0.01; ± standard error) and complete genetic and environmental correlations with 305-d production, whereas the size measure derived from the line-depicted lactation curve showed low heritability (0.09 ± 0.01) and environmental correlation (0.02 ± 0.004) and relative high genetic correlation with 305-d production (0.48 ± 0.04). This may validate both the orbital graph to depict lactation records and the use of geometric morphometrics to split variation of lactation curve into size and shape components. The maximal heritability for shape of lactation curve was 0.55 for orbital- and 0.56 for line-depicted lactation curves. The respective patterns of variations were visualized as shape changes from the mean shape in the data set. Geometric morphometrics are well grounded within the theory of shape analysis and can be paired with conventional methods in the field to characterize the patterns of phenotypic and genetic variation of shape and size of lactation curve in dairy cattle.

Standardization of milk mid-infrared spectrometers for the transfer and use of multiple models

An increasing number of models are being developed to provide information from milk Fourier transform mid-infrared (FT-MIR) spectra on fine milk composition, technological properties of milk, or even cows' physiological status. In this context, and to take advantage of these existing models, the purpose of this work was to evaluate whether a spectral standardization method can enable the use of multiple equations within a network of different FT-MIR spectrometers. The piecewise direct standardization method was used, matching "slave" instruments to a common reference, the "master." The effect of standardization on network reproducibility was assessed on 66 instruments from 3 different brands by comparing the spectral variability of the slaves and the master with and without standardization. With standardization, the global Mahalanobis distance from the slave spectra to the master spectra was reduced on average from 2,655.9 to 14.3, representing a significant reduction of noninformative spectral variability. The transfer of models from instrument to instrument was tested using 3 FT-MIR models predicting (1) the quantity of daily methane emitted by dairy cows, (2) the concentration of polyunsaturated fatty acids in milk, and (3) the fresh cheese yield. The differences, in terms of root mean squared error, between master predictions and slave predictions were reduced after standardization on average from 103 to 17 g/d, from 0.0315 to 0.0045 g/100 mL of milk, and from 2.55 to 0.49 g of curd/100 g of milk, respectively. For all the models, standard deviations of predictions among all the instruments were also reduced by 5.11 times for methane, 5.01 times for polyunsaturated fatty acids, and 7.05 times for fresh cheese yield, showing an improvement of prediction reproducibility within the network. Regarding the results obtained, spectral standardization allows the transfer and use of multiple models on all instruments as well as the improvement of spectral and prediction reproducibility within the network. The method makes the models universal, thereby offering opportunities for data exchange and the creation and use of common robust models at an international level to provide more information to the dairy sector from direct analysis of milk.