Genetic parameters for novel mastitis traits defined by combining test-day somatic cell score and differential somatic cell count in the first lactation of Japanese Holsteins.

In this study, we aimed to improve current udder health genetic evaluations by addressing the limitations of monthly sampled somatic cell score (SCS) for distinguishing cows with robust innate immunity from those susceptible to chronic infections. The objectives were to (1) establish novel somatic cell traits by integrating SCS and the differential somatic cell count (DSCC), which represents the combined proportion of polymorphonuclear leukocytes and lymphocytes in somatic cells and (2) estimate genetic parameters for the new traits, including their daily heritability and genetic correlations with milk production traits and SCS, using a random regression test-day model (RRTDM). We derived 3 traits, termed ML_SCS_DSCC, SCS_4_DSCC_65_binary, and ML_SCS_DSCC_binary, by using milk loss (ML) estimates at corresponding SCS and DSCC levels, thresholds established in previous studies, and a threshold established from milk loss estimates, respectively. Data consisted of test-day records collected during January 2021 through March 2022 from 265 herds in Hokkaido, Japan. From these records, we extracted records between 7 to 305 d in milk (DIM) in the first lactation to fit the RRTDM. The model included the random effect of herd-test-day, the fixed effect of year-month, fixed lactation curves nested with calving age groups, and random regressions with Legendre polynomials of order 3 for additive genetic and permanent environmental effects. The analysis was performed using Gibbs sampling with Gibbsf90+ software. The averages (ranges) of the daily heritability estimates over lactation were 0.086 (0.075-0.095) for SCS, 0.104 (0.073-0.127) for ML_SCS_DSCC, 0.137 (0.014-0.297) for SCS_4_DSCC_65_binary, and 0.138 (0.115-0.185) for ML_SCS_DSCC_binary; the heritability curve for SCS_4_DSCC_65_binary was erratic. Genetic correlations within the trait decreased as the DIM interval widened, especially for those integrating DSCC, indicating that these traits should be analyzed using RRTDM rather than repeatability models. The averages (ranges) of genetic correlations with milk yield over lactation were 0.01 (-0.22 to 0.28) for SCS, -0.05 (-0.40 to 0.13) for ML_SCS_DSCC, -0.08 (-0.17 to 0.09) for SCS_4_DSCC_65_binary, and -0.08 (-0.22 to 0.27) for ML_SCS_DSCC_binary. Compared with SCS, the newly defined traits exhibited slightly stronger negative genetic correlations with milk yield. Especially in late lactation stages, the genetic correlation between ML_SCS_DSCC and milk yield was significantly below zero, with a posterior median of -0.40. Furthermore, the new traits showed positive correlations with SCS, having estimates varying from 0.68 to 0.85 for ML_SCS_DSCC, 0.14 to 0.47 for SCS_4_DSCC_65_binary, and 0.61 to 0.66 for ML_SCS_DSCC_binary, depending on DIM. Considering that ML_SCS_DSCC and ML_SCS_DSCC_binary have relatively high heritability (compared with SCS) and favorable genetic correlations with milk production traits and SCS, their incorporation into breeding programs appears promising. Nevertheless, their genetic relationships with (sub)clinical mastitis require further investigation.

Image analysis traits of multiple muscles and intermuscular/subcutaneous fat influence on Japanese Black beef carcass price and genetic parameters.

Objective: The purposes of this study were to investigate the relationship between carcass unit price per 1 kg (UP) and multiple muscles and intermuscular fat (IF)/subcutaneous fat (SF) of beef carcasses using image analysis of cross-section images for Wagyu beef cattle in Japan, and to estimate their genetic parameters. Methods: The carcasses used in this study were 1,807 Japanese Black (Wagyu) cattle (1,216 steers and 591 heifers). An analysis of variance (ANOVA) was conducted with UP as the dependent variable and market date, age in months, sex, and image analysis traits (IAT) as fixed effects, and standard partial regression coefficients were calculated for each IAT on UP. Also, the heritability of each IAT that affected UP and genetic correlation among IAT vs. carcass grading traits were estimated. Results: Not only IAT related to carcass grading traits, M. trapezius dorsi, M. latissimus dorsi, and IF traits were significant differences in UP (p < 0.05). The heritability of IAT associated with UP was estimated at 0.38~0.85. The genetic correlations between the area and thickness of M. trapezius dorsi and M. latissimus dorsi vs. rib eye area (REA) were estimated to be moderately positive (0.53~0.66), while the genetic correlations between the IF area percentage vs. carcass weight (CWT), REA and yield score (YS) were estimated to be negative (-0.40, -0.56 and -0.34). Conclusion: UP was influenced by various traits, including M. trapezius dorsi, M. latissimus dorsi, and IF traits, as well as image analysis associated with carcass grading traits. Since these IAT associated with UP had hereditary and desirable genetic correlations with carcass grading traits, these traits were also important for genetic improvement.

Optimal value for the exponential term of Wilmink's function according to current Holstein lactation curves in Japan.

We compared values of Wilmink's exponential term to describe the lactation curves of Holstein cows in Japan. Data were a total of 100,971,798 test-day records from the first through fifth parities during 1991 through 2018. The lactation curve model used fourth-order Legendre polynomials and Wilmink's exponential term. In total, 810 analyses were executed to compare six values (-0.02, -0.03, …, -0.07) for the exponential term to select the one that yielded the smallest root mean square error. For all parities, daily milk yield and lactation persistency increased consistently and peak lactation days occurred later from year to year. For the years evaluated, the optimal exponential term was -0.05 for first and second parities, -0.04 for third parity, and -0.03 for fourth and fifth parities. The change in the exponential parameter with increasing year was related to delays in peak lactation.

Longitudinal genetic analysis of semen production traits in Holstein bulls according to a random regression animal models on age.

Here we used random regression animal models (RRAMs) to investigate genetic change over age in the semen volume (VOL) and sperm concentration (CON) of Holstein bulls. We used 35,294 collection records from 1284 Holstein bulls and their 4166 pedigree records. The models included year and month of collection, collection place, collection method, and number of collections attempted for each day and month of age (second-order regressions) as fixed effects; technician as a random effect; and additive genetic and permanent environment as random regressions (first-order regressions). We examined two RRAMs with homogeneous and heterogeneous residual variances (RRAM1 and RRAM2, respectively). By using RRAM1, heritability for VOL and CON increased from 0.08 to 0.61 and 0.18 to 0.57, respectively, between 10 and 126 months of age. By using RRAM2, heritability for VOL increased from 0.11 to 0.28 between 10 and 24 months of age for young bulls and increased from 0.08 to 0.48 between 25 and 126 months of age for mature bulls; heritability for CON ranged from 0.18 to 0.19 for young bulls and increased from 0.10 to 0.48 for mature bulls. Posterior genetic correlations between young ages and older ages were strongly positive for VOLs but weak for CONs.

Factors contributing to the swimmer puppy syndrome found in Labrador retrievers.

BACKGROUND: Swimmer puppy syndrome is a disease found in neonatal puppies mainly characterized by the inability to stand, but its direct cause is unknown. Since swimmer puppies were observed infrequently but continuously among the Labrador retriever colony at the Hokkaido Guide Dogs for the Blind Association in Japan, based on their birth record and pedigree, factors related to the onset of swimmer puppy syndrome in Labrador retrievers were examined. RESULTS: The total number of offspring over seven years was 436, of which 16 were swimmer puppies. Most of the affected puppies except one recovered steadily. As for the swimmer puppies, the litter size was significantly lower, and the body weights on the 10th and 28th day after delivery were significantly higher than the non-symptomatic puppies. These results suggested that the onset may be related to weight gain in the neonatal stages due to a small litter size. According to the genetic analysis, 26 ancestors common to the affected individuals were confirmed, but the causative individual could not be identified with the inbreeding coefficient. The heritability of the swimmer-puppy onset trait was 0.80, and the heritability for the the 10th-day body-weight trait was equally high at 0.78, both of which strongly suggest genetic involvement. CONCLUSIONS: In this study, the onset of swimmer puppy syndrome in the Labrador retrievers was associated with litter size and early weight gain, and result of study suggests that genetic influence might be involved.

Suitable indicator of heat stress for genetic evaluation of heat tolerance in Holstein cows in Japan.

Only a few, principal, weather stations in Japanese prefectures have the daily humidity records required to calculate the temperature-humidity index (THI) as a dairy cow heat-stress indicator. We compared three heat-stress indices: (1) THI calculated from daily average temperature and daily relative humidity at a principal weather station (PTHI); (2) daily average temperature at each herd's closest local weather station (TEMP); and (3) THI calculated from daily average temperature at each herd's closest local weather station and daily relative humidity at the principal weather station (HTHI). We used daily records from 532 provincial weather stations and test-day records of milk production from Days 6 to 305 post-first-calving in Holsteins to compare the indices as indicators of heat-stress effects on milk yield and somatic cell score (SCS). The models used the BLUPF90 package to analyze the effects of herd-year, calving age, days in milk, and PTHI, TEMP, or HTHI. We estimated each model's mean square error (MSE) and compared suitabilities among indices for each trait. TEMP heat-stress thresholds were ~18°C (milk yield) and 15-20°C (SCS). The MSE of the HTHI model was the smallest, but no significant differences were found among the indices for milk yield.

Genetic evaluation of heat tolerance in Holstein cows in Japan.

We used test-day records and daily records from provincial weather stations in Japan to evaluate heat tolerance (HT) in Holstein cows according to a random regression test-day model. Data were a total of 1,641,952 test-day records for heritability estimates and 17,245,694 test-day records for genetic evaluation of HT by using milk yield and somatic cell score (SCS) in Holstein cows that had calved for the first time in 2000 through 2015. Temperature-humidity index (THI) values were estimated by using average daily temperature and average daily relative humidity records from 60 provincial Japanese weather stations. The model contained herd-test-day, with lactation curves on days in milk within month-age group as a fixed effect. General additive genetic effect and HT of additive genetic effect were included as random effects. The threshold value of THI was set to 60. For milk yield, estimated mean heritabilities were lower during heat stress (THI = 78; 0.20 and 0.28) than when below the heat stress threshold (THI ≤ 60; 0.26 and 0.31). For SCS, heritability estimates (range 0.08-0.10) were similar under all heat stress conditions. Genetic trends of HT indicated that EBVs of HT are changing in an undesirable direction.

Length of lags in responses of milk yield and somatic cell score on test day to heat stress in Holsteins.

We used daily records from provincial Japanese weather stations and monthly test-day records of milk production to investigate the length of the lags in the responses of cows' milk yield and somatic cell score (SCS) to heat stress (HS). We also investigated the HS thresholds in milk yield and SCS. Data were a total of 17,245,709 test-day records for milk and SCS in Holstein cows that had calved for the first time between 2000 and 2015, along with weather records from 60 weather stations. Temperature-humidity index (THI) values were estimated by using average daily temperature and average daily relative humidity. Adjusted THI values were calculated by using temperature, relative humidity, wind speed, and solar radiation. The model contained herd, calving year, month of test day, age group, days in milk, and THI as a fixed effect. THIs for each day from 14 days before the test day until the test day were used to represent the HS effects. The HS occurring 3 days, and between 8 and 10 days, before the test day had the greatest effect on the milk yield and SCS, respectively. The threshold THI values for the HS effect were about 60-65 for both traits.

Development of genetic evaluation for milk production traits of Holsteins in Japan.

The procedure used for the genetic evaluation of dairy cattle in Japan has developed from a lactation sire-MGS model to a multiple-lactation random regression test-day animal model. Genetic evaluation of Holstein bulls in Japan began in 1989 with the use of field-style progeny testing; dairy herd improvement program data from all over Japan were used, along with a sire and maternal grandsire model. In 1993, an animal model was introduced to estimate breeding values for yield and type traits. A random regression test-day model was first applied in 2010. In the business of breeding dairy cattle, it is very important to users that estimated breeding values are reliable and stable among subsequent routine evaluations. With experience in the genetic evaluation of dairy cattle in Japan, Japanese researchers have found ways to improve the stability of estimated breeding values. These modifications involve changes in data editing, development of evaluation models, changes to the structures of unknown-parent groups, awareness of the problems of predicting lactation yield from partial test-day records, and adjustment for heterogeneity within herd variances. Here, I introduce developments in, and our experiences with, the genetic evaluation of yield traits of Holstein cattle in Japan.

Prediction of random-regression coefficient for daily milk yield after 305 days in milk by using the regression-coefficient estimates from the first 305 days.

OBJECTIVE: Because lactation periods in dairy cows lengthen with increasing total milk production, it is important to predict individual productivities after 305 days in milk (DIM) to determine the optimal lactation period. We therefore examined whether the random regression (RR) coefficient from 306 to 450 DIM (M2) can be predicted from those during the first 305 DIM (M1) by using a RR model. METHODS: We analyzed test-day milk records from 85,690 Holstein cows in their first lactations and 131,727 cows in their later (second to fifth) lactations. Data in M1 and M2 were analyzed separately by using different single-trait RR animal models. We then performed a multiple regression analysis of the RR coefficients of M2 on those of M1 during the first and later lactations. RESULTS: The first-order Legendre polynomials were practical covariates of RR for the milk yields of M2. All RR coefficients for the additive genetic (AG) effect and the intercept for the permanent environmental (PE) effect of M2 had moderate to strong correlations with the intercept for the AG effect of M1. The coefficients of determination for multiple regression of the combined intercepts for the AG and PE effects of M2 on the coefficients for the AG effect of M1 were moderate to high. The daily milk yields of M2 predicted by using the RR coefficients for the AG effect of M1 were highly correlated with those obtained by using the coefficients of M2. CONCLUSION: Milk production after 305 DIM can be predicted by using the RR coefficient estimates of the AG effect during the first 305 DIM.

Effects of heat stress on production, somatic cell score and conception rate in Holsteins.

We examined the effects of heat stress (HS) on production traits, somatic cell score (SCS) and conception rate at first insemination (CR) in Holsteins in Japan. We used a total of 228 242 records of milk, fat and protein yields, and SCS for the first three lactations, as well as of CR in heifers and in first- and second-lactation cows that had calved for the first time between 2000 and 2012. Records from 47 prefectural weather stations throughout Japan were used to calculate the temperature-humidity index (THI); areas were categorized into three regional groups: no HS (THI < 72), mild HS (72 ≤ THI < 79), and moderate HS (THI ≥ 79). Trait records from the three HS-region groups were treated as three different traits and trivariate animal models were used. The genetic correlations between milk yields from different HS groups were very high (0.91 to 0.99). Summer calving caused the greatest increase in SCS, and in the first and second lactations this increase became greater as THI increased. In cows, CR was affected by the interaction between HS group and insemination month: with summer and early autumn insemination, there was a reduction in CR, and it was much larger in the mild- and moderate-HS groups than in the no-HS group.

Effects of a breeding scheme combined by genomic pre-selection and progeny testing on annual genetic gain in a dairy cattle population.

The effectiveness of the incorporation of genomic pre-selection into dairy cattle progeny testing (GS-PT) was compared with that of progeny testing (PT) where the fraction of dam to breed bull (DB) selected was 0.01. When the fraction of sires to breed bulls (SB) selected without being progeny tested to produce young bulls (YB) in the next generation was 0.2, the annual genetic gain from GS-PT was 13% to 43% greater when h(2) = 0.3 and 16% to 53% greater when h(2) = 0.1 compared with that from PT. Given h(2) = 0.3, a selection accuracy of 0.8 for both YB and DB, and selected fractions of 0.117 for YB and 0.04 for DB, GS-PT produced 40% to 43% greater annual genetic gain than PT. Given h(2) = 0.1, a selection accuracy of 0.6 for both YB and DB, and selected fractions of 0.117 for YB and 0.04 for DB, annual genetic gain from GS-PT was 48% to 53% greater than that from PT. When h(2) = 0.3, progeny testing capacity had little effect on annual genetic gain from GS-PT. However, when h(2) = 0.1, annual genetic gain from GS-PT increased with increasing progeny testing capacity.

Genetic evaluation of the longevity of the Holstein population in Japan using a Weibull proportional hazard model.

The objective of this study was to confirm the stability of the genetic estimation of longevity of the Holstein population in Japan. Data on the first 10 lactation periods were obtained from the Livestock Improvement Association of Japan. Longevity was defined as the number of days from first calving until culling or censoring. DATA1 and DATA2 included the survival records for the periods 1991-2003 and 1991-2005, respectively. The proportional hazard model included the effects of the region-parity-lactation stage-milk yield class, age at first calving, the herd-year-season, and sire. The heritabilities on an original scale of DATA1 and DATA2 were 0.119 and 0.123, respectively. The estimated transmitting abilities (ETAs) of young sires in DATA1 may have been underestimated, but coefficient δ, which indicated the bias of genetic trend between DATA1 and DATA2, was not significant. The regression coefficient of ETAs between DATA1 and DATA2 was very close to 1. The proportional hazard model could steadily estimate the ETA for longevity of the sires in Japan.