Genetic relationships between weight loss in early lactation and daily milk production throughout lactation in Holstein cows.

After calving, high-yielding dairy cows mobilize body reserves for energy, sometimes to the detriment of health and fertility. This study aimed to estimate the genetic correlation between body weight loss until nadir and daily milk production (MY24) in first- (L1) and second-lactation (L2) Holstein cows. The data set included 859,020 MY24 records and 570,651 daily raw body weight (BWr) phenotypes from 3,989 L1 cows, and 665,361 MY24 records and 449,449 BWr phenotypes from 3,060 L2 cows, recorded on 36 French commercial farms equipped with milking robots that included an automatic weighing platform. To avoid any bias due to change in digestive content, BWr was adjusted for variations in feed intake, estimated from milk production and BWr. Adjusted body weight was denoted BW. The genetic parameters of BW and MY24 in L1 and L2 cows were estimated using a 4-trait random regression model. In this model, the random effects were fitted by second-order Legendre polynomials on a weekly basis from wk 1 to 44. Nadir of BW was found to be earlier than reported in the literature, at 29 d in milk, and BW loss from calving to nadir was also lower than generally assumed, close to 29 kg. To estimate genetic correlations between body weight loss and production, we defined BWL5 as the loss of weight between wk 1 and 5 after calving. Genetic correlations between BWL5 and MY24 ranged from -0.26 to 0.05 in L1 and from -0.11 to 0.10 in L2, according to days in milk. These moderate to low values suggest that it may be possible to select for milk production without increasing early body mobilization.

Short communication: Confirmation of candidate causative variants on milk composition and cheesemaking properties in Montbéliarde cows.

In a previous study, we identified candidate causative variants located in 24 functional candidate genes for milk protein and fatty acid composition in Montbéliarde, Normande, and Holstein cows. We designed these variants on the custom part of the EuroG10K BeadChip (Illumina Inc., San Diego, CA), which is routinely used for genomic selection analyses in French dairy cattle. To validate the effects of these candidate variants on milk composition and to estimate their effects on cheesemaking properties, a genome-wide association study was performed on milk protein, fatty acid and mineral composition, as well as on 9 cheesemaking traits (3 laboratory cheese yields, 5 coagulation traits, and milk pH). All the traits were predicted from midinfrared spectra in the Montbéliarde cow population of the Franche-Comté region. A total of 194 candidate variants located in 24 genes and 17 genomic regions were imputed on 19,862 cows with phenotypes and genotyped with either the BovineSNP50 (Illumina Inc.) or the EuroG10K BeadChip. We then tested the effect of each SNP in a mixed linear model including random polygenic effects estimated with a genomic relationship matrix. We confirm here the effects of candidate causative variants located in 17 functional candidate genes on both cheesemaking properties and milk composition traits. In each candidate gene, we identified the most plausible causative variant: 4 are missense in the ALPL, SLC26A4, CSN3, and SCD genes, 7 are located in 5'UTR (AGPAT6), 3' untranslated region (GPT), or upstream (CSN1S1, CSN1S2, PAEP, DGAT1, and PICALM) regions, and 6 are located in introns of the SLC37A1, MGST1, CSN2, BRI3BP, FASN, and ANKH genes.

Genetic parameters for cheese-making properties and milk composition predicted from mid-infrared spectra in a large data set of Montbéliarde cows.

Cheese-making properties of pressed cooked cheeses (PCC) and soft cheeses (SC) were predicted from mid-infrared (MIR) spectra. The traits that were best predicted by MIR spectra (as determined by comparison with reference measurements) were 3 measures of laboratory cheese yield, 5 coagulation traits, and 1 acidification trait for PCC (initial pH; pH0PPC). Coefficients of determination of these traits ranged between 0.54 and 0.89. These 9 traits as well as milk composition traits (fatty acid, protein, mineral, lactose, and citrate content) were then predicted from 1,100,238 MIR spectra from 126,873 primiparous Montbéliarde cows. Using this data set, we estimated the corresponding genetic parameters of these traits by REML procedures. A univariate or bivariate repeatability animal model was used that included the fixed effects of herd × test day × spectrometer, stage of lactation, and year × month of calving as well as the random additive genetic, permanent environmental, and residual effects. Heritability estimates varied between 0.37 and 0.48 for the 9 cheese-making property traits analyzed. Coagulation traits were the ones with the highest heritability (0.42 to 0.48), whereas cheese yields and pH0 PPC had the lowest heritability (0.37 to 0.39). Strong favorable genetic correlations, with absolute values between 0.64 and 0.97, were found between different measures of cheese yield, between coagulation traits, between cheese yields and coagulation traits, and between coagulation traits measured for PCC and SC. In contrast, the genetic correlations between milk pH0 PPC and CY or coagulation traits were weak (-0.08 to 0.09). The genetic relationships between cheese-making property traits and milk composition were moderate to high. In particular, high levels of proteins, fatty acids, Ca, P, and Mg in milk were associated with better cheese yields and improved coagulation. Proteins in milk were strongly genetically correlated with coagulation traits and, to a lesser extent, with cheese yields, whereas fatty acids in milk were more genetically correlated with cheese yields than with coagulation traits. This study, carried out on a large scale in Montbéliarde cows, shows that MIR predictions of cheese yields and milk coagulation properties are sufficiently accurate to be used for genetic analyses. Cheese-making traits, as predicted from MIR spectra, are moderately heritable and could be integrated into breeding objectives without additional phenotyping cost, thus creating an opportunity for efficient improvement via selection.

Consequences of the inscription of local breeds in protected designation of origin cow cheese specifications for the genetic management of the herds.

In the mountainous areas of Europe with a humid climate, dairy cattle production is a major agricultural activity, and the milk is often processed into cheese according to protected designation of origin (PDO) specifications. We analyzed the extent to which PDO specifications and/or a mountain environment influence the spatial distribution of estimated breeding values (EBVs) of cows and the herd-year effects (HYEs) for milk yield (kg/lactation) and protein and fat contents (g/kg), as well as lactation ranks and calving months. The study focused on the northern French Alps. A total of 37 023 lactations, recorded in 2006, in 1153 herds were analyzed. The cows belonged to the Montbéliarde (21 516 lactations), Abondance (10 346 lactations) and Tarentaise (5161 lactations) breeds. The two factors of variation considered were the status of the commune where the farm was located in relation to PDO (three categories: area with no PDO, area with a PDO with no milk yield limit, area with a PDO with a milk yield limit) and 'mountain' environment (four categories based on the European regulation: plain, piedmont, mountain and high mountain). In the Abondance breed, the average lactation rank increased with an increase in production constraints due to the PDO or to a mountain environment. In the Abondance and Tarentaise breeds, grouping of calving in winter was most marked in the 'PDO with a milk yield limit' and 'high-mountain' categories. In the Tarentaise breed, no significant effect on any trait and any variable was found in the 'PDO' or 'mountain' categories. In the other two breeds, the average EBV for milk yield decreased with an increase in the constraints due to PDO, with differences of 226 and 93 kg between extreme values in the Abondance and Montbéliarde breeds, respectively. The average HYE for milk yield was higher in the Abondance breed in the 'PDO with no milk yield limit' category than in the other categories (+740 and +1110 kg, respectively); HYE was not affected by the 'PDO' factor in the Montbéliarde breed or by the 'mountain' factor in either breed. Concerning the protein and fat contents, the effect of the 'PDO' and 'mountain' factors depended on the trait, the variable and the breed. The proportion of individual decisions (the farmer makes the decision) v. collective decisions (breed management) concerning herd dynamics in the face of existing constraints is discussed.

Estimation of genetic correlations among countries in international dairy sire evaluations with structural models.

The increase in the number of participating countries and the lack of genetic ties between some countries has lead to statistical and computational difficulties in estimating the genetic (co)variance matrix needed for international sire evaluation of milk yield and other traits. Structural models have been proposed to reduce the number of parameters to estimate by exploiting patterns in the genetic correlation matrix. Genetic correlations between countries are described as a simple function of unspecified country characteristics that can be mapped in a space of limited dimensions. Two link functions equal to the exponential of minus the Euclidian distance between the coordinates of two countries and the exponential of minus the square of this Euclidian distance were used for the study on international simulated and field data. On simulated data, it was shown that structural models might allow an easier estimation of genetic correlations close to the border of the parameter space. This is not always possible with an unstructured model. On milk yield data, genetic correlations obtained from 22 countries for structural models based on 2 and 7 dimensions, respectively, were analyzed. Only a structural model with a large number of axes gave reasonable estimates of genetic correlations compared with correlations obtained for an unstructured model: 76.7% of correlations deviated by less than 0.030. Such a model reduces the number of parameters from 231 genetic correlations to 126 coordinates. On foot angle data, large deviations were observed between genetic correlations estimated with an unstructured model and correlations estimated with a structural model, regardless of the number of axes taken into account.