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.

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 (CH4) emissions from ruminant livestock, in particular from dairy cows, have acquired global importance for sustainable, climate-smart cattle production. Based on CH4 reference measurements obtained with the SF6 tracer technique to determine ruminal CH4 production, a current equation permits evaluation of individual daily CH4 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 SF6 to measure CH4 production from cattle. This study aimed to develop an equation that allows estimating CH4 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 SF6 reference measurements regarding the potential to estimate daily CH4 emissions of dairy cows from milk FT-MIR spectra. A total of 584 RC reference CH4 measurements (mean ± standard deviation of 400 ± 72 g of CH4/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 SF6 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 SF6-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 SF6 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 CH4 production, which is 23 g/d lower than the RMSE for the SF6-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 CH4/d was estimated, which is consistent with the range of the values measured using both RC and SF6 techniques. This study confirmed that milk FT-MIR spectra could be used as a potential proxy to estimate daily CH4 emissions from dairy cows provided that the variability to predict is covered by the model.

Relationships between milk mid-IR predicted gastro-enteric methane production and the technical and financial performance of commercial dairy herds.

Considering economic and environmental issues is important in ensuring the sustainability of dairy farms. The objective of this study was to investigate univariate relationships between lactating dairy cow gastro-enteric methane (CH4) production predicted from milk mid-IR (MIR) spectra and technico-economic variables by the use of large scale and on-farm data. A total of 525 697 individual CH4 predictions from milk MIR spectra (MIR-CH4 (g/day)) of milk samples collected on 206 farms during the Walloon milk recording scheme were used to create a MIR-CH4 prediction for each herd and year (HYMIR-CH4). These predictions were merged with dairy herd accounting data. This allowed a simultaneous study of HYMIR-CH4 and 42 technical and economic variables for 1024 herd and year records from 2007 to 2014. Pearson correlation coefficients (r) were used to assess significant relationships (P<0.05). Low HYMIR-CH4 was significantly associated with, amongst others, lower fat and protein corrected milk (FPCM) yield (r=0.18), lower milk fat and protein content (r=0.38 and 0.33, respectively), lower quantity of milk produced from forages (r=0.12) and suboptimal reproduction and health performance (e.g. longer calving interval (r=-0.21) and higher culling rate (r=-0.15)). Concerning economic results, low HYMIR-CH4 was significantly associated with lower gross margin per cow (r=0.19) and per litre FPCM (r=0.09). To conclude, this study suggested that low lactating dairy cow gastro-enteric CH4 production tended to be associated with more extensive or suboptimal management practices, which could lead to lower profitability. The observed low correlations suggest complex interactions between variables due to the use of on-farm data with large variability in technical and management practices.

The feeding system impacts relationships between calving interval and economic results of dairy farms.

The calving interval (CI) can potentially impact the economic results of dairy farms. This study highlighted the most profitable CI and innovated by describing this optimum as a function of the feeding system of the farm. On-farm data were used to represent real farm conditions. A total of 1832 accounts of farms recorded from 2007 to 2014 provided economic, technical and feeding information per herd and per year. A multiple correspondence analysis created four feeding groups: extensive, low intensive, intensive and very intensive herds. The gross margin and some of its components were corrected to account for the effect of factors external to the farm, such as the market, biological status, etc. Then the corrected gross margin (cGMc) and its components were modelled by CI parameters in each feeding system by use of GLM. The relationship between cGMc and the proportion of cows with CI<380 days in each feeding group showed that keeping most of the cows in the herd with CI near to 1 year was not profitable for most farms (for the very intensive farms there was no effect of the proportion). Moreover, a low proportion of cows (0% to 20%) with a near-to-1-year CI was not profitable for the extensive and low intensive farms. Extending the proportion of cows with CI beyond 459 days until 635 days (i.e. data limitation) caused no significant economic loss for the extensive and low intensive farms, but was not profitable for the intensive and very intensive farms. Variations of the milk and feeding components explained mainly these significant differences of gross margin. A link between the feeding system and persistency, perceptible in the milk production and CI shown by the herd, could explain the different relationships observed between the extent of CI and the economic results in the feeding groups. This herd-level study tended to show different economic optima of CI as a function of the feeding system. A cow-level study would specify these tendencies to give CI objectives to dairy breeders as a function of their farm characteristics.

Bayesian single-step genomic evaluations combining local and foreign information in Walloon Holsteins.

Most dairy cattle populations found in different countries around the world are small to medium sized and use many artificial insemination bulls imported from different foreign countries. The Walloon population in the southern part of Belgium is a good example for such a small-scale population. Wallonia has also a very active community of Holstein breeders requesting high level genetic evaluation services. Single-step Genomic BLUP (ssGBLUP) methods allow the simultaneous use of genomic, pedigree and phenotypic information and could reduce potential biases in the estimation of genomically enhanced breeding values (GEBV). Therefore, in the context of implementing a Walloon genomic evaluation system for Holsteins, it was considered as the best option. However, in contrast to multi-step genomic predictions, natively ssGBLUP will only use local phenotypic information and is unable to use directly important other sources of information coming from abroad, for example Multiple Across Country Evaluation (MACE) results as provided by the Interbull Center (Uppsala, Sweden). Therefore, we developed and implemented single-step Genomic Bayesian Prediction (ssGBayes), as an alternative method for the Walloon genomic evaluations. The ssGBayes method approximated the correct system of equations directly using estimated breeding values (EBV) and associated reliabilities (REL) without any explicit deregression step. In the Walloon genomic evaluation, local information refers to Walloon EBV and REL and foreign information refers to MACE EBV and associated REL. Combining simultaneously all available genotypes, pedigree, local and foreign information in an evaluation can be achieved but adding contributions to left-hand and right-hand sides subtracting double-counted contributions. Correct propagation of external information avoiding double counting of contributions due to relationships and due to records can be achieved. This ssGBayes method computed more accurate predictions for all types of animals. For example, for genotyped animals with low Walloon REL (<0.25) without MACE results but sired by genotyped bulls with MACE results, the average increase of REL for the studied traits was 0.38 points of which 0.08 points could be traced to the inclusion of MACE information. For other categories of genotyped animals, the contribution by MACE information was also high. The Walloon genomic evaluation system passed for the first time the Interbull GEBV tests for several traits in July 2013. Recent experiences reported here refer to its use in April 2016 for the routine genomic evaluations of milk production, udder health and type traits. Results showed that the proposed methodology should also be of interest for other, similar, populations.

Genetic evaluation for birth and conformation traits in dual-purpose Belgian Blue cattle using a mixed inheritance model.

The segregation of the causal mutation () in the muscular hypertrophy gene in dual-purpose Belgian Blue (dpBB) cattle is considered to result in greater calving difficulty (dystocia). Establishing adapted genetic evaluations might overcome this situation through efficient selection. However, the heterogeneity of dpBB populations at the locus implies separating the major gene and other polygenic effects in complex modeling. The use of mixed inheritance models may be an interesting option because they simultaneously assume both influences. A genetic evaluation in dpBB based on a mixed inheritance model was developed for birth and conformation traits: gestation length (GL), calving difficulty (CD), birth weight (BiW), and body conformation score (BC). A total of 27,362 animals having records were used for analyses. The total number of animals in the pedigree used to build the numerator relationship matrix was 62,617. Genotypes at the locus were available for 2,671 animals. Missing records at this locus were replaced with genotype probabilities. A total of 13,221 (48.3%) were registered as dpBB, 1,287 (4.7%) as beef Belgian Blue, and 12,854 (47.0%) were unknown. From those 13,221 dpBB animals, 650, 849, and 534 had double or single copies or no copy, respectively, of the causal mutation () in the muscular hypertrophy gene, whereas 11,188 had missing genotypes. This heterogeneity at the locus may be the reason for high variability in the studied traits, that is, high heritability estimates of 0.33, 0.30, 0.38, and 0.43 for GL, CD, BiW, and BC, respectively. In general, additive ( < 0.05) and dominance ( < 0.001) allele substitution for calves and dams had significant impact for all traits. The moderate coefficient of genetic variation (27.80%) and high direct heritability (0.28) for CD suggested genetic variability in dpBB and possible genetic improvement through selection. This variability has allowed dpBB breeders to successfully apply mass selection in the past. Genetic trend means from 1988 to 2016 showed that sire selection for CD within genotype was progressively applied by breeders. The selection intensity was more important for CD in double-muscled lines than in segregated lines. Our study illustrated the possible confusion caused by the use of major genes in selection and the importance of fitting appropriate models such as mixed inheritance models that combine polygenic and gene content information.

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.

Assessing the effect of pregnancy stage on milk composition of dairy cows using mid-infrared spectra.

Changes in milk production traits (i.e., milk yield, fat, and protein contents) with the pregnancy stage are well documented. To our knowledge, the effect of pregnancy on the detailed milk composition has not been studied so far. The mid-infrared (MIR) spectrum reflects the detailed composition of a milk sample and is obtained by a nonexhaustive and widely used method for milk analysis. Therefore, this study aimed to investigate the effect of pregnancy on milk MIR spectrum in addition to milk production traits (milk yield, fat, and protein contents). A model including regression on the number of days pregnant was applied on milk production traits (milk yield, fat, and protein contents) and on 212 spectral points from the MIR spectra of 9,757 primiparous Holstein cows from Walloon herds. Effects of pregnancy stage were expressed on a relative scale (effect divided by the squared root of the phenotypic variance); this allowed comparisons between effects on milk traits and on 212 spectral points. Effect of pregnancy stage on production traits were in line with previous studies indicating that the model accounted well for the pregnancy effect. Trends of the relative effect of the pregnancy stage on the 212 spectral points were consistent with known and observed effect on milk traits. The highest effect of the pregnancy was observed in the MIR spectral region from 968 to 1,577 cm-1. For some specific wavenumbers, the effect was higher than for fat and protein contents in the beginning of the pregnancy (from 30 to 90 or 120 d pregnant). In conclusion, the effect of early pregnancy can be observed in the detailed milk composition through the analysis of the MIR spectrum of bovine milk. Further analyses are warranted to explore deeply the use of MIR spectra of bovine milk for breeding and management of dairy cow pregnancy.

Development of Fourier transform mid-infrared calibrations to predict acetone, ß-hydroxybutyrate, and citrate contents in bovine milk through a European dairy network.

To manage negative energy balance and ketosis in dairy farms, rapid and cost-effective detection is needed. Among the milk biomarkers that could be useful for this purpose, acetone and ß-hydroxybutyrate (BHB) have been proved as molecules of interest regarding ketosis and citrate was recently identified as an early indicator of negative energy balance. Because Fourier transform mid-infrared spectrometry can provide rapid and cost-effective predictions of milk composition, the objective of this study was to evaluate the ability of this technology to predict these biomarkers in milk. Milk samples were collected in commercial and experimental farms in Luxembourg, France, and Germany. Acetone, BHB, and citrate contents were determined by flow injection analysis. Milk mid-infrared spectra were recorded and standardized for all samples. After edits, a total of 548 samples were used in the calibration and validation data sets for acetone, 558 for BHB, and 506 for citrate. Acetone content ranged from 0.020 to 3.355mmol/L with an average of 0.103mmol/L; BHB content ranged from 0.045 to 1.596mmol/L with an average of 0.215mmol/L; and citrate content ranged from 3.88 to 16.12mmol/L with an average of 9.04mmol/L. Acetone and BHB contents were log-transformed and a part of the samples with low values was randomly excluded to approach a normal distribution. The 3 edited data sets were then randomly divided into a calibration data set (3/4 of the samples) and a validation data set (1/4 of the samples). Prediction equations were developed using partial least square regression. The coefficient of determination (R(2)) of cross-validation was 0.73 for acetone, 0.71 for BHB, and 0.90 for citrate with root mean square error of 0.248, 0.109, and 0.70mmol/L, respectively. Finally, the external validation was performed and R(2) obtained were 0.67 for acetone, 0.63 for BHB, and 0.86 for citrate, with respective root mean square error of validation of 0.196, 0.083, and 0.76mmol/L. Although the practical usefulness of the equations developed should be further verified with other field data, results from this study demonstrated the potential of Fourier transform mid-infrared spectrometry to predict citrate content with good accuracy and to supply indicative contents of BHB and acetone in milk, thereby providing rapid and cost-effective tools to manage ketosis and negative energy balance in dairy farms.

Capitalizing on fine milk composition for breeding and management of dairy cows.

The challenge of managing and breeding dairy cows is permanently adapting to changing production circumstances under socio-economic constraints. If managing and breeding address different timeframes of action, both need relevant phenotypes that allow for precise monitoring of the status of the cows, and their health, behavior, and well-being as well as their environmental impact and the quality of their products (i.e., milk and subsequently dairy products). Milk composition has been identified as an important source of information because it could reflect, at least partially, all these elements. Major conventional milk components such as fat, protein, urea, and lactose contents are routinely predicted by mid-infrared (MIR) spectrometry and have been widely used for these purposes. But, milk composition is much more complex and other nonconventional milk components, potentially predicted by MIR, might be informative. Such new milk-based phenotypes should be considered given that they are cheap, rapidly obtained, usable on a large scale, robust, and reliable. In a first approach, new phenotypes can be predicted from MIR spectra using techniques based on classical prediction equations. This method was used successfully for many novel traits (e.g., fatty acids, lactoferrin, minerals, milk technological properties, citrate) that can be then useful for management and breeding purposes. An innovation was to consider the longitudinal nature of the relationship between the trait of interest and the MIR spectra (e.g., to predict methane from MIR). By avoiding intermediate steps, prediction errors can be minimized when traits of interest (e.g., methane, energy balance, ketosis) are predicted directly from MIR spectra. In a second approach, research is ongoing to detect and exploit patterns in an innovative manner, by comparing observed with expected MIR spectra directly (e.g., pregnancy). All of these traits can then be used to define best practices, adjust feeding and health management, improve animal welfare, improve milk quality, and mitigate environmental impact. Under the condition that MIR data are available on a large scale, phenotypes for these traits will allow genetic and genomic evaluations. Introduction of novel traits into the breeding objectives will need additional research to clarify socio-economic weights and genetic correlations with other traits of interest.

Integration of external estimated breeding values and associated reliabilities using correlations among traits and effects.

Based on a Bayesian view of linear mixed models, several studies showed the possibilities to integrate estimated breeding values (EBV) and associated reliabilities (REL) provided by genetic evaluations performed outside a given evaluation system into this genetic evaluation. Hereafter, the term "internal" refers to this given genetic evaluation system, and the term "external" refers to all other genetic evaluations performed outside the internal evaluation system. Bayesian approaches integrate external information (i.e., external EBV and associated REL) by altering both the mean and (co)variance of the prior distributions of the additive genetic effects based on the knowledge of this external information. Extensions of the Bayesian approaches to multivariate settings are interesting because external information expressed on other scales, measurement units, or trait definitions, or associated with different heritabilities and genetic parameters than the internal traits, could be integrated into a multivariate genetic evaluation without the need to convert external information to the internal traits. Therefore, the aim of this study was to test the integration of external EBV and associated REL, expressed on a 305-d basis and genetically correlated with a trait of interest, into a multivariate genetic evaluation using a random regression test-day model for the trait of interest. The approach we used was a multivariate Bayesian approach. Results showed that the integration of external information led to a genetic evaluation for the trait of interest for, at least, animals associated with external information, as accurate as a bivariate evaluation including all available phenotypic information. In conclusion, the multivariate Bayesian approaches have the potential to integrate external information correlated with the internal phenotypic traits, and potentially to the different random regressions, into a multivariate genetic evaluation. This allows the use of different scales, heritabilities, variance components, measurement units, or trait definitions for external and internal traits. However, one possible issue for implementing multivariate Bayesian approaches could be the availability or estimation of genetic correlations between external and internal traits.

The use of SWOT analysis to explore and prioritize conservation and development strategies for local cattle breeds.

SWOT (Strengths, Weaknesses, Opportunities and Threats) analysis is a tool widely used to help in decision making in complex systems. It suits to exploring the issues and measures related to the conservation and development of local breeds, as it allows the integration of many driving factors influencing breed dynamics. We developed a quantified SWOT method as a decision-making tool for identification and ranking of conservation and development strategies of local breeds, and applied it to a set of 13 cattle breeds of six European countries. The method has four steps: definition of the system, identification and grouping of the driving factors, quantification of the importance of driving factors and identification and prioritization of the strategies. The factors were determined following a multi-stakeholder approach and grouped with a three-level structure. Animal genetic resources expert groups ranked the factors, and a quantification process was implemented to identify and prioritize strategies. The proposed SWOT methodology allows analyzing the dynamics of local cattle breeds in a structured and systematic way. It is a flexible tool developed to assist different stakeholders in defining the strategies and actions. The quantification process allows the comparison of the driving factors and the prioritization of the strategies for the conservation and development of local cattle breeds. We identified 99 factors across the breeds. Although the situation is very heterogeneous, the future of these breeds may be promising. The most important strengths and weaknesses were related to production systems and farmers. The most important opportunities were found in marketing new products, whereas the most relevant threats were found in selling the current products. The across-breed strategies utility decreased as they gained specificity. Therefore, the strategies at European level should focus on general aspects and be flexible enough to be adapted to the country and breed specificities.

Farmer's views and values to focus on cattle conservation policies: the case of eight European countries.

Our aim was to identify elements useful in designing policies and programmes for conservation of farm animal genetic resources, taking as case study a group of European local cattle breeds. We first investigated the implications of differences among countries in the policies and programmes to be developed. Secondly, we analysed key elements common to countries, which may affect local breed viability. We used the herd size trend expected by the farmer in the near future as an indicator of breed viability. Fifteen breeds, for a total of 355 farms, were surveyed. To take into account the multiple factors influencing breeds' demographic trends, the questionnaire included economical, technical and social aspects. Among the major differences across countries was the perception of the farmer on the value attributed to the local breed by society. Concerning the elements common to countries and their association to breed viability, the greater the collaboration among farmers and the stakeholders' appreciation as perceived by the farmer, the greater the viability of the farm. An opposite trend was observed for the age of the farmer. Older farmers generally planned to soon cease farming or decrease herd size, whereas young farmers planned to increase the size of their herds. Implications of including these elements in conservation policies are discussed.

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

Lactoferrin (LTF) is a milk glycoprotein favorably associated with the immune system of dairy cows. Somatic cell count is often used as an indicator of mastitis in dairy cows, but knowledge on the milk LTF content could aid in mastitis detection. An inexpensive, rapid and robust method to predict milk LTF is required. The aim of this study was to develop an equation to quantify the LTF content in bovine milk using mid-infrared (MIR) spectrometry. LTF was quantified by enzyme-linked immunosorbent assay (ELISA), and all milk samples were analyzed by MIR. After discarding samples with a coefficient of variation between 2 ELISA measurements of more than 5% and the spectral outliers, the calibration set consisted of 2499 samples from Belgium (n = 110), Ireland (n = 1658) and Scotland (n = 731). Six statistical methods were evaluated to develop the LTF equation. The best method yielded a cross-validation coefficient of determination for LTF of 0.71 and a cross-validation standard error of 50.55 mg/l of milk. An external validation was undertaken using an additional dataset containing 274 Walloon samples. The validation coefficient of determination was 0.60. To assess the usefulness of the MIR predicted LTF, four logistic regressions using somatic cell score (SCS) and MIR LTF were developed to predict the presence of mastitis. The dataset used to build the logistic regressions consisted of 275 mastitis records and 13 507 MIR data collected in 18 Walloon herds. The LTF and the interaction SCS × LTF effects were significant (P < 0.001 and P = 0.02, respectively). When only the predicted LTF was included in the model, the prediction of the presence of mastitis was not accurate despite a moderate correlation between SCS and LTF (r = 0.54). The specificity and the sensitivity of models were assessed using Walloon data (i.e. internal validation) and data collected from a research herd at the University of Wisconsin - Madison (i.e. 5886 Wisconsin MIR records related to 93 mastistis events - external validation). Model specificity was better when LTF was included in the regression along with SCS when compared with SCS alone. Correct classification of non-mastitis records was 95.44% and 92.05% from Wisconsin and Walloon data, respectively. The same conclusion was formulated from the Hosmer and Lemeshow test. In conclusion, this study confirms the possibility to quantify an LTF indicator from milk MIR spectra. It suggests the usefulness of this indicator associated to SCS to detect the presence of mastitis. Moreover, the knowledge of milk LTF could also improve the milk nutritional quality.

Genetic analysis of lactoferrin content in bovine milk.

Bovine lactoferrin (LF) is mainly present in milk and shows important physiological and biological functions. The aim of this study was to estimate the heritability and correlation values of LF content in bovine milk with different economic traits as milk yield (MY), fat and protein percentages, and somatic cell score (SCS). Variance components of the studied traits were estimated by REML using a multiple-trait mixed model. The obtained heritability (0.22) for LF content predicted using mid-infrared spectrometry (pLF) suggested the possibility of animal selection based on the increase of LF content in milk. The phenotypic and genetic correlation values calculated between pLF and SCS were moderate (0.31 and 0.24, respectively). Furthermore, a preliminary study of bovine LF gene polymorphism effects was performed on the same production traits. By PCR, all exons of the LF gene were amplified and then sequenced. Three new polymorphisms were detected in exon 2, exon 11, and intron 8. We examined the effects of LF gene polymorphisms of exons 2, 4, 9, 11, and 15, and intron 8 on pLF, MY, fat and protein percentages, and SCS. The different observed effects did not reach a significant level probably because of the characteristics of the studied population. However, the results were promising, and LF may be a potential indicator of mastitis. Further studies are necessary to evaluate the effect of genetic selection based on LF content on the improvement of mastitis resistance.

Genomic location of the bovine growth hormone secretagogue receptor (GHSR) gene and investigation of genetic polymorphism.

The growth hormone secretagogue receptor (GHSR) is involved in the regulation of energetic homeostasis and GH secretion. In this study, the bovine GHSR gene was mapped to BTA1 between BL26 and BMS4004. Two different bovine GHSR CDS (GHSR1a and GHSR1b) were sequenced. Six polymorphisms (five SNPs and one 3-bp indel) were also identified, three of them leading to amino acid variations L24V, D194N, and Del R242. These variations are located in the extracellular N-terminal end, the exoloop 2, and the cytoloop 3 of the receptor, respectively.

Genetic variability of lactoferrin content estimated by mid-infrared spectrometry in bovine milk.

The effects of lactoferrin (LF) on the immune system have already been shown by many studies. Unfortunately, the current methods used to measure LF levels in milk do not permit the study of the genetic variability of lactoferrin or the performance of routine genetic evaluations. The first aim of this research was to derive a calibration equation permitting the prediction of LF in milk by mid-infrared spectrometry (MIR). The calibration with partial least squares on 69 samples showed a ratio of standard error of cross-validation to standard deviation equal to 1.98. Based on this value, the calibration equation was used to establish an LF indicator trait (predicted LF; pLF) on a large number of milk samples (n = 7,690). A subsequent study of its variability was conducted, which confirmed that stage of lactation and lactation number influence the overall pLF level. Small differences in mean pLF among 7 dairy breeds were also observed. The pLF content of Jersey milk was significantly higher than that in Holstein milk. Therefore, the choice of breed could change the expected LF level. Heritability estimated for pLF was 19.7%. The genetic and phenotypic correlations between somatic cell score and pLF were 0.04 and 0.26, respectively. As somatic cell score increases in presence of mastitis, this observation seems to indicate that pLF, or a function of observed pLF, compared with expected LF might have potential as an indicator of mastitis. The negative genetic correlation (-0.36) between milk yield and pLF could indicate an undesirable effect of selection for high milk production on the overall LF level.