Short communication: genetic relationships among daily energy balance, feed intake, body condition score, and fat to protein ratio of milk in dairy cows

Postpartum energy status is critically important to health and fertility, and it remains a major task to find suitable indicator traits for energy balance. Therefore, genetic parameters for daily energy balance (EB) and dry matter intake (DMI), weekly milk fat to protein ratio (FPR), and monthly body condition score (BCS) were estimated using random regression on data collected from 682 Holstein-Friesian primiparous cows recorded between lactation d 11 to 180. Average energy-corrected milk (ECM), EB, DMI, BCS, and FPR were 32.0 kg, 9.6 MJ of NE_L, 20.3 kg, 2.95, and 1.12, respectively. Heritability estimates for EB, DMI, BCS, and FPR ranged from 0.03 to 0.13, 0.04 to 0.19, 0.34 to 0.59, and 0.20 to 0.54. Fat to protein ratio was a more valid measure for EB in early lactation than DMI, BCS, or single milk components. Correlations between FPR and EB were highest at the beginning of lactation [genetic correlation (r_g) = −0.62 at days in milk (DIM) 15] and decreased toward zero. Dry matter intake was the trait most closely correlated with EB in mid lactation (r_g = 0.73 at DIM 120 and 150). Energy balance in early lactation was negatively correlated to EB in mid lactation. The same applied to DMI. Genetic correlations between FPR across lactation stages were all positive; the lowest genetic correlation (0.55) was estimated between the beginning of lactation and early mid lactation. Hence, to improve EB at the beginning of lactation, EB and indicator traits need to be recorded in early lactation. We concluded that FPR is an adequate indicator for EB during the energy deficit phase. Genetic correlations of FPR with ECM, fat percentage, and protein percentage showed that a reduction of FPR in early lactation would have a slightly negative effect on ECM, whereas milk composition would change in a desirable manner.     

Genetic and genomic dissection of dry matter intake at different lactation stages in primiparous Holstein cows

Dry matter intake (DMI) and feed efficiency are economically relevant traits. Simultaneous selection for low DMI and high milk yield might improve feed efficiency, but bears the risk of aggravating the negative energy balance and related health problems in early lactation. Lactation stage-specific selection might provide a possibility to optimize the trajectory of DMI across days in milk (DIM), but requires in-depth knowledge about genetic parameters within and across lactation stages. Within the current study, daily heritabilities and genetic correlations between DMI records from different lactation stages were estimated using random regression models based on 910 primiparous Holstein cows. The heritability estimates from DIM 11 to 180 follow a slightly parabolic curve varying from 0.26 (DIM 121) to 0.37 (DIM 11 and 180). Genetic correlations estimated between DIM 11, 30, 80, 130, and 180 were all positive, ranging from 0.29 (DIM 11 and 180) to 0.97 (DIM 11 and 30; i.e., the correlations are inversely related to the length of the interval between compared DIM). Deregressed estimated breeding values for the same lactation days were used as phenotypes in sequential genome-wide association studies using 681 cows drawn from the study population and genotyped for the Illumina SNP50 BeadChip (Illumina Inc., San Diego, CA). A total of 21 SNP on 10 chromosomes exceeded the chromosome-wise significance threshold for at least 1 analyzed DIM, pointing to some interesting candidate genes directly involved in the regulation of feed intake. Association signals were restricted to certain lactation stages, thus supporting the genetic correlations. Partitioning the explained variance onto chromosomes revealed a large contribution of Bos taurus autosome 7 not harboring any associated marker in the current study. The results contribute to the knowledge about the genetic architecture of the complex phenotype DMI and might provide valuable information for future selection efforts.     

Analysis of the genetic architecture of energy balance and its major determinants dry matter intake and energy-corrected milk yield in primiparous Holstein cows

The focus of modern dairy cow breeding programs has shifted from being mainly yield based toward balanced goals that increasingly consider functional traits such as fertility, metabolic stability, and longevity. To improve these traits, a less pronounced energy deficit postpartum is considered a key challenge. On the other hand, feed efficiency and methane emissions are gaining importance, possibly leading to conflicts in the design of breeding goals. Dry matter intake (DMI) is one of the major determinants of energy balance (EB), and recently some efforts were undertaken to include DMI in genomic breeding programs. However, there is not yet a consensus on how this should be achieved as there are different goals in the course of lactation (i.e., reducing energy deficit postpartum vs. subsequently improving feed efficiency). Thus, the aim of this study was to gain more insight into the genetic architecture of energy metabolism across lactation by genetically dissecting EB and its major determinants DMI and energy-corrected milk (ECM) yield at different lactation stages applying random regression methodology and univariate and multivariate genomic analyses to data from 1,174 primiparous Holstein cows. Daily heritability estimates ranged from 0.29 to 0.49, 0.26 to 0.37, and 0.58 to 0.68 for EB, DMI, and ECM, respectively, across the first 180 d in milk (DIM). Genetic correlations between ECM and DMI were positive, ranging from 0.09 (DIM 11) to 0.36 (DIM 180). However, ECM and EB were negatively correlated (r_g = ?0.26 to ?0.59). The strongest relationship was found at the onset of lactation, indicating that selection for increased milk yield at this stage will result in a more severe energy deficit postpartum. The results also indicate that EB is more affected by DMI (r_g = 0.71 to 0.81) than by its other major determinant, ECM. Thus, breeding for a higher DMI in early lactation seems to be a promising strategy to improve the energy status of dairy cows. We found evidence that genetic regulation of energy homeostasis is complex, with trait- and lactation stage-specific quantitative trait loci suggesting that the trajectories of the analyzed traits can be optimized as mentioned above. Especially from the multivariate genomic analyses, we were able to draw some conclusions on the mechanisms involved and identified the genes encoding fumarate hydratase and adiponectin as highly promising candidates for EB, which will be further analyzed.     

Genetic parameters for energy balance, feed efficiency, and related traits in Holstein cattle

Objectives of the current study were to estimate genetic parameters in Holstein cows for energy balance (EB) and related traits including dry matter intake (DMI), body weight (BW), body condition score (BCS), energy-corrected milk (ECM) production, and gross feed efficiency (GFE), defined as the ratio of total ECM yield to total DMI over the first 150d of lactation. Data were recorded for the first half of lactation on 227 and 175 cows in their first or later lactation, respectively. Random regression models were fitted to longitudinal data. Also, each trait was averaged over monthly intervals and analyzed by single and multivariate animal models. Heritability estimates ranged from 0.27 to 0.63, 0.12 to 0.62, 0.12 to 0.49, 0.63 to 0.72, and 0.49 to 0.53 for DMI, ECM yield, EB, BW, and BCS, respectively, averaged over monthly intervals. Daily heritability estimates ranged from 0.18 to 0.30, 0.10 to 0.26, 0.07 to 0.22, 0.43 to 0.67, and 0.25 to 0.38 for DMI, ECM yield, EB, BW, and BCS, respectively. Estimated heritability for GFE was 0.32. The genetic correlation of EB at 10d in milk (DIM) with EB at 150 DIM was -0.19, suggesting the genetic regulation of this trait differs by stage of lactation. Positive genetic correlations were found among DMI, ECM yield, and BW averaged over monthly intervals, whereas correlations of these traits with BCS depended upon stage of lactation. Total ECM yield for the lactation was positively correlated with DMI, but a negative genetic correlation between total ECM yield and EB was found. However, the genetic correlation between total ECM yield and EB in the first month of lactation was -0.02, indicating that total production is not genetically correlated with EB during the first month of lactation, when negative EB is most closely associated with diminished fitness. The genetic correlation between GFE and EB ranged from -0.73 to -0.99, indicating that selection for more efficient cows would favor a lower energy status. However, the genetic correlation between EB in the first month of lactation and GFE calculated from 75 to 150 DIM was not significant, indicating that the unfavorable correlation between GFE and EB in early lactation may be minimized with alternative definitions of efficiency. Thus, EB, GFE and related traits will likely respond to genetic selection in Holstein cows. However, the impact of selection for improved feed efficiency on EB must be carefully considered to avoid potential negative consequences of further reductions in EB at the onset of lactation.     

Liability to diseases and their relation to dry matter intake and energy balance in German Holstein and Fleckvieh dairy cows

Dairy cow efficiency is increasingly important for future breeding decisions. The efficiency is determined mostly by dry matter intake (DMI). Reducing DMI seems to increase efficiency if milk yield remains the same, but resulting negative energy balance (EB) may cause health problems, especially in early lactation. Objectives of this study were to examine relationships between DMI and liability to diseases. Therefore, cow effects for DMI and EB were correlated with cow effects for 4 disease categories throughout lactation. Disease categories were mastitis, claw and leg diseases, metabolic diseases, and all diseases. In addition, this study presents relative percentages of diseased cows per days in milk (DIM), repeatability, and cow effect correlations for disease categories across DIM. A total of 1,370 German Holstein (GH) and 287 Fleckvieh (FV) primiparous and multiparous dairy cows from 12 dairy research farms in Germany were observed over a period of 2 yr. Farm staff and veterinarians recorded health data. We modeled health and production data with threshold random regression models and linear random regression models. From DIM 2 to 305 average daily DMI was 22.1 kg/d in GH and 20.2 kg/d in FV. Average weekly EB was 2.8 MJ of NE_L/d in GH and 0.6 MJ of NE_L/d in FV. Most diseases occurred in the first 20 DIM. Multiparous cows were more susceptible to diseases than primiparous cows. Relative percentages of diseased cows were highest for claw and leg diseases, followed by metabolic diseases and mastitis. Repeatability of disease categories and production traits was moderate to high. Cow effect correlations for disease categories were higher for adjacent lactation stages than for more distant lactation stages. Pearson correlation coefficients between cow effects for DMI, as well as EB, and disease categories were estimated from DIM 2 to 305. Almost all correlations were negative in GH, especially in early lactation. In FV, the course of correlations was similar to GH, but correlations were mostly more negative in early lactation. For the first 20 DIM, correlations ranged from ?0.31 to 0.00 in GH and from ?0.42 to ?0.01 in FV. The results illustrate that future breeding for dairy cow efficiency should focus on DMI and EB in early lactation to avoid health problems.     

Genetic correlations between the cumulative pseudo-survival rate,milk yield,and somatic cell score during lactation in Holstein cattle in Japan using a random regression model

Trends in genetic correlations between longevity, milk yield, and somatic cell score (SCS) during lactation in cows are difficult to trace. In this study, changes in the genetic correlations between milk yield, SCS, and cumulative pseudo-survival rate (PSR) during lactation were examined, and the effect of milk yield and SCS information on the reliability of estimated breeding value (EBV) of PSR were determined. Test day milk yield, SCS, and PSR records were obtained for Holstein cows in Japan from 2004 to 2013. A random subset of the data was used for the analysis (825 herds, 205,383 cows). This data set was randomly divided into 5 subsets (162–168 herds, 83,389–95,854 cows), and genetic parameters were estimated in each subset independently. Data were analyzed using multiple-trait random regression animal models including either the residual effect for the whole lactation period (H0), the residual effects for 5 lactation stages (H5), or both of these residual effects (HD). Milk yield heritability increased until 310 to 351 d in milk (DIM) and SCS heritability increased until 330 to 344 DIM. Heritability estimates for PSR increased with DIM from 0.00 to 0.05. The genetic correlation between milk yield and SCS increased negatively to under ?0.60 at 455 DIM. The genetic correlation between milk yield and PSR increased until 342 to 355 DIM (0.53–0.57). The genetic correlation between the SCS and PSR was ?0.82 to ?0.83 at around 180 DIM, and decreased to ?0.65 to ?0.71 at 455 DIM. The reliability of EBV of PSR for sires with 30 or more recorded daughters was 0.17 to 0.45 when the effects of correlated traits were ignored. The maximum reliability of EBV was observed at 257 (H0) or 322 (HD) DIM. When the correlations of PSR with milk yield and SCS were considered, the reliabilities of PSR estimates increased to 0.31–0.76. The genetic parameter estimates of H5 were the same as those for HD. The rank correlation coefficients of the EBV of PSR between H0 and H5 or HD were greater than 0.9. Additionally, the reliabilities of EBV of PSR of H0 were similar to those for H5 and HD. Therefore, the genetic parameter estimates in H0 were not substantially different from those in H5 and HD. When milk yield and SCS, which were genetically correlated with PSR, were used, the reliability of PSR increased. Estimates of the genetic correlations between PSR and milk yield and between PSR and SCS are useful for management and breeding decisions to extend the herd life of cows.     

Genetic heterogeneity of feed intake,energy-corrected milk,and body weight across lactation in primiparous Holstein,Nordic Red,and Jersey cows

In this study, we aimed to estimate and compare the genetic parameters of dry matter intake (DMI), energy-corrected milk (ECM), and body weight (BW) as 3 feed efficiency–related traits across lactation in 3 dairy cattle breeds (Holstein, Nordic Red, and Jersey). The analyses were based on weekly records of DMI, ECM, and BW per cow across lactation for 842 primiparous Holstein cows, 746 primiparous Nordic Red cows, and 378 primiparous Jersey cows. A random regression model was applied to estimate variance components and genetic parameters for DMI, ECM, and BW in each lactation week within each breed. Phenotypic means of DMI, ECM, and BW observations across lactation showed to be in very similar patterns between breeds, whereas breed differences lay in the average level of DMI, ECM, and BW. Generally, for all studied breeds, the heritability for DMI ranged from 0.2 to 0.4 across lactation and was in a range similar to the heritability for ECM. The heritability for BW ranged from 0.4 to 0.6 across lactation, higher than the heritability for DMI or ECM. Among the studied breeds, the heritability estimates for DMI shared a very similar range between breeds, whereas the heritability estimates for ECM tended to be different between breeds. For BW, the heritability estimates also tended to follow a similar range between breeds. Among the studied traits, the genetic variance and heritability for DMI varied across lactation, and the genetic correlations between DMI at different lactation stages were less than unity, indicating a genetic heterogeneity of feed intake across lactation in dairy cattle. In contrast, BW was the most genetically consistent trait across lactation, where BW among all lactation weeks was highly correlated. Genetic correlations between DMI, ECM, and BW changed across lactation, especially in early lactation. Energy-corrected milk had a low genetic correlation with both DMI and BW at the beginning of lactation, whereas ECM was highly correlated with DMI in mid and late lactation. Based on our results, genetic heterogeneity of DMI, ECM, and BW across lactation generally was observed in all studied dairy breeds, especially for DMI, which should be carefully considered for the recording strategy of these traits. The genetic correlations between DMI, ECM, and BW changed across lactation and followed similar patterns between breeds.     

Accuracies of breeding values for dry matter intake using nongenotyped animals and predictor traits in different lactations

Given the interest of including dry matter intake (DMI) in the breeding goal, accurate estimated breeding values (EBV) for DMI are needed, preferably for separate lactations. Due to the limited amount of records available on DMI, 2 main approaches have been suggested to compute those EBV: (1) the inclusion of predictor traits, such as fat- and protein-corrected milk (FPCM) and live weight (LW), and (2) the addition of genomic information of animals using what is called genomic prediction. Recently, several methodologies to estimate EBV utilizing genomic information (EBV) have become available. In this study, a new method known as single-step ridge-regression BLUP (SSRR-BLUP) is suggested. The SSRR-BLUP method does not have an imposed limit on the number of genotyped animals, as the commonly used methods do. The objective of this study was to estimate genetic parameters using a relatively large data set with DMI records, as well as compare the accuracies of the EBV for DMI. These accuracies were obtained using 4 different methods: BLUP (using pedigree for all animals with phenotypes), genomic BLUP (GBLUP; only for genotyped animals), single-step GBLUP (SS-GBLUP), and SSRR-BLUP (for genotyped and nongenotyped animals). Records from different lactations, with or without predictor traits (FPCM and LW), were used in the model. Accuracies of EBV for DMI (defined as the correlation between the EBV and pre-adjusted DMI phenotypes divided by the average accuracy of those phenotypes) ranged between 0.21 and 0.38 across methods and scenarios. Accuracies of EBV for DMI using BLUP were the lowest accuracies obtained across methods. Meanwhile, accuracies of EBV for DMI were similar in SS-GBLUP and SSRR-BLUP, and lower for the GBLUP method. Hence, SSRR-BLUP could be used when the number of genotyped animals is large, avoiding the construction of the inverse genomic relationship matrix. Adding information on DMI from different lactations in the reference population gave higher accuracies in comparison when only lactation 1 was included. Finally, no benefit was obtained by adding information on predictor traits to the reference population when DMI was already included. However, in the absence of DMI records, having records on FPCM and LW from different lactations is a good way to obtain EBV with a relatively good accuracy.     

Comparison of methane production,intensity, and yield throughout lactation in Holstein cows

Genetic selection to reduce methane (CH₄) emissions from dairy cows is an attractive means of reducing the impact of agricultural production on climate change. In this study, we investigated the feasibility of such an approach by characterizing the interactions between CH₄ and several traits of interest in dairy cows. We measured CH₄, dry matter intake (DMI), fat- and protein-corrected milk (FPCM), body weight (BW), and body condition score (BCS) from 107 first- and second-parity Holstein cows from December 2019 to November 2021. Methane emissions were measured using a GreenFeed device and expressed in terms of production (MeP, in g/d), yield (MeY, in g/kg DMI), and intensity (MeI, in g/kg FPCM). Because of the limited number of cows, only animal parameters were estimated. Both MeP and MeI were moderately repeatable (>0.45), whereas MeY presented low repeatability, especially in early lactation. Mid lactation was the most stable and representative period of CH₄ emissions throughout lactation, with animal correlations above 0.9. The average animal correlations of MeP with DMI, FPCM, and BW were 0.62, 0.48, and 0.36, respectively. The MeI was negatively correlated with FCPM (<−0.5) and DMI (>−0.25), and positively correlated with BW and BCS. The MeY presented stable and weakly positive correlations with the 4 other traits throughout lactation, with the exception of slightly negative animal correlations with FPCM and DMI after the 35th week. The MeP, MeI, and MeY were positively correlated at all lactation stages and, assuming animal and genetic correlations do not strongly differ, selection on one trait should lead to improvements in all. Overall, selection for MeI is probably not optimal as its change would result more from CH₄ dilution in increased milk yield than from real decrease in methane emission. Instead, MeY is related to rumen function and is only weakly associated with DMI, FPCM, BW, and BCS; it thus appears to be the most promising CH₄ trait for selection, provided that this would not deteriorate feed efficiency and that a system of large-scale phenotyping is developed. The MeP is easier to measure and thus may represent an acceptable alternative, although care would need to be taken to avoid undesirable changes in FPCM and BW.     

Genetic analysis of milk urea nitrogen and relationships with yield and fertility across lactation

The aim of this project was to investigate the relationship of milk urea nitrogen (MUN) with 3 milk production traits [milk yield (MY), fat yield (FY), protein yield (PY)] and 6 fertility measures (number of inseminations, calving interval, interval from calving to first insemination, interval from calving to last insemination, interval from first to last insemination, and pregnancy at first insemination). Data consisted of 635,289 test-day records of MY, FY, PY, and MUN on 76,959 first-lactation Swedish Holstein cows calving from 2001 to 2003, and corresponding lactation records for the fertility traits. Yields and MUN were analyzed with a random regression model followed by a multi-trait model in which the lactation was broken into 10 monthly periods. Heritability for MUN was stable across lactation (between 0.16 and 0.18), whereas MY, FY, and PY had low heritability at the beginning of lactation, which increased with time and stabilized after 100 d in milk, at 0.47, 0.36, and 0.44, respectively. Fertility traits had low heritabilities (0.02 to 0.05). Phenotypic correlations of MUN and milk production traits were between 0.13 (beginning of lactation) and 0.00 (end of lactation). Genetic correlations of MUN and MY, FY, and PY followed similar trends and were positive (0.22) at the beginning and negative (−0.15) at the end of lactation. Phenotypic correlations of MUN and fertility were close to zero. A surprising result was that genetic correlations of MUN and fertility traits suggest a positive relationship between the 2 traits for most of the lactation, indicating that animals with breeding values for increased MUN also had breeding values for improved fertility. This result was obtained with a random regression model as well as with a multi-trait model. The analyzed group of cows had a moderate level of MUN concentration. In such a population MUN concentration may increase slightly due to selection for improved fertility. Conversely, selection for increased MUN concentration may improve fertility slightly.     

Relationships between milk yield and somatic cell score in Canadian Holsteins from simultaneous and recursive random regression models

Multiple-trait random regression animal models with simultaneous and recursive links between phenotypes for milk yield and somatic cell score (SCS) on the same test day were fitted to Canadian Holstein data. All models included fixed herd test-day effects and fixed regressions within region-age at calving-season of calving classes, and animal additive genetic and permanent environmental regressions with random coefficients. Regressions were Legendre polynomials of order 4 on a scale from 5 to 305 d in milk (DIM). Bayesian methods via Gibbs sampling were used for the estimation of model parameters. Heterogeneity of structural coefficients was modeled across (the first 3 lactations) and within (4 DIM intervals) lactation. Model comparisons in terms of Bayes factors indicated the superiority of simultaneous models over the standard multiple-trait model and recursive parameterizations. A moderate heterogeneous (both across- and within-lactation) negative effect of SCS on milk yield (from −0.36 for 116 to 265 DIM in lactation 1 to −0.81 for 5 to 45 DIM in lactation 3) and a smaller positive reciprocal effect of SCS on milk yield (from 0.007 for 5 to 45 DIM in lactation 2 to 0.023 for 46 to 115 DIM in lactation 3) were estimated in the most plausible specification. No noticeable differences among models were detected for genetic and environmental variances and genetic parameters for the first 2 regression coefficients. The curves of genetic and permanent environmental variances, heritabilities, and genetic and phenotypic correlations between milk yield and SCS on a daily basis were different for different models. Rankings of bulls and cows for 305-d milk yield, average daily SCS, and milk lactation persistency remained the same among models. No apparent benefits are expected from fitting causal phenotypic relationships between milk yield and SCS on the same test day in the random regression test-day model for genetic evaluation purposes.     

Genetics of grass dry matter intake, energy balance, and digestibility in grazing irish dairy cows

The objective of this study was to estimate genetic parameters for grass dry matter intake (DMI), energy balance (EB), and cow internal digestibility (IDG) in grazing Holstein-Friesian dairy cows. Grass DMI was estimated up to 4 times per lactation on 1,588 lactations from 755 cows on 2 research farms in southern Ireland. Simultaneously measured milk production and BW records were used to calculate EB. Cow IDG, measured as the ratio of feed and fecal concentrations of the natural odd carbon-chain n-alkane pentatriacontane, was available on 583 lactations from 238 cows. Random regression and multitrait animal models were used to estimate residual, additive genetic and permanent environmental (co)variances across lactations. Results were similar for both models. Heritability for DMI, EB, and IDG across lactation varied from 0.10 [8 days in milk (DIM)] to 0.30 (169 DIM), from 0.06 (29 DIM) to 0.29 (305 DIM), and from 0.08 (50 DIM) to 0.45 (305 DIM), respectively, when estimated using the random regression model. Genetic correlations within each trait tended to decrease as the interval between periods compared increased for DMI and EB, whereas the correlations with IDG in early lactation were weakest when measured midlactation. The lowest correlation between any 2 periods was 0.10, −0.36, and −0.04 for DMI, EB, and IDG, respectively, suggesting the effect of different genes at different stages of lactations. Eigenvalues and associated eigenfunctions of the additive genetic covariance matrix revealed considerable genetic variation among animals in the shape of the lactation profiles for DMI, EB, and IDG. Genetic parameters presented are the first estimates from dairy cows fed predominantly grazed grass and imply that genetic improvement in DMI, EB, and IDG in Holstein-Friesian cows fed predominantly grazed grass is possible.     

Dynamics of culling for Jersey,Holstein, and Jersey × Holstein crossbred cows in large multibreed dairy herds

The objective of this observational study was to describe and compare the dynamics of reason-specific culling risk for the genetic groups Jerseys (JE), Holsteins (HO), and Jersey × Holstein crossbreds (JH), considering parity, stage of lactation, and milk yield, among other variables, in large multibreed dairy herds in Texas. The secondary objective was to analyze the association between survival and management factors, such as breeding and replacement policies, type of facilities, and use of cooling systems. After edits, available data included 202,384 lactations in 16 herds, ranging from 407 to 8,773 cows calving per year during the study period from 2007 to 2011. The distribution of lactation records by genetic group was 58, 36, and 6% for HO, JE, and JH crosses, respectively. Overall culling rates across breeds were 30.1, 32.1, and 35.0% for JH, JE, and HO, respectively. The dynamics of reason-specific culling were dependent on genetic group, parity, stage of lactation, milk yield, and herd characteristics. Early lactation was a critical period for “died” and “injury-sick” culling. The risk increased with days after calving for “breeding” and, in the case of HO, “low production” culling. Open cows had a 3.5 to 4.6 times greater risk for overall culling compared with pregnant cows. The odds of culling with reason “died” within the first 60 d in milk (DIM) were not significantly associated with genetic group. However, both JE and JH crosses had lower odds of live culling within the first 60 DIM compared with HO cows (OR = 0.72 and 0.82, respectively). Other cow variables significantly associated with the risk of dying within the first 60 DIM were cow relative 305-d mature equivalent (305ME) milk yield, parity, and season of calving. Significant herd-related variables for death included herd size and origin of replacements. In addition to genetic group, the risk of live culling within 60 DIM was associated with cow-relative 305ME milk yield, parity, and season of calving. Significant herd-related variables for live culling included herd-relative 305ME milk yield, herd size, type of facility, origin of replacement, and type of maternity. Overall, reason-specific culling followed similar patterns across DIM in the 3 genetic groups.     

Random regression test-day model for clinical mastitis: Genetic parameters,model comparison,and correlations with indicator traits

The objective was to study genetic (co)variance components for binary clinical mastitis (CM), test-day protein yield, and udder health indicator traits [test-day somatic cell score (SCS) and type traits of the udder composite] in the course of lactation with random regression models (RRM). The study used a data set from selected 15 large-scale contract herds including 26,651 Holstein cows. Test-day production and CM data were recorded from 2007 to 2012 and comprised parities 1 to 3. A longitudinal CM data structure was generated by assigning CM records to adjacent official test dates. Bivariate threshold-linear RRM were applied to estimate genetic (co)variance components between longitudinal binary CM (0 = healthy; 1 = diseased) and longitudinal Gaussian distributed protein yield and SCS test-day data. Heritabilities for liability to CM (heritability ~0.15 from 0 to 305 d after calving) were slightly higher than for SCS for corresponding days in milk (DIM) in the course of lactation. Daily genetic correlations between CM and SCS were moderate to high (genetic correlation ~0.70), but substantially decreased at the very end of lactation. Genetic correlations between CM at different test days were close to 1 for adjacent test days, but were close to zero for test days far apart. Daily genetic correlations between CM and protein yield were low to moderate. For identical DIM (e.g., DIM 20, 160, and 300), genetic correlations were −0.03, 0.11, and 0.18, respectively, and disproved pronounced genetic antagonisms between udder health and productivity. Correlations between estimated breeding values (EBV) for CM from the RRM and official EBV for linear type traits of the udder composite, including EBV from 74 influential sires (sires with >60 daughters), were −0.31 for front teat placement, −0.01 for rear teat placement, −0.31 for fore udder attachment, −0.32 for udder depth, and −0.08 for teat length. Estimated breeding values for CM from the RRM were compared with EBV from a multiple-trait model and with EBV from a repeatability model. For test days covering an identical time span and on a lactation level, correlations between EBV from RRM, multiple-trait model, and repeatability model were close to 1. Most relevant results suggest the routine application of threshold RRM to binary CM to (1) allow selection of genetically superior sires for distinct stages of lactation and (2) achieve higher selection response in CM compared with selection strategies based on indicator type traits or based on the indicator-trait SCS.     

Estimating milk, fat, and protein lactation curves with a test day model.

Test day models were used to estimate lactation curves for milk, fat, protein, fat percentage, and protein percentage and to study the influence of age, season, and herd productivity on Holstein lactation curves. Random effects of lactation within herd and fixed effects of herd test date were absorbed. Fixed effects of cow's age on test day and either DIM (57 divisions) by parity (1, 2, greater than or equal to 3) class or season of calving (winter or summer) by DIM by parity class were estimated. Lactation curves for yield traits derived from DIM solutions were flatter for first versus later lactation, even without addition of age effects. Differences between lactation curves for the two seasons were slight, suggesting that most observed seasonal differences are caused by seasonal productivity accounted for by herd test date effects. At peak, winter calving cows yielded slightly more milk of similar fat percentage but of lower protein percentage than those calving in summer. Data were also partitioned into nine subsets based on rolling herd milk and fat percentage. Lactation curves for yield traits, but not percentage traits, varied with rolling herd milk. Lactation curves for fat yield and percentage varied with rolling herd fat percentage.     

Genetic association of clinical mastitis with test-day somatic cell score and milk yield during first lactation of Finnish Ayrshire cows

In this study the genetic association during lactation of 2 clinical mastitis (CM) traits: CM1 (7 d before to 30 d after calving) and CM2 (31 to 300 d after calving) with test-day somatic cell score (SCS) and milk yield (MY) was assessed using multitrait random regression sire models. The data analyzed were from 27,557 first-lactation Finnish Ayrshire cows. Random regressions on second- and third-order Legendre polynomials were used to model the daily genetic and permanent environmental variances of test-day SCS and MY, respectively, while only the intercept term was fitted for CM. Results showed that genetic correlations between CM and the test-day traits varied during lactation. Genetic correlations between CM1 and CM2 and test-day SCS during lactation varied from 0.41 to 0.77 and from 0.34 to 0.71, respectively. Genetic correlations of test-day MY with CM1 and CM2 ranged from 0.13 to 0.51 and from 0.49 to 0.66, respectively. Correlations between CM1 and SCS were strongest during early lactation, whereas correlations between CM2 and SCS were strongest in late lactation. Genetic correlations lower than unity indicate that CM and SCS measure different aspects of the trait mastitis. Milk yield in early lactation was more strongly correlated with both CM1 and CM2 than milk yield in later lactation. This suggests that selection for higher lactation MY through selection on increased milk yield in early lactation will have a more deleterious effect on genetic resistance to mastitis than selection for higher yield in late lactation. The approach used in this study for the estimation of the genetic associations between test-day and CM traits could be used to combine information from traits with different data structures, such as test-day SCS and CM traits in a multitrait random regression model for the genetic evaluation of udder health.     

Short communication: Genetic relationships between claw disorders,protein yield,and somatic cell score by days in milk

The aim of the present study was to infer daily genetic relationships between the selected claw disorders digital dermatitis, sole ulcer (SU), and interdigital hyperplasia (IH) and protein yield and the udder health indicator somatic cell score (SCS). Data were from 26,651 Holstein cows kept in 15 selected large-scale herds located in the region of Thuringia in the eastern part of Germany. Herds are characterized by organized data recording for novel health traits, and for the present study, claw disorders from the years 2008 to 2012 were used. A longitudinal and binary health data structure was created by assigning claw disorders to adjacent official test days. No entry of a claw disorder within a given interval of approximately 30 d implied a score of 0 (healthy), and otherwise, a score of 1 (diseased). Threshold random regression models (RRM) were applied to binary health data, and linear RRM to Gaussian-distributed protein yield and SCS. Genetic correlations between protein yield and SCS for identical days in milk (DIM) only revealed a tendency for genetic antagonisms between DIM 40 and DIM 180, with a maximal genetic correlation (r_g) of 0.14 at DIM 100. With regard to protein yield and claw disorders, the largest and moderate values of r_g (~0.30), indicating a genetic antagonism between productivity and claw health, were found when correlating protein yield from DIM 300 with SU from DIM 160. Especially for SU and protein yield, time-lagged relationships were more pronounced than genetic relationships from the same test days. Genetic correlations between IH and protein yield were favorable and negative from calving to DIM 300. Generally, on the genetic scale, we found heterogeneous associations between protein yield and claw disorders (i.e., different r_g at identical test days for different claw disorders, and also an alteration of r_g for identical traits at different DIM). The SCS measured at d 20, 160, and 300 was genetically positively correlated with SU over the whole trajectory of 365 d, indicating a common genetic background for claw and udder health. A maximal value of 0.36 was found for the r_g between SCS from d 300 and SU early in lactation. Additionally, a recursive effect was observed (i.e., r_g = 0.26 between SCS from d 20 and SU from d 340). Genetic correlations between SCS and IH, and between SCS and digital dermatitis, were close to zero and partly negative during lactation. Results showed the feasibility of threshold RRM applications to binary claw health data, and a changing genetic background in the course of lactation. From a practical perspective, and with regard to the herds used in this study, continuation of breeding on productivity will have different effects on incidences of different claw disorders, with the highest susceptibility to SU.     

Genetic relationships among body condition score,body weight,milk yield,and fertility in dairy cows

Genetic (co)variances between body condition score (BCS), body weight (BW), milk production, and fertility-related traits were estimated. The data analyzed included 8591 multiparous Holstein-Friesian cows with records for BCS, BW, milk production, and/or fertility from 78 seasonal calving grass-based farms throughout southern Ireland. Of the cows included in the analysis, 4402 had repeated records across the 2 yr of the study. Genetic correlations between level of BCS at different stages of lactation and total lactation milk production were negative (-0.51 to -0.14). Genetic correlations between BW at different stages of lactation and total lactation milk production were all close to zero but became positive (0.01 to 0.39) after adjusting BW for differences in BCS. Body condition score at different stages of lactation correlated favorably with improved fertility; genetic correlations between BCS and pregnant 63 d after the start of breeding season ranged from 0.29 to 0.42. Both BW at different stages of lactation and milk production tended to exhibit negative genetic correlations with pregnant to first service and pregnant 63 d after the start of the breeding season and positive genetic correlations with number of services and the interval from first service to conception. Selection indexes investigated illustrate the possibility of continued selection for increased milk production without any deleterious effects on fertility or average BCS, albeit, genetic merit for milk production would increase at a slower rate.     

The relationship between genetic merit for yield and live weight, condition score, and energy balance of spring calving Holstein Friesian dairy cows on grass based systems of milk production

The objectives of this study were to estimate the effects of genetic merit for milk yield on energy balance, DM intake (DMI), and fertility for cows managed on three different grass-based feeding systems and to estimate possible interactions between genetic merit and feeding system. Individual animal intake estimates were obtained at pasture on 11 occasions across three grazing seasons. The data set contained 96 first lactation, 96 second lactation, and 72 third lactation cows in 1995, 1996, and 1997, respectively. Half of these cows were of high genetic merit, and half were of medium genetic merit for milk solids production. Genetic effects for the traits of interest were estimated as the contrast between the two genetic groups and by the genetic regression of phenotypic performance on the estimated breeding value for fat and protein yield, based on pedigree index. Significant effects of feeding system were observed on yields, DMI, and energy balance, with no effect on live weight, condition score, or reproductive performance. The interaction between genetic merit and feeding system was not significantly different from zero for any of the traits. Yields, grass DMI, and total DMI were all higher for cows of high genetic merit than for those of medium genetic merit and were positively correlated (P < 0.001) with pedigree index. Furthermore, condition score, conception to first and second services, and pregnancy rate were significantly negatively correlated with pedigree index. While at pasture, energy balance was positively (P < 0.01) correlated with pedigree index, although the contrast between high genetic merit and medium genetic merit was not significantly different from zero. This positive energy balance was unexpected and was probably due to the lactation stage that intake was measured. Condition score changes and energy balance measures on a small subgroup of the animals, while indoors offered a diet of silage and concentrates (n = 33), demonstrated that high genetic merit had a more negative energy balance than did medium genetic merit. The results clearly illustrate the production potential of high genetic merit cows on grass-based systems. The reduced reproductive performance questions the suitability of high yield for seasonal calving systems.     

Genetic relationships of clinical mastitis,cystic ovaries,and lameness with milk yield and somatic cell score in first-lactation Canadian Holsteins

The objective of this study was to investigate the genetic relationships of the 3 most frequently reported dairy cattle diseases (clinical mastitis, cystic ovaries, and lameness) with test-day milk yield and somatic cell score (SCS) in first-lactation Canadian Holstein cows using random regression models. Health data recorded by producers were available from the National Dairy Cattle Health System in Canada. Disease traits were defined as binary traits (0 = healthy, 1 = affected) based on whether or not the cow had at least one disease case recorded within 305 d after calving. Mean frequencies of clinical mastitis, cystic ovaries, and lameness were 12.7, 8.2, and 9.1%, respectively. For genetic analyses, a Bayesian approach using Gibbs sampling was applied. Bivariate linear sire random regression model analyses were carried out between each of the 3 disease traits and test-day milk yield or SCS. Random regressions on second-degree Legendre polynomials were used to model the daily sire additive genetic and cow effects on test-day milk yield and SCS, whereas only the intercept term was fitted for disease traits. Estimated heritabilities were 0.03, 0.03, and 0.02 for clinical mastitis, cystic ovaries, and lameness, respectively. Average heritabilities for milk yield were between 0.41 and 0.49. Average heritabilities for SCS ranged from 0.10 to 0.12. The average genetic correlations between daily milk yield and clinical mastitis, cystic ovaries, and lameness were 0.40, 0.26, and 0.23, respectively; however, the last estimate was not statistically different from zero. Cows with a high genetic merit for milk yield during the lactation were more susceptible to clinical mastitis and cystic ovaries. Estimates of genetic correlations between daily milk yield and clinical mastitis were moderate throughout the lactation. The genetic correlations between daily milk yield and cystic ovaries were near zero at the beginning of lactation and were highest at mid and end lactation. The average genetic correlation between daily SCS and clinical mastitis was 0.59 and was consistent throughout the lactation. The average genetic correlation between daily SCS and cystic ovaries was near zero (−0.01), whereas a moderate, but nonsignificant, correlation of 0.27 was observed between SCS and lameness. Unfavorable genetic associations between milk yield and diseases imply that production and health traits should be considered simultaneously in genetic selection.