Catherine Bastin

Genetic parameters of saturated and monounsaturated fatty acid content and the ratio of saturated to unsaturated fatty acids in bovine milk.

Fatty acid composition influences the nutritional quality of milk and the technological properties of butter. Using a prediction of fatty acid (FA) contents by mid-infrared (MIR) spectrometry, a large amount of data concerning the FA profile in bovine milk was collected. The large number of records permitted consideration of more complex models than those used in previous studies. The aim of the current study was to estimate the effects of season and stage of lactation as well as genetic parameters of saturated (SAT) and monounsaturated (MONO) fatty acid contents in bovine milk and milk fat, and the ratio of SAT to unsaturated fatty acids (UNSAT) that reflect the hardness of butter (SAT:UNSAT), using 7 multiple-trait, random-regression test-day models. The relationship between these FA traits with common production traits was also studied. The data set contained 100,841 test-day records of 11,626 Holstein primiparous cows. The seasonal effect was studied based on unadjusted means. These results confirmed that milk fat produced during spring and summer had greater UNSAT content compared with winter (63.13 vs. 68.94% of SAT in fat, on average). The effect of stage of lactation on FA profile was studied using the same methodology. Holstein cows in early first lactation produced a lower content of SAT in their milk fat. Variance components were estimated using a Bayesian method via Gibbs sampling. Heritability of SAT in milk (0.42) was greater than heritability of SAT in milk fat (0.24). Estimates of heritability for MONO were also different in milk and fat (0.14 vs. 0.27). Heritability of SAT:UNSAT was moderate (0.27). For all of these traits, the heritability estimates and the genetic and phenotypic correlations varied through the lactation.

Phenotypic and genetic variability of production traits and milk fatty acid contents across days in milk for Walloon Holstein first-parity cows.

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.

Genetic and environmental relationships between body condition score and milk production traits in Canadian Holsteins.

The objective of this research was to estimate genetic parameters of first-lactation body condition score (BCS), milk yield, fat percentage (Fat%), protein percentage (Prot%), somatic cell score (SCS), milk urea nitrogen (MUN), lactose percentage (Lact%), and fat to protein ratio (F:P) using multiple-trait random regression animal models. Changes in covariances between BCS and milk production traits on a daily basis have not been investigated before and could be useful for determining which BCS estimated breeding values (EBV) might be practical for selection in the future. Field staff from Valacta milk recording agency (Sainte-Anne-de-Bellevue, QC, Canada) collected BCS from Québec herds several times per cow throughout the lactation. Average daily heritabilities and genetic correlations among the various traits were similar to literature values. On an average daily basis, BCS was genetically unfavorably correlated with milk yield (i.e., increased milk yield was associated with lower body condition). The unfavorable genetic correlation between BCS and milk yield became stronger as lactation progressed, but was equivalent to zero for the first month of lactation. Favorable genetic correlations were found between BCS with Prot%, SCS, and Lact% (i.e., greater BCS was associated with greater Prot%, lower SCS, and greater Lact%). These correlations were strongest in early lactation. On an average daily basis, BCS was not genetically correlated with Fat% or MUN, but was negatively correlated with F:P. Furthermore, BCS at 5 and 50 d in milk (DIM) had the most favorable genetic correlations with milk production traits over the lactation (at 5, 50, 150, and 250 DIM). Thus, early lactation BCS EBV shows potential for selection. Regardless, this study showed that the level of association BCS has with milk production traits is not constant over the lactation. Simultaneous selection for both BCS and milk production traits should be considered, mainly due to the unfavorable genetic correlation between BCS with milk yield.

Genetic relationships between body condition score and reproduction traits in Canadian Holstein and Ayrshire first-parity cows

The objective of this study was to investigate the genetic relationship between body condition score (BCS) and reproduction traits for first-parity Canadian Ayrshire and Holstein cows. Body condition scores were collected by field staff several times over the lactation in herds from Québec, and reproduction records (including both fertility and calving traits) were extracted from the official database used for the Canadian genetic evaluation of those herds. For each breed, six 2-trait animal models were run; they included random regressions that allowed the estimation of genetic correlations between BCS over the lactation and reproduction traits that are measured as a single lactation record. Analyses were undertaken on data from 108 Ayrshire herds and 342 Holstein herds. Average daily heritabilities of BCS were close to 0.13 for both breeds; these relatively low estimates might be explained by the high variability among herds and BCS evaluators. Genetic correlations between BCS and interval fertility traits (days from calving to first service, days from first service to conception, and days open) were negative and ranged between -0.77 and -0.58 for Ayrshire and between -0.31 and -0.03 for Holstein. Genetic correlations between BCS and 56-d nonreturn rate at first insemination were positive and moderate. The trends of these genetic correlations over the lactation suggest that a genetically low BCS in early lactation would increase the number of days that the primiparous cow was not pregnant and would decrease the chances of the primiparous cow to conceive at first service. Genetic correlations between BCS and calving traits were generally the strongest at calving and decreased with increasing days in milk. The correlation between BCS at calving and maternal calving ease was 0.21 for Holstein and 0.31 for Ayrshire and emphasized the relationship between fat cows around calving and dystocia. Genetic correlations between calving traits and BCS during the subsequent lactation were moderate and favorable, indicating that primiparous cows with a genetically high BCS over the lactation would have a greater chance of producing a calf that survived (maternal calf survival) and would transmit the genes that allowed the calf to be born more easily (maternal calving ease) and to survive (direct calving ease).

Genetic parameters of milk production traits and fatty acid contents in milk for Holstein cows in parity 1 – 3

The objective of this study was to estimate genetic parameters of milk, fat, and protein yields, fat and protein contents, somatic cell count, and 17 groups and individual milk fatty acid (FA) contents predicted by mid-infrared spectrometry for first-, second- and third-parity Holstein cows. Edited data included records collected in the Walloon region of Belgium from 37,768 cows in parity 1,22,566 cows in parity 2 and 8221 in parity 3. A total of 69 (23 traits for three parities) single-trait random regression animal test-day models were run. Approximate genetic correlations among traits were inferred from pairwise regressions among estimated breeding values of cow having observations. Heritability and genetic correlation estimates from this study reflected the origins of FA: de novo synthetized or originating from the diet and the body fat mobilization. Averaged daily heritabilities of FA contents in milk ranged between 0.18 and 0.47. Average daily genetic correlations (averaged across days in milk and parities) among groups and individual FA contents in milk ranged between 0.31 and 0.99. The genetic variability of FAs in combination with the moderate to high heritabilities indicated that FA contents in milk could be changed by genetic selection; however, desirable direction of change in these traits remains unclear and should be defined with respect to all issues of importance related to milk FA.

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 correlations of days open with production traits and contents in milk of major fatty acids predicted by mid-infrared spectrometry

The objective of this study was to estimate the genetic relationships between days open (DO) and both milk production traits and fatty acid (FA) content in milk predicted by mid-infrared spectrometry. The edited data set included 143,332 FA and production test-day records and 29,792 DO records from 29,792 cows in 1,170 herds. (Co)variances were estimated using a series of 2-trait models that included a random regression for milk production and FA traits. In contrast to the genetic correlations with fat content, those between DO and FA content in milk changed considerably over the lactation. The genetic correlations with DO for unsaturated FA, monounsaturated FA, long-chain FA, C18:0, and C18:1 cis-9 were positive in early lactation but negative after 100 d in milk. For the other FA, genetic correlations with DO were negative across the whole lactation. At 5 d in milk, the genetic correlation between DO and C18:1 cis-9 was 0.39, whereas the genetic correlations between DO and C6:0 to C16:0 FA ranged from -0.37 to -0.23. These results substantiated the known relationship between fertility and energy balance status, explained by the release of long-chain FA in early lactation, from the mobilization of body fat reserves, and the consequent inhibition of de novo FA synthesis in the mammary gland. At 200 d in milk, the genetic correlations between DO and FA content ranged from -0.38 for C18:1 cis-9 to -0.03 for C6:0. This research indicates an opportunity to use FA content in milk as an indicator trait to supplement the prediction of genetic merit for fertility.

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.

Genetic analysis of heat stress effects on yield traits, udder health, and fatty acids of Walloon Holstein cows

Genetic parameters that considered tolerance for heat stress were estimated for production, udder health, and milk composition traits. Data included 202,733 test-day records for milk, fat, and protein yields, fat and protein percentages, somatic cell score (SCS), 10 individual milk fatty acids (FA) predicted by mid-infrared spectrometry, and 7 FA groups. Data were from 34,468 first-lactation Holstein cows in 862 herds in the Walloon region of Belgium and were collected between 2007 and 2010. Test-day records were merged with daily temperature-humidity index (THI) values based on meteorological records from public weather stations. The maximum distance between each farm and its corresponding weather station was 21km. Linear reaction norm models were used to estimate the intercept and slope responses of 23 traits to increasing THI values. Most yield and FA traits had phenotypic and genetic declines as THI increased, whereas SCS, C18:0, C18:1 cis-9, and 4 FA groups (unsaturated FA, monounsaturated FA, polyunsaturated FA, and long-chain FA) increased with THI. Moreover, the latter traits had the largest slope-to-intercept genetic variance ratios, which indicate that they are more affected by heat stress at high THI levels. Estimates of genetic correlations within trait between cold and hot environments were generally high (>0.80). However, lower estimates (<=0.67) were found for SCS, fat yield, and C18:1 cis-9, indicating that animals with the highest genetic merit for those traits in cold environments do not necessarily have the highest genetic merit for the same traits in hot environments. Among all traits, C18:1 cis-9 was the most sensitive to heat stress. As this trait is known to reflect body reserve mobilization, using its variations under hot conditions could be a very affordable milk biomarker of heat stress for dairy cattle expressing the equilibrium between intake and mobilization under warm conditions.

Modeling milk urea of Walloon dairy cows in management perspectives.

The aim of this study was to develop an adapted random regression test-day model for milk urea (MU) and to study the possibility of using predictions and solutions given by the model for management purposes. Data included 607,416 MU test-day records of first-lactation cows from 632 dairy herds in the Walloon Region of Belgium. Several advanced features were used. First, to detect the herd influence, the classical herd x test-day effect was split into 3 new effects: a fixed herd x year effect, a fixed herd x month-period effect, and a random herd test-day effect. A fixed time period regression was added in the model to take into account the yearly oscillations of MU on a population scale. Moreover, first autoregressive processes were introduced and allowed us to consider the link between successive test-day records. The variance component estimation indicated that large variance was associated with the random herd x test-day effect (48% of the total variance), suggesting the strong influence of herd management on the MU level. The heritability estimate was 0.13. By comparing observed and predicted MU levels at both the individual and herd levels, target ranges for MU concentrations were defined to take into account features of each cow and each herd. At the cow level, an MU record was considered as deviant if it was <200 or >400 mg/L (target range used in the field) and if the prediction error was >50 mg/L (indicating a significant deviation from the expected level). Approximately 7.5% of the MU records collected between June 2007 and May 2008 were beyond these thresholds. This combination allowed for the detection of potentially suspicious cows. At the herd level, the expected MU level was considered as the sum of the solutions for specific herd effects. A herd was considered as deviant from its target range when the prediction error was greater than the standard deviation of MU averaged by herd test day. Results showed that 6.7% of the herd test-day MU levels between June 2007 and May 2008 were considered deviant. These deviations seemed to occur more often during the grazing period. Although theoretical considerations developed in this study should be validated in the field, this research showed the potential use of a test-day model for analyzing functional traits to advise dairy farmers.