Q fever, caused by
Q fever is an occupational zoonosis caused by
A total of 3,000 DCFs housing
We used an immunofluorescence assay (IFA) (Focus Diagnostics, Cypress, CA, USA) to test serum samples for
Participating and nonparticipating farms were compared with respect to herd size; distance to nearest
Risk factors determined to be significant (p<0.20) in univariate analyses of the participant-based and farm-based data were incorporated into multivariate logistic regression and multivariate multilevel analyses, respectively. Stratified multivariate analyses for participant risk factors were performed separately for farmers and for the remaining group. Starting with a full model, manual backward elimination was performed; all variables meeting the 10.0% significance level in the likelihood ratio test were kept in the final model. Two-way interactions between biologically plausible variables in the multivariate model were investigated. Last, variables included in the final multivariate model for participant-based factors and those included in the multilevel model for farm-based factors were combined in a multivariate multilevel analysis to identify the independent risk determinants for seropositivity. The final model fit was assessed by the quasi-likelihood under the independence model criterion goodness-of -fit statistic for generalized estimation equation models. SAS version 9.2 (SAS Institute, Cary, NC, USA) was used for all analyses.
Of the 3,000 invited farms, 311 provided a BTM sample, and 755 persons from 309 (10.3%) farms participated in this study by providing a serum sample. A farm-based questionnaire was available for 736 (97.5%) persons from 301 farms, and a participant-based questionnaire was completed by 729 (96.6%) persons from 308 farms. Compared with nonparticipating farms, participating farms were a median of 1.5 km closer to small ruminant farms with
| Variable | Farms | p value | |
|---|---|---|---|
| Participating, n = 311 | Nonparticipating, n = 2,685 | ||
| Categorical, no. (%) | |||
| Farm located inside vaccination area | 83 (26.4) | 590 (21.9) | 0.08 |
| Farm region* | 0.36 | ||
| North | 80 (25.4) | 781 (29.1) | |
| East | 104 (33.7) | 911 (33.9) | |
| West | 57 (18.7) | 494 (18.3) | |
| South | 70 (22.2) | 503 (18.7) | |
| Degree of urbanization of the farm municipality | 0.77 | ||
| Moderately, strongly, or extremely (>1,000 addresses/km2) | 1 (0.3) | 17 (0.6) | |
| Hardly (500–1,000 addresses/km2) | 10 (3.2) | 94 (3.5) | |
| Not (<500 addresses/km2) | 300 (96.5) | 2,574 (95.9) | |
| Numerical, median no. | |||
| No. cows in 2008 | |||
| <1 | 35 | 35 | 0.44 |
| 1–2 | 26 | 26 | 0.65 |
| >2 | 85 | 86 | 0.16 |
| Nearest bulk tank milk positive small-ruminant farm (meters) | 9,793 | 11,301 | 0.01 |
| Goat density (animals/km2)† | |||
| Within 5-km radius | 9.2 | 6.7 | 0.27 |
| Within 10-km radius | 9.3 | 9.2 | 0.26 |
| Sheep density (animals/km2)† | |||
| Within 5-km radius | 30 | 33 | 0.04 |
| Within 10-km radius | 34 | 35 | 0.11 |
| Cattle density (animals/km2) within 5-km radius† | |||
| Including own animals | 178 | 181 | 0.29 |
| Excluding own animals | 175 | 179 | 0.27 |
| Cattle density (animals/km2) within 10-km radius† | |||
| Including own animals | 170 | 170 | 0.99 |
| Excluding own animals | 169 | 169 | 0.91 |
*North represents Groningen, Friesland, and Drenthe Provinces; East represents Gelderland, Overijssel, and Flevoland Provinces; West represents Noord–Holland, Zuid–Holland, Utrecht, and Zeeland Provinces; and South represents Limburg and Noord–Brabant Provinces. †Corrected for area in the Netherlands.
Overall
| Participant characteristic | Total no. residents/no. positive (%) | 95% CI |
|---|---|---|
| All participants | 755/544 (72.1) | 68.8–75.3 |
| Sex | ||
| M | 431/368 (85.4) | 82.0–88.7 |
| F | 323/176 (54.5) | 49.0–59.9 |
| Age, y | ||
| <35 | 169/107 (63.3) | 56.0–70.7 |
| 35–44 | 176/131 (74.4) | 67.9–80.9 |
| 45–54 | 252/185 (73.4) | 67.9–78.9 |
|
| 132/106 (80.3) | 73.4–87.2 |
| Role | ||
| Farmer | 361/315 (87.2) | 83.8–90.7 |
| Spouse | 222/121 (54.5) | 47.9–61.1 |
| Child <18 y | 52/23 (44.2) | 30.3–58.2 |
| Child | 54/40 (74.1) | 62.0–86.1 |
| Other* | 40/30 (75.0) | 61.0–89.0 |
*Represents other family members and employees.
| Variable | No. residents (% positive) | OR (95% CI) |
|---|---|---|
| Sex | ||
| M | 431 (85.4) | 4.9 (3.5–6.9) |
| F | 323 (54.5) | Reference |
| Age, y† | ||
| <35 | 169 (63.3) | Reference |
| 35–44 | 176 (74.4) | 1.7 (1.1–2.7) |
| 45–54 | 252 (73.4) | 1.6 (1.1–2.4) |
|
| 132 (80.3) | 2.4 (1.4–4.0) |
| Role | ||
| Farmer | 361 (87.3) | 5.7 (3.8–8.6) |
| Spouse | 222 (54.5) | Reference |
| Child <18 y | 52 (44.2) | 0.7 (0.4–1.2) |
| Child | 54 (74.1) | 2.4 (1.2–4.6) |
| Other‡ | 40 (75.0) | 2.5 (1.2–5.4) |
| Presence in stable | ||
| Every day | 584 (78.3) | 8.4 (2.1–32.9) |
| Every week | 117 (51.3) | 2.5 (0.6–10.0) |
| Every month | 18 (50.0) | 2.3 (0.5–12.0) |
| <1 time/month | 10 (30.0) | Reference |
| Work on the farm† | ||
| No | 42 (31.0) | Reference |
| Part time (1–39 h/wk) | 333 (60.7) | 3.3 (1.7–6.6) |
| Full time ( | 354 (88.7) | 17.5 (8.4–36.4) |
| Care for calves | ||
| Yes | 501 (78.4) | 2.5 (1.8–3.5) |
| No | 228 (59.6) | Reference |
| Milk cows | ||
| Yes | 442 (84.4) | 4.5 (3.2–6.4) |
| No | 287 (54.4) | Reference |
| Feed cattle | ||
| Yes | 515 (81.6) | 4.3 (3.0–6.0) |
| No | 214 (50.9) | Reference |
| Transport cattle | ||
| Yes | 414 (85.0) | 4.4 (3.1–6.3) |
| No | 315 (56.2) | Reference |
| Provide health care for cattle | ||
| Yes | 440 (83.0) | 3.7 (2.7–5.3) |
| No | 288 (56.6) | Reference |
| Provide birth care for cattle | ||
| Yes | 516 (80.8) | 3.8 (2.7–5.4) |
| No | 212 (52.4) | Reference |
| Remove manure | ||
| Yes | 478 (83.3) | 4.6 (3.2–6.4) |
| No | 251 (52.2) | Reference |
| Spread manure | ||
| Yes | 254 (86.6) | 3.5 (2.3–5.2) |
| No | 474 (65.2) | Reference |
| Clean stables | ||
| Yes | 505 (80.2) | 3.3 (2.3–4.6) |
| No | 222 (55.4) | Reference |
| No. of tasks (0–9) performed on farm†§ | 725 (72.7) | 1.3 (1.2–1.4) |
| Wear overalls† | ||
| Almost every time/regularly | 584 (76.7) | 2.5 (1.7–3.6) |
| Sometimes/never | 142 (57.0) | Reference |
| Wear boots† | ||
| Almost every time/regularly | 646 (74.6) | 2.1 (1.3–3.4) |
| Sometimes/never | 81 (58.0) | Reference |
| Use gloves during cattle birth care† | ||
| Fully compliant | 61 (70.5) | 0.5 (0.3–0.9) |
| Partly or noncompliant | 456 (82.2) | Reference |
| No birth care | 211 (52.4) | 0.2 (0.2–0.3) |
| Use disinfectant during cattle birth care | ||
| Fully compliant | 141 (82.3) | 1.1 (0.7–1.9) |
| Partly or noncompliant | 374 (80.5) | Reference |
| No birth care | 212 (53.3) | 0.3 (0.2–0.4) |
| Have contact with pigs at own or other farm†¶ | ||
| Yes | 83 (85.5) | 2.4 (1.3–4.6) |
| No | 646 (70.9) | Reference |
| Have contact with sheep at own or other farm†¶ | ||
| Yes | 232 (75.9) | 1.3 (0.9–1.8) |
| No | 497 (71.0) | Reference |
| Have contact with poultry at own farm†¶ | ||
| Yes | 177 (67.8) | 0.8 (0.5–1.1) |
| No | 578 (73.4) | Reference |
| Have indirect contact with rats/mice at own farm†# | ||
| Yes | 213 (82.6) | 2.2 (1.5–3.3) |
| No | 516 (68.4) | Reference |
| Have contact with dogs at own farm** | ||
| Yes | 538 (71.0) | 0.7 (0.5–1.1) |
| No | 190 (76.8) | Reference |
| Have contact with dogs at other farm†** | ||
| Yes | 273 (75.8) | 1.3 (0.9–1.8) |
| No | 456 (70.6) | Reference |
| Have contact with cows at other farm†** | ||
| Yes | 195 (82.6) | 2.2(1.4–3.3) |
| No | 532 (68.8) | Reference |
| Have indirect contact with cows at other farm# | ||
| Yes | 326 (77.6) | 1.6 (1.1–2.3) |
| No | 400 (68.3) | Reference |
| Have contact with compost† | ||
| Yes | 139 (66.2) | 0.7 (0.5–1.0) |
| No | 590 (74.1) | Reference |
| Have contact with hay† | ||
| Daily (almost) | 505 (80.2) | 3.3 (2.3–4.6) |
| Weekly or less often | 222 (55.4) | Reference |
| Have contact with cattle food† | ||
| Daily (almost) | 530 (80.9) | 4.2 (3.0–6.0) |
| Weekly or less often | 199 (50.3) | Reference |
| Have contact with raw milk† | ||
| Daily (almost) | 507 (78.9) | 2.7 (1.9–3.8) |
| Weekly or less often | 221 (57.9) | Reference |
| Have contact with sheep or goat manure† | ||
| Yes | 96 (80.2) | 1.6 (1.0–2.8) |
| No | 632 (71.4) | Reference |
| Have contact with cattle manure† | ||
| Daily (almost) | 477 (83.0) | 4.4 (3.1–6.2) |
| Weekly or less often | 252 (52.8) | Reference |
| Have contact with live–born animals | ||
| Yes | 626 (76.7) | 3.7 (2.4–5.7) |
| No | 102 (47.1) | Reference |
| Have contact with stillborn animals | ||
| Yes | 304 (83.9) | 2.9 (2.0–4.1) |
| No | 424 (64.5) | Reference |
| Have contact with placenta material or amniotic fluid† | ||
| Yes | 425 (84.2) | 4.2 (2.9–5.9) |
| No | 304 (56.3) | Reference |
| Consume Brie cheese† | ||
| Yes | 166 (63.9) | 0.6 (0.4–0.9) |
| No | 563 (75.1) | Reference |
| Lived on a farm with animals during childhood† | ||
| Yes | 604 (76.4) | 2.8 (1. 9–4.2) |
| No | 125 (53.6) | Reference |
*Analysis was performed to assess the main participant-based characteristics associated with positivity (p<0.20 in likelihood ratio test). OR, odds ratio. †Variables included in later multivariate analysis before manual backward elimination. ‡Represents other family members and employees. §No. tasks (i.e., caring for calves, milking, feeding cattle, transport of cattle, general cattle health, assistance during calving, removing manure, spreading manure and cleaning the stables) performed on the farm by 1 person (0–9 tasks). Risk increases per task. ¶See animals at <5 meters or touch animals (contact). #See animals at <5 meters (indirect contact). **Touch animals (direct contact).
Nine (1.2%) participants from 8 farms were classified as having a relatively recent infection (IgM phase II titer range 32–256). All 8 farms were within 2.5–21.2 km of the nearest
Four participants reported having had Q fever diagnosed by a physician during 2008–2010. On the basis of serum samples obtained at study entry, 3 of these participants had a serologic profile indicating past infection. These 3 participants lived in the southern or eastern region of the Netherlands on farms within a 3-km radius of the nearest small-ruminant farm with
Risk factors for seropositivity for farmers/workers and residents included age
| Variable | No. residents total (% positive) | OR (95% CI) |
|---|---|---|
| No. cows on farm in 2008†‡ | 755 (72.1) | 1.0 (1.0–1.0) |
| Nearest bulk tank milk positive small-ruminant farm† | ||
| <8 km | 331 (75.8) | 1.4 (1.0–1.9) |
|
| 424 (69.1) | Reference |
| Municipal cattle density, including beef calves§ | 755 (72.1) | 1.0 (1.0–1.0) |
| Farm location | ||
| Inside small-ruminant vaccination area | 202 (78.2) | 1.6 (1.0–2.3) |
| Outside small-ruminant vaccination area | 553 (69.8) | Reference |
| Farm region† | ||
| South | 170 (80.6) | 1.8 (1.2–2.7) |
| Other | 585 (69.6) | Reference |
| Beef cattle on the farm† | ||
| Yes | 79 (82.3) | 1.9 (1.1–3.4) |
| No | 652 (70.7) | Reference |
| Annual peak in calving | ||
| Yes | 135 (76.3) | 1.3 (0.9–2.0) |
| No | 601 (71.1) | Reference |
| Automatic milking† | ||
| Yes | 154 (65.6) | 0.7 (0.5–1.0) |
| No | 580 (73.8) | Reference |
| Use of bedding in stables | ||
| Yes | 717 (72.4) | 1.9 (1.2–2.9) |
| No | 19 (57.9) | Reference |
| Pet cat | ||
| Yes | 444 (69.1) | 0.6 (0. 5–0.9) |
| No | 285 (77.9) | Reference |
| Pet rabbit | ||
| Yes | 202 (64.4) | 0.6 (0.4–0.8) |
| No | 527 (75.7) | Reference |
| Birds in stable† | ||
| Yes | 90 (82.2) | 1.9 (1.0–3.6) |
| No | 644 (70.5) | Reference |
| Use of by-product feedstuffs† | ||
| Yes | 229 (77.3) | 1.5 (1.0–2.1) |
| No | 507 (69.6) | Reference |
| No. cows that calved in 2009‡ | 720 (71.8) | 1.0 (1.0–1.0) |
| No. live-born calves | ||
| <78 | 335 (69.0) | Reference |
|
| 344 (74.1) | 1.3 (0.9–1.8) |
| No. twin calves | ||
| 1–2 | 272 (69.9) | Reference |
|
| 313 (76.4) | 1.4 (1.0–2.0) |
| Type of farm management† | ||
| Closed herd | 515 (73.4) | Reference |
| Purchase of cattle | 213 (68.1) | 0.8 (0.6–1.1) |
| No. cattle purchase addresses in 2007† | ||
| 0 or 1 | 649 (72.7) | Reference |
|
| 76 (64.5) | 0.7 (0.4–1.0) |
| Farm boots and work clothes available for professional visitors | ||
| Yes | 662 (71.3) | 0.7 (0.4–1.1) |
| No | 74 (78.4) | Reference |
| Work clothes available for own personnel | ||
| Yes | 556 (73.6) | 1.4 (1.0–1.9) |
| No | 180 (67.2) | Reference |
*The analysis included the primary farm-based factors associated with positivity (p<0.20 in likelihood ratio test). OR, odds ratio. †Variable included in later multivariate analysis before manual backward elimination. ‡Risk increases per cow. §Risk decreases per cow.
Of the 21 variables considered in the multivariate participant model, 8 were independently associated with seropositivity: age
| Association with positivity, characteristic | OR (95% CI) |
|---|---|
| Positive association | |
| Age, y | |
| <35 | Reference |
| 35–44 | 1.4 (0.8–2.3) |
| 45–54 | 1.0 (0.6–1.6) |
|
| 1.9 (1.0–3.5) |
| Work on farm | |
| No | Reference |
| Part time (1–39 h/wk) | 2.4 (1.1–5.2) |
| Full time ( | 10.4 (4.2–25.7) |
| Contact with pigs at own or other farm | |
| Yes | 2.6 (1.2–5.4) |
| No | Reference |
| Contact with cows at other farm | |
| Yes | 1.6 (1.0–2.6) |
| No | Reference |
| Indirect contact with rats/mice at own farm | |
| Yes | 1.5 (1.0–2.4) |
| No | Reference |
| Negative association | |
| Use of gloves during cattle birth care | |
| Fully compliant | 0.4 (0.2–0.8) |
| Partly or noncompliant | Reference |
| No birth care | 0.7 (0.4–1.1) |
| Contact with poultry at own farm | |
| Yes | 0.6 (0.4–0.9) |
| No | Reference |
| Contact with compost | |
| Yes | 0.6 (0.3–0.9) |
| No | Reference |
*The analysis included the primary participant-based characteristics associated with positivity (p<0.10 in likelihood ratio test). The number of observations was 712. Model fit was assessed by use of the Hosmer–Lemeshow goodness-of-fit test (p = 0.91). OR, odds ratio.
Of the 9 variables considered in the multilevel farm model, 6 were independently associated with seropositivity; larger herd size, farm location in the southern region, beef cattle on the farm, use of food concentrate, and presence of birds in the stable were risk factors, and automatic milking was a protective factor (
| Variable | OR (95% CI) |
|---|---|
| No. cows on farm in 2008† | 1.0 (1.0–1.0) |
| Farm region | |
| South | 1.8 (1.2–2.8) |
| Other | Reference |
| Beef cattle on farm | |
| Yes | 1.7 (1.0–2.8) |
| No | Reference |
| Automatic milking | |
| Yes | 0.7 (0.4–1.0) |
| No | Reference |
| Birds in stable | |
| Yes | 2.0 (1.1–3.8) |
| No | Reference |
| Use of by-product feedstuffs | |
| Yes | 1.4 (1.0–2.0) |
| No | Reference |
*The analysis included the primary farm-based characteristics associated with positivity (p<0.10 in likelihood ratio test). The number of observations was 716; the number of levels used was 309 (quasi-likelihood under the independence model criterion 832.88). OR, odds ratio. †Risk increased per cow.
| Variable | OR (95% CI) |
|---|---|
| Age, y | |
| <35 | Reference |
| 35–44 | 1.3 (0.8–2.3) |
| 45–54 | 1.2 (0.7–2.0) |
|
| 1.9 (1.1–3.5) |
| Work on farm | |
| No | Reference |
| Part time (1–39 h/wk) | 2.5 (1.1–5.6) |
| Full time ( | 10.7 (4.2–27.0) |
| Use of gloves during cattle birth care | |
| Fully compliant | 0.4 (0.2–0.8) |
| Partly or noncompliant | Reference |
| No birth care | 0.7 (0.4–1.1) |
| Contact with pigs at own or other farm | |
| Yes | 2.4 (1.1–5.1) |
| No | Reference |
| Contact with poultry at own farm | |
| Yes | 0.5 (0.3–0.8) |
| No | Reference |
| Indirect contact with rats/mice at own farm | |
| Yes | 1.6 (1.0–2.7) |
| No | Reference |
| No. cows on farm in 2008† | 1.0 (1.0–1.0) |
| Farm region | |
| South | 1.9 (1.2–3.1) |
| Other | Reference |
| Automatic milking | |
| Yes | 0.6 (0.4–1.0) |
| No | Reference |
| Birds in stable | |
| Yes | 2.3 (1.2–4.4) |
| No | Reference |
| Contact with cows at other farm | |
| Yes | 1.8 (1.0–3.2) |
| No | Reference |
| Contact with compost | |
| Yes | 0.6 (0.3–0.9) |
| No | Reference |
| Beef cattle on the farm | |
| Yes | 1.9 (1.0–3.7) |
| No | Reference |
*The analysis included the primary participant- and farm-based characteristics associated with positivity (p<0.10 in likelihood ratio test). The number of observations was 708; the number of levels used was 309 (quasi-likelihood under the independence model criterion 695.52). OR, odds ratio. †Risk increases per cow.
The overall seroprevalence of 72.1% among DCF residents, including employees, was high, indicating a considerable lifetime risk for acquiring
Seroprevalence studies of other farmer populations, particularly dairy cattle farmers, are scarce, and, in general, it is difficult to compare international studies because of different study populations, tests, or cutoff values used. However, published seroprevalence estimates are generally lower than what we observed. A study using IFA with the same cutoff value that we used estimated a seroprevalence of 27.0% among a UK farm cohort (385 residents/workers) (
Farm residents living in the southern part of the Netherlands were more likely to be seropositive. This was not surprising because living in the south was a risk factor for dairy goat farmers (
Persons
We identified several cattle-related risk factors for seropositivity among cattle farm residents/staff: herd size, cattle contact at other farms, and presence of beef cattle on their own farm. A larger herd size could pose a risk because of an increased chance for
Protective factors included use of automatic milking and fully compliant use of gloves during birth care. Birth products of
Fully compliant use of gloves during birth care can help farmers protect themselves against
The relatively low response rate of 10.4% in this study can be explained by a general lack of motivation or awareness among cattle farmers because Q fever was mainly considered a problem among small-ruminant dairy farms. A general fear of consequences resulting from possible control measures targeting the cattle sector comparable with implemented control measures for Q fever in the small-ruminant sector might also have played a role. Study results are, however, considered representative for the Dutch dairy cattle sector because participating and nonparticipating farms were generally comparable.
The overall
These authors contributed equally to this article.
We thank all participants for their cooperation in this study; Diagnostiek Nederland for collecting blood samples at the participating farms; Jeroen Bosch Hospital, especially Jamie Meekelenkamp, for examining the blood samples; Noel Peters for sending test results to the participants’ doctors; Sanne Kelderman for mailing invitations and providing reference data; Helen Aangenend for coordinating the human data collection; and Ben Bom for generating geographic information. We also thank Jan van de Bergh, Olaf Stenvers, Rob van Oosterom, Mark Paauw, Harry Stinis, Ad de Rooij, Margo Vonk, Clementine Wijkmans, and Wim van der Hoek for their support and advice during the study and Roel Coutinho for his comments on the manuscript.
The study was funded by the Netherlands Organization for Health Research and Development (grant no. 50-50800-98-100: an integrated study on Q fever in livestock farmers and their [small] ruminants in the Netherlands) and cofinanced by the Ministry of Health, Welfare and Sport and the Ministry of Economic Affairs, Agriculture and Innovation.
Dr Schimmer is a medical epidemiologist at RIVM. Her research interests are Q fever in the Netherlands and other nonalimentary zoonoses.