Effects of breed and region on longevity traits through five years of age in Brown Swiss, Holstein, and Jersey cows in the United States

The objective of this study was to assess breed, and breed × region interactions for several longevity-related traits, measured up to 5 yr of age in Brown Swiss, Holstein, and Jersey cows in 7 regions of the United States. Data were analyzed using logistic, poisson, and linear models, and survival analyses. The traits were stayability (yes/no survived to 5 yr of age), number of completed lactations, days lived, herd-life, and days in milk (DIM) to 5 yr of age. Probable lifetime DIM were also estimated using data from the first 5 yr of age of the cows. Herd-life was defined as the days lived up to 5 yr of age minus the age at first calving. Days in milk consisted of herd-life up to 5 yr of age minus the dry periods. Three data files were analyzed: herds with one breed of cows, herds with Holstein and Brown Swiss, and herds with Holstein and Jersey cows. Breed × region interaction was usually significant, with larger effects for the southern regions. Jerseys obtained largest values for the ratio of DIM to days lived, and for the number of completed lactations to 5 yr of age. Brown Swiss had the largest probabilities of surviving to 5 yr of age (stayabilities) in all regions. For the other traits, the results for Brown Swiss were inconsistent, but usually the cows of this breed had shorter herd-life and DIM to 5 yr of age than Holsteins. Brown Swiss cows were expected to have more total DIM in their lifetime in the Southeast than Holsteins. Survival analysis gave the most readily interpretable information, although the linear, poisson, and logistic analyses answered slightly different questions. Adjustment for herd size did not modify the results.     

The association between milk urea nitrogen and DHI production variables in western commercial dairy herds

A retrospective observational study was conducted using data from Dairy Herd Improvement monthly tests to investigate the association between milk urea nitrogen (MUN) concentration and milk yield, milk protein, milk fat percentage, SCC, and parity for commercial Holstein and Jersey herds in Utah, Idaho, and Montana. Mean MUN for Holstein cows was 15.5 mg/ dl (5.5 mmol/L) MUN and 14.1 mg/dl (5.0 mmol/L) for Jersey cows. Mean MUN, categorized by 30-d increments of days in milk (DIM), paralleled changes in milk values and followed a curvilinear shape. For Holstein cows, concentrations of MUN were different among lactation groups 1, 2, and 3+ for the first 90 DIM for Holsteins. Overall, concentrations of MUN were lower during for the first 30 DIM compared with all other DIM categories for both Holstein and Jersey cows. Multivariate regression models of MUN by milk protein showed that as the milk protein percentage increased, MUN concentration decreased; however, models for Jersey cows showed that MUN did not decrease significantly until above 3.4% milk protein. Milk fat percentage also decreased as MUN increased, but by only 1 mg/dl MUN over the range of 2.2 to 5.8% milk fat. Somatic cell count showed a negative relationship with MUN. Holstein cows with milk protein percentage >3.2% had lower MUN compared with cows having milk protein <3.2% for milk yields from 27.3 to 54.5 kg/d and lower than cows having a milk protein <3.0% for milk yield of 54.5 to 63.6 kg/d. In Jersey cows, MUN concentrations were not different among milk protein percentage categorized by milk yield. This study found that MUN was inversely associated with milk protein percentage and paralleled change in milk yield over time.     

Genomic prediction of lactation curves for milk,fat, protein,and somatic cell score in Holstein cattle

Application of random regression models (RRM) in a 2-step genomic prediction might be a feasible way to select young animals based on the complete pattern of the lactation curve. In this context, the prediction reliability and bias of genomic estimated breeding value (GEBV) for milk, fat, and protein yields and somatic cell score over days in milk (DIM) using a 2-step genomic approach were investigated. In addition, the effect of including cows in the training and validation populations was investigated. Estimated breeding values for each DIM (from 5 to 305 d) from the first 3 lactations of Holstein animals were deregressed and used as pseudophenotypes in the second step. Individual additive genomic random regression coefficients for each trait were predicted using RRM and genomic best linear unbiased prediction and further used to derive GEBV for each DIM. Theoretical reliabilities of GEBV obtained by the RRM were slightly higher than theoretical reliabilities obtained by the accumulated yield up to 305 d (P305). However, validation reliabilities estimated for GEBV using P305 were higher than for GEBV using RRM. For all traits, higher theoretical and validation reliabilities were estimated when incorporating genomic information. Less biased GEBV estimates were found when using RRM compared with P305, and different validation reliability and bias patterns for GEBV over time were observed across traits and lactations. Including cows in the training population increased the theoretical reliabilities and bias of GEBV; nonetheless, the inclusion of cows in the validation population does not seem to affect the regression coefficients and the theoretical reliabilities. In summary, the use of RRM in 2-step genomic prediction produced fairly accurate GEBV over the entire lactation curve for all analyzed traits. Thus, selecting young animals based on the pattern of lactation curves seems to be a feasible alternative in genomic selection of Holstein cattle for milk production traits.     

Use of test-day records beyond three hundred five days for estimation of three hundred five-day breeding values for production traits and somatic cell score of Canadian Holsteins

The Canadian Test-Day Model includes test-day (TD) records from 5 to 305 d in milk (DIM). Because 60% of Canadian Holstein cows have at least one lactation longer than 305 d, a significant number of TD records beyond 305 DIM could be included in the genetic evaluation. The aim of this study was to investigate whether TD records beyond 305 DIM could be useful for estimation of 305-d estimated breeding value (EBV) for milk, fat, and protein yields and somatic cell score. Data were 48,638,184 TD milk, fat, and protein yields and somatic cell scores from the first 3 lactations of 2,826,456 Canadian Holstein cows. All production traits were preadjusted for the effect of pregnancy. Subsets of data were created for variance-component estimation by random sampling of 50 herds. Variance components were estimated using Gibbs sampling. Full data sets were used for estimation of breeding values. Three multiple-trait, multiple-lactation random regression models with TD records up to 305 DIM (M305), 335 DIM (M335), and 365 DIM (M365) were fitted. Two additional models (M305a and M305b) used TD records up to 305 DIM and variance components previously estimated by M335 and M365, respectively. The effects common to all models were fixed effects of herd × test-date and DIM class, fixed regression on DIM nested within region × age × season class, and random regressions for additive genetic and permanent environmental effects. Legendre polynomials of order 6 and 4 were fitted for fixed and random regressions, respectively. Rapid increase of additive genetic and permanent environmental variances at extremes of lactations was observed with all 3 models. The increase of additive genetic and permanent environmental variances was at earlier DIM with M305, resulting in greater variances at 305 DIM with M305 than with M335 and M365. Model M305 had the best ability to predict TD yields from 5 through 305 DIM and less error of prediction of 305-d EBV than M335 and M365. Model M335 had smaller change of 305-d EBV of bulls over the period of 7 yr than did M305 and M365. Model M305a had the least error of prediction and change of 305-d EBV from all models. Therefore, the use of TD records of Holstein cows from 5 through 305 DIM and variance components estimated using records up to 335 DIM is recommended for the Canadian Test-Day Model.     

Use of early lactation days open records for genetic evaluation of cow fertility

National genetic evaluations for cow fertility were introduced by the USDA in February, 2003. These evaluations, reported as daughter pregnancy rate, are based on days open. One requirement of the evaluation system is that lactations be at least 250 d in milk (DIM) to be included for analysis. The objective of this research was to develop a predictor of days open, usable in genetic evaluation, to allow for earlier predicted transmitting abilities (PTA), especially for young bulls. The final prediction equation included an overall intercept, the effects of lactation and calving ease score, the linear and quadratic effects of age at calving, and a regression on days open based on last breeding. Data used for estimation were breeding records from 4 dairy records processing centers for the years 1995 through 1998. Genetic correlations were > or =0.91 by d 130, and phenotypic means of predicted days open were in agreement with means for final days open, indicating that the predictions were phenotypically unbiased. Comparison of mean PTA based on actual and predicted days open indicated no bias in PTA, and correlations between PTA were > or =0.92 by d 130. The earlier use of data increased reliability by about 5% for sires between the ages of 4 and 5 yr. The USDA began using predicted days open for records that are at least 130 DIM in national evaluations starting November, 2003.     

Predicting hyperketonemia by logistic and linear regression using test-day milk and performance variables in early-lactation Holstein and Jersey cows

Although cowside testing strategies for diagnosing hyperketonemia (HYK) are available, many are labor intensive and costly, and some lack sufficient accuracy. Predicting milk ketone bodies by Fourier transform infrared spectrometry during routine milk sampling may offer a more practical monitoring strategy. The objectives of this study were to (1) develop linear and logistic regression models using all available test-day milk and performance variables for predicting HYK and (2) compare prediction methods (Fourier transform infrared milk ketone bodies, linear regression models, and logistic regression models) to determine which is the most predictive of HYK. Given the data available, a secondary objective was to evaluate differences in test-day milk and performance variables (continuous measurements) between Holsteins and Jerseys and between cows with or without HYK within breed. Blood samples were collected on the same day as milk sampling from 658 Holstein and 468 Jersey cows between 5 and 20 d in milk (DIM). Diagnosis of HYK was at a serum β-hydroxybutyrate (BHB) concentration ≥1.2 mmol/L. Concentrations of milk BHB and acetone were predicted by Fourier transform infrared spectrometry (Foss Analytical, Hillerød, Denmark). Thresholds of milk BHB and acetone were tested for diagnostic accuracy, and logistic models were built from continuous variables to predict HYK in primiparous and multiparous cows within breed. Linear models were constructed from continuous variables for primiparous and multiparous cows within breed that were 5 to 11 DIM or 12 to 20 DIM. Milk ketone body thresholds diagnosed HYK with 64.0 to 92.9% accuracy in Holsteins and 59.1 to 86.6% accuracy in Jerseys. Logistic models predicted HYK with 82.6 to 97.3% accuracy. Internally cross-validated multiple linear regression models diagnosed HYK of Holstein cows with 97.8% accuracy for primiparous and 83.3% accuracy for multiparous cows. Accuracy of Jersey models was 81.3% in primiparous and 83.4% in multiparous cows. These results suggest that predicting serum BHB from continuous test-day milk and performance variables could serve as a valuable diagnostic tool for monitoring HYK in Holstein and Jersey herds.     

Manure nutrient excretion by Jersey and Holstein cows

The objective of this study was to evaluate feces, urine, and N excretion by Jersey and Holstein cows. Sixteen multiparous cows (n = 8 per breed) were fed 2 experimental rations at calving in a switchback experimental design. Diets were 50% forage and based on corn meal (control) or whole cottonseed. Half the cows in each breed started on the control diet and half started on the whole cottonseed diet. Cows were switched to the other diet at 60 d in milk and switched back to their original diet at 165 d in milk. Pairs of cows were moved into open-circuit respiration chambers on d 49, 154, and 271 of lactation for 7-d measurement periods. While in the chambers, total collection of feed refusals, milk, recovered hair, feces, and urine was conducted. No effect of the interaction of diet and breed was observed for measures of nutrient digestibility and manure excretion. Total daily manure excretion was lower in Jersey cows than in Holstein cows, with reductions generally proportional to changes in feed intake. Jersey cows consumed 29% less feed and excreted 33% less wet feces and 28% less urine than Holstein cows. Intake, fecal, and urinary N were reduced by 29, 33, and 24%, respectively, in Jersey cows compared with Holstein cows. Equations from American Society of Agricultural and Biological Engineers underpredicted observed values for all manure measures evaluated (urine, manure solids, N, wet manure), and breed bias was observed in equations predicting excretion of urine, N, and wet manure. Although these equations include animal and dietary factors, intercepts of regression of observed values on predicted values differed between Holsteins and Jerseys for those 3 measures. No breed bias was observed in the prediction of manure solids excretion, however, making that equation equally appropriate for Jerseys and Holsteins. The effect of breed on manure and nutrient excretion has significant nutrient management implications.     

Effect of prostaglandin F2α on subclinical endometritis and fertility in dairy cows

The objectives were to determine the effects of PGF_2α treatment on the prevalence of subclinical endometritis (SCE) and fertility of dairy cows. A total of 406 Holstein cows (167 primiparous and 239 multiparous) from 5 herds were used. Uterine lavage for diagnosis of SCE, PGF_2α treatment, evaluation of body condition scores (BCS), and collection of blood samples for estrous cyclicity determination were performed at 21, 35, and 49 d in milk (DIM). Polymorphonuclear cells (PMN) were quantified and thresholds for diagnosing SCE were selected by receiver operating characteristics analysis. Cows classified as having SCE at 35 DIM (≥6.5% PMN) and 49 DIM (≥4.0% PMN) had increased time to pregnancy; however, cows classified as having SCE at 21 DIM (≥8.5% PMN) did not. Median days to pregnancy were delayed by 30 (151 vs. 121 d) and 40 (169 vs. 129) d for cows classified as having SCE at 35 and 49 DIM, respectively. Treatment with PGF_2α did not affect the prevalence of SCE either at 35 (37.9 vs. 38.4%) or at 49 DIM (34.0 vs. 40.4%). Treatment with PGF_2α did not affect time to first insemination (AI; median 76 DIM for cows treated with PGF_2α; 79 DIM for control. Nonetheless, PGF_2α treatment increased pregnancy to first AI in all the cows (35.5 vs. 24.1%) and hazard ratio (HR) of pregnancy in cows with BCS ≤2.5 when all of the cows were evaluated (HR = 1.5; 95% confidence interval; CI = 1.1 to 2.0) and when only cows without SCE were evaluated (HR = 1.8; 95% CI = 1.2 to 2.7). Treatment with PGF_2α did not affect the hazard of pregnancy in cows with SCE at 49 DIM (HR = 0.9; 95% CI = 0.6 to 1.3). In these farms, treatment with PGF_2α did not affect SCE or time to first insemination, but did increase first-service pregnancy per AI and decreased time to pregnancy in cows with low BCS.     

Purulent vaginal discharge in grazing dairy cows: Risk factors,reproductive performance,and prostaglandin F2α treatment

The objectives of this study were to assess the association of a 4-point scale of vaginal discharge score (VDS) with time to pregnancy to define criteria for a practical case of purulent vaginal discharge (PVD) in dairy cows, to test the risk factors for PVD, and, finally, the effect of a dose of PGF_2α on cure and reproductive performance. In experiment 1, grazing Holstein cows (n = 2,414) had their vaginal discharge scored at ～32 d in milk (DIM) on a 4-point scale, the effect of VDS on the hazard of pregnancy by 300 DIM was then assessed to derive a case definition of PVD. Risk factors for PVD and self-cure were also assessed. In experiment 2, grazing Holstein cows (n = 6,326) from 5 herds were checked for PVD at ～30 DIM. Cows with PVD were assigned to receive one dose of 500 μg of PGF_2α analog (Cloprostenol; Ciclase, Syntex SA, Buenos Aires, Argentina) per cow (odd ear tag number) or to remain untreated (even tag number). Cure was declared if cows presented clear normal vaginal discharge (VDS-0) at visit 2 (～62 DIM). Data were analyzed with Cox's regression and mixed logistic models. In experiment 1, cows with VDS ≥1 had lower hazard of pregnancy and longer calving to pregnancy interval than cows with VDS-0. This finding was not affected by the time at which the diagnosis was performed. Therefore, a cow ≥21 DIM and having VDS ≥1 was used to define a case of PVD. The odds of PVD were greater in primiparous cows compared with multiparous, in cows with abnormal calving compared with those with normal calving, and in those losing BCS peripartum. In experiment 2, PGF_2α treatment tended to slightly increase the hazard of pregnancy (adjusted hazard ratio = 1.13). Conversely, PGF_2α had no effect on the odds of cure of PVD [adjusted odds ratio (AOR) = 1.19], pregnancy at first service (AOR = 1.03), or pregnancy by 100 DIM (AOR = 0.89) or 200 DIM (AOR = 1.27). In conclusion, cows with VDS ≥1 can be considered to have PVD because of their lower hazard of pregnancy and longer calving to pregnancy interval (up to 48 d). Important risk factors are parity, calving, and body condition score loss peripartum. Optimal time of diagnosis is ≥28 to 35 DIM because cows experience a high self-cure rate. Self-cure is also affected by parity, prepartum BCS, and VDS. Finally, as treatment with one dose of PGF_2α had a small effect on the hazard of pregnancy and no effect on clinical cure, its therapeutic use in grazing dairy cows with PVD is not recommended.     

Factors Influencing Profitability of Jersey and Holstein Lactations

Dairy Herd Improvement Association production records and farm records on clinical mastitis treatment from 1617 Holstein and Jersey cows calving between July 1980 and July 1983 on one farm in north Florida were used to evaluate the profitability of 2414 lactations. The investigation demonstrates the use of a profit function to determine the relative importance of factors influencing profitability of lactations of different breed, parity and season of calving. It serves to indicate changes that would yield the greatest financial reward by identifying major inherent economic strengths and weaknesses.Breed, parity, and season of calving accounted for less than 10% of the variation in profit per lactation. Holstein lactations were more profitable than Jersey lactations except for parity 4 lactations and for cows calving between January and June when no significant breed differences were apparent. Milk yield was the single most important determinant of profitability in Jersey lactations, regardless of parity or season of calving. Days of mastitis treatment explained more of the variation in profit per lactation than did milk yield for Holstein lactations calving between January and March and for Holstein lactations of parity greater than 4. Mastitis costs were a maximum of 8.2% of total costs for Holstein lactations of parity 6.     

Dynamics of automatically generated body condition scores during early lactation and pregnancy at first artificial insemination of Holstein cows

The objective of this study was to characterize the association between body condition score (BCS) and BCS change (ΔBCS), determined by an automated camera system at multiple time points, and the subsequent pregnancy per first artificial insemination (P/AI1) of Holstein cows. A retrospective observational study was completed using data collected from 11,393 lactations in 7,928 Holstein cows calving between April 2019 and March 2021 in a commercial dairy operation located in Colorado. Cows were classified as primiparous or multiparous. Scores generated by BCS cameras at dry-off, calving, 21 days in milk (DIM), 56 DIM, and first artificial insemination were selected for the analyses and subsequently categorized as low (≤lower quartile), moderate (interquartile range), and high (≥upper quartile). Changes in BCS were calculated by periods of interest as change from dry-off to calving (multiparous cows); change from calving to 21 DIM; change from calving to 56 DIM; and change from calving to first artificial insemination and assigned into categories as large loss of BCS (top 25% of cows losing BCS); moderate loss (bottom 75% of cows losing BCS); no change (ΔBCS = 0); or gain of BCS (ΔBCS > 0). Data were examined in primiparous and multiparous cows separately using logistic regression and time-to-event analyses. Initial univariable models were followed by multivariable models that considered calving season, occurrence of disease, and milk yield up to 60 DIM as covariables. The logistic regression analyses indicated that in both parity groups the associations between BCS category and P/AI1 were more evident at 21 DIM, 56 DIM, and first artificial insemination, with lower odds of P/AI1 in cows in the low BCS category. Likewise, cows with large loss in BCS between calving and 21 DIM, calving and 56 DIM, and calving and first artificial insemination had lower odds of P/AI1 compared with other categories of ΔBCS within the same period of interest. Similarly, survival analyses evidenced that cows in the low BCS category required more time to get pregnant. In agreement, differences in the dynamics of the average daily BCS during the first 90 DIM were evident when cows were grouped by first AI outcome (pregnant vs. nonpregnant) and by their time to pregnancy category (<90 DIM; 91–150 DIM; or >150 DIM), with cows with reduced fertility showing lower BCS up to 90 DIM. Overall, low BCS and more pronounced reductions in BCS occurring closer to first artificial insemination resulted in lower odds of pregnancy per artificial insemination.     

Effect of early or late resynchronization based on different methods of pregnancy diagnosis on reproductive performance of dairy cows

The aim of this study was to compare the reproductive performance of dairy cows subjected to early (ER) or late (LR) resynchronization programs after nonpregnancy diagnoses based on either pregnancy-associated glycoproteins (PAG) ELISA or transrectal palpation, respectively. In addition, the accuracy of the PAG ELISA for early pregnancy diagnosis was assessed. Lactating Holstein cows were subjected to a Presynch-Ovsynch protocol with timed artificial insemination (AI) performed between 61 and 74 DIM. On the day of the first postpartum AI, 1,093 cows were blocked by parity and assigned randomly to treatments; however, because of attrition, 452 ER and 520 LR cows were considered for the statistical analyses. After the first postpartum AI, cows were observed daily for signs of estrus and inseminated on the same day of detected estrus. Cows from ER that were not reinseminated in estrus received the first GnRH injection of the Ovsynch protocol for resynchronization 2 d before pregnancy diagnosis. On d 28 after the previous AI (d 27 to 34), pregnancy status was determined by PAG ELISA, and nonpregnant cows continued on the Ovsynch protocol for reinsemination. Pregnant cows had pregnancy status reconfirmed on d 46 after AI (d 35 to 52) by transrectal palpation, and those that lost the pregnancies were resynchronized. Cows assigned to LR had pregnancy diagnosed by transrectal palpation on d 46 after AI (d 35 to 52) and nonpregnant cows were resynchronized with the Ovsynch protocol. Blood was sampled on d 28 after AI (d 27 to 34) from cows in both treatments that had not been reinseminated on estrus and again on d 46 after AI (d 35 to 52) for assessment of PAG ELISA to determine the accuracy of the test. Cows were subjected to treatments for 72 d after the first insemination. Pregnancy per AI (P/AI) at first postpartum timed AI did not differ between treatments and averaged 28.9%. The proportion of nonpregnant cows that were resynchronized and received timed AI was greater for ER than for LR (30.0 vs. 7.6%). Cows in ER had a shorter interval between inseminations when inseminated following spontaneous estrus (21.7 ± 1.1 vs. 27.8 ± 0.8 d) or after timed AI (35.3 ± 1.2 vs. 55.2 ± 1.4 d). Nevertheless, the ER did not affect the rate of pregnancy (adjusted hazard ratio = 1.23; 95% confidence interval = 0.94 to 1.61) or the median days postpartum to pregnancy (ER = 132 vs. LR = 140). A total of 2,129 PAG ELISA were evaluated. Overall, sensitivity, specificity, and positive and negative predictive values averaged 95.1, 89.0, 90.1, and 94.5%, respectively, and the accuracy was 92.1%. In conclusion, PAG ELISA for early diagnosis of pregnancy had acceptable accuracy, but early resynchronization after nonpregnancy diagnosis with PAG ELISA did not improve the rate of pregnancy or reduce days open in dairy cows continuously observed for estrus.     

Short communication: Jersey × Holstein crossbreds compared with pure Holsteins for body weight, body condition score, fertility, and survival during the first three lactations

Crossbred cows (n=80) resulting from the use of Jersey (JE) semen on their pure Holstein (HO) dams were compared with pure HO cows (n=77) for body weight, body condition score, fertility, and survival during their first 3 lactations. Cows were in 2 research herds of the University of Minnesota and calved from September 2003 to June 2008. The JE × HO crossbred cows had significantly less body weight during the first (-56 kg), second (-67 kg), and third (-82 kg) lactations than pure HO cows. However, JE × HO cows had significantly greater body condition score during the first (2.94 vs. 2.84), second (2.97 vs. 2.84), and third (2.99 vs. 2.87) lactations than pure HO cows. For fertility, JE × HO cows had fewer days to first breeding during the first (-10.6d), second (-8.4d), and third (-12.3d) lactations than pure HO cows. Crossbred cows were not significantly different from pure HO cows for number of services during first lactation; however, JE × HO cows had significantly fewer services (2.2) than pure HO cows (2.7) during the second lactation. Also, JE × HO cows had significantly fewer days open than pure HO cows in the first (-24 d), second (-42 d), and third (-42 d) lactations. For survival, JE × HO cows were not significantly different from pure HO cows for percentage of cows calving a second time; however, a tendency existed for a higher percentage of JE × HO cows (63.8%) than pure HO cows (49.4%) to calve a third time, and a higher percentage of JE × HO cows calved a third time within 28, 34, and 40 mo of first calving than pure HO cows.     

Phenotypic relationships of common health disorders in dairy cows to lactation persistency estimated from daily milk weights

The objective of this study was to investigate the phenotypic relationship between common health disorders in dairy cows and lactation persistency, uncorrelated with 305-d yield. The relationships with peak yield and days in milk (DIM) at peak were also studied. Daily milk weights and treatment incidence records of 991 Holstein lactations from experimental dairy herds at Virginia Tech and Pennsylvania State University were used. Persistency was calculated as a function of daily yield deviations from standard lactation curves, developed separately for first (FL) and later lactations (LL), and deviations of DIM around reference dates: 128 for FL and 125 for LL. Days in milk at peak and peak yield were computed for each lactation by using Wood's function. The disease traits studied were mastitis (MAST) only during the first 100 d (MAST1), only after 100 DIM (MAST2), both before and after 100 DIM (MAST12), and at any stage of lactation (MAST1/2), as well as metritis, displaced abomasums, lameness, and metabolic diseases. Each disease was defined as a binary trait, distinguishing between lactations with at least one incidence (1) and lactations with no incidences (0). The relationships of diseases to persistency, DIM at peak, and peak yield were investigated separately for FL and LL for all disease traits except MAST12, which was investigated across parities. The relationships of persistency to probability of the diseases in the same lactation and in the next lactation were examined using odds ratios from a logistic regression model. Metritis and displaced abomasums in FL and LL were associated with significantly higher persistencies. Metabolic diseases and MAST1 in LL were significantly related to higher persistencies. The relationships of MAST2 in both FL and LL, and MAST12 across parities were significant but negative. Cows affected by MAST tended to have less persistent lactations. Most of the diseases had a significant impact on DIM at peak in LL. In LL, metritis, metabolic diseases, and displaced abomasums tended to significantly delay DIM at peak. Mastitis only after 100 DIM was associated with significantly earlier DIM at peak in LL. Increasing persistency was associated with low MAST2 and MAST1/2 in primiparous cows. None of the diseases studied was significantly related to persistency of the previous lactation.     

Genetic parameters for female fertility traits and a fertility index in Spanish dairy cattle

Genetic correlations among female fertility traits (linear and binary) were estimated using 225,085 artificial insemination records from 120,713 lactations on 63,160 Holstein cows. Fertility traits were: calving interval, days open, a linear transformation of days open, days to first insemination, interval between first and last insemination, number of inseminations per service period, pregnancy within 56 and 90 d after first insemination, and success in first insemination. A bivariate animal model was implemented using Bayesian methods in the case of binary traits. Low heritabilities (0.02 to 0.06) were estimated for these fertility traits. Strong genetic correlations (0.89 to 0.99) were found among traits, except for days to first service, where the genetic correlation with other fertility traits ranged from −0.52 to −0.18 for binary traits, and from 0.50 to 0.82 for days to first service, calving interval, and days open. Four fertility indices were proposed utilizing information from insemination records; these indices combined one indicator of the beginning of the service period and one indicator of conception rate. Two additional indices used information from the milk-recording scheme, including calving interval and a linear transformation of days open. The fertility index composed of days to first service and pregnancy within 56 d achieved the highest genetic gain for reducing fertility cost, reducing days to first service, and reducing the number of inseminations per lactation ($8.60, −1.31 d, and −0.03 AI, respectively). This index achieved at least 15% higher genetic gain than obtained from indices with information from the milk recording scheme only (calving interval and days open).     

Effect of endocervical inflammation on days to conception in dairy cows

In contrast to endometritis, now diagnosed by cytological examination, the effect of endocervical inflammation on reproductive performance has been inadequately investigated. In this study, endocervical and endometrial cytological specimens were collected from 168 Holstein cows between 21 and 60 days in milk (DIM) to investigate the prevalence of endocervical inflammation and effect on days to conception. Statistical analyses were stratified based on DIM at examination (<35 vs. ≥35 DIM). Endocervical inflammation with ≥5% neutrophils before 35 DIM (disregarding the level of endometrial inflammation) was associated with decreased hazard of pregnancy within 300 DIM (adjusted hazard ratio, 0.4; 95% confidence interval: 0.3-0.8). A decrease in hazard of pregnancy was observed when >6% neutrophils were counted in endometrial smears (adjusted hazard ratio, 0.4; 95% confidence interval: 0.2-0.7). The study suggested an additive effect of combined endocervical and endometrial inflammation on the hazard of pregnancy within 300 DIM. Using the thresholds of 5% neutrophils for the cervix and 6% neutrophils for the uterus, 11% of the cows examined before 35 DIM presented cervicitis only, 13% were affected by endometritis only, and 32% suffered from both endometrial and endocervical inflammation. The presence (absence) of cervicitis was not indicative of the presence (absence) of endometritis. This study showed that in addition to uterine inflammation, endocervical inflammation in early lactation affects conception. Thus, the global evaluation of genital tract health may be more beneficial for reproductive performance than that of endometrial inflammation.     

Short communication: Multivariate outlier detection for routine Nordic dairy cattle genetic evaluation in the Nordic Holstein and Red population

It is of practical importance to ensure the data quality from a milk-recording system before use for genetic evaluation. A procedure was developed for detection of multivariate outliers based on an approximation for Mahalanobis distance and was implemented in the Nordic Holstein and Red population. The general target of this procedure is based on the Nordic Cattle Genetic Evaluation yield model, which is a 9-trait model for milk, protein, and fat in the first 3 lactations. The procedure is based on the phenotypic correlation structure as a function of days in milk (DIM) and on computation of trait means and standard deviations within a production year, lactation, and DIM. For each record in the data, a Mahalanobis distance value was computed based on the trait mean and the covariance matrix for the actual production year, lactation, and DIM. A set of cutoff values, ranging from 10 to 100 with steps of 10, for discarding multivariate outliers was investigated. Prediction accuracy was calculated as the Pearson correlations between estimated breeding values predicted by full data set and estimated breeding values predicted by reduced data set for cows without records in the reduced data set and with 1 or more records deleted due to the editing rules on Mahalanobis distance. The results showed that, averaged over all scenarios, gains of 0.005 to 0.048 on prediction accuracy have been obtained by deleting the multivariate outliers. The improvements were more profound for progeny of young bulls compared with progeny of proven bulls. It is easy to implement this multivariate outlier-detection procedure in the routine genetic evaluation for different dairy cattle breeds; however, an optimal cutoff value for Mahalanobis distance needs to be defined to achieve an acceptable compromise between genetic evaluation accuracy and data deletion.     

Association of blood calcium concentration in the first 3 days after parturition and energy balance metabolites at day 3 in milk with disease and production outcomes in multiparous Jersey cows

Research exploring specific associations of markers of negative energy balance and Ca in postpartum Jersey cows with lactation performance is lacking. Our objectives were to evaluate the associations of total Ca concentration (tCa) measured at 1 through 3 d in milk (DIM) and free fatty acids (FFA), β-hydroxybutyrate (BHB), and glucose measured at 3 DIM with (1) the risk of multiparous Jersey cows being diagnosed with early-lactation diseases and culling, (2) milk production in the first 9 wk of lactation, and (3) the risk of pregnancy in the first 150 DIM. A cohort study was performed in 1 dairy herd in Texas. Multivariable Poisson regression models were built to evaluate the association of the analytes of interest with the risks of early-lactation diseases and culling in the first 60 DIM (i.e., binary outcomes). Linear mixed models were used to evaluate the association of the analytes of interest with milk production within the first 9 wk of lactation, and a Cox proportional hazard model was built to assess the risk of pregnancy within 150 DIM. A total of 380 cows were used in the final analyses. Total Ca measured at 1 through 3 DIM was not associated with the risk of metritis. Cows with increased FFA and BHB had an increased risk of being diagnosed with metritis and clinical mastitis, respectively. Increased concentrations of glucose and FFA and decreased tCa at 3 DIM were associated with an increased risk of culling. Reduced tCa concentrations at 1 DIM (≤1.84 mmol/L) and 2 DIM (≤2.04 mmol/L) were associated with increased milk production across the first 9 wk of lactation compared with tCa concentrations above those thresholds. Total Ca was not associated with milk production when assessed at 3 DIM, whereas increased FFA (≥0.37 mmol/L) and decreased glucose (≤2.96 mmol/L) at 3 DIM were associated with increased milk production. None of the metabolites measured were associated with the risk of pregnancy in the first 150 DIM. Our results demonstrate that tCa concentration assessed in the first 3 DIM show temporary associations with milk production and culling in multiparous Jersey cows. Although increased concentration of FFA assessed at 3 DIM was associated with greater milk yield, it was a detrimental factor for the risk of metritis. This study attempted to better elucidate the relationship of tCa, FFA, BHB, and glucose assessed in early postpartum with health and performance of Jersey cows. Based on this study, assessments performed at 3 DIM using tCa concentration ≤1.99 mmol/L for increased risk of early-lactation culling and FFA ≥0.43 mmol/L for increased risk of metritis could be used as starting points. More studies evaluating the dynamics of energy balance markers and tCa in postpartum Jersey cows using a greater number of herds are needed to better inform dairy consultants on critical levels for exacerbated postpartum negative energy balance and subclinical hypocalcemia for the Jersey breed.     

Milk production change following clinical mastitis and reproductive performance compared among J5 vaccinated and control dairy cattle

Naturally occurring cases of bovine clinical mastitis (CM) were studied among J5 vaccinates and controls on 3 commercial dairy farms. Milk production change and reproductive performance following CM were compared between the 2 groups. Among 306 controls and 251 vaccinates, there were 221 new cases of CM affecting 120 cows; 437 lactations never had a case of CM. Environmental pathogens made up 90% (159/176) of etiologic agents isolated. Change in daily milk production following CM was associated with J5 vaccination, days in milk (DIM) at onset of CM, and herd effect as well as each 2-way interaction between the 3 factors. The adjusted daily milk for 21 d following CM was 7.6 kg greater among J5 vaccinates than controls; however, this protective effect of vaccination waned with increasing DIM at onset of CM. A mixed linear model with autoregressive order 1 [AR(1)] correlation structure estimated the daily milk production of any cow (whether or not she had CM) on a given DIM. Cows with CM caused by nonagalactiae streptococci, Staphylococcus aureus, Escherichia coli, or Klebsiella lost significant daily milk production for the entire lactation relative to nonmastitic cows. Another mixed linear model for only coliform CM cases (E. coli, Klebsiella, and Enterobacter) within the first 50 DIM showed milk loss for 21 d following coliform CM to be significantly less for J5 vaccinates than for controls, by 6 to 15 kg per day. Cows were significantly less likely to become pregnant if they had CM caused by E. coli (42% pregnant) or Streptococcus spp. (38% pregnant), whereas 78% (342/437) of cows with no mastitis conceived. Days open (number of days from calving until pregnancy) averaged 131 d for cows with no CM and 162 d for cows that had at least one case of CM. Days until conception, days until last breeding, days open, times bred, and percentage of cows pregnant by 200 DIM were not changed with J5 vaccination. Nonetheless, an important benefit of the use of J5 bacterin appears to be reduction of the loss of daily milk production following CM, whether all cases or only those caused by coliform bacteria were considered.     

Comparison of random regression models with Legendre polynomials and linear splines for production traits and somatic cell score of Canadian Holstein cows

A random regression model with both random and fixed regressions fitted by Legendre polynomials of order 4 was compared with 3 alternative models fitting linear splines with 4, 5, or 6 knots. The effects common for all models were a herd-test-date effect, fixed regressions on days in milk (DIM) nested within region-age-season of calving class, and random regressions for additive genetic and permanent environmental effects. Data were test-day milk, fat and protein yields, and SCS recorded from 5 to 365 DIM during the first 3 lactations of Canadian Holstein cows. A random sample of 50 herds consisting of 96,756 test-day records was generated to estimate variance components within a Bayesian framework via Gibbs sampling. Two sets of genetic evaluations were subsequently carried out to investigate performance of the 4 models. Models were compared by graphical inspection of variance functions, goodness of fit, error of prediction of breeding values, and stability of estimated breeding values. Models with splines gave lower estimates of variances at extremes of lactations than the model with Legendre polynomials. Differences among models in goodness of fit measured by percentages of squared bias, correlations between predicted and observed records, and residual variances were small. The deviance information criterion favored the spline model with 6 knots. Smaller error of prediction and higher stability of estimated breeding values were achieved by using spline models with 5 and 6 knots compared with the model with Legendre polynomials. In general, the spline model with 6 knots had the best overall performance based upon the considered model comparison criteria.