Integrating interdisciplinary methodologies for One Health: goat farm re-implicated as the probable source of an urban Q fever outbreak, the Netherlands, 2009
BMC Infectious Diseases volume 15, Article number: 372 (2015)
In spring 2008, a goat farm experiencing Q fever abortions (“Farm A”) was identified as the probable source of a human Q fever outbreak in a Dutch town. In 2009, a larger outbreak with 347 cases occurred in the town, despite no clinical Q fever being reported from any local farm.
Our study aimed to identify the source of the 2009 outbreak by applying a combination of interdisciplinary methods, using data from several sources and sectors, to investigate seventeen farms in the area: namely, descriptive epidemiology of notified cases; collation of veterinary data regarding the seventeen farms; spatial attack rate and relative risk analyses; and GIS mapping of farms and smooth incidence of cases. We conducted further spatio-temporal analyses that integrated temporal data regarding date of onset with spatial data from an atmospheric dispersion model with the most highly suspected source at the centre.
Our analyses indicated that Farm A was again the most likely source of infection, with persons living within 1 km of the farm at a 46 times larger risk of being a case compared to those living within 5-10 km. The spatio-temporal analyses demonstrated that about 60 – 65 % of the cases could be explained by aerosol transmission from Farm A assuming emission from week 9; these explained cases lived significantly closer to the farm than the unexplained cases (p = 0.004). A visit to Farm A revealed that there had been no particular changes in management during the spring/summer of 2009, nor any animal health problems around the time of parturition or at any other time during the year.
We conclude that the probable source of the 2009 outbreak was the same farm implicated in 2008, despite animal health indicators being absent. Veterinary and public health professionals should consider farms with past as well as current history of Q fever as potential sources of human outbreaks.
Q fever is a zoonosis caused by the bacterium Coxiella burnetii and is present worldwide with the exception of New Zealand. Clinical disease in humans can range from an acute influenza-like illness to a chronic infection manifesting mainly as endocarditis or vascular infection , though 60 % of infections are asymptomatic. There is also evidence of a “post Q-fever fatigue syndrome” manifesting as chronic fatigue following acute infection . The animal reservoir for Q fever is large as a wide range of mammals, birds, reptiles and fish can become infected with C. burnetii, but cattle and small ruminants (sheep and goats) are the most affected . Clinical presentation of Q fever in domestic ruminants is most commonly associated with late abortions, stillbirth and delivery of weak offspring [3–5]. The most frequent route of infection for humans is from inhalation of aerosols contaminated by infected ruminants at slaughter or parturition . Windborne spread of such aerosols can cause infections several kilometres from their origin [1, 6–9]. Furthermore, environmental contamination related to parturition or slaughter can last potentially from months to years, making inhalation of contaminated dust a risk .
Between 2007 and 2010, the Netherlands experienced the largest outbreak of Q fever ever reported, with over 4,000 notified cases [10, 11] and mainly associated with dairy goat farms [10–13]. Since that time, a number of methodologies have been developed and applied to investigate the epidemic. These methods incorporate traditional epidemiological approaches alongside Geographic Information Systems (GIS) [12, 13], mathematical modelling , and atmospheric dispersion models ; and use data from a range of sectors and sources, including animal, human, environmental and meteorological. In this sense, the Q-fever epidemic has offered a unique opportunity to develop interdisciplinary, “One-Health” methodologies to investigate airborne and zoonotic disease.
In the spring/summer of 2008, a Q fever outbreak occurred in an urban area (approximately 88,000 inhabitants) in the south of the Netherlands with 96 cases notified to the local Municipal Health Service (MHS), Brabant-Zuidoost. Through use of a Geographic Information System, this outbreak was subsequently linked to a commercial dairy goat farm (Farm A) that had experienced a Q fever-related abortion wave in the weeks preceding the human outbreak . Persons living within 2 km of the farm had over a thirty times higher risk for Q fever than those living more than 5 km away, with risk declining with increasing distance from the farm. This spatial relationship, together with the fact that no other farms in the locality had reported Q-fever problems, supported the hypothesis that this single goat farm was the source of the human outbreak.
In 2009, an even larger human outbreak occurred in the town, with 347 cases notified to the same MHS between April and July. However, in this year no farm in the locality reported any clinical Q fever or related animal health problems. The reasons for the upsurge in human cases were unclear. Hypotheses included: increased reporting of Q fever owing to greater awareness of the disease amongst medical professionals and the community; the existence of a new source or sources of infection; or Farm A being again the source due to sub-clinical shedding or environmental contamination. Our study aimed to integrate the interdisciplinary methodologies and inter-sectoral data developed and gathered since the start of the national epidemic to identify the most likely origin of the 2009 outbreak in the area.
We defined cases as residents of the Dutch Municipal Health Service region ‘Brabant-Zuidoost’ who were notified to the national surveillance system for infectious diseases as having Q fever with date of onset of between weeks 11 and 36 (9 March-6 September), 2009. In 2009, the notification criteria for Q fever were (1) fever, pneumonia or hepatitis, and (2) evidence of C. burnetii infection through a fourfold rise in IgG antibody titre in paired sera, presence of IgM phase II antibodies, or detection of nucleic acid by polymerase chain reaction (PCR) in blood or respiratory tract samples . We excluded 48 cases who were known to have attended farm viewing days at a Q fever positive local petting farm in 2009, which had previously been implicated in transmission of the disease to visitors .
Comparison of 2008 to 2009 cases
We compared the 2009 cases to cases notified in 2008 to the same MHS using descriptive analysis and standard statistics (Student’s T-test or chi-square test as appropriate). Specifically, we examined the cases in terms of time, place and person (epidemic curve, age, gender ratio, and postcode of residence accurate to neighbourhood level), together with hospitalisation rates and time from date of onset to laboratory confirmation to see if there were any significant differences in these parameters in 2009 as compared to 2008. For 2008 cases, we used the same case definition that Schimmer et al.  employed in their investigation of the 2008 outbreak. Data were analysed using Stata 12.
Suspect source farm identification
Van Leuken et al.  identified 56 farms holding small ruminants in the area of the 2009 outbreak, of which seventeen were indicated to be potential sources using an exponential model to fit Q fever incidence against increasing distance from each farm. We used three approaches to investigate these seventeen farms further:
Collation of veterinary and farm data: we obtained data from the Animal Health Service (GD) regarding the farm location and number of goats and sheep held in 2009 for each of the seventeen farms. We also collated data on Q-fever history of each farm from three sources: (1) clinical Q-fever notifications in 2008 and 2009; (2) results of the non-mandatory bulk tank milk testing for C. burnetii in 2008 on large dairy goat farms (for which the take-up rate was ~75 % of dairy goat farms with >50 animals) ; (3) results of the mandatory, systematic bulk tank milk testing for C. burnetii, introduced in October 2009 and conducted at all dairy farms with >50 goats/sheep ; and (4) results of previous animal and environmental Q-fever investigations in the study area [19–21].
Concentric ring attack rate analyses: following the approach used by Schimmer et al.  to investigate the 2008 outbreak in the area, we used distance of a case’s residence to a farm as a proxy for exposure together with 2009 population data to calculate attack rates (AR) for residents in concentric 1 km rings around each potential source farm, up to a 5 km radius (i.e. 0-1 km, 1-2 km, 2-3 km, 3-4 km and 4-5 km). We also calculated the attack rate in a reference ring of 5-10 km around each potential source farm and used this to compute relative risks (RR) for each 1 km concentric ring. Data were analysed using R version 3.0.3.
Mapping: using ArcGIS 9 software (ESRI, Redlands, CA, USA), we geo-referenced the residential addresses of each case according to six figure postcode (discriminatory down to street level) and used 2009 population data to produce smooth incidence maps, as developed by Van der Hoek et al.  Briefly, this method imposes a 500 x 500 m grid over the area of interest; counts the number of cases and the population number in each cell; and models the number of cases as an inhomogenous Poisson process whereby the underlying incidence is estimated by 2-dimensional P-spline smoothing. We also mapped the seventeen farms identified as potential sources, together with farm size and animal species. As several farms were in close proximity to each other, the ARs and RRs of a particular farm calculated in the concentric ring analyses could have been an artefact of the C. burnetii emissions from a nearby farm; therefore mapping served to visualise whether this phenomenon might be occurring.
Further investigation of the most strongly suspected farm
We then performed the following additional investigations for the farm most strongly associated with occurrence of cases according to the previous methods (Farm A): (1) further spatio-temporal analyses using this farm as the “emission point” for C. burnetii, and (2) an on-site farm visit:
1) Spatio-temporal analyses using an atmospheric dispersion model: we defined a “period of probable infection (PPI)” denoting the period during which a case was likely to have been infected, i.e. 8–33 days (95 % CI) prior to the date of onset of disease . We hypothesised that relatively high modelled concentrations of C. burnetii originating from Farm A during the PPI preceding a particular case’s date of onset may indicate that the case in question was infected via airborne transmission from the farm. That is, if a case was infected by exposure to C. burnetii from Farm A, it would have been exposed mainly during its PPI rather than during the period prior to that.
To test this hypothesis, we applied an atmospheric dispersion model previously described by Van Leuken et al.  to model the airborne transmission of C. burnetii. The model calculates the physical dispersion of C. burnetii bacteria from a particular emission point (in this case a farm) by taking into account hourly-based meteorological data, including wind velocity and direction. The model output consists of hourly averaged C. burnetii concentration matrices which are then converted to period-specific averaged concentration maps.
We inserted hourly averaged local meteorological conditions from the Royal Netherlands Meteorological Institute into the model and performed sensitivity analyses on the threshold wind velocity (0, 2, 4 and 6 m/s) and the emission profile: conYear (constant emission in time) or lNormEpi (log-normal emission curve) starting from the first day of the kidding season (for details see Van Leuken et al. ). For a case i, we determined the cumulative modelled concentration at its street-level postcode during its PPI (CPPI,i), and also prior to the PPI (CPRE,i). Subsequently, we divided all cases into two groups: “group 1”, containing all cases who had a higher cumulative concentration during their PPI compared to prior to it (CPPI > CPRE); and “group 2”’ containing all cases with a higher cumulative concentration prior to their PPI compared to during the PPI (CPPI ≤ CPRE). In other words, a case with a cumulative exposure ratio CPPI: CPRE of 90:10 was clearly attributed to group 1, whereas a case with a ratio of 10:90 was attributed to group 2. The cut-off was thus set at 50:50. Figure 1 gives an example of the classification of cases into groups 1 and 2.
We created histograms of both groups as a function of distance from Farm A, and, in addition, we determined the relative abundance of cases (ϕ) per group per circular ring r of 1 km:
where nr,1 and nr,2 are the number of cases in ring r in group 1 and 2 respectively, and n1 and n2 are the total number of cases per group. Thus, the denominator and the numerator give the proportion of cases in each group as a function of distance. Variable ϕ then gives the ratio of these proportions and thus determines which group is over-represented as a function of distance from Farm A.
For the purposes of the model, we assumed that emission of C. burnetii from the farm started in week 9 in 2009, i.e. two weeks before the first cases were notified (to correspond approximately with the incubation period of 8 – 33 days ). However, to test the dependence of ϕ on the week that the emissions commenced, we additionally performed a sensitivity analysis regarding week of emission.
Finally, we performed a Kruskal-Wallis test (a non-parametric test) to compare the median distance of both groups to the suspected source farm (since the cases are spatially clustered and therefore the incidence-distance relationship is exponentially curved , these medians (or means) cannot be compared using a standard t-test, which assumes normality). We compared the resulting T-statistic to a χ2-distribution with one degree of freedom, assuming a threshold for significance of p < 0.05.
2) Farm visit: on 17 August 2012, we visited Farm A and interviewed the farmer using a semi-structured questionnaire to gain insight into the farm characteristics, husbandry methods, breeding periods, and animal health status in spring-summer 2009.
Q fever is a notifiable disease in the Netherlands and according to Dutch legislation, written consent from patients for the use of anonymised surveillance data is not necessary. Therefore, ethical review was not required. The Municipal Health Service Brabant-Zuidoost made the case information available to the researchers. However, publication of this case information is not permitted.
Comparison of 2008 to 2009 cases
1) Time and Place: a total of 320 cases met the 2009 case definition, compared to 96 in 2008. All cases identified in 2009 were different from those identified in 2008, ie were new cases. Figure 2 shows the distribution of cases over time from 2008 to 2009. The 2008 and 2009 outbreaks both started and ended at broadly similar times, with date of onset of cases between weeks 16 and 32 in 2008 compared to weeks and 11 and 36 in 2009. The peak number of cases occurred in week 23 in 2008 (n = 19) compared to week 15 in 2009 (n = 54). Cases lived in similar postcode areas to the 2008 outbreak , with the same two neighbourhoods being most affected both in 2008 and 2009 (65 % of cases and 53 % of cases respectively); these two neighbourhoods are contiguous.
2) Demographics: the age and gender of cases did not differ significantly between 2008 and 2009. In 2008, median age affected was 50, compared to 49 in 2009 (p = 0.79). In 2008, 68.8 % of cases were male compared to 61.6 % in 2009 (p = 0.20).
3) Hospitalisation rate and laboratory confirmation delay: hospitalisation rate in 2009 was higher than in 2008 (23.7 % and 18.8 % respectively), but this difference was not significant (p = 0.56). Laboratory confirmation delay (time between date of symptom onset and date of laboratory confirmation) was significantly shorter in 2009 compared to 2008 (p = 0.0001), with a median delay of 20 days in 2009 compared to 29 days in 2008.
Suspect source farm identification
Collation of veterinary and farm data
Table 1 describes the farms in terms of their size, species and Q fever history. Of the seventeen farms, eight had only goats, seven only sheep and two were mixed. Farm A was the largest farm with 791 goats (no sheep); Farm H was the second largest with 175 sheep (no goats), and the remaining farms all had fewer than 40 animals. Three of the seventeen farms had a history of Q fever. Farm A had clinical Q fever (abortion wave) in 2008, had tested positive for C. burnetii during bulk tank milk (BTM) testing in 2008 and 2009; and had almost 100 % of various samples testing positive in previous investigations since the original 2008 outbreak (vaginal and stable environmental swabs in 2009 ; and aerosol samples in the stable and the outdoor environment, together with swabs from the stable environment, in 2010 . Farm H had 5/20 positive samples , and Farm C had 3/12 positive swabs in 2008 . The remaining fourteen farms were too small (<50 animals) to have qualified for bulk tank milk testing, and had neither been included in other previous animal/environmental investigations; however, none had reported any clinical Q fever nor abortion problems to the Animal Health Service in 2008 or 2009.
Concentric ring attack rate analyses
Figure 3 shows the Relative Risk for Q fever with increasing distance from each of the seventeen farms (data are presented in Additional file 1: Table S1). Farms A, B and E had the most striking distance-response relationship, followed by G, and then D, C, and F to a lesser extent. The remaining ten farms showed no clear relationship between distance from farm and RR. Farm A had the highest overall attack rate in the 0–1 km ring at 2156 cases per 100,000 population, although the 0-1 km attack rates for Farms B (1802 cases per 100,000 population), E (1818 cases per 100,000 population), F (1031 cases per 100,000 population) and G (1132 cases per 100,000 population) were not statistically significantly different to this figure. Farm A also had the highest relative risk in the 0–1 km ring compared to the 5–10 km reference ring (RR = 46; 95 % CIs 22–97).
Figure 4 shows the map of the residential addresses of notified cases, the location and size of the seventeen potential source farms, and the smooth incidence of human cases in the locality over the outbreak period (weeks 11–36). Reflecting the results of the concentric rings analysis, Farms B and E were in closest proximity to Farm A (502 m and 196 m respectively), followed by Farm G (801 m).
Further investigation of the most strongly suspected farm
Taking all the above information together, Farm A was identified as the most strongly suspected potential source farm. Our reasoning for this is as follows: it was by far the largest farm of the seventeen; it had a confirmed history of clinical Q fever (2008) plus positive BTM investigations subsequently in 2008 and 2009; it was a goat farm rather than a sheep farm (goat farms being more frequently indicated in the Dutch Q fever epidemic than sheep farms ); and the results of the concentric rings analysis combined with the mapping pointed most strongly towards it (although Farms B and E also had very strong distance-RR relationships in the concentric rings analysis, the mapping showed that they were both very close to Farm A. Therefore we considered that these relationships were likely to be an artefact of proximity, particularly given that Farm B held only five goats and Farm E, two sheep).
Spatio-temporal analyses using an atmospheric dispersion model
Figure 5 (upper two rows) shows the histograms of the number of cases in groups 1 and 2 as a function of distance for emission profile conYear, stratified by the threshold wind velocity (see Additional file 2: Figure S1 for emission profile lNormEpi). These results show that approximately 60 – 65 % of the cases are attributed to group 1. Also, the emission profile and threshold wind speed have little influence on the shape of the distributions. The third row shows the ratio of the proportions (ϕ) from equation 1: group 1 is over-represented at distances up to 5 km (i.e. the points lie above 1.0). At distances of 5 – 7 km, in contrast, group 2 seems to be over-represented (i.e. the points lie below 1.0), although the absolute number of cases is smaller. Note that the first point should be interpreted carefully as the number of cases up to 1000 m in group 1 is equal to zero. Additional file 3: Figure S2 shows the actual cumulative concentrations (CPPI and CPRI) per case.
Figure 6 shows the results of the sensitivity analysis of the number of cases in group 1 as a function of the starting week of emissions (assuming emission profile conYear and a threshold wind velocity of 2 m/s). In our model, we assumed that the emissions commenced in week 9, corresponding with about 61 % of the cases classified to group 1; however, if we had set the first week of emissions to be week 12 or 13, then this number would have been slightly higher at approximately 69 %.
Table 2 shows the median distances of a case’s residence to Farm A per group as a function of the threshold wind velocity. The 95 % confidence intervals were determined by bootstrapping 10,000 times. According to the Kruskal-Wallis test the median distance of group 1 is always significantly lower (p ≈ 0.004) than that of group 2, indicating that cases of group 1 are, on average, living closer to Farm A than cases of group 2. The fraction of cases classified to group 1 has little association with emission profile or threshold wind velocity.
Finally, the bottom row in Fig. 5 shows the logs of the cumulative concentrations that cases were exposed to during their PPI (CPPI) for emission profile conYear and threshold wind velocities of 0, 2, 4, and 6 m/s (V0, V2, V4, V6 respectively). The results show that the differences between the groups were very small, but for higher threshold wind speeds the cumulative concentration during the PPI was higher in group 1. Note that the data become more scattered for high threshold wind speeds: this is caused by a more plume-shaped exposure pattern (with higher contrast) for the cumulative concentration maps (see Additional file 3: Figure S2).
Visit to Farm A: in 2009, the goats were kept on a deep litter system with fresh straw added daily. The shed was mechanically ventilated. The litter was removed three to four times during the year. No further cleaning was done before replenishing with fresh litter. Birthing products were not cleared from the pens during the kidding season. The farmer could not recall any particular changes in management in the spring/summer of 2009, nor any animal health problems around the time of parturition or at any other time during the year. The goats had not been vaccinated against Q fever in 2008, but were vaccinated against the disease in May/June 2009. Although the farmer had collected data regarding the period of kidding, numbers of kids born and number of live/still births in 2009 for submission to the Animal Health Service, he no longer had these data available. We were unable to access these data from the Animal Health Service independently.
This study aimed to identify the most likely origin of an outbreak of Q fever in an urban area in 2009 that occurred in the absence of any clinical Q-fever notifications from surrounding farms that year. Our results suggest that the most likely source of this outbreak was the same dairy goat farm, Farm A, as had been implicated in a smaller outbreak in the area the previous year. We gathered data from a wide range of sources and sectors (animal, human, environmental and meteorological) and applied a range of interdisciplinary methods incorporating both traditional and innovative techniques in order to reach this conclusion. Integrated approaches such as these are increasingly recognised as key to One Health research. The mechanisms underlying zoonotic disease emergence and transmission are invariably complex and multi-factorial, such that understanding them requires interdisciplinary work that applies a holistic set of methodologies and data from the human health, veterinary, and environmental sectors [24, 25]. For the purposes of this study, we also developed a new method that enabled the integration of spatial data generated through an atmospheric dispersion model with temporal data concerning time of probable infection. This method could prove valuable to future Q fever outbreak investigations as it can serve to disentangle acute aerosol transmission from a point source from diffuse environmental transmission.
While comparison of the area’s 2008 outbreak data to its 2009 outbreak data suggests that there was a higher awareness of Q fever among the public and medical community in the area in 2009 compared to 2008 (demonstrated by the significantly shorter period between date of onset of disease and date of laboratory confirmation in 2009), the similar demographics of cases in terms of time, place and person and the lack of a reduction in hospitalisation rate do not support the hypothesis that the increase in cases of Q fever seen in 2009 was caused solely by changes in reporting behaviour. Instead, our results suggest that this was a true outbreak rather than a reporting artefact, with Farm A as the probable source. The timing of the 2009 outbreak (in the spring around the kidding season) suggests that the outbreak was due to goats sub-clinically shedding C. burnetii during this period, rather than due to environmental contamination of Farm A and its surrounds during 2008.
This is supported by the results of spatio-temporal analyses using an atmospheric distribution model. By means of this method, we calculated the proportion of cases that were mainly exposed during their probable period of infection (PPI), namely 60 – 65 % of the cases. Thus, these cases could be explained by airborne transmission of C. burnetii from Farm A. The other cases were mainly exposed prior to their PPI; explanations for this could be that (1) they were absent from home during high-exposure hours/days, or (2) they did not happen to inhale an infected particle during the PPI by chance, despite high concentrations around their area of residence at that time, or (3) they were infected by a different source (including re-aerosolised bacteria from contaminated environments).
It has been suggested that C. burnetii may still be excreted via the placenta during normal parturition [3, 4, 26, 27] even in vaccinated animals, albeit to a lesser extent . A recent meta-analysis of studies investigating the effect of vaccination on shedding suggested that while vaccination significantly reduced the risk of shedding form uterine secretions in previously exposed goats, shedding through all other routes (milk, vaginal secretions and faeces) was not significantly reduced.  Furthermore, farms that have experienced Q fever-related abortions in the past have been shown to have higher proportions of C. burnetii-positive animal and environmental samples, and higher levels of C. burnetii DNA within those positive samples compared to both farms negative for Q fever and farms positive on bulk milk sampling but with no history of abortion. Unfortunately detailed veterinary/husbandry data were unavailable regarding the start, end and peak kidding periods or other on-farm events in 2009 (e.g. removal of litter from the pen) that might have contributed to the rise in human cases. For the purposes of the model, we assumed that C. burnetti emission started in week 9, two weeks before the first notification. Although the proportion of explained cases would have been larger had we assumed an emission start in week 12 or 13 (Fig. 6), cases with weeks of onset between weeks 11 and 14 would subsequently have been missed by the model.
The model also demonstrated a strong relationship between distance from cases’ residence to Farm A and incidence of Q fever stratified by group. Cases belonging to group 1 (i.e. mainly exposed during their PPI) were over-represented at distances of up to 4 to 5 kilometres from Farm A, while cases in group 2 were living at relatively further distances of 5 to 8 kilometres (there were too few cases at 9 km and 10 km to draw conclusions about these distances). Also, cases in group 1 were generally exposed to higher doses during their PPI (Fig. 5) compared to cases in group 2. Note that we did not look at the difference in absolute exposure between cases, but instead at the exposure per case. After all, the probability of infection is always 100 % at the end of a case’s PPI, regardless of the actual dose. Nevertheless, other sources could have increased the actual dose and thus the probability of infection. Furthermore, while choice of emission profile (i.e. conYear and lNormEpi) and threshold wind velocity had only a minor effect on the proportion of explained cases (Fig. 5 and Additional file 2: Figure S1), exposure at short distances up to 4-5 km was an important indicator for disease. These results complement the findings from concentric rings attack rate analyses in both this study and previous studies [12, 13] which have indicated that persons living within 5 km of an infected farm are at particular risk of contracting Q fever.
Although seventeen farms were included for investigation at the outset, we focussed our farm visit and spatio-temporal analyses to include Farm A only. It could be argued that, despite the weight of evidence pointing to Farm A, similar investigations should have been performed at all farms for the purposes of fair comparison. One argument for focussing on Farm A was its size; with 791 animals it was by far the largest farm of the seventeen: indeed, fifteen farms had under 40 animals, twelve of which had under ten. The only other large farm, Farm H (175 sheep), was not further investigated because the concentric rings analysis and smooth incidence mapping did not point towards it; however, five out of 20 environmental and aerosol samples taken at and around the farm in 2009/2010 tested positive for C. burnetii . Furthermore, although it may be considered unlikely that the small farms would be the source of such a large human outbreak, it is theoretically possible. Unfortunately, because these farms were so small, they did not fall under statutory regulations for bulk tank milk sampling; neither had they been included in previous investigations, so we had no data to rule out Q fever in these farm animals. Ideally, all seventeen farms would have been individually visited and animals sampled as part of our investigations; however, our resources did not permit this and we focussed on the most strongly suspected farm only. With regards the spatio-temporal modelling: the concentric rings analyses suggested most strongly that farms A, B and E, and possibly G, warranted further investigation. However, as these farms were so close to Farm A, any results yielded from atmospheric dispersion model centred on these farms would have been almost identical to that for Farm A and therefore would have yielded almost identical information.
Caution is called for in the interpretation of the spatio-temporal analyses. The method requires that cases are assigned into one of two distinct groups: those who were exposed to a higher modelled concentration of C. burnetii in the 8–33 days before date of onset (the period of probable infection, PPI) compared to prior to it (group 1), and those who had a higher modelled concentration prior to their PPI compared to during the it (group 2). However, this was not always a straightforward distinction (Fig. 5) and some cases may have been mis-assigned. It is difficult to extrapolate from modelled concentration at a particular postcode area to the actual level of exposure of a person in reality: for example, a case may have been absent from their residence during their PPI; there may have been a relatively low number of contaminated particles in their particular street; or they may have been exposed to other sources in addition to Farm A.
Our study found that most likely source of a large urban outbreak of Q fever in 2009 was the same dairy goat farm implicated in a smaller outbreak in the area the previous year. The results of our investigations highlight that farms with a history of Q fever can still pose a risk to public health even when there are no animal health indicators such as abortion storms or other fertility problems. We recommend that veterinary and public health professionals consider farms with past as well as current history of Q fever as potential sources of human outbreaks. We used a variety of methods to reach our conclusions, including a newly developed method to investigate such outbreaks in space and time using an atmospheric dispersion model. Using meteorological forecast data, this method could potentially be used prospectively to visualise the spread of Coxiella bacteria and inform infectious disease prevention and response measures.
Bulk Tank Milk
Animal Health Service
Geographic Information System
Municipal Health Service
Period of Probable Infection
Polymerase Chain Reaction
Parker NR, Barralet JH, Bell AM. Seminar Q fever. Lancet. 2006;367:679–88.
Marmion BP, Shannon M, Maddocks I, Storm P, Penttila I. Protracted debility and fatigue after acute Q fever. Lancet. 1996;347:977–8.
Rousset E, Berri M, Durand B, Dufour P, Prigent M, Delcroix T, et al. Coxiella burnetii shedding routes and antibody response after outbreaks of Q fever-induced abortion in dairy goat herds. Appl Environ Microbiol. 2009;75:428–33.
Berri M, Rousset E, Champion JL, Russo P, Rodolakis A. Goats may experience reproductive failures and shed Coxiella burnetii at two successive parturitions after a Q fever infection. Res Vet Sci. 2007;83:47–52.
Roest H-J, Van Gelderen B, Dinkla A, Frangoulidis D, Van Zijderveld F, Rebel J, et al. Q fever in pregnant goats: pathogenesis and excretion of Coxiella burnetii. PLoS ONE. 2012;7.
Hawker JI, Ayres JG, Blair I, Evans MR, Smith DL, Smith EG, et al. A large outbreak of Q fever in the West Midlands: windborne spread into a metropolitan area? Commun Dis Public Health. 1998;1:180–7.
Tissot-Dupont H, Amadei M, Nezri M, Raoult D. Wind in November, Q fever in December. Emerg Infect Dis. 2004;10:1264–9.
EFSA. Scientific Opinion on Q fever. EFSA Journal. 2010;8(5):1595.
Wallensten A, Moore P, Webster H, Johnson C, Van der Burgt G, Pritchard G, et al. Q fever outbreak in Cheltenham, United Kingdom, in 2007 and the use of dispersion modelling to investigate the possibility of airborne spread. Eurosurveillance. 2010;15:1–7.
Roest HIJ, Tilburg JJHC, Van der Hoek W, Vellema P, Van Zijderveld FG, Klaassen CHW, et al. The Q fever epidemic in The Netherlands: history, onset, response and reflection. Epidemiol Infect. 2011;139:1–12.
Dijkstra F, Van der Hoek W, Wijers N, Schimmer B, Rietveld A, Wijkmans CJ, et al. The 2007–2010 Q fever epidemic in The Netherlands: characteristics of notified acute Q fever patients and the association with dairy goat farming. FEMS Immunol Med Microbiol. 2012;64:3–12.
Schimmer B, Ter Schegget R, Wegdam M, Züchner L, De Bruin A, Schneeberger PM, et al. The use of a geographic information system to identify a dairy goat farm as the most likely source of an urban Q-fever outbreak. BMC Infect Dis. 2010;10.
Hackert VH, Van der Hoek W, Dukers-Muijrers N, De Bruin A, Al Dahouk S, Neubauer H, et al. Q Fever: Single-point source outbreak with high attack rates and massive number of undetected infections across an entire region. Clin Infect Dis. 2012.
Van Leuken JPG, Havelaar AH, Van der Hoek W, Ladbury GAF, Hackert VH, Swart AN. A model for the early identification of sources of airborne pathogens in an outdoor environment. PLoS ONE. 2013;8, e80412.
Van Leuken JPG, van de Kassteele J, Sauter FJ, Van der Hoek W, Heederik D, Havelaar AH, et al. Improved correlation of human Q fever incidence to modelled C. burnetii concentrations by means of an atmospheric dispersion model. Int J Health Geogr. 2015;14:14.
Van der Hoek W, Morroy G, Renders NH, Wever PC, Hermans MH, Leenders AC, et al. Epidemic Q Fever in Humans in the Netherlands. In: Coxiella burnetii: Recent Advances and New Perspectives in Research of the Q Fever Bacterium. Advances in Experimental Medicine and Biology. 2012. p. 329–64.
Whelan J, Schimmer B, De Bruin A, Robert-Du Ry Van Beest Holle M, Van der Hoek W, Ter Schegget R. Visits on “lamb-viewing days” at a sheep farm open to the public was a risk factor for Q fever in 2009. Epidemiol Infect. 2012;140:858–64.
Van den Brom R, Van Engelen E, Luttikholt S, Moll L, Van Maanen K, Vellema P. Coxiella burnetii in bulk tank milk samples from dairy goat and dairy sheep farms in The Netherlands in 2008. Vet Rec. 2012;170:310.
De Bruin A, De Groot A, De Heer L, Bok J, Wielinga PR, Hamans M, et al. Detection of Coxiella burnetii in complex matrices by using multiplex quantitative PCR during a major Q fever outbreak in The Netherlands. Appl Environ Microbiol. 2011;77:6516–23.
De Bruin A, Van der Plaats RQJ, De Heer L, Paauwe R, Schimmer B, Vellema P, et al. Detection of Coxiella burnetii DNA on small-ruminant farms during a Q fever outbreak in the Netherlands. Appl Environ Microbiol. 2012;78:1652–7.
De Bruin A, Janse I, Koning M, De Heer L, Van der Plaats RQJ, Van Leuken JPG, et al. Detection of Coxiella burnetii DNA in the environment during and after a large Q fever epidemic in the Netherlands. J Appl Microbiol. 2013;114:1395–404.
Van der Hoek W, Van de Kassteele J, Bom B, De Bruin A, Dijkstra F, Schimmer B, et al. Smooth incidence maps give valuable insight into Q fever outbreaks in The Netherlands. Geospatial Health. 2012;7:127–34.
Porten K, Rissland J, Tigges A, Broll S, Hopp W, Lunemann M, et al. A super-spreading ewe infects hundreds with Q fever at a farmers’ market in Germany. BMC Infect Dis. 2006;6:147.
Daszak P, Zambrana-Torrelio C, Bogich TL, Fernandez M, Epstein JH, Murray K, et al. Interdisciplinary approaches to understanding disease emergence: The past, present, and future drivers of Nipah virus emergence. Proc Natl Acad Sci. 2012;110:3681. doi:3688.
Zinsstag J, Schelling E, Waltner-Toews D, Tanner M. From “one medicine” to “one health” and systemic approaches to health and well-being. Preventive veterinary medicine. 2011;101:148–56.
Hatchette TF, Hudson RC, Schlech WF, Campbell N, Hatchette JE, Ratnam S, et al. Goat-associated Q fever: a new disease in Newfoundland. Emerg Infect Dis. 2001;7:413–9.
Roest HIJ, Ruuls RC, Tilburg JJHC, Nabuurs-Franssen MH, Klaassen CHW, Vellema P, et al. Molecular epidemiology of Coxiella burnetii from ruminants in Q fever outbreak, the Netherlands. Emerg Infect Dis. 2011;17:668–75.
Hogerwerf L, Van den Brom R, Roest HIJ, Bouma A, Vellema P, Pieterse M, et al. Reduction of Coxiella burnetii prevalence by vaccination of goats and sheep, The Netherlands. Emerg Infect Dis. 2011;17:379–86.
O’Neill TJ, Sargeant JM, Poljak Z. A systematic review and meta-analysis of phase I inactivated vaccines to reduce shedding of Coxiella burnetii from sheep and goats from routes of public health importance. Zoonoses and public health. 2014;61:519–33.
We thank Ben Bom for producing the maps; Maaike van Veen for critical review of the study proposal; and Ioannis Karagiannis and Marianne van der Sande for critical review of the study protocol and manuscript. This project was funded by the core funding of the RIVM/EPIET programme and the RIVM strategic research facility (SOR project S/210206/01/RQ).
The authors declare that they have no competing interests.
GL designed and implemented the study methods, except for those involving the application of the atmospheric distribution model and spatio-temporal analysis, which were designed and implemented by JvL. AS supervised the design, implementation and interpretation of the methods conducted by JvL. PV supervised and facilitated all parts of the study that related to animal data, and co-conducted the on-farm visit together with GL. BS contributed to the design of the methods. RTS provided data and background on the 2008 and 2009 outbreaks in the town. WvdH conceived of the project and provided overall supervision throughout. GL drafted the manuscript, which was reviewed and approved by all authors.
Joint first author: Georgia Ladbury and Jeroen van Leuken.
Georgia A. F. Ladbury and Jeroen P.G. Van Leuken contributed equally to this work.
Attack rates and relative risks per potential source farm. Number of cases, number of total inhabitants, attack rates (AR) per 100,000 inhabitants, and relative risks (RR) per concentric 1 km ring around each potential source farm A – Q. (XLS 29 kb)
Results of spatio-temporal analysis using an alternative emission profile. As for Fig. 5, but for emission profile lNormEpi. (PDF 564 kb)
Actual cumulative concentrations (CPPI and CPRI) per case. Cumulative concentration during the PPI (CPPI) as a function of the cumulative concentration prior to the PPI (CPRE) for each case (emission profile conYear, threshold wind velocity 2 m/s). Cases in group 1 are denoted by red closed circles [red bullet]; cases in group 2 by blue open circles [blue bullet]. (PDF 7 kb)
About this article
Cite this article
Ladbury, G.A.F., Van Leuken, J.P., Swart, A. et al. Integrating interdisciplinary methodologies for One Health: goat farm re-implicated as the probable source of an urban Q fever outbreak, the Netherlands, 2009. BMC Infect Dis 15, 372 (2015). https://doi.org/10.1186/s12879-015-1083-9