# Quantification of the increase in the frequency of early calving associated with late exposure to bluetongue virus serotype 8 in dairy cows: implications for syndromic surveillance

Simon Nusinovici, Aurélien Madouasse, Christine Fourichon

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

## To cite this version:

Simon Nusinovici, Aurélien Madouasse, Christine Fourichon. Quantification of the increase in the frequency of early calving associated with late exposure to bluetongue virus serotype 8 in dairy cows: implications for syndromic surveillance. Veterinary Research, 2016, 47 (1), pp.18. 10.1186/s13567-015-0296-7. http://hal-01341470

HAL Id: hal-01341470

https://hal.science/hal-01341470

Submitted on 4 Jul 2016

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

Open Access

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

CrossMark

# Quantification of the increase in the frequency of early calving associated with late exposure to bluetongue virus serotype 8 in dairy cows: implications for syndromic surveillance

Simon Nusinovici¹,<fnref n="2" />, Aurélien Madouasse¹,<fnref n="2" />* and Christine Fourichon¹,<fnref n="2" />

## Abstract

A recent study evaluating whether reproductive data could be used for syndromic surveillance found an increased frequency of early calving (calving occurring a few days earlier than expected) in areas exposed to the Bluetongue virus serotype 8 (BTV-8) in northern Europe. A high proportion of herds infected during the 2006–2009 European outbreak were not reported through the surveillance system. The objectives of this study were (1) to quantify the increase in the frequency of early calving associated with the exposure to BTV-8 in late gestation and (2) to determine whether this association could be found in populations exposed to BTV-8 but without reported clinical signs. Increases in frequency of early calving were quantified for cows in herds located in the 2007 outbreak area in France, reported or not as cases. Increases were detected for cows in both categories of herds with a larger effect in herds reported after clinical signs. Moreover, the largest effect was found for exposures occurring during the latest stage of pregnancy, suggesting that BTV infection could trigger calving in cows in late gestation, a few days earlier than expected. This is the first study quantifying the association between a viral infection and a shortened pregnancy length (still within a normal range). The high magnitude of the increase in frequency of early calving, their occurrence in herds from infected areas but not reported, and the short time interval between exposure and the occurrence of the event confirm the interest of using early calving as an indicator for syndromic surveillance.

## Introduction

The number of infectious diseases to emerge in humans has been increasing since the 1940s [1]. In Europe, climate change could favor the emergence of vector-borne emerging diseases by expanding the geographical distribution of vectors further north [2, 3] while international trade increases the chance of introducing vectors and pathogens. In animals, the emergence in northern Europe of the Bluetongue virus serotype 8 (BTV-8) in 2006 and of the Schmallenberg virus in 2011 illustrate this increasing risk in cattle,

associated with the introduction of known pathogens in free areas or with the emergence of an unknown pathogen. Precise description of the disease effects and of their frequencies in exposed populations is needed to quantify diseases consequences and to refine surveillance. In the case of BTV8, clinical signs [4–6] as well as increased mortality [7–9] were described. The impact of BTV-8 on milk production [10, 11] and reproductive performance were quantified. BTV-8 was associated with both a decrease in fertility [12, 13] and an increased risk of abortion [14]. This impact on reproduction was due to the infection in pregnant cows (with a tropism of the virus in the genital tract) possibly followed by the infection of the fetus.

A recent study evaluated whether reproduction data could be used to build indicators suitable for syndromic

<footnote>*Correspondence: aurelien.madouasse@oniris-nantes.fr</footnote>

<footnote n="2">² LUNAM Université, Oniris, UMR1300 Biology, Epidemiology and Risk Analysis in Animal Health, CS 40706, F-44307 Nantes, France
Full list of author information is available at the end of the article</footnote>
surveillance [15]. This study showed an increased frequency of early calving (after a pregnancy length still within a normal range) in cows located in BTV infected areas, concomitant with the BTV notifications. Overall, an indicator based on the frequency of early calving was able to detect BTV outbreaks at an early stage. Given the short time interval between the notification of the first case and the occurrence of a detectable effect at a regional level, the delay between infection and the occurrence of an early calving should be short. Furthermore, because infections during early stages of gestation induced the most severe malformations [16], an event such as early calving in the normal range of gestation would be more likely associated with an exposure at late stages of gestation. Consequently, it is hypothesized that a viral infection at a late stage of pregnancy could trigger calving. No study has investigated the effect of BTV-8 infection for cows in late gestation on the frequency of early calving.

The surveillance system implemented in France was able to show the geographical progression of the BTV-8 epizootic but did not allow a precise estimation of its extent in terms of prevalence of infection [17]. This cross-sectional serologic study conducted in 2007 showed a high under-reporting rate. Moreover, a decrease in fertility was found for cows in herds located in the 2007 outbreak area but not reported as cases [18]. This decreased fertility should be at least partly due to the undetected or unreported viral circulation.

The objectives of this study were (1) to quantify the possible increase in frequency of early calving after a normal length of gestation associated with exposure to BTV-8 in late gestation and (2) to determine whether

this association could be found in populations exposed to BTV-8 but without reported clinical signs.

## Materials and methods

### General study design and available data

Increases in the frequency of early calving were quantified for cows in herds reported as cases, based on the detection of clinical signs, as well as for cows in herds located in the 2007 outbreak area in France but not reported as cases. These frequencies were compared to frequencies of early calving in cows in herds not exposed to the virus.

Information about notification of herds positive to BTV-8 during 2007 (from July to December) was obtained from the official veterinary surveillance system. Among herds reported during 2007, only herds with a confirmed detection reported after clinical suspicion were included. Information about BTV-8 exposure was available at the herd level only. Thus, all the cows of a herd were considered exposed if at least one animal with clinical signs had tested positive for BTV-8. The proportion of infected animals in reported herds was unknown. Herds selected for this study will be referred to as “case herds”. Case herds located in 23 departments (geographical and administrative unit) were selected (n = 8494), corresponding to 96% of all cattle herds reported with clinical signs during the 2007 epizootic (the non-selected case herds were located in departments with very low prevalence of BTV-positive herds). Cattle herds that were not reported during 2007 and located in these departments were considered as likely to have been exposed to BTV-8 (n = 46 569) (Figure 1). To see a description of the BTV notifications in France in 2007, see Durand et al.

<figure><img src="image_4.png" /><figcaption>**Figure 1** Geographical location of cattle herds in the 23 departments located in the outbreak area. A 8494 case herds reported after clinical signs of Bluetongue virus serotype 8 (BTV-8) in 2007, B 46 569 non-reported herds (e.g., herds with an interpolated date of BTV-8 exposure), 2007; France.</figcaption></figure>
[17]. The geographical coordinates of herds were available at the municipality level.

The reproductive data were obtained from dairy herds where artificial insemination was used between October 2005 and September 2009. For each cow, obtained data were dates of AIs, date of culling (if it had occurred during the study period), calving date and data used to adjust for factors known to influence the length of gestation: cow and bull breeds and parity.

A single date of exposure to BTV-8 was estimated from recorded data for all the cows of a case herd. For herds located in the 2007 outbreak area but not reported as cases, a date of exposure was interpolated. This interpolation was based on the spatio-temporal dynamics of detection of confirmed case herds that reported clinical signs in 2007. Increase in early calving frequency was quantified in both case herds and herds located in the 2007 outbreak area but not reported as cases.

## Estimated dates of exposure for reported case herds

For each case herd, available data included the date at which clinical signs of disease were first suspected and the date at which disease was confirmed via diagnostic tests. Some of the cows were exposed earlier than the date of first suspicion since the BTV incubation period is between 2 and 18 days in cattle [19]. But not all the cows from a herd would be sick on the same day and clinical signs were considered to be associated with viremia

in sick animals and could therefore serve as a source of virus for midges. Therefore, the estimated date of exposure for reported herds was defined as the recorded date of suspicion which corresponded to the first detection of clinical signs in the herd. The same date of exposure was assigned to all cows in a herd. For 6.1% of the case herds, the date of clinical suspicion was missing but the date of confirmation by a diagnostic test was known. In order to assign a date of suspicion, an imputation procedure based on the distribution of the time intervals between dates of suspicion and dates of confirmation was applied (values selected at random around a median interval of 4 days). Moreover, 181 case herds that had a non-plausible time interval between dates of suspicion and confirmation (interval >30 days or date of suspicion posterior to the date of confirmation) were excluded.

## Interpolation of dates of exposure for herds located in exposed areas but not reported as cases

A date of exposure to BTV-8 was interpolated for each herd located in the 2007 outbreak area but not reported as case (Figure 2). Kriging, a geostatistical interpolation method, was used to estimate a date of exposure for these herds. Details on the method and on the steps followed can be found in Nusinovici et al. [18]. Dates were converted to numbers of days since the first case herd reported in 2007. Kriging uses a data sample (case herds) to predict values at unsampled locations (herd located in

<figure><img src="image_5.png" /><figcaption>**Figure 2** Kriging map of the dates of exposure to Bluetongue virus serotype 8. The dates were expressed as a number of days since the first clinical case herd during the 2007 epizootic in France (31st July 2007), and location of reported case herds (black crosses). The hatched areas correspond to regions with no data (from [15]).</figcaption></figure>
the 2007 outbreak area but not reported). All the cattle case herds (dairy and beef) were included as data samples because they could all have played a role in the virus spread and give a more precise picture of the epizootic wave diffusion than dairy herds alone. Kriging is based on assumptions regarding the form of the trend of the sample data, its variance and spatial correlation. The first step consisted in modelling the spatial correlation of the data. Two models were compared using a cross-validation process with observed data to determine each model's goodness of fit and to compare their predictions. The final model was based on a Gaussian spatial component for filtering the random local component. Finally, to account for the non-stationarity of the BTV-8 spreading process, the gradient of the viral diffusion was also included in the model by the use of universal Kriging in place of ordinary Kriging.

## Selection of unexposed herds and cows

To test the robustness of results and to limit possible selection biases, two references unexposed populations were selected. The first reference population was composed of cows belonging to herds exposed in 2007 that had calved in 2005 and 2006 (before virus introduction in the area). The second was composed of cows belonging to herds situated in areas unexposed to BTV 1 and 8 in 2007, Britany and the Southwest of France (no other serotype was detected in continental France that year). These areas were chosen based on the geographical distribution of detected case herds in France in 2007 [17].

The month of calving is a known variation factor of gestation length [14, 20]. Cows in unexposed herds were therefore selected according to their calving date so that the cows in exposed (reported or not) and unexposed populations calved during the same period of the year.

## Definition of early calving and data selection

The last AI before calving was considered to be the insemination that led to conception. The gestation length was calculated as the interval between that AI and the date of calving. Early calving was defined as calving occurring between expected percentile 1 (LB for lower bound) and percentile 25 (HB for higher bound) of the normal gestation length. Expected percentiles were calculated with the Gaussian distributions of normal gestation lengths estimated using a sample of gestations between 260 and 320 days over a period without any major epidemics (from 2003 to 2005) [15]. Breed and parity influence the length of pregnancies. Therefore, LB and HB were calculated for each breed and each parity group (Additional file 1). Gestation lengths shorter than P1 (corresponding to 4.3% of the total number of gestations) were excluded because BTV-8 increases the risk of

abortion which can be evidenced by very short gestation lengths (e.g., between 175 and 270 days [14]).

Data from the most common breeds in France, i.e., Holstein, Montbéliarde and Normande cows, were selected, representing 94.5% of the dataset. Cows inseminated with the semen of Holstein, Montbéliarde and Normande bulls were selected assuming that the gestation length could differ according to the bull's breed (due to an effect of the calf weight). Cows with extreme or aberrant reproduction data were excluded: calving-first AI interval <35 or >180 days (corresponding to 1.8% of the total number of cows), assumed conceiving AI-calving interval >320 days (corresponding to 3.8% of the total number of cows). Herds under special management were also excluded (corresponding to 5.9% of the total number of herds): special demographic structure (% of primiparous cows within herds <10 or >75%), late breeding (calving-AI1 interval >120 days >75%) and herds likely to have used a bull (%AI2/AI1 < 10%). After these selections, populations of case herds, herds located in outbreak areas but not reported as cases and the two reference populations were composed of 2616, 4008, 11038 (unexposed herds in 2007) and 8015 (herds in 2005/2006 located in the 2007 outbreak area) herds respectively.

## Selection and classification of cows according to the date of exposure within gestation

The quantification of the BTV-8 exposure effect on the frequency of early calving took into account the date of exposure (observed or interpolated) within gestation. Only pregnancies exposed at the end of gestation were selected, from 63 days prior to LB to HB (exposures at HB were not included). Within this time interval, 8 categories of exposure were considered, corresponding to 9-day time intervals (except the last category LB-HB of 8-day length). The corresponding intervals in days of pregnancies for Holstein's multiparous cows are detailed in the Figure 3 (for all intervals, see Additional file 1).

## Statistical models

The relationship between exposure to BTV-8 and occurrence of early calving was assessed using multivariable statistical models. To assess the frequency of early calving, a mixed-logistic regression model was used. To account for factors likely to influence the frequency of early calving, the association was adjusted for parity and month of calving as described by the following equation:

$$Y_{ij} \sim \text{Bernoulli}(p_{ij})$$

$$\log \left( \frac{p_{ij}}{1 - p_{ij}} \right) = \alpha + \beta_1 \text{EXP}_{ij} + \beta_2 \text{PAR}_{ij} + \beta_3 \text{MTH}_{ij} + \text{RANEF}_j$$

$$\text{RANEF}_j \sim \text{Normal} (0, \sigma_j^2)$$
<figure><img src="image_6.png" /><figcaption>**Figure 3** Categories of exposure during gestation depending on the stage of pregnancy at the date of exposure (observed or inter-polated). The values considered corresponded to multiparous Holstein cows (time intervals for other breeds and parity in Additional file 1). The histogram represents the distribution of gestation lengths for multiparous Holstein cows, with the gray area corresponding to early calving (defined as the interval between the expected percentile 1 and 25 of the Gaussian distributions of normal gestation lengths).</figcaption></figure>

where the outcome Yij was a binary variable denoting
the occurrence of an early calving in cow i from herd j,
with a probability of occurrence pij; α was the intercept;
EXPij was the exposure category. The exposure variable
was categorized in 17 classes, depending on the time of
exposure during pregnancy: eight classes corresponded to
time intervals considered for cows with an observed date
of exposure, eight classes corresponded to cows with an
interpolated date of exposure and one class corresponded
to the unexposed population (i.e., the reference popu-
lation); PARj was the parity (4 classes: 1, 2, 3 and ≥4);
MTHij was the month of calving (8 classes from August to
March) and RANEFj was a random variable correspond-
ing to the herd number. The random variable allowed
adjusting for herd clustering. As odds-ratio (OR) overes-
timate the true relative risk (RR) when the incidence of
the study event is high (frequency of early calving), OR
were converted into RR using Beauveau and Fourichon's
method [21]. The effect in percentage points of early calv-
ing frequency were calculated from RR estimates. Statisti-
cal analyses were performed using R software [22].

Results

Unadjusted frequencies of early calving

Raw frequencies of early calving in the two reference pop-
ulations were comparable (Table 1). Without adjustment

for parity and month of calving, frequencies were highest
both for cows in case herds and cows in herds located in
the outbreak area but not reported as cases.

Except during the 2007 outbreak period, the frequency
of early calving was quite stable over the years with a high
seasonal effect (Figure 4). Frequencies of early calving of
cows located in Brittany and Southwest of France showed
similar patterns. When selecting only cows exposed
between LB-63 days and HB, the increased frequency of
early calving was much more pronounced (Figure 4, red
curve). In the latter population, the highest increases
were observed for calvings occurring in September and
October 2007.

Increased of early calving frequencies associated with BTV
infection

Exposure to BTV-8 was associated with an increased
frequency of early calving in both case herds and herds
located in the outbreak area but not reported as cases
(Figure 5). This increase was higher for exposures at the
latest stage of pregnancy. The increase was higher for cows
in case herds compared to cows in herds located in the
outbreak area but not reported as cases. Results were simi-
lar whatever the reference population considered (Table 2).
Exposure between day 270 and 278 yielded the highest
OR (2.09 [1.91–2.29]), which corresponds to a RR of 1.73
<table><caption>Table 1 Number of herds, number of gestations and raw frequency of early calvings according to exposure statuses of the herd to Bluetongue virus serotype 8 (BTV-8)</caption><thead><tr><th></th><th>Number of herds</th><th>Number of gestations</th><th>Frequency of early calving (%)</th></tr></thead><tbody><tr><td>Reported case herds with clinical signs in 2007</td><td>2616</td><td>17 139<sup>a</sup></td><td>24.4</td></tr><tr><td>Non-reported herds located in the 2007 outbreak area</td><td>4008</td><td>20 139<sup>a</sup></td><td>21.8</td></tr><tr><td>Unexposed herds in 2007</td><td>11 038</td><td>148 857<sup>b</sup></td><td>19.5</td></tr><tr><td>Herds in 2005/2006 that were located in the 2007 outbreak area</td><td>8015</td><td>189 359</td><td>19.4</td></tr></tbody></table>

Early calving was defined as a calving after a pregnancy length between expected percentiles 1 and 25 of the Gaussian distribution of normal gestation length. 2005–2006: unexposed population, 2007: exposed population, France

<table><caption>Table 2 Increase in frequencies of early calving, expressed as odds-ratios (OR) for cows in herds reported as cases of Bluetongue virus serotype 8 during the 2007 outbreak and cows in herds located in the 2007 outbreak area but not reported as cases</caption><thead><tr><th rowspan="3">Type of herd</th><th rowspan="3">Category of exposure<sup>a</sup></th><th colspan="3">Reference population</th><th colspan="3">Cows in unexposed herds in 2007</th></tr><tr><th colspan="3">Cows that calved in 2005 or 2006 belonging to herds located in the 2007 outbreak area</th><th colspan="3">Cows in unexposed herds in 2007</th></tr><tr><th>Odds ratio</th><th>95% CI</th><th>p value</th><th>Odds ratio</th><th>95% CI</th><th>p value</th></tr></thead><tbody><tr><td rowspan="10">Cows in case herds reported in 2007</td><td>(207–216)</td><td>1.06</td><td>0.94–1.21</td><td>0.34</td><td>1.10</td><td>0.97–1.26</td><td>0.14</td></tr><tr><td>(216–225)</td><td>1.18</td><td>1.05–1.33</td><td>0.01</td><td>1.20</td><td>1.06–1.35</td><td>0.004</td></tr><tr><td>(225–234)</td><td>1.25</td><td>1.11–1.40</td><td>&lt;0.001</td><td>1.25</td><td>1.11–1.40</td><td>&lt;0.001</td></tr><tr><td>(234–243)</td><td>1.37</td><td>1.23–1.52</td><td>&lt;0.001</td><td>1.37</td><td>1.23–1.53</td><td>&lt;0.001</td></tr><tr><td>(243–252)</td><td>1.51</td><td>1.36–1.66</td><td>&lt;0.001</td><td>1.51</td><td>1.36–1.67</td><td>&lt;0.001</td></tr><tr><td>(252–261)</td><td>1.73</td><td>1.57–1.90</td><td>&lt;0.001</td><td>1.70</td><td>1.54–1.88</td><td>&lt;0.001</td></tr><tr><td>(261–270)</td><td>1.81</td><td>1.65–1.98</td><td>&lt;0.001</td><td>1.74</td><td>1.58–1.92</td><td>&lt;0.001</td></tr><tr><td>(270–278)</td><td>2.09</td><td>1.91–2.29</td><td>&lt;0.001</td><td>1.97</td><td>1.79–2.16</td><td>&lt;0.001</td></tr><tr><td>(207–216)</td><td>1.18</td><td>1.05–1.33</td><td>&lt;0.001</td><td>1.30</td><td>1.15–1.46</td><td>&lt;0.001</td></tr><tr><td>(216–225)</td><td>1.27</td><td>1.14–1.42</td><td>&lt;0.001</td><td>1.45</td><td>1.30–1.62</td><td>&lt;0.001</td></tr><tr><td rowspan="7">Cows in herds located in the 2007 outbreak area but not reported as cases</td><td>(225–234)</td><td>1.25</td><td>1.12–1.38</td><td>&lt;0.001</td><td>1.43</td><td>1.28–1.60</td><td>&lt;0.001</td></tr><tr><td>(234–243)</td><td>1.30</td><td>1.18–1.44</td><td>&lt;0.001</td><td>1.49</td><td>1.35–1.66</td><td>&lt;0.001</td></tr><tr><td>(243–252)</td><td>1.25</td><td>1.14–1.38</td><td>&lt;0.001</td><td>1.39</td><td>1.25–1.53</td><td>&lt;0.001</td></tr><tr><td>(252–261)</td><td>1.23</td><td>1.12–1.35</td><td>&lt;0.001</td><td>1.31</td><td>1.19–1.45</td><td>&lt;0.001</td></tr><tr><td>(261–270)</td><td>1.42</td><td>1.30–1.56</td><td>&lt;0.001</td><td>1.47</td><td>1.34–1.61</td><td>&lt;0.001</td></tr><tr><td>(270–278)</td><td>1.34</td><td>1.22–1.47</td><td>&lt;0.001</td><td>1.40</td><td>1.27–1.54</td><td>&lt;0.001</td></tr></tbody></table>

In both populations, cows with dates of exposure (reported or interpolated) within the last 10 weeks of gestation were selected. The different OR values correspond to each category of exposure, defined according to the dates of exposure during pregnancy. Two reference populations were considered: cows that calved in 2005 or 2006 belonging to herds located in the 2007 outbreak area, and cows in unexposed herds on 2007 (Brittany and Southwest of France)

[1.62–1.83]. Consequently, for 100 multiparous Holstein cows exposed between 270 and of 278 days of pregnancy, the frequency of early calving would be on average 33.6% [31.4–35.5], compared to 19.4% without exposure (corresponding to an average increase of 14.2%). The effects of parity and month of calving were significant (Table 3).

Discussion

During the outbreak of 2007 in France, infection by BTV-8 in late gestation resulted in gestations that,

although they were in the normal range of gestation lengths, were a few days shorter than expected. Increases in the frequency of early calvings were detected for both cows in herds reported as cases in 2007 and cows in herds located in the outbreak area but not reported as cases. For cows in case herds, increases in the frequency of early calving reached high values (OR > 2). To our knowledge, this is the first study quantifying the association of disease and triggering of an early calving. Such an effect would have been very difficult to detect with an

<footnote>a Gestations for which the date of herd exposure (observed or interpolated) occurred during the last 10 weeks of pregnancy (see "Materials and methods" section for more details).</footnote>

<footnote>b Gestations with calving occurring during the same period in the year than herds located in the outbreak area.</footnote>

<footnote>a The values considered corresponded to multiparous Holstein cows (for time intervals for other breeds and parity, see Additional file 1).</footnote>
<figure><img src="image_7.png" /><figcaption>**Figure 4** Raw frequencies of early calving per calving month from October 2005 to September 2009. (A) 34 248 cows in 2616 reported case herds with clinical signs in 2007; (B) 209 231 cows in 9540 unexposed herds located in Brittany; (C) 26 701 cows in 1498 unexposed herds located in the Southwest of France. B and C composed the unexposed population BTV-8 in 2007; (D) 17 139 cows in 2616 herds that reported clinical signs in 2007 exposed during the last 10 weeks of pregnancies. Population D is a subset of the population A. (E) 20 139 cows in 4008 herds for which a date of exposure was interpolated during the last 10 weeks of pregnancies. These herds were non-reported during the 2007 outbreak but located in the outbreak area. The two red vertical lines correspond to the 2007 outbreak period in France.</figcaption></figure>

experimental study for several reasons. Firstly, the event of interest lies within the normal/physiological range of variation. Indeed, the early calvings that were investigated in the present study would not qualify as abortions or any other pathological event to be recorded. Secondly, the time window between exposure to the virus and the observed increase in the frequency of early calvings was very narrow and the magnitude of the effect increased as the expected time to calving decreased. Thirdly, given the high influence of the season (in the same range of magnitude as the effect of exposure), an unbiased estimation of the effect would have required to stratify the observations by season, which would require a very large number of cows. The present study was carried out on a very large cattle population allowing a high statistical power. Moreover, the multivariable approach allowed controlling for known factors that influence gestation length. Finally, two reference populations were considered to limit possible bias because of variations in livestock management either between regions or over time. The results were similar regardless of the reference population. This

robustness allows concluding that there was no major bias regarding variations in livestock management that would have not been accounted in the models.

The estimated magnitude of increased early calving frequency corresponds to mean effects in a herd under natural conditions of exposure in late gestation. In these herds, the proportion of infected cows was unknown. Therefore, these values likely underestimate the effect of infection for a single infected cow. BTV-8 is known to affect reproductive performances (fertility and abortion) due to infection in pregnant cows (with a tropism of the virus for the genital tract) possibly followed by the infection of the fetus. Interestingly, the effect of BTV-8 on early calving, for which no physiopathological mechanism is yet known, is higher in terms of magnitude of effect at a population level. However, lengths of pregnancies were comprised in a normal range and were expected to result in the birth of live calves. It can therefore be expected that the economic impact at the herd level may be low.

The highest effect of BTV-8 infection was found for exposures occurring during the latest stage of
<figure><img src="image_8.png" /><figcaption>**Figure 5** Increase in frequencies of early calving. The increases were expressed as odds-ratios (OR), for cows in herds reported as cases of Blue-tongue virus serotype 8 during the 2007 outbreak (black curve) and cows in herds located in the 2007 outbreak area but not reported as cases (gray curve). In both populations, cows with dates of exposure (reported or interpolated) within the last 10 weeks of gestation were selected. The different OR values correspond to each category of exposure, defined according to the dates of exposure during pregnancy. The reference population considered was composed of cows that calved in 2005 or 2006 from herds located in the 2007 outbreak area. Vertical bars around ORs correspond to their 95% confidence intervals. The histogram represents the distribution of gestation lengths for multiparous Holstein cows, with the gray area corresponding to early calving (defined as the interval between the expected percentile 1 and 25 of the Gaussian distributions of normal gestation lengths).</figcaption></figure>

<table><caption>**Table 3** Effect of adjustment variables on the frequency of early calving estimated with mixed logistic models</caption><thead><tr><th>Variable and class</th><th>Number of pregnancies</th><th>Odds Ratio</th><th>95% CI</th><th>p value</th></tr></thead><tbody><tr><td>Parity</td><td></td><td></td><td></td><td></td></tr><tr><td>1</td><td>96 048</td><td>1</td><td>NA</td><td></td></tr><tr><td>2</td><td>62 470</td><td>1.18</td><td>1.15-1.21</td><td>&lt;0.001</td></tr><tr><td>3</td><td>37 036</td><td>1.11</td><td>1.08-1.15</td><td>&lt;0.001</td></tr><tr><td>4 or more</td><td>37 293</td><td>1.10</td><td>1.07-1.14</td><td>&lt;0.001</td></tr><tr><td>Month of calving</td><td></td><td></td><td></td><td></td></tr><tr><td>August</td><td>17 525</td><td>1</td><td>NA</td><td>&lt;0.001</td></tr><tr><td>September</td><td>24 525</td><td>0.59</td><td>0.57-0.62</td><td>&lt;0.001</td></tr><tr><td>October</td><td>41 507</td><td>0.63</td><td>0.61-0.66</td><td>&lt;0.001</td></tr><tr><td>November</td><td>52 305</td><td>0.49</td><td>0.47-0.51</td><td>&lt;0.001</td></tr><tr><td>December</td><td>44 618</td><td>0.52</td><td>0.5-0.55</td><td>&lt;0.001</td></tr><tr><td>January</td><td>26 490</td><td>0.51</td><td>0.48-0.53</td><td>&lt;0.001</td></tr><tr><td>February</td><td>14 341</td><td>0.53</td><td>0.5-0.56</td><td>&lt;0.001</td></tr><tr><td>March</td><td>11 536</td><td>0.57</td><td>0.54-0.61</td><td>&lt;0.001</td></tr></tbody></table>

The unexposed population considered in this analysis corresponded to cows belonging to herds exposed in 2007 that calved in 2005 and 2006 (17 139 cows in 2616 case herds, 20 139 cows in 4008 herds located in the 2007 outbreak area but not reported and 189 359 cows in 8015 reference herds).

NA not applicable.

pregnancy. The impact of exposure to BTV on the mechanism triggering calving could therefore be higher in the very last stage of gestation compared to earlier exposures. This finding is in accordance with Marceau et al. [15] that indicated that delay between infection and the occurrence of an early calving should be short. Consequently, cows infected at the latest stage of pregnancy calve a few days earlier than expected i.e., that the infection by BTV-8 triggers calving. Increases in early calving were also observed for exposures from 200 days of gestation. Because the dates of infections of cows within these herds are unknown, it is difficult to conclude about effects associated with these early exposures. On the one hand, this could be related to the time needed for the virus to spread within herds: due to this infection dynamic, some cows within infected herds could have been infected weeks after the date of first infection. On the other hand, this could be due a delayed effect of infection, but no results are available to support or infirm this hypothesis.

An increase in the frequency of early calving was quantified in herds located in the outbreak area but not
reported as cases. This effect corresponded to approximately half of the increase quantified in herds reported with clinical signs. A similar proportion was found for the quantification of fertility decrease in herds located in the outbreak area but not reported compared to the effect of BTV-8 in case herds [18]. Because of the high under-reporting rate, some of the herds located in the outbreak area but not reported have probably been infected during the outbreak. Infected herds could have not been reported as cases either because no animal showed any clinical signs (asymptomatic BTV-8 infections in cattle were quite frequent [4, 23]) or because the clinical signs were not specific and thus not attributed by farmers to BTV-8 infection.

It has been demonstrated before that an indicator based on short gestation was able to detect the 2007 BTV outbreak at an early stage with a large number of abnormal elevations in infected areas [15]. The present study strengthens the interest of considering such an indicator for the following reasons. Firstly, the magnitude of the increased early calving frequency associated with BTV-8 infection was high. It is known that the ability to detect a disease using syndromic surveillance is highly related to the magnitude of its effect on the indicator used [24]. The related biological mechanism behind early calving is unknown. This event is different from an abortion because it occurs within the normal range of gestation lengths and is expected to result in the birth of live calves. It would have been interesting to investigate whether other diseases are associated with such increases frequency of early calving. If so, that would confirm the interest of considering early calving as an indicator for syndromic surveillance. Secondly, increased frequency of early calvings was detected in areas with infected herds that were not reported through a clinical surveillance system. Infec-tions with very moderate clinical signs or clinical signs not specific to BTV-8 were frequent in cattle during the epizootic in North-west Europe [4, 6]. That suggests that early calving could be associated with non-specific biological phenomena such as fever in infected cows, as suggested by [15]. An indicator associated with non-specific symptoms would allow detecting diseases that are particularly difficult to detect through clinical surveillance. Moreover, such symptoms might be associ-ated with a wide range of infectious diseases. Finally, on the contrary to other effects of BTV-8 on reproductive performance (increase embryonic mortality and abor-tion), the increase in early calving should occur quickly after infection. Using this indicator for syndromic sur-veillance could therefore contribute to the early detec-tion of pathogens with such an effect on reproductive performance.

**Additional file**

**Additional file 1. Time intervals considered for each category of exposure depending on breed and parity.** This additional file contained the time intervals considered for each category of exposure depending on breed and parity. These time intervals varied depend-ing on breed and parity because these factors influence the length of pregnancies.

**Competing interests**

The authors declare that they have no competing interests.

**Authors' contributions**

SN carried out the statistical analysis and drafted the manuscript. AM and CF participated in the design of the study and helped to draft the manuscript. All authors read and approved the final manuscript.

**Acknowledgements**

We would like to thank the INRA Centre de Traitement de l'Information Géné-tique (CTIG) for providing the reproductive data and the Direction Générale de l'Alimentation (DGAL) of the French Ministry of Agriculture for providing the BTV notification data.

**Author details**

¹ INRA, UMR1300 Biology, Epidemiology and Risk Analysis in Animal Health, CS 40706, F-44307 Nantes, France. ² LUNAM Université, Oniris, UMR1300 Biology, Epidemiology and Risk Analysis in Animal Health, CS 40706, F-44307 Nantes, France.

Received: 13 April 2015 Accepted: 17 December 2015

Published online: 13 January 2016

**References**

1. Jones KE, Patel NG, Levy MA, Storeygard A, Balk D, Gittleman JL, Daszak P (2008) Global trends in emerging infectious diseases. Nature 451:990–993

2. Caminade C, Medlock JM, Ducheyne E, McIntyre KM, Leach S, Baylis M, Morse AP (2012) Suitability of European climate for the Asian tiger mosquito Aedes albopictus: recent trends and future scenarios. J R Soc Interface 9:2708–2717

3. Guis H, Caminade C, Calvete C, Morse AP, Tran A, Baylis M (2011) Modeling the effects of past and future climate on the risk of bluetongue emergence in Europe. J R Soc Interface 9:339–350

4. Elbers AR, Backx A, Meroc E, Gerbier G, Staubach C, Hendrickx G, van der Spek A, Mintiens K (2008) Field observations during the bluetongue sero-type 8 epidemic in 2006. I. Detection of first outbreaks and clinical signs in sheep and cattle in Belgium, France and the Netherlands. Prev Vet Med 87:21–30

5. Elbers AR, van der Spek AN, van Rijn PA (2009) Epidemiologic characteristics of bluetongue virus serotype 8 laboratory-confirmed outbreaks in The Netherlands in 2007 and a comparison with the situation in 2006. Prev Vet Med 92:1–8

6. Le Gal MC, Dufour B, Geoffroy E, Zanella G, Moutou F, Millemann Y, Rieffel JN, Pouilly F (2008) Bluetongue virus serotype 8 in the Ardennes in 2007. Vet Rec 163:668

7. Elbers AR, Backx A, Mintiens K, Gerbier G, Staubach C, Hendrickx G, van der Spek A (2008) Field observations during the Bluetongue serotype 8 epidemic in 2006. II. Morbidity and mortality rate, case fatality and clinical recovery in sheep and cattle in the Netherlands. Prev Vet Med 87:31–40

8. Perrin JB, Ducrot C, Vinard JL, Morignat E, Gauffier A, Calavas D, Hendrikx P (2010) Using the National Cattle Register to estimate the excess mortality during an epidemic: application to an outbreak of Bluetongue serotype 8. Epidemics 2:207–214

9. Santman-Berends IM, van Schaik G, Bartels C J, Stegeman JA, Vellema P (2011) Mortality attributable to bluetongue virus serotype 8 infection in Dutch dairy cows. Vet Microbiol 148:183–188
10. Santman-Berends IM, Hage JJ, Lam TJ, Sampimon OC, van Schaik G (2011) The effect of bluetongue virus serotype 8 on milk production and somatic cell count in Dutch dairy cows in 2008. J Dairy Sci 94:1347–1354

11. Nusinovici S, Souty C, Seegers H, Beauveau F, Fourichon C (2013) Decrease in milk yield associated with exposure to bluetongue virus serotype 8 in cattle herds. J Dairy Sci 96:877–888

12. Santman-Berends IM, Hage JJ, Rijn PA, Stegeman JA, Schaik GV (2010) Bluetongue virus serotype 8 (BTV-8) infection reduces fertility of Dutch dairy cattle and is vertically transmitted to offspring. Theriogenology 74:1377–1384

13. Nusinovici S, Seegers H, Joly A, Beauveau F, Fourichon C (2011) Quantification and at-risk period of decreased fertility associated with exposure to Bluetongue virus serotype 8 in naïve dairy cattle herds. J Dairy Sci 95:3008–3020

14. Nusinovici S, Seegers H, Joly A, Beauveau F, Fourichon C (2012) Increase in the occurrence of abortions associated with exposure to the Bluetongue virus serotype 8 in naïve dairy herds. Theriogenology 78:1140–1151

15. Marceau A, Madouasse A, Lehébel A, van Schaik G, Veldhuis A, Van der Stede Y, Fourichon C (2014) Can routinely recorded reproductive events be used as indicators of disease emergence in dairy cattle? An evaluation of 5 indicators during the emergence of bluetongue virus in France in 2007 and 2008. J Dairy Sci 97:6135–6150

16. MacLachlan N, Conley A, Kennedy P (2000) Bluetongue and equine viral arteritis viruses as models of virus-induced fetal injury and abortion. Anim Reprod Sci 60–61:643–651

17. Durand B, Zanella G, Biteau-Coroller F, Locatelli C, Baurier F, Simon C, Le Dréan E, Delaval J, Prengère E, Beauté V, Guis H (2010) Anatomy of blue-tongue virus serotype 8 epizootic wave, France, 2007–2008. Emerg Infect Dis 16:1861–1868

18. Nusinovici S, Monestiez P, Seegers H, Beauveau F, Fourichon C (2014) Using animal performance data to evidence the under-reporting of case herds during an epizootic: application to an outbreak of bluetongue in cattle. PLoS One 9:e100137

19. Thiry E, Zimmer JY, Haubruge E (2008) In: Saegerman C, Reviriego-Gordejo F, Pastoret PP. Fièvre catarrhale ovine en Europe du nord. 13-25

20. Norman HD, Wright JR, Kuhn MT, Hubbard SM, Cole JB, VanRaden PM (2009) Genetic and environmental factors that affect gestation length in dairy cattle. J Dairy Sci 92:2259–2269

21. Beauveau F, Fourichon C (1998) Estimating relative risk of disease from outputs of logistic regression when the disease is not rare. Prev Vet Med 36:243–256

22. R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2012 (http://www.R-project.org/)

23. Dal Pozzo F, Saegerman C, Thiry E (2009) Bovine infection with bluetongue virus with special emphasis on European serotype 8. Vet J 182:142–151

24. Reis BY, Mandl KD (2003) Time series modeling for syndromic surveillance. BMC Med Inform Decis Mak 3:2

Submit your next manuscript to BioMed Central and we will help you at every step:

* We accept pre-submission inquiries

* Our selector tool helps you to find the most relevant journal

* We provide round the clock customer support

* Convenient online submission

* Thorough peer review

* Inclusion in PubMed and all major indexing services

* Maximum visibility for your research

Submit your manuscript at
www.biomedcentral.com/submit

BioMed Central