Quantifying cost of disease in livestock: a new metric for the Global Burden of Animal Diseases

BACKGROUND

Increasing awareness of the environmental and public health impacts of expanding and intensifying animal-based food and farming systems creates discord, with the reliance of much of the world's population on animals for livelihoods and essential nutrition. Increasing the efficiency of food production through improved animal health has been identified as a step towards minimising these negative effects without compromising global food security. The Global Burden of Animal Diseases (GBADs) programme aims to provide data and analytical methods to support positive change in animal health across all livestock and aquaculture animal populations.

METHODS

In this study, we present a metric that begins the process of disease burden estimation by converting the physical consequences of disease on animal performance to farm-level costs of disease, and calculates a metric termed the Animal Health Loss Envelope (AHLE) via comparison between the status quo and a disease-free ideal. An example calculation of the AHLE metric for meat production from broiler chickens is provided.

FINDINGS

The AHLE presents the direct financial costs of disease at farm-level for all causes by estimating losses and expenditure in a given farming system. The general specification of the model measures productivity change at farm-level and provides an upper bound on productivity change in the absence of disease. On its own, it gives an indication of the scale of total disease cost at farm-level.

INTERPRETATION

The AHLE is an essential stepping stone within the GBADs programme because it connects the physical performance of animals in farming systems under different environmental and management conditions and different health states to farm economics. Moving forward, AHLE results will be an important step in calculating the wider monetary consequences of changes in animal health as part of the GBADs programme.

FUNDING

Bill & Melinda Gates Foundation, the UK Foreign, Commonwealth and Development Office, EU Horizon 2020 Research and Innovation Programme.

Association between days post conception and lactation persistency in dairy cattle

Determining the optimal insemination moment for individual cows is complex, particularly when considering the impact of pregnancy on milk production. The effect of pregnancy on the absolute milk yield has already been reported in several studies. Currently, there is limited quantitative knowledge about the association between days post conception (DPC) and lactation persistency, based on a lactation curve model, and, specifically, how persistency changes during pregnancy and relates to the days in milk at conception (DIMc). Understanding this association might provide valuable insights to determine the optimal insemination moment. This study, therefore, aimed to investigate the association between DPC and lactation persistency, with an additional focus on the influence of DIMc. Available milk production data from 2005 to 2022 were available for 23,908 cows from 87 herds located throughout the Netherlands and Belgium. Persistency was measured by a lactation curve characteristic decay, representing the time taken to halve milk production after peak yield. Decay was calculated for 8 DPC (0, 30, 60, 90, 120, 150, 180 and 210 d after DIMc) and served as the dependent variable. Independent variables included DPC, DIMc (< = 60, 61-90, 91-120, 121-150, 151-180, 181-210, > 210), parity group, DPC × parity group, DPC × DIMc and variables from 30 d before DIMc as covariates. The results showed an increase in decay, i.e., a decrease in persistency, during pregnancy for both parity groups, albeit in different ways. Specifically, from DPC 150 to DPC 210, multiparous cows showed a higher decline in persistency compared with primiparous cows. Furthermore, a later DIMc (cows conceiving later) was associated with higher persistency. Except for the early DIMc groups (DIMc < 90), DIMc does not impact the change in persistency by gestation. The findings from this study contribute to a better understanding of how DPC and DIMc during lactation influence lactation persistency, enabling more informed decision-making by farmers who wish to take persistency into account in their reproduction management.

The association between time-series milk β-hydroxybutyrate dynamics and early reproductive performance of dairy cows

The DeLaval Herd Navigator is an on-farm sensor system that measures on a frequent basis milk progesterone (P4) and β-hydroxybutyrate (BHB) in individual cows to closely monitor reproductive performance and energy balance. This information provides the opportunity to investigate the dynamics of BHB measured in milk (mBHB) and study the association between mBHB and early reproductive performance. The objectives of the study were (1) to describe mBHB dynamics within the first 20 d in milk (DIM), and (2) to evaluate the association between mBHB dynamics and early reproductive performance at cow-level. Two-year time-series data from 4,133 dairy cows in 38 Dutch dairy farms were available for analysis. Data included information on mBHB, daily milk yield and the indicators of early reproductive performance, days from calving to resumption of cyclicity, days from calving to first estrus, and days from calving to first insemination. The following mBHB dynamic parameters were defined based on the first 20 DIM for each individual cow: average mBHB (AvgBHB), DIM when mBHB was for the first time ≥80 μmol/L (OnsetKeto), the total number of consecutive days a cow had mBHB concentration ≥80 μmol/L, and the number of measurements mBHB concentration was ≥80 μmol/L. Three Cox proportional hazard regression models with random herd effect were developed to evaluate the association between cow level mBHB dynamics and days from calving to resumption of cyclicity, first estrus, and first insemination. Results showed that the mean AvgBHB within 20 DIM among all cows was 73 μmol/L. The mean OnsetKeto within 20 DIM, was 8 DIM. Among all cows having hyperketolactia, 55.8% (1,350/2,419) had OnsetKeto in the first week of lactation. In total, 41.5% (1,714/4,133) of the cows did not have OnsetKeto in the first 20 DIM. An early onset of hyperketolactia was associated with delayed fertility events. Cows with higher AvgBHB have a prolonged time interval from calving to resumption of cyclicity and first estrus. Information on mBHB dynamics and the association with early reproductive performance provides insights that might be helpful to improve reproductive performance of individual dairy cows.

Prediction of persistency for day 305 of lactation at the moment of the insemination decision

When deciding on the voluntary waiting period of an individual cow, it might be useful to have insight into the persistency for the remainder of that lactation at the moment of the insemination decision, especially for farmers who consider persistency in their reproduction management. Currently, breeding values for persistency are calculated for dairy cows but, to our knowledge, prediction models to accurately predict persistency at different moments of insemination are lacking. This study aimed to predict lactation persistency for DIM 305 at different insemination moments (DIM 50, 75, 100, and 125). Available cow and herd level data from 2005 to 2022 were collected for a total of 20,508 cows from 85 herds located in the Netherlands and Belgium. Lactation curve characteristics were estimated for every daily record using the data up to and including that day. Persistency was defined as the number of days it takes for the milk production to decrease by half during the declining stage of lactation, and calculated from the estimated lactation curve characteristic 'decay'. Four linear regression models for each of the selected insemination moment were built separately to predict decay at DIM 305 (decay-305). Independent variables included the lactation curve characteristics at the selected insemination moment, daily milk yield, age, calving season, parity group and other herd variables. The average decay-305 of primiparous cows was lower than that of multiparous cows (1.55 *10-3 vs. 2.41*10-3, equivalent to a persistency of 447 vs. 288 days, respectively). Results showed that our models had limitations in accurately predicting persistency, although predictions improved slightly at later insemination moments, with R2 values ranging between 0.27 and 0.41. It can thus be concluded that, based only on cow and herd milk production information, accurate prediction of persistency for DIM 305 is not feasible.

Milk β-hydroxybutyrate metrics and its consequences for surveillance of hyperketonaemia on commercial dairy farms

Dairy cows that are unable to adapt to a change in their metabolic status are at risk for hyperketonaemia (HK). Reported HK herd level prevalences range a lot and we hypothesized that this is partly due to differences in used tests and monitoring protocols. Insights in milk β-hydroxybutyrate (BHB) metrics can potentially explain why the reported incidences or prevalences vary between test strategies. Automated collection and repeated analyses of individual milk samples with the DeLaval Herd Navigator™ (HN) provides real-time data on milk BHB concentrations. We aimed to use that information to gain insight in BHB metrics measured in milk from 3 to 60 days in milk (DIM). Using different cut-offs (0.08, 0.10 and 0.15 mmol/L), 5 BHB metrics were determined. Furthermore, the impact of 4 arbitrary test protocols on the detected incidence of HK was assessed. We used HN data of 3,133 cows from 35 herds. The cumulative incidence of HK between 3 and 60 DIM varied between 30.5 and 76.7% for different cut-off values. We found a higher HK incidence for higher parity cows. The first elevated BHB concentrations were roughly found between one and two weeks after calving. For higher parity cows the maximum BHB concentrations were higher, the onset of HK was earlier after calving, and the number of episodes of HK was higher. It appeared that the sensitivity of a HK test protocol can be increased by increasing the testing frequency from once to twice a week. Also extending the number of days of the test window from 4-14 to 4-21 days enhances the chance to find cows experiencing HK. In conclusion, HN data provided useful insights in milk BHB metrics. The chosen cut-off value had a large effect on the reported metrics which explains why earlier reported incidences or prevalences vary such a lot. Differences in test period and sample selection also had a large impact on the observed HK incidence. We suggest to take this in consideration while evaluating whether HK is an issue on farm level and use a uniform protocol for benchmarking of HK between farms.

Stochastic model to assess bioeconomic impact of PRRS on pig farms in Costa Rica

Despite the economic importance of PRRS and its high prevalence in Costa Rica, there are no studies on the bioeconomic impact of the disease in the country or, even, in Central America. Such studies are essential in finding cost-effective preventive measures tailored for different production circumstances. Therefore, the objective of this study was to evaluate economic and production parameters of a PRRSV-infection for a medium-sized farrow-to-finish pig farm system in Costa Rica with a farm-level stochastic Monte Carlo simulation model. The effect of PRRS was assessed by scenario analysis, in which a baseline PRRS-free situation was compared against three alternative scenarios that assumed low, medium and high PRRS effects. The PRRS effects were based on data from local farms, scientific literature and expert opinion. Sensitivity analyses were performed to assess the impact of key input parameters on output variables. Results show that at the animal level, changes between the baseline and the PRRS-high scenario were estimated as: + 25 d in age to slaughter, - 9.9 pigs to slaughter (per breeding sow/yr), + 6% annual replacement rate, - 255 d in sow productive lifetime, - 6.9 mo in age at culling of sows, and + 24 non- productive days. For a medium size local farm (n = 588 sows), a reduction of 5826 fat pigs to slaughter per farm/yr from baseline compared to PRRS-high scenario was observed. PRRS-induced loss per farm per year was estimated at -US $142,542, US $180,109 and -US $524,719 for PRRS-low, medium and high scenarios, respectively. Revenues/costs ratio changed from 1.12 in the baseline to 0.89 in the PRRS-high scenario. The production cost per kg carcass weight increased from US $2.63 for the baseline to US $3.35 in the PRRS-high scenario. PRRS-induced loss was estimated at US $77.1 per slaughtered pig/yr and US $892 per breeding sow/yr for the PRRS-high scenario. Results from the model indicate that pig farms with medium to high prevalence of PRRS will require optimal market conditions in order to have positive economic outcomes. These results can be helpful in the design of better control strategies for PRRS.

The economic impact of endemic respiratory disease in pigs and related interventions - a systematic review

BACKGROUND

Understanding the financial consequences of endemically prevalent pathogens within the porcine respiratory disease complex (PRDC) and the effects of interventions assists decision-making regarding disease prevention and control. The aim of this systematic review was to identify what economic studies have been carried out on infectious endemic respiratory disease in pigs, what methods are being used, and, when feasible, to identify the economic impacts of PRDC pathogens and the costs and benefits of interventions.

RESULTS

By following the PRISMA method, a total of 58 studies were deemed eligible for the purpose of this systematic review. Twenty-six studies used data derived from European countries, 18 from the US, 6 from Asia, 4 from Oceania, and 4 from other countries, i.e., Canada, Mexico, and Brazil. Main findings from selected publications were: (1) The studies mainly considered endemic scenarios on commercial fattening farms; (2) The porcine reproductive and respiratory syndrome virus was by far the most studied pathogen, followed by Mycoplasma hyopneumoniae, but the absence or presence of other endemic respiratory pathogens was often not verified or accounted for; (3) Most studies calculated the economic impact using primary production data, whereas twelve studies modelled the impact using secondary data only; (4) Seven different economic methods were applied across studies; (5) A large variation exists in the cost and revenue components considered in calculations, with feed costs and reduced carcass value included the most often; (6) The reported median economic impact of one or several co-existing respiratory pathogen(s) ranged from €1.70 to €8.90 per nursery pig, €2.30 to €15.35 per fattening pig, and €100 to €323 per sow per year; and (7) Vaccination was the most studied intervention, and the outcomes of all but three intervention-focused studies were neutral or positive.

CONCLUSION

The outcomes and discussion from this systematic review provide insight into the studies, their methods, the advantages and limitations of the existing research, and the reported impacts from the endemic respiratory disease complex for pig production systems worldwide. Future research should improve the consistency and comparability of economic assessments by ensuring the inclusion of high impact cost and revenue components and expressing results similarly.

Dutch dairy farmers' perspectives on culling reasons and strategies

Since the abolishment of the milk quota system in Europe in 2014 and the introduction of environmental policies such as the phosphate rights system in the Netherlands, the reasons for culling dairy cows might have changed. The aim of this study was to determine the culling reasons for dairy cattle and to identify farmers' culling strategies and their intentions regarding the alteration of indicated culling strategies. To this end, an online questionnaire was distributed among dairy farmers nationally that resulted in 207 responses. Results showed that the most frequent culling reasons were related to problems with reproduction, udder, and hoof health. Primiparous cows were primarily culled for miscellaneous reasons such as injury, reproduction failure, and low milk yield. Multiparous cows were culled predominantly for reproduction failure, udder health and hoof health reasons. Most respondents indicated that they consider formulating a culling strategy, based on certain rules of thumb regarding the most common reasons for culling. Most farmers also reported that culling decisions on their farms were perceived to be unavoidable, though reproductive culling decisions are primarily voluntary. Most respondents stated that they intended to reduce the culling rate for better economic gain did not intend to alter the amount of replacement stock reared. The applied rules of thumb regarding culling strategies do not seem to have changed since the policy changes in dairy farming. The question remains whether farmers' rules of thumb might have made them unaware of the actual economic consequences of their culling strategies under the altered situation.

The DECIDE project: from surveillance data to decision-support for farmers and veterinarians

Farmers, veterinarians and other animal health managers in the livestock sector are currently missing sufficient information on prevalence and burden of contagious endemic animal diseases. They need adequate tools for risk assessment and prioritization of control measures for these diseases. The DECIDE project develops data-driven decision-support tools, which present (i) robust and early signals of disease emergence and options for diagnostic confirmation; and (ii) options for controlling the disease along with their implications in terms of disease spread, economic burden and animal welfare. DECIDE focuses on respiratory and gastro-intestinal syndromes in the three most important terrestrial livestock species (pigs, poultry, cattle) and on reduced growth and mortality in two of the most important aquaculture species (salmon and trout). For each of these, we (i) identify the stakeholder needs; (ii) determine the burden of disease and costs of control measures; (iii) develop data sharing frameworks based on federated data access and meta-information sharing; (iv) build multivariate and multi-level models for creating early warning systems; and (v) rank interventions based on multiple criteria. Together, all of this forms decision-support tools to be integrated in existing farm management systems wherever possible and to be evaluated in several pilot implementations in farms across Europe. The results of DECIDE lead to improved use of surveillance data and evidence-based decisions on disease control. Improved disease control is essential for a sustainable food chain in Europe with increased animal health and welfare and that protects human health.

Hormone use for reproductive diseases and heat induction in relation to herd-level reproductive performance in Dutch dairy farms

This ecological study aimed to associate hormone use for reproductive diseases and heat induction with reproductive performance at herd level. Hormone use, herd characteristics, and test-day recording data were obtained from 754 representative Dutch dairy farms belonging to five large veterinary practices from 2017 to 2019 (1679 observations in total). Hormone use was classified into prostaglandin, gonadotropin-releasing hormone (GnRH), and progesterone, and was expressed at herd level as the annual number of hormone doses per 100 adult dairy cows. Hormone use was categorized into four levels (no usage, low, medium, and high use), following the 33rd and 66th percentiles of herds that applied them. Three herd-level reproductive performance indicators (calving interval, calving-to-1st insemination interval, number of inseminations per cow) were analyzed using multivariable General Estimating Equations models. The median annual total hormone use was 36.1 (mean=43.1; min=0.0; max=248.2) doses per 100 adult dairy cows in all herds while the median was 39.2 (mean=46.8; min=0.4; max=248.2) doses per 100 adult dairy cows among the user-herds. The median annual group-specific hormone use was 21.3 (mean=26.1; min 0.0; max=180.0), 11.0 (mean=15.3; min=0.0; max=127.0) and 0.0 (mean=1.8; min=0.0; max=40.3) doses per 100 adult dairy cows for prostaglandin, GnRH, and progesterone, respectively. The final statistical models identified that herds with a high hormone use had a calving interval and a calving-to-1st insemination interval that was 9.3 ± 2.6 and 16.4 ± 2.1 days shorter than that of non-user herds (424.0 ± 2.7 and 114.0 ± 2.1 days), respectively. Furthermore, high-user herds needed on average 0.3 ± 0.04 inseminations more to get their cows pregnant compared to non-user herds (1.83 ± 0.04 no. of inseminations per cow). Medium-user herds had a 6.5 ± 2.6 days shorter calving interval and a 12.0 ± 2.1 days shorter calving-to-1st insemination interval with 0.2 ± 0.04 additional inseminations per cow compared to non-user herds. Low-user herds had a 6.2 ± 2.7 days shorter calving interval and a 7.9 ± 2.2 days shorter calving-to-1st insemination interval compared to non-user herds. The model produced the same trend for prostaglandin and GnRH use, with the higher use being associated with a shorter calving interval, a shorter calving-to-1st insemination interval, and a higher insemination per cow number. For progesterone use the opposite effect was observed. In conclusion, using a large representative herd-level dataset, hormone use was associated with a better reproductive performance in terms of calving interval and calving-to-1st insemination interval but gave extra average number of inseminations per cow. It should be monitored how reproduction performance changes when striving for a more prudent hormone use.

Associations between dairy farm performance indicators and culling rates under policy-driven herd size constraints

Introduction

This article aimed to study cross-sectional associations between the performance of dairy farms and their corresponding culling proportions under the herd size constraint as imposed in 2018 by the new phosphate regulation in the Netherlands.

Methods

To this end, production data from 10,540 Dutch dairy farms were analyzed to capture the inflow and outflow of both primiparous and multiparous cows. Farm performance was measured by 10 indicators structured in four areas of longevity, production, reproduction, and udder health. Farm culling proportions were represented by the overall culling (OC) and the number of culled primiparous cows in relation to (i) the total number of producing cows (PC), (ii) the number of producing primiparous cows (PPC), and (iii) the number of culled producing cows (POC). Spearman's rank correlation and weighted logistic regression were adopted to study associations.

Results

In 2018, on average, 28% of producing cows were culled (OC). The number of primiparous cows culled represented 4.5% of the total number of producing cows (PC) and the mean proportion of culled primiparous cows was 18.8% of the total number of producing primiparous cows (PPC), and, of the total number of producing culled cows, 15% were primiparous cows (POC). However, the variance around the mean, and among individual farms, was high (SD 4-15% for all four culling proportions). Results from rank correlation showed very low-rank conformity (<12%) between the areas of production, reproduction, and udder health to the culling proportions. Results from logistic regression showed that higher farm levels of production and higher percentages of cows with poor udder health were associated with more overall culling but with less primiparous culling. For reproduction indicators, the associations were similar for overall and primiparous culling. However, except for the average age of culled animals, the odds ratios for indicators were close to 1 (range: 0.92-1.07 and 0.68-1.07 for OC and PPC, respectively), indicating only weak associations to culling proportions.

Discussion

In conclusion, although the introduction of phosphate regulation resulted in an increased outflow of cattle, corresponding culling proportions were not associated with the level of farm performance measured in terms of production, reproduction, or udder health.

The association of herd performance indicators with dairy cow longevity: An empirical study

The associations between reproductive performance, milk yield and health status with the risk of culling, and thus with a cow's longevity, have been well documented at the individual cow level. Associations at individual cow level may, however, not be valid at herd level due to interrelated herd management aspects and/or policy restrictions. The objective of this study was to explore the association of herd performance indicators with herd-level dairy cow longevity under Dutch production conditions. Longevity was expressed by three different measures, viz. age at culling, lifetime milk production of culled cows and culling rate. The evaluated herd performance indicators included factors on milk production, youngstock rearing, reproduction and health performance as registered on 10 719 Dutch commercial dairy herds during the period 2007-2016. Averaged over herds and the evaluated period, the age of culled milking cows was 2 139 days (5.8 years, SD±298 days), the lifetime milk production of culled cows was 31 238 kg (SD±7,494 kg), and the culling rate was 0.24 (SD±0.08). A mixed linear regression modelling approach was applied to evaluate the association of each of the three longevity measures with the selected herd performance indicators. The results indicated that only four herd performance indictors (herd size, herd expansion, heifer ratio and the proportion of cows with potential subclinical ketosis) shared significant associations with all three longevity variables. Generally, the strength of the associations between each of the evaluated longevity measures and herd performance indicators was only limited. The absence of strong associations between the longevity measures and herd performance indicators reveal that there is potential of extending cattle longevity without affecting the herd performance in terms of milk production, reproduction and health. Moreover, only part of the observed variance in longevity among the herds over time was explained by the herd performance variables, indicating that differences in longevity at herd level may predominantly be determined by other factors, like farmers' attitude and strategic management.

Linear Mixed-Effects Model to Quantify the Association between Somatic Cell Count and Milk Production in Italian Dairy Herds

Milk production loss due to mastitis in dairy herds is economically important. Before estimating the economic impacts of mastitis, it is crucial to quantify the association between mastitis and milk production. The objective of this study was to estimate the association between somatic cell count (SCC, as an indicator of intramammary infection due to mastitis) and milk production for dairy cows in Lombardy, Italy. The test-day (TD) records data of 3816 dairy herds located in three different geographical areas of Lombardy from January 2016 to December 2018 were used. After data editing, the final dataset comprised 10,445,464 TD records from 2970 farms and 826,831 cows. The analysis was carried out by using a mixed-effects model with six fixed effects (geographical Area, Breed, Days in Milk, Parity, Season and Year) and nested random effects for each cow and herd. The results confirmed that the SCC had a negative association with milk production. On average, this study found that any two-fold increase of SCC resulted in a milk production loss of 0.830 (95% CI: -0.832, -0.828) kg/cow/day in the whole of Lombardy. These results can be used for economic calculations on the costs of mastitis.

Herd level economic comparison between the shape of the lactation curve and 305 d milk production

Herd milk production performance is generally evaluated using the herd's average 305-day milk production (HM305). Economic comparisons between herds are also often made using HM305. Comparing herds is thus based on summarized milk production, and not on the form of the lactation curves of the cows within the herd. Cow lactation curve characteristics can be aggregated on a calendar year basis to herd lactation curve characteristics (HLCC) (herd magnitude, herd time to peak yield and herd persistency). Thus far, no literature has evaluated whether the shape of the lactation curve (described by HLCC) is better able to explain the economic variation of herds than summarized milk production such as HM305 does. This study aims to determine whether HM305 or HLCC is better able to explain the variation in economic performance between herds. To do so, we evaluated 8 years of Dutch longitudinal data on milk production and the financial accounts of 1,664 herds. Cow lactation curve characteristics were calculated through lactation curve modeling and aggregated to HLCC on a calendar year basis for two parity groups (primiparous cows and multiparous cows). Using income over feed cost per cow (IOFC-cow) or per 100 kg milk (IOFC-milk) as the dependent variable separately, we developed four linear mixed models. Two models were used to analyse the association between herd economic performance and HLCC; the other two models were used to analyse the association between herd economic performance and HM305. A Cox test and J test were used to compare two non-nested models to investigate whether HM305 or HLCC better explain IOFC. The average IOFC-cow was €2,305 (SD = 408) per year, while the average IOFC-milk was €32.1 (SD = 4.6). Results showed that HLCC and HM305 explain the same amount of variance of IOFC-cow or IOFC-milk. IOFC-cow was associated with HM305 and HLCC (except herd time to peak yield for primiparous cows). Herd magnitude was most strongly associated with IOFC-cow, followed by herd persistency and herd time to peak yield of multiparous cows. IOFC-milk was not associated with HM305 or HLCC (except for a weak negative association with herd persistency for primiparous cows). IOFC-cow and IOFC-milk were driven most by time effects. In conclusion, HLCC and HM305 explain the same amount of variance in IOFC-cow or IOFC-milk. HLCC is more computationally expensive, while HM305 is more readily available.

Estimating the Effect of a Bovine Viral Diarrhea Virus Control Program: An Empirical Study on the Performance of Dutch Dairy Herds

More and more European countries have implemented a bovine viral diarrhea virus (BVDV) control program. The economic effects of such programs have been evaluated in simulations, but empirical studies are lacking, especially in the final stage of the program. We investigated the economic (gross margin) and production effects (milk yield, somatic cell count, and calving interval) of the herds obtaining BVDV-free certification based on longitudinal annual accounting and herd performance data from Dutch dairy herds between 2014 and 2019, the final stages of the Dutch national BVDV-free program. This study was designed as a case-control study: two types of case herds were defined for two analyses. The case herds in the first analysis are herds where the BVDV status changed from “BVDV not free” to “BVDV free” during the study period. The not-free status refers to a herd that participated in the BVDV-free program but had not yet obtained the BVDV-free certification. In the second analysis, the case herds started participating in the Dutch BVDV-free program during the study period and obtained the BVDV-free certification. Control herds in both analyses were BVDV-free during the entire study period. Potential bias between the covariates of the two herd groups was reduced by matching case and control herds using the propensity score matching method. To compare the differences between case and control herds before and after BVDV-free certification, we used the time-varying Difference-in-Differences estimation (DID) methodology. The results indicate that there was no significant change in milk yield, somatic cell count, calving interval, and gross margin upon BVDV-free certification. There are several possible explanations for the non-significant effects observed in our study, such as the final stage of the BVDV control program, not knowing the true BVDV infection situation in case herds and not knowing if control measures were implemented in case herds prior to participating in the BVDV-free program. In our study, the effects of BVDV-free certification might have been underestimated, given that the Dutch BVDV control program became mandatory during the study period, and some of the case herds might have never experienced any BVDV infection. The results of this study suggest that in the final stage of the BVDV control program, the program may no longer have a clear benefit to the herd performance of participating dairy herds. When designing national programs to eradicate BVDV, it is therefore important to include incentives for such farms to motivate them to join the program.

The effect of new bovine viral diarrhea virus introduction on somatic cell count, calving interval, culling, and calf mortality of dairy herds in the Dutch bovine viral diarrhea virus-free program

Bovine viral diarrhea virus (BVDV) infection has a major effect on the health of cows and consequently on herd performance. Many countries have implemented control or eradication programs to mitigate BVDV infection and its negative effects. These negative effects of BVDV infection on dairy herds are well documented, but there is much less information about the effects of new introduction of BVDV on dairy herds already participating in a BVDV control program. The objective of our study was to investigate the effect of a new BVDV introduction in BVDV-free herds participating in the Dutch BVDV-free program on herd performance. Longitudinal herd-level surveillance data were combined with herd information data to create 4 unique data sets, including a monthly test-day somatic cell count (SCC) data set, annual calving interval (CIV) and culling risk (CR) data sets, and a quarterly calf mortality rate (CMR) data set. Each database contained 2 types of herds: herds that remained BVDV free during the whole study period (defined as free herds), and herds that lost their BVDV-free status during the study period (defined as breakdown herds). The date of losing the BVDV-free status was defined as breakdown date. To compare breakdown herds with free herds, a random breakdown date was artificially generated for free herds by simple random sampling from the distribution of the breakdown month of the breakdown herds. The SCC and CIV before and after a new introduction of BVDV were compared through linear mixed-effects models with a Gaussian distribution, and the CR and CMR were modeled using a negative binomial distribution in generalized linear mixed-effects models. The explanatory variables for all models included herd type, BVDV status, year, and a random herd effect. Herd size was included as an explanatory variable in the SCC, CIV, and CMR model. Season was included as an explanatory variable in the SCC and CMR model. Results showed that free herds have lower SCC, CR, CMR, and shorter CIV than the breakdown herds. Within the breakdown herds, the new BVDV introduction affected the SCC and CMR. In the year after BVDV introduction, the SCC was higher than that in the year before BVDV introduction, with a factor of 1.011 [2.5th to 97.5th percentile (95% PCTL): 1.002, 1.020]. Compared with the year before BVDV breakdown, the CMR in the year of breakdown and the year after breakdown was higher, with factors of 1.170 (95% PCTL: 1.120; 1.218) and 1.096 (95% PCTL: 1.048; 1.153), respectively. This study reveals that a new introduction of BVDV had a negative but on average relatively small effect on herd performance in herds participating in a BVDV control program.

An Empirical Analysis on the Longevity of Dairy Cows in Relation to Economic Herd Performance

Several studies have stated the various effects of an increased dairy cow longevity on economic herd performance, but empirical studies are lacking. This study aimed to investigate the association between longevity of dairy cows and the economic performance of dairy herds based on longitudinal Dutch accounting data. Herd and farm accounting data (n = 855 herds) over the years 2007–2016 were analyzed. Herd data contained yearly averages on longevity features, herd size and several production variables. Longevity was defined as the age of cows at culling and by lifetime milk production of culled cows. Farm accounting data contained yearly averages on revenues, fixed and variable costs of the herds, by which gross margins were defined. Data was analyzed using generalized linear mixed modeling, with gross margin as dependent variable. The independent variables consisted of average age of culled cows, average lifetime production of culled cows, year, herd size, herd intensity (milk production per ha), herd expansion rate, soil type, milking system, successor availability, total full-time equivalent, heifer ratio (% of heifers per cow) and use of outsourced heifer rearing. Herd was included as a random effect to account for the heterogeneity among herds. Descriptive statistics showed that the average age of culled cows was 5.87 (STD = 0.78) years and the average lifetime milk production of culled cows was 31.87 (STD = 7.56) tons per cow with an average herd size of 89 cows (STD = 38.85). The average age of culled cows was stable over the 10 years (variation between 5.79 AND 5.90 years). The gross margin was on average €24.80/100 kg milk (STD = 4.67), with the lowest value in year 2009 and the highest value in year 2013. Gross margin was not significantly associated with age of culled cows and lifetime milk production of culled cows. Variance in longevity between herds was large (STD = 0.78 years) but herds with a higher longevity did not perform economically better nor worse than herds resulting in lower longevity. This indicates that, within current practice, there is potential for improving longevity in order to meet society's concerns on animal welfare and environmental pollution, without affecting the economic performance of the herd.

The effect of bovine viral diarrhea virus introduction on milk production of Dutch dairy herds

Dairy cows are negatively affected by the introduction of bovine viral diarrhea virus (BVDV), and consequently, produce less milk. Existing literature on potential milk production losses is based on relatively outdated data and hardly evaluates milk production loss in relation to a new BVDV infection in a surveillance system. This study determined the annual and quarterly loss in milk production of BVDV introduction in 3,126 dairy herds participating in the Dutch BVDV-free program between 2007 and 2017. Among these herds, 640 were "breakdown-herds" that obtained and subsequently lost their BVDV-free status during the study period, and 2,486 herds obtained and retained their BVDV-free status during the study period. Milk yields before and after BVDV introduction were compared through annual and quarterly linear mixed models. The fixed variables for both models included herd type (breakdown-herd or free-herd), bovine viral diarrhea status (on an annual and quarterly basis), year, season, and a random herd effect. The dependent variable was the average daily milk yield on the test day. To define the possible BVDV-introduction dates, 4 scenarios were developed. In the default scenario, the date of breakdown (i.e., loss of the BVDV-free status) was assumed as the BVDV-introduction date. For the other 3 scenarios, the BVDV-introduction dates were set at 4, 6, and 9 mo before the date of breakdown, based on the estimated birth date of a persistently infected calf. In the default scenario, the loss in milk yield due to BVDV introduction occurred mainly in the first year after breakdown, with a reduction in yield of 0.08 kg/cow per day compared with the last year before breakdown. For the other 3 scenarios, the greatest yield reduction occurred in the second year after BVDV introduction, with a loss of 0.09, 0.09, and 0.1 kg/cow per day, respectively. For the first 4 quarters after BVDV introduction in the default scenario, milk yield loss was 0.14, 0.09, 0.02, and 0.08 kg/cow per day, respectively. These quarterly results indicated that milk yield loss was greatest in the first quarter after BVDV introduction. Overall, BVDV introduction had a negative, but on average a relatively small, effect on milk yield for herds participating in the BVDV-free program. This study will enable dairy farmers and policymakers to have a clearer understanding of the quantitative milk production effect of BVDV on dairy farms in a control program.

The economic impact of drying off cows with a dry-off facilitator (cabergoline) compared with 2 methods of gradual cessation of lactation for European dairy farms

An abrupt method to dry off cows has disadvantages and is considered inappropriate for current dairy cows due to welfare issues and risks for intramammary infections (IMI). A gradual cessation of lactation (by feeding or milking frequency reduction) has been the generally recommended method for drying off cows to prevent these adverse effects. However, a new alternative to the gradual approach is to abruptly stop milking at the same time as using cabergoline (CAB), a prolactin inhibitor. The aim of the study was to compare the net costs of 3 different methods of drying off cows [gradual reduction in feed (referred to as gradual feeding), gradual reduction in milking frequency (referred to as gradual milking), and abrupt cessation of milking with CAB]. A stochastic Monte Carlo simulation model, at cow level, was developed to calculate the net costs of applying these methods. All inputs for the model were based on literature information, authors' expertise, and expert knowledge. The net costs were determined by only including costs and benefits, which varied between the 3 methods. The model simulated a cow from 7 d before the day of drying off until the end of the next lactation. The likelihood of whether a cow was leaking milk early in the dry period was determined. Subsequently, it was determined whether or not the cow will get an IMI during the dry period, where the probability of getting an IMI was higher for cows leaking milk than for cows not leaking milk. If the IMI was not cured during the dry period, the cow had an IMI at calving. Also, milk production and feed requirements were modeled, and labor for applying the drying off method was included. For all methods, the net costs were calculated as the sum of costs for feed during the gradual feed reduction period, costs for applying the gradual-milking method, and the IMI costs during the dry period and lactation, minus the milk revenues during the transition from lactation to the dry period. Under default conditions, the average net cost of abrupt cessation of milking with CAB was €49.6/cow. The data showed that 90% of the net costs ranged from -€13.7 to €307.8/cow. The average net costs for gradual feeding and gradual milking were €99.1 and €71.5/cow, respectively. In conclusion, abrupt cessation of milking with CAB saved €49.5 and €21.9/cow on average compared with gradual feeding and gradual milking, respectively. This difference was mainly due to more milk returns and lower labor and IMI costs during lactation.

Economic and epidemiological impact of different intervention strategies for subclinical and clinical mastitis

The objective of this study was to evaluate and compare different combinations of intervention strategies for contagious or opportunistic subclinical and clinical intramammary infections (IMI). We simulated two different Danish dairy cattle herds with ten different intervention strategies focusing on cow-specific treatment or culling, including three baseline strategies without subclinical interventions. In one herd, the main causative pathogen of IMI was Staphylococcus (S.) aureus. In the other herd, Streptococcus (St.) agalactiae was the main causative agent. For both herds, we investigated costs and effectiveness of all ten intervention strategies. Intervention strategies consisted of measures against clinical and subclinical IMI, with baselines given by purely clinical intervention strategies. Our results showed that strategies including subclinical interventions were more cost-effective than the respective baseline strategies. Increase in income and reduction of IMI cases came at the cost of increased antibiotic usage and an increased culling rate in relation to IMI. However, there were differences between the herds. In the St. agalactiae herd, the clinical intervention strategy did not seem to have a big impact on income and number of cases. However, intervention strategies which included cow-specific clinical interventions led to a higher income and lower number of cases in the S. aureus herd. The results show that intervention strategies including interventions against contagious or opportunistic clinical and subclinical IMI can be highly cost-effective, but should be herd-specific.

Cow characteristics and their association with udder health after different dry period lengths

Shortening or omitting the dry period (DP) in dairy cows is of interest because of potential beneficial effects on energy balance and metabolic health. Reported effects of a short or omitted dry period on udder health are ambiguous. This study aimed to evaluate the effect of no DP (0d), a short DP (30d), or a conventional DP (60 d) on the occurrence of intramammary infections (IMI) during the precalving period and on somatic cell counts (SCC), elevations of SCC (SCC≥200,000 cells/mL), and clinical mastitis in the subsequent lactation. The study also aimed to analyze which prepartum cow characteristics are associated with udder health after different DP lengths. Holstein-Friesian dairy cows (n=167) were randomly assigned to a DP length (0, 30, or 60 d). Cows with a 0-d DP had a greater occurrence of chronic IMI and a lower occurrence of cured IMI during the precalving period than cows with a 30-d or 60-d DP. Postpartum average SCC for lactation was greater in cows with a 0-d DP than in cows with a 30-d or 60-d DP. The number of cows with at least 1 elevation of SCC, the number of elevations of SCC per affected cow, the number of cows treated for clinical mastitis, and the number of cases of mastitis per affected cow did not differ among DP lengths. Cow characteristics related to postpartum average SCC for lactation were DP length, parity, and the following interactions: DP length with prepartum elevation of SCC, DP length with fat- and protein-corrected milk (FPCM) reduction between 150 and 67d prepartum, DP length with parity and with average SCC for lactation, and last FPCM before the conventional drying-off day with average SCC for lactation. Cows with prepartum parity 1 had a lower occurrence of at least 1 elevation of SCC in subsequent lactation compared with cows with parity >2. Last SCC before the conventional drying-off day was positively associated with occurrence of clinical mastitis in the subsequent lactation. In this study, DP length was not a risk factor for either elevation of SCC or occurrence of clinical mastitis in the subsequent lactation. The identified cow characteristics could be used in a decision support model to optimize DP length for individual cows.

The total cost of rearing a heifer on Dutch dairy farms: calculated versus perceived cost

Background

As farmers do not often keep a record of the expenditures for rearing, an economic tool that provides insight into the cost of rearing is useful. In the Netherlands, an economic tool (Jonkos) has been developed that can be used by farmers to obtain insight into the cost of rearing on their farm. The first objective of this study is to calculate the total cost of rearing young stock in Dutch dairy herds using Jonkos. The second objective is to compare the calculated total cost of rearing with the farmers’ own estimation of the cost of rearing (the perceived cost).

Findings

Information was available for 75 herds that reared their own young stock and who had used the Jonkos tool. The perceived cost of rearing young stock was only available for 36 herds. In the 75 herds, the average herd size was 100 dairy cows. The average calculated total cost of rearing a heifer was €1,790. The average perceived total cost of rearing a heifer (including labour and housing costs) was €1,030.

Conclusion

Most Dutch farmers in the study underestimated the total cost of rearing. The Jonkos economic tool has the advantage that herd-specific information can be entered as input values. The output of the tool can improve the awareness of farmers about the total costs of rearing. This awareness can lead to a higher priority of young stock rearing and consequently to an improved quality of young stock rearing.

Effect of sensor systems for cow management on milk production, somatic cell count, and reproduction

To improve management on dairy herds, sensor systems have been developed that can measure physiological, behavioral, and production indicators on individual cows. It is not known whether using sensor systems also improves measures of health and production in dairy herds. The objective of this study was to investigate the effect of using sensor systems on measures of health and production in dairy herds. Data of 414 Dutch dairy farms with (n=152) and without (n=262) sensor systems were available. For these herds, information on milk production per cow, days to first service, first calving age, and somatic cell count (SCC) was provided for the years 2003 to 2013. Moreover, year of investment in sensor systems was available. For every farm year, we determined whether that year was before or after the year of investment in sensor systems on farms with an automatic milking system (AMS) or a conventional milking system (CMS), or whether it was a year on a farm that never invested in sensor systems. Separate statistical analyses were performed to determine the effect of sensor systems for mastitis detection (color, SCC, electrical conductivity, and lactate dehydrogenase sensors), estrus detection for dairy cows, estrus detection for young stock, and other sensor systems (weighing platform, rumination time sensor, fat and protein sensor, temperature sensor, milk temperature sensor, urea sensor, β-hydroxybutyrate sensor, and other sensor systems). The AMS farms had a higher average SCC (by 12,000 cells/mL) after sensor investment, and CMS farms with a mastitis detection system had a lower average SCC (by 10,000 cells/mL) in the years after sensor investment. Having sensor systems was associated with a higher average production per cow on AMS farms, and with a lower average production per cow on CMS farms in the years after investment. The most likely reason for this lower milk production after investment was that on 96% of CMS farms, the sensor system investment occurred together with another major change at the farm, such as a new barn or a new milking system. Most likely, these other changes had led to a decrease in milk production that could not be compensated for by the use of sensor systems. Having estrus detection sensor systems did not improve reproduction performance. Labor reduction was an important reason for investing in sensor systems. Therefore, economic benefits from investments in sensor systems can be expected more from the reduction in labor costs than from improvements in measures of health and production in dairy herds.

The optimal number of heifer calves to be reared as dairy replacements

Dairy farmers often keep almost all their newborn heifer calves despite the high cost of rearing. By rearing all heifer calves, farmers have more security and retain flexibility to cope with the uncertainty in the availability of replacement heifers in time. This uncertainty is due to mortality or infertility during the rearing period and the variation in culling rate of lactating cows. The objective of this study is to provide insight in the economically optimal number of heifer calves to be reared as replacements. A herd-level stochastic simulation model was developed specific for this purpose with a herd of 100 dairy cows; the biological part of the model consisted of a dairy herd unit and rearing unit for replacement heifers. The dairy herd unit included variation in the number of culled dairy cows. The rearing unit incorporated variation in the number of heifers present in the herd by including uncertainty in mortality and variation in fertility. The dairy herd unit and rearing unit were linked by the number of replacement heifers and culled dairy cows. When not enough replacement heifers were available to replace culled dairy cows, the herd size was temporarily reduced, resulting in an additional cost for the empty slots. When the herd size reached 100 dairy cows, the available replacement heifers that were not needed were sold. It was assumed that no purchase of cows and calves occurred. The optimal percentage of 2-wk-old heifer calves to be retained was defined as the percentage of heifer calves that minimized the average net costs of rearing replacement heifers. In the default scenario, the optimal retention was 73% and the total net cost of rearing was estimated at €40,939 per herd per year. This total net cost was 6.5% lower than when all heifer calves were kept. An earlier first-calving age resulted in an optimal retention of 75%, and the net costs of rearing were €581 per herd per year lower than in the default scenario. For herds with a lower or higher culling rate of dairy cows (10 or 40% instead of 25% in the default scenario), it was optimal to retain 35 or 100% of the heifer calves per year. Herds that had a lower or higher cost of empty slots (€50 or €120 per month instead of €82 in the default scenario) had an optimal retention of 49 or 83% per year; the optimal retention percentage was dependent on farm and herd characteristics. For Dutch dairy farming conditions, it was not optimal to keep all heifer calves.

Invited review: sensors to support health management on dairy farms

Since the 1980s, efforts have been made to develop sensors that measure a parameter from an individual cow. The development started with individual cow recognition and was followed by sensors that measure the electrical conductivity of milk and pedometers that measure activity. The aim of this review is to provide a structured overview of the published sensor systems for dairy health management. The development of sensor systems can be described by the following 4 levels: (I) techniques that measure something about the cow (e.g., activity); (II) interpretations that summarize changes in the sensor data (e.g., increase in activity) to produce information about the cow's status (e.g., estrus); (III) integration of information where sensor information is supplemented with other information (e.g., economic information) to produce advice (e.g., whether to inseminate a cow or not); and (IV) the farmer makes a decision or the sensor system makes the decision autonomously (e.g., the inseminator is called). This review has structured a total of 126 publications describing 139 sensor systems and compared them based on the 4 levels. The publications were published in the Thomson Reuters (formerly ISI) Web of Science database from January 2002 until June 2012 or in the proceedings of 3 conferences on precision (dairy) farming in 2009, 2010, and 2011. Most studies concerned the detection of mastitis (25%), fertility (33%), and locomotion problems (30%), with fewer studies (16%) related to the detection of metabolic problems. Many studies presented sensor systems at levels I and II, but none did so at levels III and IV. Most of the work for mastitis (92%) and fertility (75%) is done at level II. For locomotion (53%) and metabolism (69%), more than half of the work is done at level I. The performance of sensor systems varies based on the choice of gold standards, algorithms, and test sizes (number of farms and cows). Studies on sensor systems for mastitis and estrus have shown that sensor systems are brought to a higher level; however, the need to improve detection performance still exists. Studies on sensor systems for locomotion problems have shown that the search continues for the most appropriate indicators, sensor techniques, and gold standards. Studies on metabolic problems show that it is still unclear which indicator reflects best the metabolic problems that should be detected. No systems with integrated decision support models have been found.

First-calving age and first-lactation milk production on Dutch dairy farms

Farmers attempting to reduce first-calving age (FCA) need to understand which rearing management factors influence FCA and first-lactation milk production (FLP). Reduced FCA might be associated with lower FLP. This study describes the association between herd FCA, FLP, and several herd-level health and rearing management variables and describes the association between FCA and FLP at the cow level. It uses data from a 2010 survey of 100 Dutch dairy farms about general management, colostrum and milk feeding, housing, cleanliness, healthcare, disease, and breeding. It also used available data on FCA and 305-d FLP at both cow and herd level. The associations between median FCA and median FLP of the herd and herd-level health and rearing management variables were determined using multivariate regression analysis. The median FCA was associated with minimum age of first insemination, feeding of waste milk, and the amount of milk given preweaning. The median FLP was associated with median FCA and vaccination status for bovine respiratory syncytial virus. The association between FCA and FLP (based on 8,454 heifers) was analyzed with a single-effect linear mixed model, where the dependent variable was either FCA or relative FCA (defined as the difference between FCA of the heifer and median FCA of the herd to which they belonged). Heifers having an FCA of 24 mo produced, on average, 7,164 kg of milk per 305 d, and calving 1 mo earlier gave 143 kg less milk per 305 d. When FCA did not deviate from the median herd FCA, heifers produced, on average, 7,272 kg of milk per 305 d. From the median FCA of the herd, heifers calving 1 mo earlier produced 90 kg of milk per 305 d less, and heifers calving 1 mo later produced 86 kg per 305 d more. This is the first study that explained FLP using relative FCA. It assumes that heifers raised within the same farm have similar development because they are similarly managed. Similar management is reflected by the median FCA of the herd, with a deviation of the heifer's FCA from median FCA reflecting the heifer's development relative to the herd's average. The advantage of using relative FCA was that it accounts for between-farm differences in rearing management. It showed that earlier insemination without adjusting management to ensure sufficient development lowers FLP. An economic optimum exists between rearing costs, FCA, and FLP and, as a consequence, decisions with regard to young stock management should be made with care.

Sensors and clinical mastitis--the quest for the perfect alert

When cows on dairy farms are milked with an automatic milking system or in high capacity milking parlors, clinical mastitis (CM) cannot be adequately detected without sensors. The objective of this paper is to describe the performance demands of sensor systems to detect CM and evaluats the current performance of these sensor systems. Several detection models based on different sensors were studied in the past. When evaluating these models, three factors are important: performance (in terms of sensitivity and specificity), the time window and the similarity of the study data with real farm data. A CM detection system should offer at least a sensitivity of 80% and a specificity of 99%. The time window should not be longer than 48 hours and study circumstances should be as similar to practical farm circumstances as possible. The study design should comprise more than one farm for data collection. Since 1992, 16 peer-reviewed papers have been published with a description and evaluation of CM detection models. There is a large variation in the use of sensors and algorithms. All this makes these results not very comparable. There is a also large difference in performance between the detection models and also a large variation in time windows used and little similarity between study data. Therefore, it is difficult to compare the overall performance of the different CM detection models. The sensitivity and specificity found in the different studies could, for a large part, be explained in differences in the used time window. None of the described studies satisfied the demands for CM detection models.

The influence of cow factors on the incidence of clinical mastitis in dairy cows

Many cow-specific risk factors for clinical mastitis (CM) are known. Other studies have analyzed these risk factors separately or only analyzed a limited number of risk factors simultaneously. The goal of this study was to determine the influence of cow factors on the incidence rate of CM (IRCM) with all cow factors in one multivariate model. Also, using a similar approach, the probability of whether a CM case is caused by gram-positive or gram-negative pathogens was calculated. Data were used from 274 Dutch dairy herds that recorded CM over an 18-mo period. The final dataset contained information on 28,137 lactations of 22,860 cows of different parities. In total 5,363 CM cases were recorded, but only 2,525 CM cases could be classified as gram-positive or gram-negative. The cow factors parity, lactation stage, season of the year, information on SCC from monthly test-day records, and CM history were included in the logistic regression analysis. Separate analyses were performed for heifers and multiparous cows in both the first month of lactation and from the second month of lactation onward. For investigating whether CM was caused by gram-positive or gram-negative pathogens, quarter position was included in the logistic regression analysis as well. The IRCM differed considerably among cows, ranging between 0.0002 and 0.0074 per cow-day at risk for specific cows depending on cow factors. In particular, previous CM cases, SCC in the previous month, and mean SCC in the previous lactation increased the IRCM in the current month of lactation. Results indicate that it is difficult to distinguish between gram-positive and gram-negative CM cases based on cow factors alone.