A structural equation model for describing relationships between somatic cell score and milk yield in first-lactation dairy cows

Relationships between production and diseases may involve recursive or simultaneous effects between traits. Four structural equation models (SEqM) for somatic cell score and milk yield, with varying specifications for the effects relating the 2 traits, were compared. Data consisted of repeated records of milk yield and somatic cell score of 33,453 first-lactation daughters of 245 Norwegian Red sires that had their first progeny test in 1991 and 1992. All models included random effects of the sire and of the cow and were fitted using the LISREL software. The Bayesian information criterion clearly favored a model with a recursive effect from somatic cell score on milk yield over the 3 other models fitted (absence of recursive effects; an effect from milk yield on somatic cell score; simultaneity of effects between the 2 traits). This provides evidence that the negative association between milk yield and somatic cell score is more likely due to an effect of infection (measured indirectly by the somatic cell score) on production than to a dilution effect. Estimates indicated that a mastitis event would reduce milk yield in the following 15 d by about 900 g/d. The estimated genetic (co)variances did not change sizably when the specification of recursive or simultaneous effects was varied. However, estimates of the phenotypic covariance were altered when a recursive effect from somatic cell score on milk yield was included in the model.     

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

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

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.     

Negative influence of high maternal milk production before and after conception on offspring survival and milk production in dairy cattle

There is a paucity of studies on the effect of intrauterine conditions on subsequent progeny performance in dairy cows. Using a large national data set on Irish Holstein-Friesian dairy cows, the objective of this study was to determine if intrauterine conditions, quantified by a maternal genetic variance component, significantly affected milk production, age at first calving, calving interval, somatic cell score (natural logarithm of somatic cell count) and survival in first-, second-, and third-parity female offspring. Maternal genetic variance for each trait in each parity was estimated in a linear mixed model which included, other than fixed effects, direct additive genetic, maternal genetic, cytoplasmic and permanent environmental effect of the dam, and residual component. A covariance was also estimated between the direct additive and maternal genetic components where possible. Because calves in Irish dairy herds are separated from dams at birth, a significant maternal genetic variance (with all other random effects in the model) indicates a prepartum maternal effect. A significant maternal genetic variance was estimated for 305-d milk yield in first and third lactation, somatic cell score in first lactation, and survival to second lactation from 188,144 lactations on 80,881 animals. Where estimated, a negative correlation existed between the direct additive and maternal genetic components. Regression of maternal mixed model solutions on dam milk production at different stages relative to conception revealed that greater milk yield preconception and during gestation was associated with reduced survival and milk yield and greater somatic cell count in the progeny. This study suggests that offspring survival and performance are affected by prepartum conditions that offspring experience as an oocyte, embryo, or fetus, one of which is mediated through milk production (or factors related to milk production) of the dam.     

Effects of postpartum oral calcium supplementation on milk yield,milk composition,and reproduction in multiparous Jersey and Jersey × Holstein crossbreed cows

The objective of the present study was to evaluate the effects of postpartum oral calcium supplementation on milk yield, energy-corrected milk yield, milk fat concentration, milk protein concentration, and somatic cell count linear score across the first 3 monthly tests postpartum, peak milk yield, risk of pregnancy at first service, and hazard of pregnancy by 150 d in milk on 1,129 multiparous Jersey and Jersey × Holstein crossbreed cows from 2 commercial dairies. After calving, cows were systematically assigned to control (no oral calcium supplementation; n = 567) or oral calcium supplementation at 0 and 1 d in milk (oral Ca; 50 to 60 g of calcium as boluses; n = 562). Monthly test milk yield, composition, and somatic cell count information was obtained from the Dairy Herd Improvement Association. Herd records were used for reproductive data. Statistical analysis was conducted using generalized multiple linear, Poisson, and Cox's hazard regressions. Treatment effects were evaluated considering cow-level information available at parturition (parity, breed, previous lactation milk yield, previous lactation length, dry period length, gestation length, body condition, and locomotion score at calving, calving ease, and calf sex). In addition, for a subset of cows serum calcium concentration before treatment administration was evaluated (n = 756). Overall, oral calcium supplementation did not affect the evaluated productive and reproductive variables. However, effects conditional to previous lactation length and calving locomotion score were observed. Milk yield and energy-corrected milk yield across the first 3 monthly tests were 1.8 kg/d higher for supplemented cows with a previous lactation length within the fourth quartile, compared with control cows on the same quartile. Energy-corrected milk yield tended to be 1.1 kg/d lower for supplemented cows with a previous lactation length within the first quartile, compared with control counterparts. Peak milk yield tended to be 1.6 kg higher for supplemented cows with a calving locomotion score ≥2, compared with control cows with the same locomotion score. Treatment effects were not conditional to serum calcium concentration before treatment administration. Our results suggest that postpartum oral calcium supplementation effects are conditional to cow-level factors such as previous lactation length and calving locomotion score in multiparous Jersey and Jersey × Holstein crossbreed cows.     

Association between somatic cell count during the first lactation and the cumulative milk yield of cows in Irish dairy herds

Reduced potential milk yield is an important component of mastitis costs in dairy cows. The first aim of this study was to assess associations between somatic cell count (SCC) during the first lactation, and cumulative milk yield over the first lactation and subsequent lifetime of cows in Irish dairy herds. The second aim was to assess the association between SCC at 5 to 30 d in milk during parity 1 (SCC1), and SCC over the entire first lactation for cows in Irish dairy herds. The data set studied included records from 51,483 cows in 5,900 herds. Somatic cell count throughout the first lactation was summarized using the geometric mean and variance of SCC. Data were analyzed using linear models that included random effects to account for the lack of independence between observations, and herd-level variation in coefficients. Models were developed in a Bayesian framework and parameters were estimated from 10,000 Markov chain Monte Carlo simulations. The final models were a good fit to the data. A 1-unit increase in mean natural logarithm SCC over the first lactation was associated with a median decrease in first lactation and lifetime milk yield of 135 and 1,663 kg, respectively. A 1-unit increase in the variance of natural logarithm SCC over the first lactation was associated with a median decrease in lifetime milk yield of 719 kg. To demonstrate the context of lifetime milk yield results, microsimulation was used to model the trajectory of individual cows and evaluate the expected outcomes for particular changes in herd-level geometric mean SCC over the first lactation. A 75% certainty of savings of at least €199/heifer in the herd was detected if herd-level geometric mean SCC over the first lactation was reduced from ≥120,000 to ≤72,000 cells/mL. The association between SCC1 and SCC over the remainder of the first lactation was highly herd dependent, indicating that control measures for heifer mastitis should be preferentially targeted on an individual-herd basis toward either the pre- and peripartum period, or the lactating period, to optimize the lifetime milk yield of dairy cows.     

Effect of herd environment on the genetic and phenotypic relationships among milk yield, conception rate, and somatic cell score in Holstein cattle

A total of 248,230 primiparous records of Holstein cows calving from 1987 to 1994 (daughters of 588 sires in 3042 herds) was used to evaluate potential genotype by environment interactions among mature equivalent milk yield, lactation mean somatic cell score, and conception rate at first service. Herds were classified into low and high environmental groups using three different criteria: standard deviation of herd mature equivalent milk yield, a combination of herd mature equivalent milk yield mean and standard deviation, and the herd mean of body weight at first calving divided by age at first calving. Genetic parameters were modeled by using multiple-trait linear mixed models and were fitted using the multiple-trait derivative-free software. Heritabilities for mature equivalent milk yield, lactation mean somatic cell score, and conception rate at first service were 0.221, 0.106, and 0.015 in low environment herds and 0.300, 0.093, and 0.009 in high environment herds, respectively. Genetic (and phenotypic) correlations between mature equivalent milk yield and lactation mean somatic cell score, mature equivalent milk yield and conception rate at first service, and lactation mean somatic cell score and conception rate at first service were 0.277, -0.417, and -0.209, (-0.049, -0.180, and -0.040) and 0.173, -0.318, and -0.144, (-0.087, -0.166, and -0.035) in low and high environment herds, respectively. The genetic correlations between pairs of traits were consistently smaller in high environment herds, suggesting that differences in management between the two environment levels lessened the antagonistic genetic association between the traits studied. A long-range plan for low environment herds should focus on improving the level of management, which would greatly reduce the unfavorable correlated changes in lactation mean somatic cell score and conception rate at first service associated with the genetic improvement of mature equivalent milk yield.     

Effect of gradual or abrupt cessation of milking at dry off on milk yield and somatic cell score in the subsequent lactation

The objective of this study was to assess the effect of milk cessation method (abrupt or gradual) at dry off on milk yield and somatic cell score (SCS) up to 120 d in milk during the subsequent lactation. Data from 428 cows from 8 dairy herds in Ohio were analyzed. Abrupt cessation cows kept the farm's regular milking schedule (2 or 3 times) through dry off and gradual cessation cows were milked once daily for the final week of lactation. Milk yield and SCS were collected using Dairy Herd Improvement Association test-day records. Aseptic quarter milk samples were collected approximately 1 wk before dry off, at dry off, and within 1 wk after calving for bacterial culture to determine the presence of intramammary infections. Overall, milk cessation method was not significantly associated with either milk yield or SCS in early lactation; however, interaction between the milk cessation method and herd was highly significant. Cows producing greater amounts of milk around dry off had significantly higher SCS in the following lactation. Shorter dry periods were significantly associated with decreased milk yield in the following lactation, especially among abruptly dried off cows. Additionally, as expected, several other factors, such as parity of cows and stage of lactation, were significantly associated with both outcomes. No interactions between the milk cessation method and the other explanatory variables in the final models were significant. The results of the current study suggest that higher milk yield at dry off was associated with higher SCS in the following lactation, even though milk cessation method at the end of lactation had a varying effect on test-day milk yield and SCS in different herds during the first 120 d in milk in the following lactation. The specific herd characteristics influencing this could not be identified within this study, warranting further research.     

Hot topic: effects of frequent milking in early lactation on milk yield and udder health

A field study was conducted to evaluate the influence of milking frequency (3 or 6 times/d [3x or 6x, respectively]) during the initial 21 d of lactation on milk and milk component yield and mammary gland health as indicated by somatic cell count. During 2 seasons, spring and fall, multiparous cows were milked 6 times/d until d 21 of lactation and then returned to the 3 times/d frequency for the remainder of lactation (6x; n = 9). Multiparous cows milked 3 times/d from the beginning of lactation served as a control group (3x; n = 17). With the exception of milking frequency, all other aspects of management, including housing, milk harvesting, and feeding, were identical between the groups and were consistent with industry norms. Retrospective analysis of Dairy Herd Improvement Association records was used to evaluate milk yield, milk component yield, and somatic cell scores. As expected, 6x cows produced more milk on the first test day than 3x cows. Compared with 3x cows, higher milk yields persisted for 6x cows from test day 2 through 6, indicating a persistent effect of early lactation milking frequency on milk yield potential for that lactation. Milk component yield followed a similar pattern: 6x cows produced significantly more protein, fat, and total solids than did control cows throughout the study. With regard to udder health, 6x cows had lower somatic cell counts at the first test day relative to 3x cows and had reduced somatic cell scores for the first 3 mo of lactation, which suggests that early lactation milking frequency influences the mammary gland capacity to resist infection in addition to improving milk production efficiency.     

Hot topic: an association between a leptin single nucleotide polymorphism and milk and protein yield

Allelic variation (C to T transition that results in an Arg25Cys) in the leptin gene has been associated with increased fat deposition in beef cattle. We report that this same genetic variant is also present in dairy breeds. Body fat reserves play an important role in sustaining high milk production in early lactation, when energy intake is limited. To test for an association between the leptin single nucleotide polymorphism and milk productivity, we genotyped 416 Holstein cows and compared lactation performance data using a mixed model. Animals homozygous for the T allele produced more milk (1.5 kg/d vs. CC animals) and had higher somatic cell count linear scores, without significantly affecting milk fat or protein percent over the entire lactation. The increase in milk yield is most prominent in the first 100 d of lactation (2.44 kg/d), declining to 1.74 kg/d between 101 and 200 d in lactation. The milk yield advantage, observed in cows homozygous for the T allele, could represent a major economic advantage to dairy producers.     

Genetic parameters for milk somatic cell scores and relationships with production traits in French Lacaune dairy sheep

Records from 94,445 and 45,499 French Lacaune dairy ewes in first and second lactations, respectively, were used to estimate genetic parameters for somatic cell scores. Somatic cell count data came from an extensive recording scheme and sample testing that began in 1999 using the flocks enrolled in the official milk recording system. Somatic cell count data were from 2 to 4 test days per lactation. Lactation average and single test-day somatic cell scores were considered in multitrait sire models. The heritability estimate of lactation somatic cell score was close to 0.13 and similar for first and second parity. Heritabilities of somatic cell scores increased from first to fourth test day (from 0.07 to 0.11 in first lactation and from 0.05 to 0.13 in second lactation). Genetic correlations between somatic cell scores were high, usually more than 0.91, but lower between first test day and later test days in first lactation (0.64 to 0.88). The genetic correlations between lactation somatic cell score and milk yield, between lactation somatic cell score and fat content, and between lactation somatic cell score and protein content were 0.18, 0.04, and 0.03 in first lactation, respectively. The genetic antagonism between test day somatic cell score and milk yield measured in first lactation increased from beginning to the end of the lactation (0.05 to 0.23). This antagonism was slightly lower for somatic cell score in second lactation (from 0.09 to 0.14, and 0.08 for lactation mean). Environmental correlations in first lactation between lactation somatic cell score and milk yield, between lactation somatic cell score and fat content, and between lactation somatic cell score and protein content were -0.18, 0.13, and 0.30, respectively.     

Genetic analyses of mastitis data using animal threshold models and genetic correlation with production traits

In the present study, 6 different mastitis data sets of 3 dairy herds with an overall herd size of 3200 German Holstein cows were analyzed. Data collection periods included the first 50, 100, or 300 d of lactation. The 3 data collection periods were analyzed with a lactation model and a test-day model. All models were animal threshold models. Mastitis frequencies in the lactation model data sets varied between 29 and 45%, and varied between 3 and 6% in the test-day model data sets. Depending on the period of data collection, heritabilities of liability to mastitis in the lactation models were 0.05 (50 d), 0.06 (100 d), and 0.07 (300 d). In the test-day models, heritabilities were slightly higher with values of 0.09 (50 and 100 d), and 0.06 (300 d). Between lactation models, the rank correlations between the relative breeding values were high and varied between 0.86 and 0.94. Rank correlations between the relative breeding values of the test-day models ranged from 0.68 to 0.87. The rank correlations between the relative breeding values of lactation models and test-day models varied from 0.51 and 0.80. Genetic correlations between mastitis and milk production traits were estimated with a linear animal test-day model. The correlations with mastitis were 0.29 (milk yield), 0.30 (fat yield), 0.20 (fat content), 0.34 (protein yield), and 0.20 (protein content). The estimated genetic correlation between mastitis and somatic cell score was 0.84.     

Effect of subclinical mastitis detected in the first month of lactation on somatic cell count linear scores,milk yield,fertility, and culling of dairy cows in certified organic herds

It is well established that subclinical mastitis (SCM), characterized by somatic cell count (SCC) >200,000 cells/mL, has a negative effect on the productivity, reproductive performance, and survivability of cows from conventional dairy herds. However, in organic herds, where the use of antimicrobial drugs is restricted for the treatment and control of intramammary infections (IMI) in dairy cows, little is known about the effect of SCM on performance and survivability. The objective of this study was to evaluate whether SCM diagnosed during the first month of lactation was associated with SCC linear score dynamics, milk production, fertility, and culling of dairy cows in USDA-certified organic herds. We collected data from 2 organic herds in New Mexico and Texas. A total of 1,511 cows that calved between June 2018 and May 2019 were included in the study and were followed until month 10 of the current lactation. Cows with SCC >200,000 cells/mL in the first month of lactation were considered to have SCM. We used mixed linear regression models accounting for repeated measures to assess the effect of SCM on monthly milk production and SCC linear scores. We used Cox proportional hazards models to evaluate the effect of SCM on the risk of pregnancy and culling. We considered parity, farm, previous gestation length, stillbirth, twinning, dystocia, and 2- and 3-way interactions as potential confounders. Cows diagnosed with SCM during the first month of lactation produced less milk than cows without SCM. Cows with SCM had elevated SCC linear scores during their previous lactation and throughout the subsequent months of lactation compared to cows without SCM. The effect of SCM on SCC linear scores was more pronounced in multiparous than primiparous cows. Subclinical mastitis during the first month of lactation did not affect the likelihood of pregnancy during the first 300 d in milk. Cows with SCM in the first month were more likely to die or be culled during the 300 d of lactation than cows without SCM. We observed that elevated SCC in the first month of lactation had detrimental effects on the milk yield and survivability of dairy cows in USDA organic herds, but it did not affect reproductive performance. We demonstrated that cows with SCM diagnosed in the first month of lactation continued to have elevated SCC linear scores throughout their entire lactation, and that elevated SCC was carried over from the previous lactation.     

Relationships among severity and duration of clinical mastitis and sire transmitting abilities for somatic cell score,udder type traits,productive life,and protein yield

The objective of this study was to determine the relationships among severity and duration of clinical mastitis during first and second lactation and sire transmitting abilities for somatic cell score, udder type traits, productive life, and protein yield. Recording of clinical episodes began at first parturition for 1704 Holstein cows (in six Pennsylvania herds and one Nebraska herd) and continued into second lactation for 1055 of these cows. A total of 456 cows (sired by 168 bulls) had at least one clinical episode during first lactation, and 230 cows (sired by 100 bulls) had at least one clinical episode during second lactation. A severity code from 1 (normal milk) to 5 (acute systemic mastitis) was assigned daily (for up to 30 d after detection) to all quarters that had clinical mastitis. Only the severity codes for the first clinical episode to occur during first and second lactation are considered here. The initial and maximum severity codes, as well as the natural logarithms of both the sum of severity codes that were above normal (> 1) and the total days severity codes were above normal were regressed on herd (a classification variable), age at first calving, days in milk at clinical detection, and sire transmitting abilities taken one at a time. Linear and nonlinear effects were estimated for sire transmitting abilities. Separate analyses were conducted on dependent variables that considered severity and duration of clinical mastitis from: all organisms, coagulase-negative staphylococci, coliform species, streptococci other than Streptococcus agalactiae, and the most common environmental organisms (coliform species and streptococci other than Streptococcus agalactiae). Daughters of sires that transmit the lowest somatic cell score had the least severe and shortest clinical episodes from environmental organisms during first lactation. Selection for lower somatic cell score may reduce the severity and duration of clinical episodes from environmental organisms during first lactation.     

Dry period length in US Jerseys: characterization and effects on performance

The objectives of this research were to characterize dry period lengths for US Jerseys, determine the effects of days dry (DD) on subsequent lactation actual milk, fat, and protein yields, fat and protein percentages, somatic cell score (SCS), and days open (DO), and to determine the dry period length that maximizes yield across lactations. Field data, collected through the Dairy Herd Improvement Association, on US Jersey cows first calving between January 1997 and November 2004 were used. Characterization of DD included a frequency distribution of dry period lengths as well as factors affecting US Jersey DD. Of the factors considered in this research, the primary ones affecting dry period length were DO, milk yield, and SCS. Cows with longer DO, lower milk yield, and higher SCS received longer dry periods. The model for analyses included herd-year of calving, year-state-month of calving, parity of calving, previous lactation record, age at calving, and DD as a categorical variable; records were preadjusted for cow effects. A total of 123,032 records from 73,797 cows in 808 herds were used for estimation of DD effects on subsequent lactation actual milk yield. Jersey milk, fat, and protein yields in the subsequent lactation were maximized with 61 to 65 DD. Dry periods of 30 d or fewer resulted in large reductions in subsequent lactation production. A short dry period was beneficial for fat and protein percentages in the subsequent lactation. Short dry periods also resulted in fewer DO in the subsequent lactation; however, this was entirely due to the lower milk yield associated with shortened dry periods. The biggest difference between Jerseys and Holsteins was a much larger detrimental effect on SCS in Jerseys for dry periods of 30 d or less. Jersey SCS increased 10%, relative to the overall mean, for dry periods of 20 d or less and 4.6% for DD between 21 and 30 d. Dry periods of 45 to 70 d maximized yields across adjacent lactations. A dry period length, after first lactation, of 45 to 70 d also maximized actual milk yield across lactations 1, 2, and 3. The final recommendation to Jersey producers is to avoid dry periods of <45 d. Long dry periods (>70 d) should also be avoided because these are even more costly to total yield than dry periods <30 d.     

Use of repeated measures analysis for evaluation of genetic background of dairy cattle behavior in automatic milking systems

Milking frequencies measured at official test days were used with repeated measurement analysis to reveal the environmental and genetic impact on the milking frequency of cows in automatic milking systems. Repeated measurements were 3 test-day observations per cow within days in milk (DIM) classes, with 1,216 cows in DIM class 1 (d 0 to 99), from 1,112 cows in DIM class 2 (d 100 to 199), and from 1,004 cows in DIM class 3 (d 200 to 299) kept in 15 farms. Selection criteria for models analyzing repeated measurements were Akaike and Schwarz Bayesian values, which favored the autoregressive [AR(1)] covariance structure over the compound symmetry model. Results from the AR(1) model indicated a significant impact of fixed herd and parity effects. Milking frequencies decreased with increasing parities and were greatest for first-parity cows. High daily milk yield was associated with higher milking frequencies. Heritabilities for milking frequency were 0.16, 0.19, and 0.22 in DIM classes 1, 2, and 3, respectively, from the AR(1) model. Higher heritabilities in the later stage of lactation were due to a substantial reduction of the residual variance. Genetic correlations between test-day milk yield and daily milking frequency were in the range of 0.46 to 0.57 for all DIM classes and between milking frequency and somatic cell score were near zero. For verification of results, milking frequencies of the same cows obtained from herd management programs were averaged within DIM classes. Heritabilities were slightly above the values from the AR(1) model. In conclusion, heritabilities for milking frequency in automatic milking systems are moderate enough to incorporate this behavioral trait in a combined breeding goal. The inevitable improvement of labor efficiency in dairy cattle farming demands such cows going easily and voluntarily in automatic milking systems.     

Piecewise modeling of the associations between dry period length and milk,fat, and protein yield changes in the subsequent lactation

Our objective was to develop predictive models of 305-d mature-equivalent milk, fat, and protein yields in the subsequent lactation as continuous functions of the number of days dry (DD) in the current lactation. In this retrospective cohort study with field data, we obtained DHIA milk recording lactation records with the last DD in 2014 or 2015. Cows included had DD from 21 to 100 d. After editing, 1,030,141 records from cows in 7,044 herds remained. Three parity groups of adjacent (current, subsequent) lactations were constructed. We conducted all analyses by parity group and yield component. We first applied control models to pre-adjust the yields in the subsequent lactation for potentially confounding effects. Control models included the covariates mature-equivalent yield, days open, somatic cell score at 180 d pregnant, daily yield at 180 d pregnant, and a herd-season random effect, all observed in the current lactation. Days dry was not included. Second, we modeled residuals from control models with smooth piecewise regression models consisting of a simple linear, quadratic, and another simple linear equation depending on DD. Yield deviations were calculated as differences from predicted mature-equivalent yield at 50 DD. For validation, predictions of yield deviations from piecewise models by DD were compared with predictions from local regression for the DHIA field records and yield deviations reported in 38 experimental and field studies found in the literature. Control models reduced the average root mean squared prediction error by approximately 21%. Yield deviations were increasingly more negative for DD shorter than 50 d, indicating lower yields in the subsequent lactation. For short DD, the decrease in 305-d mature-equivalent milk yield ranged from 43 to 53 kg per DD. For mature-equivalent fat and protein yields, decreases were between 1.28 and 1.71 kg per DD, and 1.06 and 1.50 kg per DD, respectively. Yield deviations often were marginally positive and increasing for DD >50, so that the highest yield in the subsequent lactation was predicted for 100 DD. For long DD, the 305-d mature-equivalent milk yield increased at most 4.18 kg per DD. Patterns in deviations for fat and protein yield were similar to those for milk yield deviations. Predictions from piecewise models and local regressions were very similar, which supports the chosen functional form of the piecewise models. Yield deviations from field studies in the literature typically were decreasing when DD were longer, likely because of insufficient control for confounding effects. In conclusion, piecewise models of mature-equivalent milk, fat, and protein yield deviations as continuous functions of DD fit the observed data well and may be useful for decision support on the optimal dry period length for individual cows.     

Genetic parameters of milking frequency and milk production traits in Canadian Holsteins milked by an automated milking system

Twice-a-day milking is currently the most frequently used milking schedule in Canadian dairy cattle. However, with an automated milking system (AMS), dairy cows can be milked more frequently. The objective of this study was to estimate genetic parameters for milking frequency and for production traits of cows milked within an AMS. Data were 141,927 daily records of 953 primiparous Holstein cows from 14 farms in Ontario and Quebec. Most cows visited the AMS 2 (46%) or 3 (37%) times a day. A 2-trait [daily (24-h) milking frequency and daily (24-h) milk yield] random regression daily animal model and a multiple-trait (milk, fat, protein yields, somatic cell score, and milking frequency) random regression test-day animal model were used for the estimation of (co)variance components. Both models included fixed effect of herd × test-date, fixed regressions on days in milk (DIM) nested within age at calving by season of calving, and random regressions for additive genetic and permanent environmental effects. Both fixed and random regressions were fitted with fourth-order Legendre polynomials on DIM. The number of cows in the multiple-trait test-day model was smaller compared with the daily animal model. Heritabilities from the daily model for daily (24-h) milking frequency and daily (24-h) milk yield ranged between 0.02 and 0.08 and 0.14 and 0.20, respectively. Genetic correlations between daily (24-h) milk yield and daily (24-h) milking frequency were largest at the end of lactation (0.80) and smallest in mid-lactation (0.27). Heritabilities from the test-day model for test-day milking frequency, milk, fat and protein yield, and somatic cell score were 0.14, 0.26, 0.20, 0.21, and 0.20, respectively. The genetic correlation was positive between test-day milking frequency and official test-day milk, fat, and protein yields, and negative between official test-day somatic cell score and test-day milking frequency.     

Blood calcium concentration and performance in periparturient and early lactating dairy cows is influenced by plant bioactive lipid compounds

Previous studies ex vivo suggested that plant bioactive lipid compounds (PBLC) can increase ruminal calcium absorption. Therefore, we hypothesized that PBLC feeding around calving may potentially counteract hypocalcemia and support performance in postpartum dairy cows. The corresponding aim of the study was to investigate the effect of PBLC feeding on blood minerals in Brown Swiss (BS) and hypocalcemia-susceptible Holstein Friesian (HF) cows during the period from d −2 to 28 relative to calving and on milk performance until d 80 of lactation. A total of 29 BS cows and 41 HF cows were divided each into a control (CON) and PBLC treatment group. The latter was supplemented with 1.7 g/d menthol-rich PBLC from 8 d before expected calving to 80 d postpartum. Milk yield and composition, body condition score and blood minerals were measured. Feeding PBLC induced a significant breed × treatment interaction for iCa, supporting that PBLC increased iCa exclusively in HF cows; the increase was 0.03 mM over the whole period and 0.05 mM from d 1 to 3 after calving. Subclinical hypocalcemia was seen in one BS-CON and 8 HF-CON cows and 2 BS-PBLC and 4 HF-PBLC cows. Clinical milk fever was detected only in HF cows (2 HF-CON and one HF-PBLC). Other tested blood minerals, such as sodium, chloride, and potassium, as well as blood glucose, were neither affected by PBLC feeding nor breed, nor were their 2-way interactions, except for higher sodium levels in PBLC cows on d 21. Body condition score showed no effect of treatment, except for a lower body condition score in BS-PBLC compared with BS-CON at d 14. Dietary PBLC increased milk yield, milk fat yield, and milk protein yield at 2 consecutive dairy herd improvement test days. As indicated by treatment × day interactions, energy-corrected milk yield and milk lactose yield were increased by PBLC on the first test day only, and milk protein concentration decreased from test d 1 to test d 2 in CON only. The concentrations of fat, lactose, and urea, as well as somatic cell count, were not affected by treatment. The weekly milk yield over the first 11 wk of lactation was 29.5 kg/wk higher for PBLC versus CON across breeds. It is concluded that the applied PBLC induced a small but measurable improvement of calcium status in HF cows in the study period and had additional positive effects on milk performance in both breeds.     

Genetic correlations of mid-infrared-predicted milk fatty acid groups with milk production traits

The objective of this research was to estimate the genetic correlations between milk mid-infrared-predicted fatty acid groups and production traits in first-parity Canadian Holsteins. Contents of short-chain, medium-chain, long-chain, saturated, and unsaturated fatty acid groupings in milk samples can be predicted using mid-infrared spectral data for cows enrolled in milk recording programs. Predicted fatty acid group contents were obtained for 49,127 test-day milk samples from 10,029 first-parity Holstein cows in 810 herds. Milk yield, fat and protein yield, fat and protein percentage, fat-to-protein ratio, and somatic cell score were also available for these test days. Genetic parameters were estimated for the fatty acid groups and production traits using multiple-trait random regression test day models by Bayesian methods via Gibbs sampling. Three separate 8- or 9-trait analyses were performed, including the 5 fatty acid groups with different combinations of the production traits. Posterior standard deviations ranged from <0.001 to 0.01. Average daily genetic correlations were negative and similar to each other for the fatty acid groups with milk yield (?0.62 to ?0.59) and with protein yield (?0.32 to ?0.25). Weak and positive average daily genetic correlations were found between somatic cell score and the fatty acid groups (from 0.25 to 0.36). Stronger genetic correlations with fat yield, fat and protein percentage, and fat-to-protein ratio were found with medium-chain and saturated fatty acid groups compared with those with long-chain and unsaturated fatty acid groups. Genetic correlations were very strong between the fatty acid groups and fat percentage, ranging between 0.88 for unsaturated and 0.99 for saturated fatty acids. Daily genetic correlations from 5 to 305 d in milk with milk, protein yield and percentage, and somatic cell score traits showed similar patterns for all fatty acid groups. The daily genetic correlations with fat yield at the beginning of lactation were decreasing for long-chain and unsaturated fatty acid groups and increasing for short-chain fatty acids. Genetic correlations between fat percentage and fatty acids were increasing at the beginning of lactation for short- and medium-chain and saturated fatty acids, but slightly decreasing for long-chain and unsaturated fatty acid groups. These results can be used in defining fatty acid traits and breeding objectives.