# The Use of Meta-Analysis for the Measurement of Animal Disease Burden: Losses Due to Clinical Mastitis as an Example

Didier Raboisson, Ahmed Ferchiou, Beate Pinior, Thomas Gautier, Pierre

Sans, Guillaume Lhermie

► To cite this version:

Didier Raboisson, Ahmed Ferchiou, Beate Pinior, Thomas Gautier, Pierre Sans, et al.. The Use of Meta-Analysis for the Measurement of Animal Disease Burden: Losses Due to Clinical Mastitis as an Example. Frontiers in Veterinary Science, 2020, 7, pp.149. 10.3389/fvets.2020.00149. hal-03206995

HAL Id: hal-03206995

https://hal.inrae.fr/hal-03206995

Submitted on 23 Apr 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# The Use of Meta-Analysis for the Measurement of Animal Disease Burden: Losses Due to Clinical Mastitis as an Example

Didier Raboisson¹*, Ahmed Ferchiou¹, Beate Pinior², Thomas Gautier¹, Pierre Sans³ and Guillaume Lhermie¹

¹ IHAP, Université de Toulouse, INRAE, ENVT, Toulouse, France, ² Food Technology and Veterinary Public Health, Institute of Food Safety, University of Veterinary Medicine Vienna, Vienna, Austria, ³ ALISS, INRAE, Ivry-sur-Seine, ENVT, Toulouse, France

The literature contains an extensive panel of studies focusing on the costs of animal diseases. The losses of an agriculture holding can be influenced by many factors since farming is a complex system and diseases are closely interrelated. Meta-analysis can be used to detect effects (i.e., change in clinical mastitis losses here) across studies and to identify factors that may influence those effects. This includes the external validity of the published study results with regard to the input parameters and the internal validity of the study, particularly how other diseases related to the target disease were accounted for. Mixed-effect meta-regressions were performed to estimate the mean clinical mastitis losses per case across the literature and to elucidate to what extent clinical mastitis losses are influenced by (i) general factors, such as etiology; (ii) the types of losses that contribute to the total mastitis losses; and (iii) prices. In total, 82 observations from nine studies were included in the meta-analysis to assess mean clinical mastitis losses per case. The multivariate meta-regression showed that etiology significantly influenced the clinical mastitis loss per case. The mean loss was determined to be €224 per case for all published etiologies. In detail, mean losses equalled €457 and €101 per case of clinical mastitis due to gram-negative and gram-positive bacteria, respectively, and €428 and €74 per case of clinical mastitis due to *Escherichia coli* and *Staphylococcus aureus*, respectively. Additionally, the mean loss obtained depended on whether diagnostic costs and reduced feed intake in cases of mastitis were included in the clinical mastitis loss calculation. The monetary values of labor cost, drug cost and culling cost, as well as treatment price (all included), significantly influenced the clinical mastitis losses per case. All other tested moderators were not associated with mastitis losses, highlighting the need for more standardized economic studies, for both methods and ways results are presented, and suggesting that the mastitis losses assessed in the literature cannot be extrapolated (limited external validity). Although meta-analyses are useful to overview the burden of diseases across studies, their ability to summarize extensive literature with various economic assessments is limited. These limitations in loss assessments also suggest the need to focus on management strategies rather than on pure monetary estimations of disease costs, at least for production diseases at the farm level.

**OPEN ACCESS**

**Edited by:**
Tariq Halasa,
University of Copenhagen, Denmark

**Reviewed by:**
Anna-Maija Heikkilä,
Natural Resources Institute Finland
(Luke), Finland
Marcos Veiga,
University of São Paulo, Brazil

*Correspondence:
Didier Raboisson
didier.raboisson@envt.fr

**Specialty section:**
This article was submitted to Veterinary Epidemiology and Economics, a section of the journal Frontiers in Veterinary Science

Received: 10 December 2019

Accepted: 28 February 2020

Published: 18 March 2020

**Citation:**
Raboisson D, Ferchiou A, Pinior B, Gautier T, Sans P and Lhermie G (2020) The Use of Meta-Analysis for the Measurement of Animal Disease Burden: Losses Due to Clinical Mastitis as an Example. *Front. Vet. Sci.* 7:149.
doi: 10.3389/fvets.2020.00149

**Keywords:** clinical mastitis, economics, etiology, meta-analysis, dairy cows
# INTRODUCTION

Some reviews show large variations in the calculated impact of animal diseases, such as subclinical ketosis bovine viral diarrhea (1-3). There is increasing concerns of cost evaluation in the context of animal diseases because of difficulties in assigning disease cost to an individual disease due to the coexistence of simultaneous diseases (4). Consequently, the risk of overestimating by including the same contributors in the costs of different diseases is high (5-7). Further, key questions in disease impact evaluation is whether and how the results can be extrapolated, particularly considering the high price volatility of input and output parameters in economic assessments. Trends that focus on the whole economic strategy to manage disease rather than on the cost of disease to address this concern are increasing (5, 8). We hypothesize that meta-analysis may be an adequate approach to define how factors such as the type of incorporated losses and the associated prices may influence the value of the economic burden of the disease across the studies.

Mastitis is one of the most important diseases in dairy farms worldwide (9, 10); it is related to economic, environmental and societal stakes through losses, increased carbon and nitrogen outputs from the production process, and increased antimicrobial use (11, 12). The cost of mastitis differs across studies (1, 13, 14), particularly regarding etiology, clinical degree, types of losses included in the mastitis costs, treatment costs, level of prices, and economic assessments methods used. Mastitis is a complex disease, and its diagnosis can be clinical, bacteriologic, and cytologic. Clinical mastitis includes local and general clinical signs, and subclinical mastitis is diagnosed when no clinical signs are observed. A bacteriological diagnosis includes the identification of the etiology of mastitis and antibiotic sensitivity of the pathogen. A cytologic diagnosis is based on milk somatic cell counts (SCCs), which is a proxy generally used to measure subclinical mastitis, despite clinical mastitis lead to high SCCs.

The present work focuses on the factors that may influence losses due to animal diseases, using clinical mastitis as an example, to define whether the present state of the economics of this disease can be adapted to elucidate the (i) internal validity of the study (what is accounted for during the economic assessment) and (ii) the external validity of published study results with regard to the input parameters used. It aims to describe the usefulness of meta-analysis to evaluate which factors may influence the estimated losses due to mastitis infections in the dairy population according to the literature.

# MATERIALS AND METHODS

## Literature Search and Criteria for the Inclusion or Exclusion of Studies

Publications on the losses due to clinical mastitis were selected from English-language literature up to June 2019. The literature search was conducted in PubMed, Science Direct, and Google Scholar. The keywords were applied separately or in different combinations for the literature search in the three databases. Subsequently, the reference lists of the identified studies were also screened. All the studies were analyzed according to the

following inclusion criteria: (i) the publication included clinical mastitis and presented results for clinical mastitis separately from subclinical mastitis if both were included; (ii) the publication included the monetary losses presented per clinical mastitis case or per year and cow; (iii) the publication results were obtained in temperate-climate countries; and (iv) the publication data was obtained after 1990 to represent the modern livestock system. No restrictions were set on the intensification level (milk production), the level of monetary losses due to clinical mastitis, or the currency used. Publications focusing only on the preventive costs of clinical mastitis and/or on specific breeds (Simmental) and/or determinist methods with no variance associated with the mean losses were excluded. In the present work, production losses and curative extra costs were eligible to be considered as clinical mastitis losses. All expressions and proxies of the variance were accepted [i.e., standard deviation (SD), standard error (SE), min-max, confidence intervals] and transformed into a unique unit (i.e., SE) to compare the results across the included studies.

The total number of identified publications and the applied two-step selection process for eligible studies, which was performed in accordance with the PRISMA guidelines (Preferred Reporting Items for Systematic Reviews and Meta-Analysis), are illustrated in Figure 1. All articles were screened in full by two reviewers (TG, DR) and eligible studies, i.e., those which met the inclusion criteria, were then reviewed in full by one reviewer (TG) in accordance with the predefined variables shown in Table 1. All relevant data from the eligible studies were entered into a Microsoft Excel spreadsheet (Table S1). A publication was further divided into different observation sets if the study considered different variables, according to Table 1, into account and thus published different monetary losses per animal. Consequently, the total number of publications included in the presented study was not identical to the total number of observations. The details of the four incorporated groups of variables (called moderators) are reported in Table 1. In brief, they refer to (A) general variables (year, etiology), (B) the type of losses (i.e., nature of the contributors of the losses) considered either included or excluded from the published raw models, (C) the monetary level of losses for each contributor (in Euros), and (D) the prices used as input parameters in the published raw models.

The clinical mastitis losses reported in the literature were standardized per case of clinical mastitis. The mastitis losses were published in different national currencies and years. A standardization to the Euro (€) and the year 2018 for each respective country was performed as follows:

$$Y\_DL(€; 2018) = \frac{Y\_DL_X^i}{\tau_{conv^X(i \to €)}} * \left( \frac{I\_OCDE_{2018}}{I\_OCDE_X} \right) \quad (1)$$

where $Y\_DL(€; 2018)$ represents the mean clinical mastitis losses per animal in € in 2018, and $i$ indicates the national currency
<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Flow diagram of the selection process for the systematic literature review. The final regressions included eight articles, since the influential analysis has led to exclude one article.</figcaption></figure>

of the respective country for which losses were determined in the year X. The nominal exchange rate ($\tau_{conv}^{x(i \to e)}$) was distinguished between the Eurozone (i.e., exchange rate of the national currency *i* into the currency € in 2002) and non-Eurozone (i.e., exchange rate of the national currency *i* into the currency € in the year of publication). The index I_OCDEx includes the economic annual growth rate of the respective country and incorporates the inflation rate based on the consumer price index. The same procedure of standardization

to the Euro (€) and the year 2018 was applied for all monetary values in the dataset (Table S1).

## Meta-Analysis

The meta-analyses were implemented in R (Version 3.5.1 R Foundation for Statistical Computing, Vienna, Austria) using the Metafor package (15). Random-effects models were first used to estimate the log-effect size and its 95% confidence interval (CI) and statistical significance level. The statistical
<table><caption>TABLE 1 | Factors selected from the systematic review and considered in the meta-regression analysis.</caption><thead><tr><th>Moderators</th><th>Class</th><th>Number of observations</th><th>Mean (± SE) value of the clinical mastitis losses (in €)</th><th>Meta-regression set in which the data is included</th></tr></thead><tbody><tr><td colspan="5"><strong>(A) GENERAL</strong></td></tr><tr><td>Study type</td><td>Modeling</td><td>42</td><td>293 ± 105</td><td>1, 2, and 3</td></tr><tr><td></td><td>Descriptive</td><td>40</td><td>227 ± 159</td><td></td></tr><tr><td>Publication year</td><td>Numeric</td><td>82</td><td>262 ± 137</td><td>1, 2, and 3</td></tr><tr><td>Country</td><td>Nominal</td><td>82</td><td>262 ± 137</td><td>1, 2, and 3</td></tr><tr><td>Number of herds</td><td>Numeric</td><td>51</td><td>267 ± 145</td><td>1</td></tr><tr><td>Average herd size</td><td>Numeric</td><td>49</td><td>261 ± 146</td><td>1</td></tr><tr><td>Number of clinical cases</td><td>Numeric</td><td>20</td><td>291 ± 129</td><td>1</td></tr><tr><td>Average milk yield<sup>a</sup></td><td>Numeric</td><td>35</td><td>305 ± 110</td><td>1</td></tr><tr><td>Parity</td><td>All</td><td>68</td><td>245 ± 144</td><td>1</td></tr><tr><td></td><td>Primiparous</td><td>6</td><td>287 ± 19</td><td></td></tr><tr><td></td><td>Multiparous</td><td>8</td><td>379 ± 17</td><td></td></tr><tr><td>Incidence (%)</td><td>0.12</td><td>6</td><td>287 ± 19</td><td>1</td></tr><tr><td></td><td>0.20-0.24</td><td>8</td><td>379 ± 17</td><td></td></tr><tr><td></td><td>0.35</td><td>1</td><td>87</td><td></td></tr><tr><td>Prevalence (%)</td><td></td><td>48</td><td>225 ± 55</td><td>1</td></tr><tr><td>Etiology 1</td><td>All pathogens</td><td>48</td><td>282 ± 112</td><td>1</td></tr><tr><td></td><td>Gram positive</td><td>16</td><td>155 ± 68</td><td></td></tr><tr><td></td><td>Gram negative</td><td>6</td><td>477 ± 119</td><td></td></tr><tr><td></td><td>Other (no growth, two pathogens)</td><td>12</td><td>235 ± 173</td><td></td></tr><tr><td>Etiology 2</td><td>All</td><td>48</td><td>282 ± 112</td><td>1</td></tr><tr><td></td><td>S. aureus<sup>e</sup></td><td>4</td><td>123 ± 30</td><td></td></tr><tr><td></td><td>S. coag.<sup>e</sup></td><td>4</td><td>168 ± 104</td><td></td></tr><tr><td></td><td>S. spp.<sup>f</sup></td><td>8</td><td>165 ± 64</td><td></td></tr><tr><td></td><td>Gram negative</td><td>6</td><td>444 ± 108</td><td></td></tr><tr><td></td><td>(Other no growth, two pathogens)</td><td>12</td><td>264 ± 195</td><td></td></tr><tr><td>Etiology 3</td><td>All</td><td>48</td><td>282 ± 112</td><td>1</td></tr><tr><td></td><td>S. aureus<sup>d</sup></td><td>4</td><td>123 ± 30</td><td></td></tr><tr><td></td><td>S. coag.<sup>e</sup></td><td>4</td><td>168 ± 104</td><td></td></tr><tr><td></td><td>Streptococcus Esculine +</td><td>4</td><td>152 ± 15</td><td></td></tr><tr><td></td><td>Streptococcus Esculine -</td><td>4</td><td>178 ± 95</td><td></td></tr><tr><td></td><td>Gram negative</td><td>6</td><td>444 ± 108</td><td></td></tr><tr><td></td><td>Other (no growth, two pathogens)</td><td>12</td><td>264 ± 195</td><td></td></tr><tr><td colspan="5"><strong>(B) TYPE OF IOSSES: CONTRIBUTORS TO MASTITIS LOSSES (ACCOUNTED FOR OR NOT)</strong></td></tr><tr><td>Diagnosis (before treatment)<sup>b</sup></td><td>No</td><td>76</td><td>261 ± 140</td><td>1</td></tr><tr><td></td><td>Yes</td><td>4</td><td>281 ± 15</td><td></td></tr><tr><td>Feed intake (saved if mastitis)<sup>c</sup></td><td>No</td><td>46</td><td>221 ± 151</td><td>1</td></tr><tr><td></td><td>Yes</td><td>36</td><td>311 ± 97</td><td></td></tr><tr><td>Milk withdrawal</td><td>No</td><td>0</td><td></td><td>No</td></tr><tr><td></td><td>Yes</td><td>82</td><td>262 ± 137</td><td></td></tr><tr><td>Milk not produced</td><td>No</td><td>4</td><td>87 ± 17</td><td>1</td></tr><tr><td></td><td>Yes</td><td>78</td><td>271 ± 134</td><td></td></tr><tr><td>Veterinary cost</td><td>No</td><td>20</td><td>291 ± 129</td><td>1</td></tr><tr><td></td><td>Yes</td><td>62</td><td>252 ± 139</td><td></td></tr><tr><td>Drug cost</td><td>No</td><td>0</td><td></td><td>No</td></tr><tr><td></td><td>Yes</td><td>82</td><td>262 ± 137</td><td></td></tr><tr><td>Extra labor</td><td>No</td><td>3</td><td>87 ± 21</td><td>1</td></tr><tr><td></td><td>Yes</td><td>77</td><td>268 ± 135</td><td></td></tr></tbody></table>

(Continued)
<table><caption>TABLE 1 | Continued</caption><thead><tr><th>Moderators</th><th>Class</th><th>Number of observations</th><th>Mean (± SE) value of the clinical mastitis losses (in €)</th><th>Meta-regression set in which the data is included</th></tr></thead><tbody><tr><td rowspan="2">Culling</td><td>No</td><td>4</td><td>87 ± 17</td><td rowspan="2">1</td></tr><tr><td>Yes</td><td>76</td><td>271 ± 134</td></tr><tr><td rowspan="2">Extended day open</td><td>No</td><td>66</td><td>245 ± 144</td><td rowspan="2">1</td></tr><tr><td>Yes</td><td>14</td><td>340 ± 50</td></tr><tr><td rowspan="2">Cow Mortality</td><td>No</td><td>19</td><td>175 ± 74</td><td rowspan="2">1</td></tr><tr><td>Yes</td><td>63</td><td>289 ± 141</td></tr><tr><td rowspan="2">Carcass disposal</td><td>No</td><td>19</td><td>175 ± 74</td><td rowspan="2">1</td></tr><tr><td>Yes</td><td>63</td><td>289 ± 141</td></tr><tr><td rowspan="2">Milk replacer used</td><td>No</td><td>32</td><td>257 ± 125</td><td rowspan="2">1</td></tr><tr><td>Yes</td><td>48</td><td>265 ± 145</td></tr><tr><td colspan="5" style="text-align: center; font-weight: bold;">(C) MONETARY LEVEL OF LOSSES (VALUE OF EACH CONTRIBUTOR)</td></tr><tr><td>Milk withdrawal</td><td></td><td>13</td><td>261 ± 78</td><td>2</td></tr><tr><td>Milk not produced</td><td></td><td>24</td><td>245 ± 83</td><td>2</td></tr><tr><td>Veterinary cost</td><td></td><td>16</td><td>221 ± 100</td><td>2</td></tr><tr><td>Drug cost</td><td></td><td>24</td><td>328 ± 78</td><td>2</td></tr><tr><td>Extra labor</td><td></td><td>26</td><td>239 ± 87</td><td>2</td></tr><tr><td>Culling</td><td></td><td>25</td><td>245 ± 83</td><td>2</td></tr><tr><td>Extended day open</td><td></td><td>2</td><td>328 ± 61</td><td>No</td></tr><tr><td>Cow Mortality</td><td></td><td>3</td><td>306 ± 57</td><td>No</td></tr><tr><td colspan="5" style="text-align: center; font-weight: bold;">(D) PRICES (OF INPUT PARAMETERS OF RAW MODELS)</td></tr><tr><td rowspan="2">Cow culled (€/kg carcass)</td><td>1.69</td><td>12</td><td>341 ± 50</td><td rowspan="2">3</td></tr><tr><td>1.94</td><td>2</td><td>334 ± 70</td></tr><tr><td rowspan="2">Replacement heifer (€/head)</td><td>1502</td><td>12</td><td>341 ± 51</td><td rowspan="2">3</td></tr><tr><td>1684</td><td>2</td><td>330 ± 59</td></tr><tr><td rowspan="3">Milk (€/kg)</td><td>0.31</td><td>12</td><td>375 ± 82</td><td rowspan="3">3</td></tr><tr><td>0.33-0.37</td><td>15</td><td>316 ± 83</td></tr><tr><td>0.41-0.49</td><td>5</td><td>366 ± 62</td></tr><tr><td rowspan="2">Feed (€/kg) dry matter</td><td>0.16</td><td>12</td><td>345 ± 48</td><td rowspan="2">3</td></tr><tr><td>0.18</td><td>2</td><td>307 ± 64</td></tr><tr><td rowspan="4">Labor (€/h)</td><td>5.89-10.58</td><td>13</td><td>353 ± 112</td><td rowspan="4">3</td></tr><tr><td>19.42-23.88</td><td>16</td><td>199 ± 64</td></tr><tr><td>28.88-30.36</td><td>14</td><td>340 ± 50</td></tr><tr><td>36.76</td><td>35</td><td>224 ± 153</td></tr><tr><td>Treatment (€/treatment, all included)</td><td>Numeric</td><td>32</td><td>253 ± 124</td><td>3</td></tr></tbody></table>

heterogeneity between and within studies was assessed using the Cochran Q statistic and the $I^2$ statistic, respectively (16). For response variables with high $I^2$, univariate meta-regression was then performed to explore the sources of heterogeneity. The meta-regression was conducted by screening for the moderators, as described in Table 1. A moderator was a variable that resulted in reduced heterogeneity when introduced in the meta-regression (i.e., factor). In the first step, the meta-regression was performed for all factors together and then separately for (A) the general factors, (B) the type of losses (i.e., contributor), (C) the monetary level of losses for each contributor, and (D) the prices used for the economic input parameters. The variable “Publication” was kept as a random effect. The inclusion of the factors in the meta-regression analysis

was conducted as follows. Univariate meta-regressions were first performed to identify factors according to Table 1 that may have had a significant association with the clinical mastitis loss per case. A reference class for each factor was chosen to allow a comparison of the effect size. Any significant factors in the univariate test were selected as a potential influencing factor for the multivariate analysis, which aimed to reduce the heterogeneity between the included studies in the meta-analysis. The $\tau^2$ (residual heterogeneity variance) denoted the amount of heterogeneity that may not be explained through the inclusion of the factors in the meta-analysis. Publication bias was identified by performing the Egger test, a regression test for funnel plot asymmetry and inspection of the associated funnel plots (Figure 2). Outliers were also identified by conducting

<footnote>a 305-days average milk production; b mastitis diagnosis before treatment; c adjustment for reduced feed intake in cases of mastitis; d Staphylococcus aureus; e coagulase negative-Staphylococcus; f Escherichia coli.</footnote>
<figure><img src="image_4.png" /><figcaption>**FIGURE 2** | Funnel plot of the random meta-analysis of studies, without the incorporating moderators in **Table 1**.</figcaption></figure>

influential case diagnostics (i.e., DFFITS value, Cook's distances, covariance ratios, estimates of τ² and test statistics for (residual) heterogeneity). Because the dataset contained moderators with different numbers of missing data, three subdatasets were proposed, as described in **Table 1**, and meta-regression was performed accordingly.

## RESULTS

In total, 82 observations from nine studies were included in the meta-analysis (Table S2). Influential case diagnostics indicated two observations and one study as sources of asymmetry. These observations were considerably higher regarding clinical mastitis loss per case (with a mean of €1,000) than other observations (with a mean ± SD of €262 ± 137) and thus highly influenced the results of the meta-regression. Consequently, two observations and one study were excluded in the presented meta-analysis. The estimated pooled effect size obtained with the random-effects model with no moderator (Table 2) was €195 (Se = 37, *p* < 0.001). No publication bias was determined with the Egger test (*t* = 269, *p* = 0.21). The heterogeneity between the studies was very high (*I*² = 99.9%; AIC = 99,711; Q-Test: *x*² = 99,136; df = 79; *p* < 0.001).

The meta-regression with all moderators, which was performed on the dataset including the moderators without any missing data (denoted subdataset 1 in **Table 1**), showed that mastitis losses were associated with the etiology of mastitis as well as with the inclusion diagnostic costs and feed intake decreases

in cases of mastitis in the raw model (Table 2). The observed decrease in the heterogeneity through the inclusion of these factors was 33%. A diagnosis before treatment was associated with an extra loss of €155 per case, and the adjustment of the cost assessment by the diet saved with a reduced the loss per case by €35. The average losses of gram positive and gram negative clinical mastitis were €101 (€224–€123) and €457 (€224+€233), respectively. The losses were €74 (€224–€150), €79 (€224–€145), €121 (€224–€103) and €428 (€224+€204) for mastitis due to *Staphylococcus aureus*, coagulase-*Staphylococcus*, *Streptococcus* spp., and *Escherichia coli*, respectively (Table 2).

The dataset focusing on the monetary levels of losses includes 26 observations (denoted subdataset 2, see **Table 1**). All moderators summarized in the term monetary levels of losses in **Table 1** were significantly associated with the cost of mastitis, but correlation was observed between the moderators (Table 3), leading to the final regression proposed in **Table 2**. The marginal values of drug cost, labor cost and culling cost in clinical mastitis losses were €0.8, €2.9, and €1.04, respectively. This means, for instance, that one extra euro for the drug cost was associated with an extra clinical mastitis losses equal to €0.8.

The dataset focusing on prices of inputs included 12 to 76 observations, depending on the moderators (denoted subdataset 3, see **Table 1**). The price of treatment (all included) was associated with a marginal value in the clinical mastitis-related loss of €0.76 (Se = 0.04, *p* < 0.0001). All other moderators, including the price of milk, were not associated with the losses of clinical mastitis (*p* > 0.7 for the three classes compared to the reference class, **Table 1**).

The influential case diagnostics of the three meta-regressions shown in **Table 2** indicated outliers in the incorporated influencing factors on mastitis losses (Figures S1–S3). The removal of the outliers did not change the coefficients of the meta-regressions, and the results shown in **Table 2** were considered final. Final forest plots are reported in **Figures 3–5**.

## DISCUSSION

The meta-regression was performed according to usual recommendations (17, 18). The final choice for the models was made considering the decrease in heterogeneity. More than one model was reported for the same outcome because the authors judged that all the models had biological significance and would be of interest to the scientific community. The multivariate models provided in **Table 2** show close coefficients compared to the univariate models, and the addition of a new moderator reduced the heterogeneity. Unfortunately, many studies did not report any estimation of the variance (determinist method) and could not be included in the present meta-regression, leading to only nine included publications, although extensive literature is available. This issue has been highlighted in a previous review (13). The present meta-analysis and the previous review (13) both highlight the high heterogeneity within the method used to assess mastitis losses, the nature of the included losses and the limits of comparing
<table><caption>**TABLE 2** | Final factors considered in the meta-regression analysis.</caption><thead><tr><th>Group of moderators<sup>a</sup></th><th>Moderator and class</th><th>Estimate (SE)</th><th>95% CI</th><th>P-value</th></tr></thead><tbody><tr><td>Without</td><td>Intercept</td><td>195 (37)</td><td>122/267</td><td>&lt;0.0001</td></tr><tr><td rowspan="6">A and B: General and type of mastitis losses</td><td>Intercept</td><td>224 (43)</td><td>139/308</td><td>&lt;0.0001</td></tr><tr><td>Gram positive (ref=All)</td><td>-123 (7.3)</td><td>-108/-137</td><td>&lt;0.0001</td></tr><tr><td>Gram negative (ref=All)</td><td>233 (16.3)</td><td>201/264</td><td>&lt;0.0001</td></tr><tr><td>Other (ref=All)</td><td>-133 (7.0)</td><td>-119/-146148/161</td><td>&lt;0.0001</td></tr><tr><td>Diagnosis (ref= No)</td><td>155 (3.5)</td><td>148/161</td><td>&lt;0.0001</td></tr><tr><td>Feed intake (ref = No)</td><td>-29 (2.5)</td><td>-24/-34 139/308</td><td>&lt;0.0001</td></tr><tr><td rowspan="7">A and B: General and type of mastitis losses<sup>b</sup></td><td>Intercept</td><td>224 (43)</td><td>139/308</td><td>&lt;0.0001</td></tr><tr><td>Staphylococcus aureus (ref=All)</td><td>-150 (9)</td><td>-132/-167</td><td>&lt;0.0001</td></tr><tr><td>Staphylococcus spp. (ref=All)</td><td>-145 (11)</td><td>123/166 166166</td><td>&lt;0.0001</td></tr><tr><td>Streptococcus spp. (ref=All)</td><td>-103 (8)</td><td>-87/-118</td><td>&lt;0.0001</td></tr><tr><td>E. coli (ref=All)</td><td>204 (18)</td><td>168/239</td><td>0.0097</td></tr><tr><td>Other (ref=All)</td><td>-131 (7.0)</td><td>-117/-144</td><td>&lt;0.0001</td></tr><tr><td>Diagnosis (ref= No)</td><td>155 (3.5)</td><td>148/161</td><td>&lt;0.0001</td></tr><tr><td rowspan="4">C- Monetary level of losses<sup>c</sup></td><td>Feed intake (ref = No)</td><td>-29 (2.5)</td><td>-24/-34</td><td>&lt;0.0001</td></tr><tr><td>Intercept</td><td>124 (43)</td><td>39/208</td><td>0.0039</td></tr><tr><td>Labour<sup>e</sup></td><td>2.9 (0.24)</td><td>2.4/3.4</td><td>&lt;0.0001</td></tr><tr><td>Drug cost<sup>e</sup></td><td>0.8 (0.04)</td><td>0.72/0.87</td><td>&lt;0.0001</td></tr><tr><td rowspan="2">D- Prices<sup>d</sup></td><td>Culling<sup>e</sup></td><td>1.04 (0.02)</td><td>1.00/1.08</td><td>&lt;0.0001</td></tr><tr><td>Intercept</td><td>150 (39)</td><td>73/226</td><td>&lt;0.0001</td></tr><tr><td></td><td>Treatment (all included)<sup>f</sup></td><td>0.76 (0.04)</td><td>0.68/0.83</td><td>&lt;0.0001</td></tr></tbody></table>

<table><caption>**TABLE 3** | Correlations (and P-values) of the values of the moderators in the category "monetary level of losses" (moderator group C, see **Table 1**).</caption><thead><tr><th></th><th>Milk not produced</th><th>Milk withdrawal</th><th>Veterinary cost</th><th>Drug cost</th><th>Extra labor</th><th>Culling</th></tr></thead><tbody><tr><td>Milk withdrawal</td><td>0.92/&lt;0.001<sup>a</sup></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Veterinary cost</td><td>0.92/&lt;0.001</td><td>0.94/&lt;0.001</td><td></td><td></td><td></td><td></td></tr><tr><td>Drug cost</td><td>-0.02/0.89</td><td>0.36/0.06</td><td>0.43/0.06</td><td></td><td></td><td></td></tr><tr><td>Extra labor</td><td>0.77/&lt;0.001</td><td>0.9/&lt;0.001</td><td>0.85/&lt;0.001</td><td>0.04/0.8</td><td></td><td></td></tr><tr><td>Culling</td><td>0.68/&lt;0.001</td><td>0.6/&lt;0.001</td><td>0.53/0.005</td><td>-0.34/0.9</td><td>0.39/0.01</td><td></td></tr><tr><td>Mortality</td><td>-0.2/0.38</td><td>0.32/0.26</td><td>-0.61/&lt;0.001</td><td>0.61/0.01</td><td>0.67/0.09</td><td>0.78/0.06</td></tr></tbody></table>

or summarizing results. In spite that one aim of the present work was to adjust the estimation of clinical mastitis losses by the occurrence of other diseases, data available did not permit to do for most of the regressions since papers included in the meta-analysis only scarcely reported other diseases that may interact with mastitis. Moreover, even if many non-significant associations were found in the present study, it helps (i) to determine factors influencing mastitis losses, (ii) to quantify heterogeneity, and (iii) to highlight the concerns faced when aiming at reducing heterogeneity in the context of clinical mastitis losses.

The present study should remind the reader that focusing on the cost or loss of disease may be of limited value for

some diseases, such as mastitis. The present results show that the losses of clinical mastitis were higher for gram negative mastitis than for gram positive mastitis, although opposite perception has been reported from the field. For instance, the management of gram negative mastitis is easier than gram positive mastitis, although the long-term consequences of udder contamination by gram positive are greater than those for gram negative mastitis (treatment efficacy, chronic infection). This result is in accordance with the fact that medium- or long-term SCCs arising from clinical mastitis are poorly accounted for in the present works since they are scarcely reported in the literature. Similarly, the way culling is integrated in the studies remains unclear. This result demonstrates that any economic

<footnote><sup>a</sup>as defined in **Table 1**.</footnote>

<footnote><sup>b,c,d</sup>the corresponding influential case diagnostics are indicated in **Figures S1, S2, and S3**, respectively.</footnote>

<footnote><sup>e</sup>the moderator is a continuous variable that equals the monetary value of the contributor of the losses due to clinical mastitis. The coefficient is then expressed as the marginal value of loss (for one extra euro of the total value of the contributors "labor cost," "treatment cost," and "culling cost").</footnote>

<footnote><sup>f</sup>the moderator is a continuous variable that equals the price of treatment considered in the raw model. The coefficient is then expressed as the marginal value of loss (for one extra euro of treatment price).</footnote>

<footnote><sup>a</sup>Expressed as correlation value/P-value.</footnote>
<figure><img src="image_5.png" /><figcaption>**FIGURE 3** | Forest graph of the meta-regressions including moderators in groups A and B. The column on the right refers to the mean loss per case with the corresponding confidence interval (in brackets). The two single letters in the left column represent the moderators diagnosis and feed intake, as defined in **Table 1**. The moderator etiology (Table 1) is located to the left of the authors. The gray diamonds represent the effect size adjusted for the moderator and are included in the meta-regression.</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**FIGURE 4** | Forest graph of the meta-regression including moderators in group C (see **Table 1**). The column on the right refers to the mean loss per case with the corresponding confidence interval (in brackets). The gray diamonds represent the effect size adjusted for the moderators included in the meta-regression.</figcaption></figure>

reasoning focusing on the losses of mastitis is inappropriate and should be avoided. The present meta-analysis does not provide additional information since the economic reasoning in most of the raw models used in the present meta-analysis was biased. Recent literature on economics of udder health has increasingly focused on the strategic management of mastitis instead of its losses (19–21). This appears to be an appropriate trend since it accounts for herd dynamics, short- and long-term issues and farmers' behaviors, for instance, through

farmers' satisfaction not only relying on income optimization (utility). This is also a reason why some recent publications cannot be included in the present meta-analysis that focus on clinical mastitis losses. The exclusion criteria have also led to focus on the period from 1990 to 2019. The milk production environment has changed a lot worldwide during this period. The period was included in the present analysis as moderator but it was not significantly associated with the clinical mastitis losses. Altogether, the PRISMA procedure, the outcome
<figure><img src="image_7.png" /><figcaption>**FIGURE 5** | Forest graph of the meta-regression including moderators in group D (see **Table 1**). The column on the right refers to the mean loss per case with the corresponding confidence interval (in brackets). The gray diamonds represent the effect size adjusted for the moderator included in the meta-regression.</figcaption></figure>

studied, the period focused and the geographical restrictions have led to exclude many studies that deal with economics of mastitis.

The present work tried to explain the losses of clinical mastitis considering different groups of moderators (**Table 1**). The moderator groups B, C, and D (see **Table 1**) referred to the questions “is the moderator accounted for? (yes/no),” “can the size of the contributor (total monetary value) be linked to clinical mastitis losses?” i.e., “is the share of the contributor almost always

the same?” and “can we summarize how the economic input parameters influence the total losses?,” respectively. Moderator group B showed that factors were systematically included, scarcely included (adjustment for reduced feed intake), or almost never included (reproduction impact, pre-treatment diagnosis). These differences contributed to the large heterogeneity observed in the outcome of the present study. Unfortunately, moderators in group C were scarcely reported. Publications precisely reported the main contributors included in the economic
assessment (moderators in group B), but most of them failed to clearly describe the value of this contributor (moderators in group C). This is a key limit to evaluating whether the share of the contributors substantially changed between publications, and it does not allow here to draw a conclusion on the external validity of the literature, even if it is extensive. Such a standardization would require precise lists of items to be reported and clear procedures to be followed to perform and evaluate the analysis, as it exists for other scientific approaches (see PRISMA as the reference method for meta-analyses). In a companion paper focusing on BVD, the same limits were highlighted, and the contributors of the total losses of BVDV infection could not be defined precisely due to imprecision in the included contributors and a lack of clear reporting of the values of each contributor (4). Last, the values of the economic parameters had almost no association with the clinical mastitis losses. This is also due to unclear results in the publications and the low variability within the values of the moderators. In addition, no association was reported between the input values and the outcome for moderators such as price of milk or labor (Table 1), although the number of trials included was moderate to high and the range of the values appeared to be appropriate for regression. The wide range in the values of milk price without any significant association was in opposition with the statement that the price of milk influences the losses of the disease and may be an item to consider to adapt the disease management strategy to the market context, which is often reported in sensitivity analyses in publications and in the field. The present results do not support such a relationship. Based on the present results and due to the heterogeneity, the conclusion is that the scientific community should be very careful to use the monetary values from one study out of its context or for their own study design, because the mastitis losses depend on many other factors which cause the heterogeneity that

could not be explained by the factors considered in the present study.

## CONCLUSIONS

The present work proposes an estimation of the losses of clinical mastitis for different etiologies. It failed to elucidate which contributor mainly influenced the losses of clinical mastitis and did not highlight any relationship between the price of milk or labor and losses due to clinical mastitis. This supports the avoidance any economic reasoning focused on the losses of each case of clinical mastitis since such a reasoning is inappropriate from an economic point of view. The internal and external validity of the losses evaluation is highly questionable in the case of clinical mastitis. These results also highlight the need for standardization on how economic assessments of losses due to animal diseases should be performed. This includes both methods and ways results are presented.

## DATA AVAILABILITY STATEMENT

All datasets generated for this study are included in the article/Supplementary Material.

## AUTHOR CONTRIBUTIONS

DR and TG designed and performed the analysis. BP, AF, PS, and GL contributed to the analysis, reviewed, and revised the manuscript.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2020.00149/full#supplementary-material

## REFERENCES

1. Raboisson D, Mounié M, Maigné E. Diseases, reproductive performance, and changes in milk production associated with subclinical ketosis in dairy cows: a meta-analysis and review. *J Dairy Sci.* (2014) 97:7547-63. doi: 10.3168/jds.2014-8237

2. Pinior B, Firth CL, Richter V, Lebl K, Trauffler M, Dzieciol M, et al. A systematic review of financial and economic assessments of bovine viral diarrhea virus (BVDV) prevention and mitigation activities worldwide. *Prev Vet Med.* (2017) 137:77-92. doi: 10.1016/j.prevetmed.2016.12.014

3. Richter V, Lebl K, Baumgartner W, Obritzauser W, Käsbohrer A, Pinior B. A systematic worldwide review of the direct monetary losses in cattle due to bovine viral diarrhoea virus infection. *Vet J.* (2017) 220:80-7. doi: 10.1016/j.tvjl.2017.01.005

4. Pinior B, Garcia S, Minviel JJ, Raboisson D. Epidemiological factors and mitigation measures influencing production losses in cattle due to bovine viral diarrhoea virus infection: a meta-analysis. *Transbound Emerg Dis.* (2019) 66:2426-39. doi: 10.1111/tbed.13300

5. Raboisson D, Ferchiou A, Lhermie G, Sans P. *Multi-criteria Optimisation to Fix the Limits of Present Standards in Microeconomics of Animal Health: The Example of Dairy Production.* Atlanta, GA: The International Society for Economics and Social Sciences of Animal Health (Atlanta: ISESSAH).

6. GBADs. GBADS - Global Burden of Animal Diseases. (2019). Available online at: https://animalhealthmetrics.org/about/ (accessed November 27, 2019).

7. Ferchiou A, Raboisson D, Minviel JJ. Bio-economic modelling of antimicrobial use and health management in French dairy production. In: *The Agricultural Economics Society Annual Meeting (AES)*. Coventry (2018). doi: 10.22004/ag.econ.273484

8. Marschik T, Obritzauser W, Wagner P, Richter V, Mayerhofer M, Egger-Danner C, et al. A cost-benefit analysis and the potential trade effects of the bovine viral diarrhoea eradication programme in Styria, Austria. *Vet J.* (2018) 231:19-29. doi: 10.1016/j.tvjl.2017.11.010

9. Hillerton JE, Berry EA. Treating mastitis in the cow-a tradition or an archaism. *J Appl Microbiol.* (2005) 98:1250-5. doi: 10.1111/j.1365-2672.2005.02649.x

10. Seegers H, Fourichon C, Beauveau F. Production effects related to mastitis and mastitis economics in dairy cattle herds. *Vet Res.* (2003) 34:475-91. doi: 10.1051/vetres:2003027

11. Schabauer A, Pinior B, Gruber C-M, Firth CL, Käsbohrer A, Wagner M, et al. The relationship between clinical signs and microbiological species, spa type, and antimicrobial resistance in bovine mastitis cases in Austria. *Vet microbiol.* (2018) 227:52-60. doi: 10.1016/j.vetmic.2018.10.024
12. Firth CL, Käsbohrer A, Egger-Danner C, Fuchs K, Pinior B, Roch F-F, et al. Comparison of defined course doses (DCDvet) for blanket and selective antimicrobial dry cow therapy on conventional and organic farms. *Animals*. (2019) 9:707. doi: 10.3390/ani9100707

13. Halasa T, Huijps K, Østerås O, Hogeveen H. Economic effects of bovine mastitis and mastitis management: a review. *Vet Q*. (2007) 29:18–31. doi: 10.1080/01652176.2007.9695224

14. Hogeveen H, Steeneveld W, Wolf CA. Production diseases reduce the efficiency of dairy production: a review of the results, methods, and approaches regarding the economics of mastitis. *Annu Rev Resour Econ*. (2019) 11:289–312. doi: 10.1146/annurev-resource-100518-093954

15. Viechtbauer W. Conducting meta-analyses in R with the metafor package. *J Stat Softw*. (2010) 36:1–48. doi: 10.18637/jss.v036.i03

16. Higgins JPT, Thompson SG. Quantifying heterogeneity in a meta-analysis. *Stat Med*. (2002) 21:1539–58. doi: 10.1002/sim.1186

17. Schwarzer G, Carpenter JR, Rücker G. *Meta-analysis with R*. Cham; Freiburg; London: Springer (2015) doi: 10.1007/978-3-319-21416-0

18. Endo R, Ishii A, Nakanishi A, Nabenishi H, Ashizawa K, Tsuzuki Y. Effect of the addition of β-hydroxybutyrate to chemically defined maturation medium on the nuclear maturation, sperm penetration and embryonic development of porcine oocytes *in vitro*. *Asian-Australas J Anim Sci*. (2010) 23:1421–6. doi: 10.5713/ajas.2010.10073

19. Gussmann M, Steeneveld W, Kirkeby C, Hogeveen H, Nielen M, Farre M, et al. Economic and epidemiological impact of different intervention strategies for clinical contagious mastitis. *J Dairy Sci*. (2019) 102:1483–93. doi: 10.3168/jds.2018-14939

20. Gussmann M, Denwood M, Kirkeby C, Farre M, Halasa T. Associations between udder health and culling in dairy cows. *Prev Vet Med*. (2019) 171:104751. doi: 10.1016/j.prevetmed.2019.104751

21. Getaneh AM, Mekonnen SA, Hogeveen H. Stochastic bio—economic modeling of mastitis in Ethiopian dairy farms. *Prev Vet Med*. (2017) 138:94–103. doi: 10.1016/j.prevetmed.2017.01.014

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2020 Raboisson, Ferchiou, Pinior, Gautier, Sans and Lhermie. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.