Increasing consumer concern exists over the relationship between food composition and human health. Because of the known effects of fatty acids on human health, the development of a quick, inexpensive, and accurate method to directly quantify the fatty acid (FA) composition in milk would be valuable for milk processors to develop a payment system for milk pertinent to their customer requirements and for farmers to adapt their feeding systems and breeding strategies accordingly. The aim of this study was (1) to confirm the ability of mid-infrared spectrometry (MIR) to quantify individual FA content in milk by using an innovative procedure of sampling (i.e., samples were collected from cows belonging to different breeds, different countries, and in different production systems); (2) to compare 6 mathematical methods to develop robust calibration equations for predicting the contents of individual FA in milk; and (3) to test interest in using the FA equations developed in milk as basis to predict FA content in fat without corrections for the slope and the bias of the developed equations. In total, 517 samples selected based on their spectral variability in 3 countries (Belgium, Ireland, and United Kingdom) from various breeds, cows, and production systems were analyzed by gas chromatography (GC). The samples presenting the largest spectral variability were used to calibrate the prediction of FA by MIR. The remaining samples were used to externally validate the 28 FA equations developed. The 6 methods were (1) partial least squares regression (PLS); (2) PLS+repeatability file (REP); (3) first derivative of spectral data+PLS; (4) first derivative+REP+PLS; (5) second derivative of spectral data+PLS; and (6) second derivative+REP+PLS. Methods were compared on the basis of the cross-validation coefficient of determination (R2cv), the ratio of standard deviation of GC values to the standard error of cross-validation (RPD), and the validation coefficient of determination (R2v). The third and fourth methods had, on average, the highest R2cv, RPD, and R2v. The final equations were built using all GC and the best accuracy was observed for the infrared predictions of C4:0, C6:0, C8:0, C10:0, C12:0, C14:0, C16:0, C18:0, C18:1 trans, C18:1 cis-9, C18:1 cis, and for some groups of FA studied in milk (saturated, monounsaturated, unsaturated, short-chain, medium-chain, and long-chain FA). These equations showed R2cv greater than 0.95. With R2cv equal to 0.85, the MIR prediction of polyunsaturated FA could be used to screen the cow population. As previously published, infrared predictions of FA in fat are less accurate than those developed from FA content in milk (g/dL of milk) and no better results were obtained by using milk FA predictions if no corrections for bias and slope based on reference milk samples with known contents of FA were used. These results indicate the usefulness of equations with R2cv greater than 95% in milk payment systems and the usefulness of equations with R2cv greater than 75% for animal breeding purposes.

Interest in the fatty acid composition of dairy products is increasing; however, the measurement of fatty acids requires using gas-liquid chromatography. Although this method is suitable, it involves a time-consuming procedure, expensive reagents, and qualified staff. By comparison, the mid-infrared (MIR) spectrometry method could be a good alternative for assessing the fatty acid profile of dairy products. The objective of this study was to explore the calibration of MIR spectrometry for estimating fatty acid concentrations in milk and milk fat. Estimated concentrations in milk fat were less reliable than those for the same fatty acids in milk. Results also showed that when the fatty acid concentrations in milk increased, the efficiency of the infrared analysis method in predicting these values simultaneously increased. Selected prediction equations must have a high cross-validation coefficient of determination, a high ratio of standard error of cross-validation to standard deviation, and good repeatability of chromatographic data. Results from this study showed that the calibration equations predicting 12:0, 14:0, 16:0, 16:1cis-9, 18:1, and saturated and monounsaturated fatty acids in milk could be used. Thus, with its potential for use in regular milk recording, this infrared analysis method offers the possibility of assessing and improving the quality of milk produced. Indeed, it enables the fatty acid composition in milk to be estimated for each cow and the estimates to be used as indicator traits to determine the genetic values of underlying fatty acid concentrations. The knowledge of these genetic values

The current cattle selection program for dairy cattle in the Walloon region of Belgium does not consider the relative content of the different fatty acids (FA) in milk. However, interest by the local dairy industry in differentiated milk products is increasing. Therefore, farmers may be interested in selecting their animals based on the fat composition. The aim of this study was to evaluate the feasibility of genetic selection to improve the nutritional quality of bovine milk fat. The heritabilities and correlations among milk yield, fat, protein, and major FA contents in milk were estimated. Heritabilities for FA in milk and fat ranged from 5 to 38%. The genetic correlations estimated among FA reflected the common origin of several groups of FA. Given these results, an index including FA contents with the similar metabolic process of production in the mammary gland could be used, for example, to increase the monounsaturated and conjugated fatty acids in milk. Moreover, the genetic correlations between the percentage of fat and the content of C14:0, C12:0, C16:0, and C18:0 in fat were −0.06, 0.55, 0.60, and 0.84, respectively. This result demonstrates that an increase in fat content is not directly correlated with undesirable changes in FA profile in milk for human health. Based on the obtained genetic parameters, a future selection program to improve the FA composition of milk fat could be initiated.

Gene content is the number of copies of a particular allele in a genotype of an animal. Gene content can be used to study additive gene action of candidate gene. Usually genotype data are available only for a part of population and for the rest gene contents have to be calculated based on typed relatives. Methods to calculate expected gene content for animals on large complex pedigrees are relatively complex. In this paper we proposed a practical method to calculate gene content using a linear regression. The method does not estimate genotype probabilities but these can be approximated from gene content assuming Hardy-Weinberg proportions. The approach was compared with other methods on multiple simulated data sets for real bovine pedigrees of 1 082 and 907 903 animals. Different allelic frequencies (0.4 and 0.2) and proportions of the missing genotypes (90, 70, and 50%) were considered in simulation. The simulation showed that the proposed method has similar capability to predict gene content as the iterative peeling method, however it requires less time and can be more practical for large pedigrees. The method was also applied to real data on the bovine myostatin locus on a large dual-purpose Belgian Blue pedigree of 235 133 animals. It was demonstrated that the proposed method can be easily adapted for particular pedigrees.

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.

This study was aimed to evaluate the degree of thermal stress exhibited by Holsteins under a continental temperate climate. Milk, fat, protein, and somatic cell count test-day records collected between 2000 and 2011 from 23,963 cows in 604 herds were combined with meteorological data from 14 public weather stations in Luxembourg. Daily values of 6 different thermal indices (TI) weighted in term of temperature, relative humidity, solar radiation, and wind speed were calculated by averaging hourly TI over 24h. Heat stress thresholds were first identified by a broken-line regression model. Regression models were thereafter applied to quantify milk production losses due to heat stress. The tipping points at which milk and protein yields declined were effectively identified. For fat yield, no valid threshold was identified for any of the studied TI. Daily fat yields tended to decrease steadily with increasing values of TI. Daily somatic cell score patterns were marked by increased values at both lowest and highest TI ranges, with a more pronounced reaction to cold stress for apparent temperature indices. Thresholds differed between TI and traits. For production traits, they ranged from 62 (TI(1)) to 80 (TI(3)) for temperature-humidity indices (THI) and from 16 (TI(5)) to 20 (TI(6)) for apparent temperature indices. Corresponding somatic cell score thresholds were higher and ranged from 66 (TI(1)) to 82 (TI(3)) and from 20 (TI(5)) to 23 (TI(6)), respectively. The largest milk decline per unit of mild, moderate, and extreme heat stress levels of 0.164, 0.356, and 0.955 kg, respectively, was observed when using the conventional THI (TI(1)). The highest yearly milk, fat, and protein losses of 54, 5.7, and 4.2 kg, respectively, were detected by TI(2), the THI index that is adjusted for wind speed and solar radiation. The latter index could be considered as the best indicator of heat stress to be used for forecast and herd management in a first step in temperate regions under anticipated climate changes.

The objective of this study was to assess the phenotypic and genetic variability of production traits and milk fatty acid (FA) contents throughout lactation. Genetic parameters for milk, fat, and protein yields, fat and protein contents, and 19 groups and individual FA contents in milk were estimated for first-parity Holstein cows in the Walloon Region of Belgium using single-trait, test-day animal models and random regressions. Data included 130,285 records from 26,166 cows in 531 herds. Heritabilities indicated that de novo synthesized FA were under stronger genetic control than FA originating from the diet and from body fat mobilization. Estimates for saturated short- and medium-chain individual FA ranged from 0.35 for C4:0 to 0.44 for C8:0, whereas those for monounsaturated long-chain individual FA were lower (around 0.18). Moreover, de novo synthesized FA were more heritable in mid to late lactation. Approximate daily genetic correlations among traits were calculated as correlations between daily breeding values for days in milk between 5 and 305. Averaged daily genetic correlations between milk yield and FA contents did not vary strongly among FA (around -0.35) but they varied strongly across days in milk, especially in early lactation. Results indicate that cows selected for high milk yield in early lactation would have lower de novo synthesized FA contents in milk but a slightly higher content of C18:1 cis-9, indicating that such cows might mobilize body fat reserves. Genetic correlations among FA emphasized the combination of FA according to their origin: contents in milk of de novo FA were highly correlated with each other (from 0.64 to 0.99). Results also showed that genetic correlations between C18:1 cis-9 and other FA varied strongly during the first 100 d in milk and reinforced the statement that the release of long-chain FA inhibits FA synthesis in the mammary gland while the cow is in negative energy balance. Finally, results showed that the FA profile in milk changed during the lactation phenotypically and genetically, emphasizing the relationship between the physiological status of cow and milk composition.

Milk and dairy products are a major source of minerals, particularly calcium, involved in several metabolic functions in humans. Currently, several dairy products are enriched with calcium to prevent osteoporosis. The development of an inexpensive and fast quantitative analysis for minerals is required to offer dairy farmers an opportunity to improve the added value of the produced milk. The aim of this study was to develop 5 equations to measure Ca, K, Mg, Na, and P contents directly in bovine milk using mid-infrared (MIR) spectrometry. A total of 1,543 milk samples were collected between March 2005 and May 2006 from 478 cows during the Walloon milk recording and analyzed by MIR spectrometry. Using a principal component approach, 62 milk samples were selected by their spectral variability and separated in 2 calibration sets. Five outliers were detected and deleted. The mineral contents of the selected samples were measured by inductively coupled plasma atomic emission spectrometry. Using partial least squares combined with a repeatability file, 5 calibration equations were built to estimate the contents of Ca, K, Mg, Na, and P in milk. To assess the accuracy of the developed equations, a full cross-validation and an external validation were performed. The cross-validation coefficients of determination (R 2 cv) were 0.80, 0.70, and 0.79 for Ca, Na, and P, respectively (n = 57), and 0.23 and 0.50 for K and Mg, respectively (n = 31). Only Ca, Na, and P equations showed sufficient R 2 cv for a potential application. These equations were validated using 30 new milk samples. The validation coefficients of determination were 0.97, 0.14, and 0.88 for Ca, Na, and P, respectively, suggesting the potential to use the Ca and P calibration equations. The last 30 samples were added to the initial milk samples and the calibration equations were rebuilt. The R 2 cv for Ca, K, Mg, Na, and P were 0.87, 0.36, 0.65, 0.65, and 0.85, respectively, confirming the potential utilization of the Ca and P equations. Even if new samples should be added in the calibration set, the first results of this study showed the feasibility to quantify the calcium and phosphorus directly in bovine milk using MIR spectrometry.