Skip to main content

Modelling effectiveness of herd level vaccination against Q fever in dairy cattle


Q fever is a worldwide zoonosis caused by the bacterium Coxiella burnetii. The control of this infection in cattle is crucial: infected ruminants can indeed encounter reproductive disorders and represent the most important source of human infection. In the field, vaccination is currently advised in infected herds but the comparative effectiveness of different vaccination protocols has never been explored: the duration of the vaccination programme and the category of animals to be vaccinated have to be determined. Our objective was to compare, by simulation, the effectiveness over 10 years of three different vaccination strategies in a recently infected dairy cattle herd.

A stochastic individual-based epidemic model coupled with a model of herd demography was developed to simulate three temporal outputs (shedder prevalence, environmental bacterial load and number of abortions) and to calculate the extinction rate of the infection. For all strategies, the temporal outputs were predicted to strongly decrease with time at least in the first years of vaccination. However, vaccinating only three years was predicted inadequate to stabilize these dynamic outputs at a low level. Vaccination of both cows and heifers was predicted as being slightly more effective than vaccinating heifers only. Although the simulated extinction rate of the infection was high for both scenarios, the outputs decreased slower when only heifers were vaccinated.

Our findings shed new light on vaccination effectiveness related to Q fever. Moreover, the model can be further modified for simulating and assessing various Q fever control strategies such as environmental and hygienic measures.


Q fever is a zoonotic disease caused by Coxiella burnetii, a bacterium found worldwide in a wide range of animals. In ruminants, the infection may cause abortions, infertility, metritis or chronic mastitis [14], which can lead to non negligible economic losses for the infected herds. Furthermore, since 2007, Q fever has become an important public health problem in several parts of Europe [57]. Although Q fever is asymptomatic in 60% of human cases, it can lead to acute or chronic infections and cause flu-like syndrome, hepatitis, pneumonia, endocarditis or abortions [8]. In the Netherlands, where a steep increase in the number of human cases was observed in 2007, 2008, and 2009, a link has been established between some human cases and farms of small ruminants where abortions due to Q fever were detected [9]. Ruminants are indeed recognized as the main source of human infection [10, 11]. Infected animals shed large quantities of bacteria into the environment through faeces, vaginal mucus, urine, milk and especially parturition products [4, 12, 13]. C. burnetii survives very well in the environment and contaminates aerosols and dust [14]. These infected particles are the main route of infection for both animals and humans. Due to its importance in both animal and public health, the control of this infection is crucial. Therefore, any control measure leading to a decrease in the prevalence of shedders and in the environmental bacterial load seems a key point to limit both the spread of the infection in ruminants and the zoonotic risk.

Nowadays, in infected cattle herds in France, control measures against Q fever consist of environmental measures such as destruction of placentas or disinfection of birth locations, antibiotic treatment like oxytetracycline injections during the last month of gestation, and vaccination [15]. According to Rodolakis et al. [15], vaccination would be an efficient tool to control the disease. A phase I vaccine was indeed shown to prevent abortions and dramatically reduce the frequency of bacterial shedding in the milk, vaginal mucus and faeces [16]. Besides, according to Guatteo et al. [17], susceptible cattle that were vaccinated when non pregnant had a five times lower probability to become a shedder than an animal receiving placebo.

Thus, in the field, vaccination is often recommended in infected herds after the occurrence of abortions due to Q fever. However, the studies assessing the vaccination efficacy in ruminants were carried out in experimental conditions or for a limited period of time and they evaluated the effect of the vaccine mostly at the individual level. Therefore, it is difficult to extrapolate these results to the case of a whole herd vaccination over several years. Furthermore, different vaccination strategies can be implemented: the duration of the vaccination programme as well as the category of vaccinated animals (e.g., the whole herd or the heifers only) have to be determined. To assess the long run effectiveness of these different strategies in reducing the infection prevalence or the environmental bacterial load, field studies are not optimal: no reference situation (without control strategy) is generally available, and long-term observations must be performed, making these studies very costly and even unfeasible. Modelling is therefore a convenient approach as it provides means to compare the effectiveness of different potential management strategies [18].

The objective of this model study was to assess the comparative effectiveness of several vaccination strategies against C. burnetii in an already infected dairy cattle herd. The criteria considered for efficacy evaluation were changes in the prevalence of shedders, the environmental bacterial load, the number of abortions, as well as the extinction rate of infection.

Materials and methods

A model representing the C. burnetii infection dynamics in a standard French dairy cattle herd and different vaccination strategies was elaborated based on a previous variant model not including interventions. First of all, the epidemic model representing the natural course of infection (i.e. without any control strategy) will be briefly described, then the inclusion of vaccination will be presented and finally, the different vaccination scenarios that we tested will be explained in detail. Some elements of a sensitivity analysis performed to identify the most influent parameters and structural characteristics of both the model variants without and with vaccination are given in Additional file 1.

General description of the epidemic model of the natural course of infection

The model represents the spread of the bacterium in a dairy herd of lactating and dry cows (diagram flow in Figure 1 and parameters in Additional file 1: Table S1). It is a stochastic individual-based model in discrete time with a time step of one week. The model also represents the heterogeneity of shedding described for C. burnetii infections [13, 19, 20]: in field data [20], differences in shedding routes, levels (i.e. quantities of bacteria shed) and duration were observed between animals. Some animals can also become chronically infected [12]. Therefore, in the model, three types of shedders with different shedding characteristics and epidemiological behaviours were distinguished: (i) shedders of type 1 which have the possibility to clear the infection (I 1 ), (ii) shedders of type 2 which are chronically infected (I 2 ), and (iii) shedders of type 3 which are chronically infected and always shed in milk, for a longer period of time and at higher levels than shedders of type 2 (I 3 ). By inhaling bacteria contained in the environment, a susceptible cow (S) can become infectious and start shedding. Either it manages to eliminate the bacterium and becomes apparently susceptible again (i.e. non-shedder without antibodies) or it becomes a chronically infected shedder (I 2 or I 3 ). Since the shedding is intermittent [12], these types of shedders can stop shedding and start again. We assume that a former I 2 or I 3 also has a chance to clear the infection. However, if this cow is infected again, it becomes a chronically infected animal again. Therefore, each cow is in one of the six mutually exclusive health states at a given time: S (non-shedder apparently susceptible cow), I 1 (shedder of type 1), I 2 (shedder of type 2), I 3 (shedder of type 3), C 1 (non-shedder but still infected cow), C 2 (non-shedder which was C 1 in the past but cleared the infection). Besides, sub-categories are defined for the shedder cows with respect to the shedding route. Thus, an I 1 or I 2 cow can shed in (1) milk only (denoted by or respectively), (2) vaginal mucus and/or faeces ( or respectively), or (3) milk and either vaginal mucus or faeces or both ( or respectively). In the same way, an sheds in milk only and an sheds in milk and vaginal mucus and/or faeces (by definition, an I 3 animal always sheds in milk and can not be in the state). The possible transitions between health states are represented in the top part of Figure 1. Shedders (I 1 , I 2 and I 3 ) fill the environment compartment (E) with bacteria. Three categories of shedding levels are represented: low, moderate and high level, corresponding respectively to a quantity of bacteria shed Qty of 1/3000, 1/30 and 1 unit of environment per week. The probability for a cow to shed at one of these levels depends on her infectious state (I 1 , I 2 or I 3 ) and on the shedding route (milk or mucus/faeces). This probability is governed by the probability distributions Q, described in Additional file 1: Table S1. Only a fraction ρ of the quantity of bacteria shed by a shedder is assumed to reach the environment. Thus, the quantity of bacteria arriving into the environment during a time step t is the sum, for all the shedding routes of all the shedders, of , the quantity of bacteria shed by the shedder i through the shedding route (route {milk, mucus/faeces}), times ρroutethe impact of this shedding on the environment (see Additional file 1: section 1.2. for details). The probability of infection or re-infection, p (transition from S to I 1 or from C 2 to I 2 ) is expressed at each time step as where E t is the quantity of bacteria in the herd environment at time t (one unit of E t corresponding to a probability of transition from S to I 1 of (1 - 1/e)). The mortality rate of C. burnetii in the environment, μ, includes the natural mortality of the bacterium and its removal in relation to the periodic cleaning of the cattle housing carried out by the farmer.

Figure 1
figure 1

Flow diagram describing the modelled spread of C. burnetii within a cattle herd. The health states are the following: S, non-shedder apparently susceptible cow, I 1 , shedder which still has the possibility to eliminate the bacterium and to become S again, I 2 , shedder which no longer has the possibility to become S again, I 3 , shedder which no longer has the possibility to become S again and sheds in milk in a persistent way, C 1 , non-shedder but still infected individual and C 2 , non-shedder which was C 1 in the past but eliminated the bacterium. The V e states (SV e , I 1 V e , I 2 V e , I 3 V e , C 1 V e and C 2 V e ) are defined in the same way as S, I 1 , I 2 , I 3 , C 1 and C 2 respectively, except that these animals have been vaccinated when susceptible and non pregnant and are then assumed "vaccinated in an effective way" (V e ). I and IV e cows are in the subcategory m if they shed in milk only, mf if they shed in vaginal mucus/faeces only and mmf if they shed in milk and vaginal mucus/faeces. E represents the environmental bacterial load and p, the probability of infection or reinfection for non V e individuals, is equal to . p v is the probability of infection or reinfection for V e individuals, which is a fraction of p. The other model parameters are presented in Additional file 1: Table S1. ε 1 , ε 2 , ε 3 , ε 1 V e , ε 2 V e and ε 3 V e are the quantities of bacteria shed during a time step by an individual I 1 , I 2 , I 3 , I 1 V e , I 2 V e and I 3 V e respectively and contaminating the environment. For a any shedder, ε represents the sum, for each shedding route, of the quantity of bacteria released, Qty, times ρ, its fraction reaching the herd environment.

Since abortions are the main clinical signs attributable to C. burnetii infections [14, 20], they are also represented in the model: a cow is assumed to have a risk to abort after her first infection, her reinfection or the resumption of shedding (i.e. after a transition from S to I 1 , from C 1 to I 2 or from C 2 to I 2 ). According to Arricau-Bouvery et al. [21], aborting females in cattle generally do not abort during the following gestations. Based on [21], it was assumed in our model that a cow can abort only once in her life. If the cow aborts in the first or second third of gestation, she sheds a moderate quantity of bacteria Qty through the mucus/faeces, whereas if the abortion occurs in the last third of gestation, a high quantity of bacteria is released through this shedding route.

The epidemic model was also coupled to a model of population dynamics in order to represent the gestation and lactation cycles of each cow. In short, for each cow the lactation number is represented, as well as the stage of lactation, the stage of gestation, the abortion history, the health state and the shedding characteristics (if the cow is shedding).

Representation of the vaccination

Based on Guatteo et al. [16], possibly due to downregulated Th1-type immune responses during pregnancy, we assumed that the vaccine is effective when applied to non pregnant uninfected individuals. Thus, in the epidemic model, non pregnant S and C 2 individuals become protected when vaccinated and move into the "vaccinated in an effective way" (V e ) state (bottom of Figure 1). Pregnant S and C 2 , as well as all I 1 , I 2 and C 1 are what we defined the uselessly vaccinated: the vaccine has no effect on the infection dynamics in these animals, and they keep moving between the states S, I 1 , I 2 , I 3 , C 1 and C 2 (top of Figure 1). Six additional health states are defined for the V e individuals. SV e and C 2 V e individuals can get infected and become I 1 V e or I 2 V e respectively with a decreased transition rate p v (equal to a fraction of p). Except for this difference between p and p v , the V e animals can evolve through the same health states with identical transition rates as the non V e animals.

Regarding the shedding levels, according to Guatteo et al. [16], the only quantified bacterium load of a V e shedder was lower than the lowest bacterium load of the placebo cows. Besides, in Rousset et al. [22], the bacterial loads in vaginal swabs were lower in vaccinated than in non vaccinated animals. Therefore, we assumed that no high level shedding is possible for V e animals and that the probability to shed at a low level is increased (expressed through probability distributions QV e in Additional file 1: Table S1). Finally, based on Arricau-Bouvery et al. [16], it was assumed that the V e cows cannot abort.

Vaccination scenarios

Scenario 1: vaccination over the whole simulation period (10 years)

At the start of the simulation, all the cows are vaccinated and all the heifers entering the herd of cows are assumed to be SV e (susceptible and vaccinated when non pregnant). In addition, all the animals are boosted every year: there is no loss of immunity and no possible transition from the V e states to the non V e states.

Scenario 2: vaccination for a limited period of time (3 years)

The assumptions are the same as those of scenario 1 except for the vaccination duration. Here, the herd is supposed to be vaccinated for three years. At the end of this 3 year period, two assumptions regarding the evolution of immunity were explored.

▪ Scenario 2A: immunity lasts for one year. One year after the end of the vaccination period, the V e animals lose their immunity and move to the non V e equivalent states (e.g. an I 2 V e cow becomes an I 2 cow).

▪ Scenario 2B: lifelong immunity. After the vaccination period, the V e animals do not lose their immunity and keep moving within the V e states until the end of their life.

Scenario 3: vaccination of the heifers only over the whole simulation period (10 years)

At the start of the simulation the cows are not vaccinated. They stay in the non V e states and progress through infection states. Only the heifers arriving thereafter are assumed to be vaccinated in an effective way. These animals are in the SV e state when entering the dairy herd. Afterwards, all the previously vaccinated animals are boosted every year: there is no loss of immunity and no transition from the V e states to the corresponding non V e states.

Negative control

No control programme was implemented and all the animals progress through the non V e states.

Parameters and initial conditions

The values of all epidemiologic parameters are displayed in Additional file 1: Table S1. Parameters m, q, r 1 , s and μ were fixed at their values estimated using data from five French chronically infected dairy cattle herds [23]; probability distributions of shedding related parameters, α, β, β calv , γ, γ calv , (governing the partition in subcategories according to the shedding route) and Q1, Q2, Q3, Q4 and Q5 (characterizing the shedding levels), were qualitatively calibrated to match field data, as well as pIp (the proportion of new chronically infected shedders becoming I3) and probav (the probability of abortion due to Q fever). The parameters governing the demography and herd management (Additional file 1: Table S2) were chosen to represent a standard French dairy cattle herd.

The transition rate p v was parameterized using the hazard ratio for a vaccinated initially susceptible animal to become a shedder in Guatteo et al. [17], which is equal to 0.21 with a 95% confidence interval of 0.05-0.90. Thus, we performed the simulations with p v = 0.21p. However, given the large confidence interval of [16], in scenario 1 two additional values were also tested (p v = 0.05p and p v = 0.90p) in order to determine the influence of this parameter value on the model output.

We simulated 100 repetitions of the introduction of a primiparous I 2 cow that has just calved into a fully susceptible herd of 50 cows to generate infected herds. We let the model run until three abortions had occurred during a period of 12 months to initiate reactive vaccination. This limit was motivated by the fact that testing for a large panel of abortive pathogens (including C. burnetii) is usually performed in France from the 3rd abortion within the calving period. Thus, we obtained 100 so called "initial herds", different from each other. Then, for each initial herd, the three vaccination scenarios and the negative control scenario were run once over a 10-year simulation period.

Outputs of the model

The mean prevalence of shedders, the number of abortions per herd per year and the environmental bacterial load were the model's dynamic outputs of interest. In addition, for each scenario, the rate of extinction over the 10 year simulation period was calculated as the ratio between the number of extinct trajectories and the total number of repetitions. The infection was assumed to be extinct when there were no more I, IV e , C 1 and C 1 V e cows in the herd at the end of the simulation time.


Description of the herds at the start of the vaccination strategy

At the start of simulations, the mean prevalence of shedders (over 100 initial herds) was equal to 28.2% (min: 0.0%, max: 59.2%) and the mean prevalence of milk shedders amounted to 12.9% (min: 0.0%, max: 32.7%). In a herd, 92.1% of the cows on average had been shedders for at least one time step (min: 58.0%, max: 100%). The mean environmental bacterial load was 0.34 units (min: 0.01, max: 1.17) and the herds consisted of 49.7 cows on average (min: 41, max: 58). The herds were rather recently infected: the mean time from the beginning of the infection was 58 weeks (min: 8, max: 154).

Influence of the vaccination scenarios on the temporal model outputs

If no control strategy was implemented, the mean prevalence of shedders, the mean environmental bacterial load and the mean number of abortions increased to a steady state around respectively 45% shedders, 0.8 unit of environment and four abortions per herd per year. On the contrary, for any vaccination scenario, all these outputs decreased with time at least for the first years of vaccination (Figure 2a, b, c). In scenario 1 (vaccination of heifers and cows during 10 years), the decrease covered the whole period. In scenario 3 (vaccination of heifers only for the whole simulation time), the decrease was much slower in the first years of vaccination than in scenario 1: the latter reached a mean prevalence of shedders of 5% and a mean environmental load of 0.05 respectively 9.2 and 9.9 months sooner than scenario 3. At the end of the vaccination period, the mean prevalence of shedders, environmental bacterial load and number of abortions were close to 0 in scenario 1 as well as in scenario 3. In scenario 2, there was an increase in the mean prevalence of shedders, the yearly number of abortions and the environmental bacterial load, after the vaccination ceased. For scenario 2A, this increase occurred immediately after the loss of immunity, whereas for scenario 2B (lifelong immunity), the increase was almost zero in the first year without vaccination and more progressive afterwards. Thus, the mean prevalence of shedders was around 9.5% for both scenarios 2A and 2B three years after the simulation started and increased to respectively 40.4% and 19.2% eight years after the simulation started. The mean number of yearly abortions increased from 0.4 and 0.54 the third year after the start of vaccination to 4.2 and 2.0 abortions per herd respectively during the eighth year after vaccination started.

Figure 2
figure 2

Temporal dynamics of the mean prevalence of shedders. (a), the mean environmental bacterial load (b) and the mean number of abortions (c) for the 4 vaccination scenarios. Scenario 1: vaccination of heifers and cows for a 10-year period (black line); scenario 2: vaccination of heifers and cows for a 3-year period with (scenario 2A - grey line) or without (scenario 2B - grey dotted line) loss of immunity one year after at the last vaccination; scenario 3: vaccination of heifers for a 10-year period (black dotted line); control: no vaccination (black thick line). Temporal dynamics of the mean prevalence of shedders (d) and mean environmental bacterial load (e) in scenario 1 with different values of p v (transition rate from SV e to I 1 V e ).

Influence of the ratio p v /p values on the model dynamics

As shown in Figure 2d and 2e for scenario 1, the mean prevalence of shedders was highly influenced by the values of the p v /p ratio (the ratio between the transition rate SV e => I 1 V e and the transition rate S => I 1 ), whereas the mean yearly number of abortions (results not shown) and the mean dynamics of environmental bacterial load were not affected by this parameter. For p v = 0.9p, the mean prevalence of shedders was almost stable within the first two years of vaccination and decreased afterwards to reach 8.3% five years after the start of the simulation. On the contrary, when considering p v = 0.05p, the decrease was much faster and the mean prevalence of shedders was less than 1% five years after the start of the simulation. In both cases, the mean prevalence of shedders was less than 1% at the end of the simulation time.

Influence of the vaccination scenarios and p v /p ratio on the extinction rate

Whereas the extinction rate was nil when no control programme was implemented, it varied from 8% to 97% between the vaccination scenarios and the values of p v /p (Table 1). It appears that most of the extinctions occurred late: as shown in Figure 3 for scenario 1, 82.5% of the extinctions occurred between the 5th and the 7th year of the vaccination programme.

Table 1 Extinction rate and mean time to extinction for each of the vaccination scenarios
Figure 3
figure 3

Distribution of times to extinction of the 97 extinct trajectories of scenario 1.


In this study, we explored the long term effectiveness of three different vaccination strategies in a recently infected dairy cattle herd, through a modelling approach. Mathematical models are nowadays one of the most effective tools to compare control measures for both human [2426] and animal infectious diseases [2730]. Here, we focused on vaccination since it is widely used in the field and was identified as a long-term control option for C. burnetii infections by the EFSA [31]. Vaccination with a phase I vaccine in cattle was indeed shown to suppress shedding in milk, placenta, uterine fluid, vagina and colostrum [32, 33]. In Rousset et al. [22], the vaccine appeared neither able to prevent infection in exposed kids, nor to clear infection in infected goats, but effectively reduced the level of shedding in a heavily infected herd. Hogerwerf et al. [34] also found that both the prevalence of shedders as well as the bacterial load in uterine fluid, vaginal swabs and milk were reduced in vaccinated dairy goats.

The vaccination scenarios tested in our study were (1) vaccination of the whole herd for 10 years, (2) vaccination of the whole herd for 3 years, and (3) vaccination of the heifers for only 10 years. Scenario 1 was predicted the most effective control strategy. In fact, all three vaccination strategies reduced the prevalence of shedders, the environmental bacterial load and the number of abortions. However, their effectivenesses were not equivalent. Since the infection was seldom eradicated in the first years of vaccination, an early cessation of vaccination (scenario 2) would be ineffective on the long run. Its short-term effect on infection dynamics depends on the duration of immunity for effectively vaccinated cows. According to Rodolakis et al. [35], in infected herds, more than 80% of the vaccinated cows still had immune markers one year after vaccination. However, at the same time, less than 60% of the vaccinated heifers were still skin-test positive. In the field, this means that immunity should last between one year (scenario 2A) and life long (scenario 2B). In this context, the increase in the prevalence of shedders, the environmental bacterial load and the number of abortions should not be observable in the first months following the cessation of vaccination. Nevertheless, the infection is spreading again. Thus, before stopping a vaccination programme on a farm, it seems essential to determine the presence or absence of C. burnetii in the herd. Diagnostic tests at a herd level (e.g. PCR in bulk tank milk) can probably be helpful [36], although they are imperfect.

According to our simulations, when only the heifers are vaccinated yearly (scenario 3), the decrease in the prevalence of shedders, the environmental bacterial load and the number of abortions is slower than when all the animals are vaccinated (scenario 1): it takes between 9 and 10 additional months to reach the same level of prevalence of shedders and environmental load, although the two strategies only differ in the initial action of the control programme. The extinction rate is high for both scenarios. Thus, although scenario 1 seems the best strategy from an epidemiological point of view, the difference between scenarios 1 and 3 is not so marked and a cost-benefit analysis would be useful to better compare the relative interest of these two strategies. It has to be highlighted that the numerical results of our study partially depend on the model structure and parameter values. The model represented the heterogeneity of shedding which is known to affect infection dynamics and hence the intervention efficacies in many diseases [37]. Indeed, sensitivity analysis show that model parameters governing the shedding levels, the characteristics of the bacterium in the environment as well as some of the physiological parameters strongly influenced the C. burnetii dynamics (see Additional File 1: section 2.1.). Here, parameter values were inferred or calibrated from field data of naturally infected dairy cattle herds [23]. Thus, we took into account the latest knowledge on C. burnetii infections. Although numerical values of the most influential parameters strongly influenced the output of numerical values, they did not change the relative ranking of the vaccination strategies. The only impact was that, for some combinations of parameter values, the differences between simulated effectiveness for the different vaccination scenarios became less marked, especially between scenarios 1 and 3.

The probability of infection for an effectively vaccinated susceptible cow p v was quantified based on the hazard ratio of the probability of shedding for vaccinated using non pregnant cows provided in Guatteo et al. [17]. Since the confidence interval of this parameter was wide, we studied the influence of this parameter value on the model outputs. Although the mean shedder prevalence was highly influenced by the p v /p ratio, the mean environmental bacterial load (which indirectly represents the infection risk for both animals and humans) decreased by roughly the same rate regardless of the ratio value. This is likely because the effectively vaccinated animals shed in decreased quantities. Therefore, irrespective of whether the mean prevalence of vaccinated shedders remains high, the prevalence of high shedders was reduced, with a major impact on the environmental load. This result has also been described by Lu et al. [38] who showed that to reduce the Salmonella prevalence in the long term, highly effective vaccines lowering the infectiousness would be a better choice than highly effective vaccines reducing susceptibility. Interestingly, the environmental bacterial load was hardly sensitive to the p v /p ratio (infection probability for effectively vaccinated cows), whereas the extinction rate was sensitive. Therefore, if the vaccine is to be used for eradication of C. burnetii from infected farms, both susceptibility and infectiousness of vaccinated animals have to be determined more accurately in order for the model to be used for prediction purposes or decision support. According to Rousset et al. [22], the lowest shedding level in vaginal swabs was shown to be more frequent in vaccinated than non vaccinated goats. However, further studies are needed to determine if a decrease of infectiousness is observed for all vaccinated animals or only for the animals vaccinated when non pregnant and still uninfected, and to quantify this decrease in all the shedding routes.

It should be noted that the extinction rate is highly influenced by the effect of vaccination on the susceptibility, the level of shedding and the mortality rate of the bacterium in the environment (see Additional file 1: section 2.1.2.), which are all not accurately documented variables. This extinction rate should then be interpreted with caution and used to compare different control strategies within the model. However, the behavior of the extinction rate suggests that it takes time to get free from C. burnetii within a herd.

This model was developed for dairy herds and the results presented here can not be straightforwardly generalized to small ruminants. A major difference is flock management. The typical size of a dairy flock is often much higher than the typical size of a dairy herd and kidding is usually synchronized. Therefore, the way to represent demography in the model should be adapted when representing C. burnetii spread in small ruminants. Besides, shedding characteristics and clinical manifestations may be different between species: according to Rodolakis et al. [15], ewes were found to shed mostly in faeces and vaginal mucus while goats seem to shed mostly in milk. Arricau-Bouvery et al. [21] reported that high abortion rates were rare, except in some caprine herds. However, all considered, our model represents an extensively documented basis for further development.

In conclusion, although an additional cost-benefit analysis considering the economic aspects of control programmes is needed to design an optimal control strategy, our modelling approach showed that a long term yearly vaccination would reduce infection risk in vaccinated herds.


  1. Plommet M, Capponi M, Gestin J, Renoux G: Fièvre Q expérimentale des bovins. Ann Rech Vet. 1973, 4: 325-346.

    Google Scholar 

  2. Aitken ID: Clinical aspects and prevention of Q fever in animals. Europ J Epidemiol. 1989, 5: 420-424. 10.1007/BF00140132.

    Article  CAS  Google Scholar 

  3. Bildfell RJ, Thomson GW, Haines DM, McEwen BJ, Smart N: Coxiella burnetii infection is associated with placentitis in cases of bovine abortion. J Vet Diagn Invest. 2000, 12: 419-425. 10.1177/104063870001200505.

    Article  CAS  PubMed  Google Scholar 

  4. Berri M, Souriau A, Crosby M, Crochet D, Lechopier P, Rodolakis A: Relationships between the shedding of Coxiella burnetii, clinical signs and serological responses of 34 sheep. Vet Rec. 2001, 148: 502-505. 10.1136/vr.148.16.502.

    Article  CAS  PubMed  Google Scholar 

  5. Karagiannis I, Schimmer B, Van Lier A, Timen A, Schneeberger P, Van Rotterdam B, De Bruin A, Wijkmans C, Rietveld A, Van Duynhoven Y: Investigation of a Q fever outbreak in a rural area of The Netherlands. Epidemiol Infect. 2009, 137: 1283-1294. 10.1017/S0950268808001908.

    Article  CAS  PubMed  Google Scholar 

  6. Panaiotov S, Ciccozzi M, Brankova N, Levterova V, Mitova-Tiholova M, Amicosante M, Rezza G, Kantardjiev T: An outbreak of Q fever in Bulgaria. Ann Ist Super Sanita. 2009, 45: 83-86.

    PubMed  Google Scholar 

  7. van der Hoek W, Dijkstra F, Schimmer B, Schneeberger PM, Vellema P, Wijkmans C, ter Schegget R, Hackert V, van Duynhoven Y: Q fever in the Netherlands: an update on the epidemiology and control measures. Euro Surveill. 2010, 15: pii:19520-

    Google Scholar 

  8. ECDC: Risk assessment on Q fever, Stockholm. 2010

    Google Scholar 

  9. Schimmer B, Dijkstra F, Vellema P, Schneeberger PM, Hackert V, ter Schegget R, Wijkmans C, van Duynhoven Y, van der Hoek W: Sustained intensive transmission of Q fever in the south of the Netherlands, 2009. Euro Surveill. 2009, 14: pii:19210-

    Google Scholar 

  10. McQuiston JH, Childs JE: Q fever in humans and animals in the United States. Vector Borne Zoonotic Dis. 2002, 2: 179-191. 10.1089/15303660260613747.

    Article  PubMed  Google Scholar 

  11. Gilsdorf A, Kroh C, Grimm S, Jensen E, Wagner-Wiening C, Alpers K: Large Q fever outbreak due to sheep farming near residential areas, Germany, 2005. Epidemiol Infect. 2008, 136: 1084-1087.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  12. Guatteo R, Beaudeau F, Joly A, Seegers H: Coxiella burnetii shedding by dairy cows. Vet Res. 2007, 38: 849-860. 10.1051/vetres:2007038.

    Article  PubMed  Google Scholar 

  13. Arricau-Bouvery N, Souriau A, Lechopier P, Rodolakis A: Experimental Coxiella burnetii infection in pregnant goats: excretion routes. Vet Res. 2003, 34: 423-433. 10.1051/vetres:2003017.

    Article  PubMed  Google Scholar 

  14. Welsh HH, Lennette EH, Abinanti FR, Winn JF: Air-borne transmission of Q fever: the role of parturition in the generation of infective aerosols. Ann N Y Acad Sci. 1958, 70: 528-540. 10.1111/j.1749-6632.1958.tb35409.x.

    Article  CAS  PubMed  Google Scholar 

  15. Rodolakis A: Q Fever in dairy animals. Ann N Y Acad Sci. 2009, 1166: 90-93. 10.1111/j.1749-6632.2009.04511.x.

    Article  PubMed  Google Scholar 

  16. Arricau-Bouvery N, Souriau A, Bodier C, Dufour P, Rousset E, Rodolakis A: Effect of vaccination with phase I and phase II Coxiella burnetii vaccines in pregnant goats. Vaccine. 2005, 23: 4392-4402. 10.1016/j.vaccine.2005.04.010.

    Article  CAS  PubMed  Google Scholar 

  17. Guatteo R, Seegers H, Joly A, Beaudeau F: Prevention of Coxiella burnetii shedding in infected dairy herds using a phase I C. burnetii inactivated vaccine. Vaccine. 2008, 26: 4320-4328. 10.1016/j.vaccine.2008.06.023.

    Article  CAS  PubMed  Google Scholar 

  18. Keeling MJ, Rohani P: Modeling infectious diseases in Humans and animals. 2008, Princeton University Press, Princeton

    Google Scholar 

  19. Durand MP: Lacteal and placental excretion of Coxiella burnetti, agent of Q fever, in the cow. Importance and prevention. Bull Acad Natl Med. 1993, 177: 935-945. discussion 945-946 (in French)

    CAS  PubMed  Google Scholar 

  20. Guatteo R, Beaudeau F, Berri M, Rodolakis A, Joly A, Seegers H: Shedding routes of Coxiella burnetii in dairy cows: implications for detection and control. Vet Res. 2006, 37: 827-833. 10.1051/vetres:2006038.

    Article  PubMed  Google Scholar 

  21. Arricau-Bouvery N, Rodolakis A: Is Q fever an emerging or re-emerging zoonosis?. Vet Res. 2005, 36: 327-349. 10.1051/vetres:2005010.

    Article  CAS  PubMed  Google Scholar 

  22. Rousset E, Durand B, Champion JL, Prigent M, Dufour P, Forfait C, Marois M, Gasnier T, Duquesne V, Thiery R, Aubert MF: Efficiency of a phase 1 vaccine for the reduction of vaginal Coxiella burnetii shedding in a clinically affected goat herd. Clin Microbiol Infect. 2009, 15 (Suppl 2): 188-189.

    Article  CAS  PubMed  Google Scholar 

  23. Courcoul A, Vergu E, Denis JB, Beaudeau F: Spread of Q fever within dairy cattle herds: key parameters inferred using a Bayesian approach. Proc Biol Sci. 2010, 277: 2857-2865. 10.1098/rspb.2010.0575.

    Article  PubMed Central  PubMed  Google Scholar 

  24. Acedo L, Diez-Domingo J, Morano JA, Villanueva RJ: Mathematical modelling of respiratory syncytial virus (RSV): vaccination strategies and budget applications. Epidemiol Infect. 2010, 138: 853-860. 10.1017/S0950268809991373.

    Article  CAS  PubMed  Google Scholar 

  25. Sharomi O, Podder CN, Gumel AB, Mahmud SM, Rubinstein E: Modelling the Transmission Dynamics and Control of the Novel 2009 Swine Influenza (H1N1) Pandemic. Bull Math Biol. 2011, 73: 515-548. 10.1007/s11538-010-9538-z.

    Article  CAS  PubMed  Google Scholar 

  26. Van Effelterre T, Moore MR, Fierens F, Whitney CG, White L, Pelton SI, Hausdorff WP: A dynamic model of pneumococcal infection in the United States: implications for prevention through vaccination. Vaccine. 2010, 28: 3650-3660. 10.1016/j.vaccine.2010.03.030.

    Article  PubMed  Google Scholar 

  27. Le Menach A, Vergu E, Grais RF, Smith DL, Flahault A: Key strategies for reducing spread of avian influenza among commercial poultry holdings: lessons for transmission to humans. Proc Biol Sci. 2006, 273: 2467-2475. 10.1098/rspb.2006.3609.

    Article  PubMed Central  PubMed  Google Scholar 

  28. Andraud M, Rose N, Grasland B, Pierre JS, Jestin A, Madec F: Influence of husbandry and control measures on porcine circovirus type 2 (PCV-2) dynamics within a farrow-to-finish pig farm: a modelling approach. Prev Vet Med. 2009, 92: 38-51. 10.1016/j.prevetmed.2009.07.009.

    Article  CAS  PubMed  Google Scholar 

  29. Backer JA, Hagenaars TJ, van Roermund HJ, de Jong MC: Modelling the effectiveness and risks of vaccination strategies to control classical swine fever epidemics. J R Soc Interface. 2009, 6: 849-861.

    Article  PubMed Central  PubMed  Google Scholar 

  30. Garner MG, Cowled B, East IJ, Moloney BJ, Kung NY: Evaluating the effectiveness of early vaccination in the control and eradication of equine influenza-A modelling approach. Prev Vet Med. 2011, 99: 15-27. 10.1016/j.prevetmed.2010.02.007.

    Article  CAS  PubMed  Google Scholar 

  31. EFSA Panel on Animal Health and Welfare (AHAW): Scientific Opinion on Q Fever. EFSA Journal. 2010, 8 (5): 1595-

    Google Scholar 

  32. Sadecky E, Brezina R, Kazar J, Urvolgyi J: Immunization against Q-fever of naturally infected dairy cows. Acta Virol. 1975, 19: 486-488.

    CAS  PubMed  Google Scholar 

  33. Biberstein EL, Riemann HP, Franti CE, Behymer DE, Ruppanner R, Bushnell R, Crenshaw G: Vaccination of dairy cattle against Q fever (Coxiella burnetii): results of field trials. Am J Vet Res. 1977, 38: 189-193.

    CAS  PubMed  Google Scholar 

  34. Hogerwerf L, van den Brom R, Roest HIJ, Bouma A, Vellema P, Pieterse M, Dercksen D, Nielen M: Reduction of Coxiella burnetii prevalence by vaccination of goats and sheep, The Netherlands. Emerg Infect Dis. 2011, 17: 379-386.

    Article  PubMed Central  PubMed  Google Scholar 

  35. Rodolakis A, Clement P, Cochonneau D, Beaudeau F, Sarradin P, Guatteo R: Investigation of humoral and cellular immunity of dairy cattle after one or two year of vaccination with a phase I Coxiella vaccine. Procedia in Vaccinology. 2009, 1: 85-88. 10.1016/j.provac.2009.07.015.

    Article  Google Scholar 

  36. Guatteo R, Beaudeau F, Joly A, Seegers H: Assessing the within-herd prevalence of Coxiella burnetii milk-shedder cows using a real-time PCR applied to bulk tank milk. Zoonoses Public Health. 2007, 54: 191-194. 10.1111/j.1863-2378.2007.01043.x.

    Article  CAS  PubMed  Google Scholar 

  37. Matthews L, Woolhouse M: New approaches to quantifying the spread of infection. Nat Rev Microbiol. 2005, 3: 529-536. 10.1038/nrmicro1178.

    Article  CAS  PubMed  Google Scholar 

  38. Lu Z, Grohn YT, Smith RL, Wolfgang DR, Van Kessel JA, Schukken YH: Assessing the potential impact of Salmonella vaccines in an endemically infected dairy herd. J Theor Biol. 2009, 259: 770-784. 10.1016/j.jtbi.2009.04.028.

    Article  PubMed  Google Scholar 

Download references


The authors would like to thank Alain Joly for fruitful discussions on vaccination strategies. The collaboration between the Faculty of Veterinary Medicine of Utrecht and the French National Institute for Agricultural Research (INRA) was financially supported by the Netherlands Organisation for Scientific Research (NWO) and the French Ministry of Foreign and European Affairs through the Van Gogh Programme.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Aurélie Courcoul.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

AC participated in the design of the study, carried out the model development and analysis, and drafted the manuscript. LH, DK and MN participated in the design of the study and helped to draft the manuscript. FB and EV participated in the design and coordination of the study, and helped to draft the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Courcoul, A., Hogerwerf, L., Klinkenberg, D. et al. Modelling effectiveness of herd level vaccination against Q fever in dairy cattle. Vet Res 42, 68 (2011).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: