Genetic analysis of somatic cell scores in US Holsteins with a Bayesian mixture model

The objective of this study was to apply finite mixture models to field data for somatic cell scores (SCS) for estimation of genetic parameters. Data were approximately 170,000 test-day records for SCS from first-parity Holstein cows in Wisconsin. Five different models of increasing level of complexity were fitted. Model 1 was the standard single-component model, and the others were 2-component Gaussian mixtures consisting of similar but distinct linear models. All mixture models (i.e., 2 to 5) included separate means for the 2 components. Model 2 assumed entirely homogeneous variances for both components. Models 3 and 4 assumed heterogeneous variances for either residual (model 3) or genetic and permanent environmental variances (model 4). Model 5 was the most complex, in which variances of all random effects were allowed to vary across components. A Bayesian approach was applied and Gibbs sampling was used to obtain posterior estimates. Five chains of 205,000 cycles were generated for each model. Estimates of variance components were based on posterior means. Models were compared by use of the deviance information criterion. Based on the deviance information criterion, all mixture models were superior to the linear model for analysis of SCS. The best model was one in which genetic and PE variances were heterogeneous, but residual variances were homogeneous. The genetic analysis suggested that SCS in healthy and infected cattle are different traits, because the genetic correlation between SCS in the 2 components of 0.13 was significantly different from unity.     

Application of a finite mixture model to somatic cell scores of Italian goats

The objectives of this study were to apply a finite mixture model (FMM) to data for somatic cell count in goats and to compare the fit of the FMM with that of a standard linear mixed effects model. Bacteriological information was used to assess the ability of the model to classify records from healthy or infected goats. Data were 4518 observations of somatic cell score (SCS) and bacterial infection from both udder halves of 310 goats from 5 herds in Northern Italy. The records were from a complete production season, and were taken monthly from February to November 2000. Explanatory factors in both models included a 3-parameter regression on days in milk (DIM); fixed class effects of herd-test-day, parity group, and udder side (left or right); and random effects of goat and udder half within goat. In addition, the 2-component FMM included a fixed mean for the second component of the model (theoretically corresponding to infected udder halves), as well as an unknown probability of membership to a given putative infection status. A Bayesian statistical approach was used for the analysis with Gibbs sampling used to obtain draws from posterior distributions of parameters of interest. Two sampling chains of 200,000 cycles each were generated for each model. The FMM yielded a much lower estimate of residual variance than the standard model (1.28 vs. 3.02 SCS²), and a slightly higher estimate for the between-goat variance (1.79 vs. 1.48). The deviance information criterion (DIC) was used to compare the fit of the 2 models. The DIC was much lower for the FMM, indicating a better fit to the data. The FMM was able to classify correctly 60 and 48% of the healthy and infected observations, respectively. This was slightly higher than what would be expected from random classification, but not high enough for useful mastitis diagnosis. Nevertheless, increased precision of genetic evaluation is the goal of applying the FMM, rather than timely and accurate mastitis diagnosis. The results suggest that more research on FMM for SCS is merited and necessary for proper application.     

Detection of mastitis in dairy cattle by use of mixture models for repeated somatic cell scores: a Bayesian approach via Gibbs sampling

The distribution of somatic cell scores could be regarded as a mixture of at least two components depending on a cow's udder health status. A heteroscedastic two-component Bayesian normal mixture model with random effects was developed and implemented via Gibbs sampling. The model was evaluated using datasets consisting of simulated somatic cell score records. Somatic cell score was simulated as a mixture representing two alternative udder health statuses ("healthy" or "diseased"). Animals were assigned randomly to the two components according to the probability of group membership (Pm). Random effects (additive genetic and permanent environment), when included, had identical distributions across mixture components. Posterior probabilities of putative mastitis were estimated for all observations, and model adequacy was evaluated using measures of sensitivity, specificity, and posterior probability of misclassification. Fitting different residual variances in the two mixture components caused some bias in estimation of parameters. When the components were difficult to disentangle, so were their residual variances, causing bias in estimation of Pm and of location parameters of the two underlying distributions. When all variance components were identical across mixture components, the mixture model analyses returned parameter estimates essentially without bias and with a high degree of precision. Including random effects in the model increased the probability of correct classification substantially. No sizable differences in probability of correct classification were found between models in which a single cow effect (ignoring relationships) was fitted and models where this effect was split into genetic and permanent environmental components, utilizing relationship information. When genetic and permanent environmental effects were fitted, the between-replicate variance of estimates of posterior means was smaller because the model accounted for random genetic drift.     

Genetic analysis of somatic cell score in Norwegian cattle using random regression test-day models

The dataset used in this analysis contained a total of 341,736 test-day observations of somatic cell scores from 77,110 primiparous daughters of 1965 Norwegian Cattle sires. Initial analyses, using simple random regression models without genetic effects, indicated that use of homogeneous residual variance was appropriate. Further analyses were carried out by use of a repeatability model and 12 random regression sire models. Legendre polynomials of varying order were used to model both permanent environmental and sire effects, as did the Wilmink function, the Lidauer-M?ntysaari function, and the Ali-Schaeffer function. For all these models, heritability estimates were lowest at the beginning (0.05 to 0.07) and higher at the end (0.09 to 0.12) of lactation. Genetic correlations between somatic cell scores early and late in lactation were moderate to high (0.38 to 0.71), whereas genetic correlations for adjacent DIM were near unity. Models were compared based on likelihood ratio tests, Bayesian information criterion, Akaike information criterion, residual variance, and predictive ability. Based on prediction of randomly excluded observations, models with 4 coefficients for permanent environmental effect were preferred over simpler models. More highly parameterized models did not substantially increase predictive ability. Evaluation of the different model selection criteria indicated that a reduced order of fit for sire effects was desireable. Models with zeroth- or first-order of fit for sire effects and higher order of fit for permanent environmental effects probably underestimated sire variance. The chosen model had Legendre polynomials with 3 coefficients for sire, and 4 coefficients for permanent environmental effects. For this model, trajectories of sire variance and heritability were similar assuming either homogeneous or heterogeneous residual variance structure.     

Genetic associations between clinical mastitis and somatic cell score in early first-lactation cows

The objectives of this study were to examine genetic associations between clinical mastitis and somatic cell score (SCS) in early first-lactation cows, to estimate genetic correlations between SCS of cows with and without clinical mastitis, and to compare genetic evaluations of sires based on SCS or clinical mastitis. Clinical mastitis records from 15 d before to 30 d after calving and first test-day SCS records (from 6 to 30 d after calving) from 499,878 first-lactation daughters of 2,043 sires were analyzed. Results from a bivariate linear sire model analysis of SCS in cows with and without clinical mastitis suggest that SCS is a heterogeneous trait. Heritability of SCS was 0.03 for mastitic cows and 0.08 for healthy cows, and the genetic correlation between the 2 traits was 0.78. The difference in rank between sire evaluations based on SCS of cows with and without clinical mastitis varied from −994 to 1,125, with mean 0. A bivariate analysis with a threshold-liability model for clinical mastitis and a linear Gaussian model for SCS indicated that heritability of liability to clinical mastitis is at least as large as that of SCS in early lactation. The mean (standard deviation) of the posterior distribution of heritability was 0.085 (0.006) for liability to clinical mastitis and 0.070 (0.003) for SCS. The posterior mean (standard deviation) of the genetic correlation between liability to clinical mastitis and SCS was 0.62 (0.03). A comparison of sire evaluations showed that genetic evaluation based on SCS was not able to identify the best sires for liability to clinical mastitis. The association between sire posterior means for liability to clinical mastitis and sire predicted transmitting ability for SCS was far from perfect.     

Recursive relationships between milk yield and somatic cell score of Canadian Holsteins from finite mixture random regression models

Finite mixture, multiple-trait, random regression animal models with recursive links between phenotypes for milk yield and somatic cell score (SCS) on the same test-day were applied to first lactation 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. Causal links between phenotypes for milk yield and SCS were fitted separately for records from healthy cows and cows with a putative, subclinical form of mastitis. Bayesian methods via Gibbs sampling were used for the estimation of model parameters. Bayes factors indicated superiority of the model with recursive link from milk to SCS over the reciprocal recursive model and the standard multiple-trait model. Differences between models measured by other, single-trait model comparison criteria (i.e., weighted mean squared error, squared bias, and correlation between observed and expected data) were negligible. Approximately 20% of test-day records were classified as originating from cows with mastitis in recursive mixture models. The proportion of records from cows infected with mastitis was largest at the beginning of lactation. Recursive mixture models exhibited different distributions of data from healthy and infected cows in different parts of lactation. A negative effect of milk to SCS (up to −0.15 score points for every kilogram of milk for healthy cows from 5 to 45 d in milk) was estimated for both mixture components (healthy and infected) in all stages of lactation for the most plausible model. The magnitude of this effect was stronger for healthy cows than for cows infected with mastitis. Different patterns of genetic and environmental correlations between milk and SCS for healthy and infected records were revealed, due to heterogeneity of structural coefficients between mixture components. Estimated breeding values for SCS from the best fitting model for sires of infected daughters were more related to estimated breeding values for the same trait from the regular multiple-trait model than evaluations for sires of mastitis-free cows.     

Genetic improvement of mastitis resistance: validation of somatic cell score and clinical mastitis as selection criteria

Mean daughter deviations for clinical mastitis among second-crop daughters were regressed on predicted transmitting abilities for clinical mastitis and lactation mean somatic cell score in first-crop daughters to validate the predictive ability of these traits as selection criteria for reduced incidence of clinical mastitis. A total of 321 sires had 684,897 second-crop daughters, while predicted transmitting abilities were calculated for 2159 sires, based on 495,681 records of first-crop daughters. Predictive ability, as a measure of efficiency of selection, was 23 to 43% higher for clinical mastitis than for lactation mean somatic cell score. Compared to single-trait selection, predictive ability improved 8 to 13% from utilizing information on both traits. The relative weight that should be assigned to standardized predicted transmitting abilities from univariate genetic analyses were 60 to 67% for clinical mastitis and 33 to 40% for lactation mean somatic cell score. No significant nonlinear genetic relationship between the two traits was found.     

Genetic improvement in mastitis resistance: comparison of selection criteria from cross-sectional and random regression sire models for somatic cell score

Several selection criteria for reducing incidence of mastitis were developed from a random regression sire model for test-day somatic cell score (SCS). For comparison, sire transmitting abilities were also predicted based on a cross-sectional model for lactation mean SCS. Only first-crop daughters were used in genetic evaluation of SCS, and the different selection criteria were compared based on their correlation with incidence of clinical mastitis in second-crop daughters (measured as mean daughter deviations). Selection criteria were predicted based on both complete and reduced first-crop daughter groups (261 or 65 daughters per sire, respectively). For complete daughter groups, predicted transmitting abilities at around 30 d in milk showed the best predictive ability for incidence of clinical mastitis, closely followed by average predicted transmitting abilities over the entire lactation. Both of these criteria were derived from the random regression model. These selection criteria improved accuracy of selection by approximately 2% relative to a cross-sectional model. However, for reduced daughter groups, the cross-sectional model yielded increased predictive ability compared with the selection criteria based on the random regression model. This result may be explained by the cross-sectional model being more robust, i.e., less sensitive to precision of (co)variance components estimates and effects of data structure.     

Dairy management practices associated with incidence rate of clinical mastitis in low somatic cell score herds in France

An epidemiological prospective study was carried out in French dairy herds with Holstein, Montbéliarde, or Normande cows and with low herd somatic cell scores. The objective was to identify dairy management practices associated with herd incidence rate of clinical mastitis. The studied herds were selected on a national basis, clinical cases were recorded through a standardized system, and a stable dairy management system existed. In the surveyed herds, mean milk yield was 7420 kg/cow per yr and mean milk somatic cell score was 2.04 (132,000 cells/mL). Overdispersion Poisson models were performed to investigate risk factors for mastitis incidence rate. From the final model, the herds with the following characteristics had lower incidence rates of clinical mastitis: 1) culling of cows with more than 3 cases of clinical mastitis within a lactation; 2) more than 2 person-years assigned to dairy herd management; 3) balanced concentrate in the cow basal diet. Moreover, herds with the following characteristics had higher incidence rates of clinical mastitis: 1) milking cows loose-housed in a straw yard; 2) no mastitis therapy performed when a single clot was observed in the milk; 3) clusters rinsed using water or soapy water after milking a cow with high somatic cell count; 4) 305-d milk yield >7435 kg; 5) herd located in the South region; 6) herd located in the North region; 7) cows with at least 1 nonfunctional quarter; and 8) premilking holding area with a slippery surface. The underlying mechanisms of some highlighted risk factors, such as milk production level and dietary management practices, should be investigated more thoroughly through international collaboration.     

Relationship between intramammary infection prevalence and somatic cell score in commercial dairy herds

We examined consistency of the relationship between intramammary infection (IMI) and somatic cell score (SCS) across several classes of cow, herd, and sampling time variables. Microbial cultures of composite milk samples were performed by New York Quality Milk Production Services from 1992 to 2004. SCS was from the most recent Dairy Herd Improvement test before IMI sampling. Records were analyzed from 79,308 cows in 1,124 commercial dairy herds representing a broad range of production systems. Three binary dependent variables were presence or absence of contagious IMI, environmental IMI, and all IMI. Independent variables in the initial models were SCS, SCS², lactation number, days in milk, sample day milk yield, use of coliform mastitis vaccine, participant type (required by regulation or voluntary), production system (type of housing, milking system, and herd size), season of sampling, year of sampling, and herd; also the initial models included interactions of SCS and SCS² with other independent variables, except herd and milk yield. Interaction terms characterize differences in the IMI-SCS relationship across classes of the independent variables. Models were derived using the Glimmix macro in SAS (SAS Institute Inc., Cary, NC) with a logistic link function and employing backward elimination. The final model for each dependent variable included all significant independent variables and interactions. Simplified models omitted SCS² and all interactions with SCS. Interactions of SCS with days in milk, use of coliform mastitis vaccine, participant type, season, and year were not significant in any of the models. Interaction of SCS with production system was significant for the all IMI model, whereas interaction of SCS with lactation number was significant for the environmental and all IMI models. Each 1-point increase in SCS (or doubling of somatic cell count) was associated with a 2.3, 5.5, and 9.1% increase in prevalence of contagious, environmental, and all IMI, respectively. Empirical receiver operator characteristic curves and areas under the curve were derived for final and simplified models. The areas under the curve for simplified and final models within each type of IMI differed by 0.009 or less. We concluded that the relationship of IMI with SCS was generally stable over time and consistent across seasons, production systems, and cow factors.     

Joint analysis of quantitative trait loci for clinical mastitis and somatic cell score on five chromosomes in three Nordic dairy cattle breeds

Five chromosomes were selected for joint quantitative trait loci (QTL) analyses for clinical mastitis (CM) and somatic cell score (SCS) in 3 breeds: Finnish Ayrshire (FA), Swedish Red and White (SRB), and Danish Red (DR). In total, 19 grandsires and 672 sons in FA, 19 grandsires and 499 sons in SRB, and 8 grandsires and 258 sons in DR were used in the study. These individuals were genotyped with the 61 microsatellite markers used in any of the previous QTL scans on the selected chromosomes. Within-family QTL analyses based on linear regression models were carried out for CM and SCS to identify the segregating sires for each region. On the segregating families, joint single-trait and 2-trait analyses were performed using variance components models. The analyses confirmed that QTL affecting CM or SCS, or both, segregate on Bos taurus autosomes (BTA) 9, 11, 14, and 18, whereas a QTL on BTA29 could not be confirmed. Our results indicate that there may be at least 2 linked QTL on BTA9, one that primarily affects CM and a second that primarily affects SCS. On chromosomes BTA11, 14, and 18, the joint analyses were only significant for SCS.     

Associations among body condition score, body weight, somatic cell count, and clinical mastitis in seasonally calving dairy cattle

The objective of this study was to determine if an association existed among body condition score (BCS), body weight (BW), and udder health, as indicated by somatic cell score (SCS) and cases of clinical mastitis (CM). The data consisted of 2,635 lactations from Holstein-Friesian (n = 523) and Jersey (n = 374) cows in a seasonal calving pasture-based research herd between the years 1986 and 2000, inclusive. Increased BCS at calving was associated with reduced SCS in first- and second-parity cows, and greater SCS in cows of third parity or greater. This relationship persisted for most BCS traits throughout lactation. Body weight was positively associated with SCS, although the effect was greater in Jersey cows than in Holstein-Friesians. Increased BCS and BW loss in early lactation were associated with lower SCS and a reduced probability of a high test-day SCC. Body condition score was not significantly related to CM with the exception of a curvilinear relationship between the daily rate of BCS change to nadir and CM in early lactation. Several BW variables were positively associated with a greater likelihood of CM. Nevertheless, most associations with udder health lacked biological significance within the ranges of BCS and BW generally observed on-farm. Results are important in assuring the public that modern dairy systems, where cows are subjected to substantial amounts of BCS mobilization in early lactation, do not unduly compromise cow udder health.     

Spatial and temporal patterns of herd somatic cell score in france

Spatial and temporal patterns of annual milk somatic cell score (ASCS) were explored in French dairy herds between 1996 and 2000 to detect regional singularities for risk of mastitis. A new cluster detection method was used, which was adapted to continuous variables and which allowed ASCS variation factors to be taken into account. The statistical unit was the herd-year. A linear regression model for each year allowed adjustment for breed, mean parity, number of calvings for each season, herd size, and farm altitude. Cluster detection was performed on raw data and on residuals of the model through a method based on the Hellinger distance between spatial distributions. The Hellinger distance between farm distributions was computed at different levels of ASCS (or residuals). Temporal ASCS patterns were explored using a computation of correlations and comparisons between spatial structures of the different years. The general ASCS trend over the study period was a decrease. The global Hellinger distance, which was higher than what could have been randomly expected for each of the 5 yr, indicated a significant spatial cluster formation. Cluster mapping over the 5 yr identified several areas, which sometimes differed between detection using raw data and that using ASCS residuals. Temporal correlations between ASCS residuals for each year were positive and decreasing, and 1996 and 2000 appeared spatially different from the other years. The more affected areas were regions that were not specialized in dairy production. During the study period, 2 progressive movements were detected, corresponding to a disappearance of clusters in the northwest and an increase of clusters in the southwest. Cluster detection could aid in the identification of new risk factors that are relevant at different spatial scales, and could help local organizations to supervise the risk of mastitis, and improve udder health management.     

Refined positioning of a quantitative trait locus affecting somatic cell score on chromosome 18 in the German Holstein using linkage disequilibrium

Combined linkage and linkage disequilibrium analysis (LALD) was conducted to more accurately map a previously reported quantitative trait locus (QTL) affecting somatic cell score on bovine chromosome 18. A granddaughter design consisting of 6 German Holstein grandsire families with 1,054 progeny-tested genotyped sons was used in this study. Twenty microsatellite markers, 5 single nucleotide polymorphisms, and an erythrocyte antigen marker with an average marker spacing of 1.95 cM were analyzed along a chromosomal segment of 50.80 cM. Variance components were estimated and restricted maximum likelihood test statistics were calculated at the midpoint of each marker interval. The test statistics calculated in single-QTL linkage analysis exceeded the genome-wide significance threshold at several putative QTL positions. Using LALD, we were successful in assigning a genome-wide significant QTL to a confidence interval of 10.8 cM between the markers ILSTS002 and BMS833. The QTL in this marker interval was estimated to be responsible for between 5.89 and 13.86% of the genetic variation in somatic cell score. In contrast to the single-QTL linkage analysis model, LALD analyses with a 2-QTL model confirmed the position of one QTL, but gave no conclusive evidence for the existence or position of a second QTL. Ultimately, the QTL position was narrowed down considerably compared with previous results with a refined confidence interval of less than 11 cM.     

Genetic correlation patterns between somatic cell score and protein yield in the Italian Holstein-Friesian population

Genetic parameters for somatic cell score (SCS) in the Italian Holstein-Friesian population were estimated addressing the pattern of genetic correlation with protein yield in different parities (first, second, and third) and on different days in milk within each parity. Three approaches for parameter estimation were applied using random samples of herds from the national database of the Italian Holstein Association. Genetic correlations for lactation measures (305-d protein yield and lactation SCS) were positive in the first parity (0.31) and close to zero in the second (0.01) and third (0.09) parities. These results indicated that larger values of SCS were genetically associated with increased production. The second and third sets of estimates were based on random regression test-day models, modeling the shape of lactation curve with the Wilmink function and fourth-order Legendre polynomials, respectively. Genetic correlations from both random regression models showed a specific pattern associated with days in milk within and across parities. Estimates varied from positive to negative in the first and second parity, and from null to negative in the third parity. Patterns were similar for both random regression models. The average overall correlation between SCS and protein yield was zero or slightly positive in the first lactation and ranged from zero to negative in later lactations. Correlation estimates differed by parity and stage of lactation. They also demonstrated the dubiousness of applying a single genetic correlation measure between SCS and protein in setting selection strategies. Differences in magnitude and the sign of genetic correlations between SCS and yields across and within parities should be accounted for in selection schemes.     

Relationship of test-day somatic cell score with test-day and lactation milk yields

To determine the relationship of test-day (TD) somatic cell score (SCS) to TD and lactation milk yields, 1,320,590 records from Holstein first and second calvings from 1995 through 2002 were examined. All lactations had recorded yield and SCS for at least the first 4 TD. Least square analyses were conducted for yields on TD 2 through 10 within herd and cow. The model included regressions on current TD SCS and mean SCS of all previous TD with separate estimates by parity; effects for parity and calving year were included as well as regression on days in milk on TD 1. Corresponding analyses were conducted without regression on current SCS. An analysis of lactation yield was performed with a similar model and regression on all TD SCS. The SCS was highest most often on TD 1 for parity 1 (22.5%) and on TD 10 for parity 2 (18.5%). Regression of TD milk yield on mean of previous TD SCS was highest during the latter half of lactation (maximum of -0.346 kg/SCS unit on TD 9) for parity 1 and during TD through 7 (maximum of -0.366 kg/SCS unit on TD 4) for parity 2. Regression of TD yield on current TD SCS tended to be larger for later lactation. Regression of lactation yield on TD SCS was negative and important for TD 1 through 6 for parity 1 and for all TD for parity 2. To minimize milk loss, mastitis control is most important immediately pre- and postcalving for parity 1 and throughout lactation for parity 2.     

Management practices from questionnaire surveys in herds with very low somatic cell score through a national mastitis program in France

French dairy herds (n = 534) were enrolled in the National 'Zero Mastitis Objective' Program to highlight management practices characterizing very low somatic cell score (SCS) herds. The herds studied were stratified into 2 groups. The first group (LOW) included herds within the first 5 percentiles and the second group (MED) herds within the 50 to 55 percentiles of herds on the basis of mean SCS for the 36 mo preceding the program. Potential explanatory variables, collected through questionnaire surveys, were analyzed using multistep logistic regression models. Twenty-six variables were significant factors in the final models, in which 18 were considered as primary factors for very low SCS. The probability for a herd belonging to the LOW group was associated with: (1) regular use of teat spraying; (2) herdsman precise in his techniques; (3) less than 1 person-year used at activities other than dairy herd; (4) teat dipping after mammary infusion at dry off; (5) heifers kept in a calving pen around parturition; (6) cows locked in feed-line lockups after milking; (7) dry cows with prepartum Ca restriction; (8) heifers on a nondamp pasture; (9) cows culled when at least one damaged teat; (10) heifers at pasture not drinking water from a river; and (11) disinfecting teat ends with alcohol before intramammary infusion at dry off. The probability for a herd belonging to the MED group was associated with: (1) milking cows housed in a straw yard; (2) checking heifers for mastitis only beginning at 2-wk prepartum; (3) no mastitis treatment when at least one clot was observed in milk at successive milkings; (4) distance of herdsman's house to cowshed >300 m; (5) only dirty teats washed before milking; (6) free access of cows from pasture to cowshed during bad weather; and (7) more than 18% of spring calvings. The variables associated with very low SCS should be applied as part of a thorough mastitis-control program adapted to each herd.     

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.     

Cow, farm, and herd management factors in the dry period associated with raised somatic cell counts in early lactation

This study investigated cow characteristics, farm facilities, and herd management strategies during the dry period to examine their joint influence on somatic cell counts (SCC) in early lactation. Data from 52 commercial dairy farms throughout England and Wales were collected over a 2-yr period. For the purpose of analysis, cows were separated into those housed for the dry period (6,419 cow-dry periods) and those at pasture (7,425 cow-dry periods). Bayesian multilevel models were specified with 2 response variables: ln SCC (continuous) and SCC >199,000 cells/mL (binary), both within 30 d of calving. Cow factors associated with an increased SCC after calving were parity, an SCC >199,000 cells/mL in the 60 d before drying off, increasing milk yield 0 to 30 d before drying off, and reduced DIM after calving at the time of SCC estimation. Herd management factors associated with an increased SCC after calving included procedures at drying off, aspects of bedding management, stocking density, and method of pasture grazing. Posterior predictions were used for model assessment, and these indicated that model fit was generally good. The research demonstrated that specific dry-period management strategies have an important influence on SCC in early lactation.